我的意思是吧sys用双线性变换离散化为sys1,然后在同一个figure里面画出两者的伯德图并比较离散化之后的差异。但是,不知道我编的程序和原程序有什么区别,画出来的图为什么不一样呢?我感觉两者表达的数学式是一样的
我编的程序:
num = [1];
den = [1 0.8 1];
sys = tf(num,den)
T = 1;
sysd1 = c2d(sys,T,'tustin')
w = 0:0.01:10;
dbode(sysd1,T,w);
书上的程序:
w=0:0.01:10;
num1=[0.1515 0.3030 0.1515];den1=[1 -0.9090 0.5132];
[m1,p1]=dbode(num1,den1,1,w);
num0=[1];den0=[1 0.8 1];
[m0,p0]=bode(num0,den0,w);
figure(1)
plot(w,m0,'-',w,m1,'--'),title('BodeDiagrames'),xlabel('Frequency(rad/s)'),
ylabel('Magnitude(dB)'),grid on
标准程序参见高金源《计算机控制系统》p114。