如何在Scipy中定义带延时的LTI系统？

具有时间延迟的LTI系统的传递函数具有分子项exp（-Td * s），其中Td是时间延迟。 在Matlab中，人们可以创建在许多方面这样的LTI系统（例如，使用“S”运营商和直接设置的指数项或通过设置inputdelay outputdelay的性质tf的对象。）但是，我找不到任何方式做到这一点Scipy Signal LTI对象。 我也检查了Python控制系统库，但仍然找不到方法。

我不想使用时间延迟的Pade逼近，并且想要为LTI系统设置精确的时间延迟。

有谁知道如何在Scipy或任何其他外部Python库中实现？

---

我在github上检出了ltisys模块，并尝试创建一个延时的LTI类。 我认为，如果用BU（t-Td）代替BU（t），其中Td是时间延迟，那么在状态方程中引入输入时间延迟应该是直接的。 以下方法适用于单输入单输出系统。 可能不是没有错误，但它解决了我的目的。

#Inherit the parent LTI class  to create LTI class with time delay 


class ltidelay(lti):
    def __init__(self,inputdelay,*args,**kwargs):
        super(ltidelay,self).__init__(*args,**kwargs)    
        self.d =inputdelay

#define a method to simulate LTI with time delay . just copied lsim2 and made 2 changes. 1. passed the delay from the `ltidelay` object and 2. modified the state equation.


def lsim3(system , U=None, T=None,X0=None, **kwargs):
    if isinstance(system,lti):
        sys = system
    else:
        sys = lti(*system)
    delay = sys.d
    if X0 is  None:
        X0 = zeros(sys.B.shape[0],sys.A.dtype)        
    if T is None:
        T = linspace(0,10,101)
    T = atleast_1d(T)
    if len(T.shape) != 1:
        raise ValueError("T must be a rank1 array")
    if U is not None:
        U = atleast_1d(U)
        if len(U.shape)==1:
            U=U.reshape(-1,1)
        sU = U.shape
        if sU[0] != len(T):
            raise ValueError("U must have the same number of rows as elements in T")
        if sU[1] != sys.inputs:
            raise ValueError("The number of inputs in U is not compatible")
        ufunc = interpolate.interp1d(T, U, kind ='linear',axis =0,bounds_error =False)
        def fprime(x,t,sys,ufunc):
            return  dot(sys.A,x)+squeeze(dot(sys.B,nan_to_num(ufunc([t-delay]))))
        xout = odeint(fprime,X0,T,args=(sys,ufunc),**kwargs)
        yout = dot(sys.C,transpose(xout))
    else:
        def fprime(x,t,sys):
            return dot(sys.A,x)
        xout = odeint(fprime,X0,T,args=(sys,),**kwargs)
        yout = dot(sys.C, transpose(xout))
    return T , squeeze(transpose(yout)),xout   

#create an LTI system with delay 10

 tf = ltidelay(10,2,[4,1])

#create a step signal and time vector to simulate the LTI and check


u = linspace(0,0,100)

u[50:100] = 1

 t = linspace(1,100,100)

#check the simulation
y = lsim3(tf,u,t,X0 =0)

plot(y[1])

# compare with  LTI without time delay
y1 =lsim2(tf, u,t, X0=0)

plot(y1[1])

#delay works

链接地址: http://www.djcxy.com/p/63745.html

上一篇: How to define LTI systems with Time delay in Scipy?

下一篇: Paypal integration to accept credit card payment directly