得到传递函数Gtf
syms t s a
f1 = expm(a*t)
F1 = laplace(f1)
sys1=tf([-2 3.6 -32 1.7],[1 3.3 16.8 25.4 40 13 ]);
sys2=tf([-2 3.6],[1 3.3 16.8 25.4]);
sys3=tf([-2],[1 3.3 16.8]);
mol1 = 5;
mol2 = [1 0];
mol3 = [1 2];
mol = conv(mol1, conv(mol2, mol3));
den1 = [1 0];
den2 = [0.1 1];
den3 = [1 2 4];
den = conv(den1, conv(den2, den3));
Gtf = tf(mol, den)
Eigenvalue = eig(Gtf)
Zeros = zero(Gtf)
Poles = pole(Gtf)
pzmap(Gtf);
Key = isstable(Gtf)
GG = feedback(Gtf,1)
用传递来转成伯德图(幅相特性曲线)
bode(mol, den)
nyquist(mol, den)
title('Bode Diagram RGD')
margin(Gtf);
[Gm,Pm,Wcg,Wcp] = margin(Gtf)
grid on
转到离散域,方便DSP实现(差分方程的形式)
Gtf = tf(mol, den)
Ztf = c2d(Gtf,0.001,'tustin')