PadeDelay.hpp

Go to the documentation of this file.

#pragma once
 
#include "math.hpp"
#include "ContinuousTransferFunction.hpp"
 
namespace controlpp
{
 
    template<class T, int NumOrder, int DenOrder=NumOrder>
    requires(NumOrder <= DenOrder)
    ContinuousTransferFunction<T, NumOrder, DenOrder> pade_delay(T delay){
        ContinuousTransferFunction<T, NumOrder, DenOrder> result;
        
        for(int k = 0; k < (NumOrder+1); ++k){
            result.num().at(k) = pade_num_param(NumOrder, DenOrder, k) * controlpp::pow(-delay, k);
        }
 
        for(int k = 0; k < (DenOrder+1); ++k){
            result.den().at(k) = pade_den_param(NumOrder, DenOrder, k) * controlpp::pow(delay, k);
        }
 
        return result;
    }
 
} // namespace controlpp