如何在倍频程中划分两个传递函数?我想计算整体传递函数:A = CPTd / 1+CPTe+CP,其中 C、P、Td 和 Te 本身都是传递函数。
我的代码:
# not a function file:
Wn = 4000;
DampingRatio = 2^-0.5;
P=tf([1],[0.001 1]);
C=tf([Wn^2],[1 2*DampingRatio*Wn 0]) / P;
Td=tf([-0.0001 1],[0.0001 1]);
Te=tf([-0.00001 1],[0.00001 1]);
CP=sysmult(C,P);
A_num=sysmult(CP,Td);
A_denum=sysadd(tf([1],[1]), sysadd(sysmult(CP,Te),CP));
A = Anum / A_denum;
sysout(A, "tf");
提前致谢。