Tustin变换，连续传递函数离散化

举例

上述说明了，不要盲目相信ChatGPT

f_vehspd = 12;
phase_vehspd = 120;
f_res = f_vehspd / tan((360-phase_vehspd) * pi/(4*180) );%连续传递函数
num = [1 -2*f_res f_res^2];
den = [1 2*f_res f_res^2];
sys = tf(num,den);
bode(sys);%离散传递函数
Ts = 0.002;
B0z = 4-4*f_res*Ts+(f_res^2)*(Ts^2);
B1z = 2*(f_res^2)*(Ts^2)-8;
B2z = 4+4*f_res*Ts+(f_res^2)*(Ts^2);
A0z = 4+4*f_res*Ts+(f_res^2)*(Ts^2);
A1z = 2*(f_res^2)*(Ts^2)-8;
A2z = 4-4*f_res*Ts+(f_res^2)*(Ts^2);
num_z = [B0z B1z B2z];
den_z = [A0z A1z A2z];
sysz = tf(num_z,den_z,Ts);
bode(sysz);

连续传递函数bode图

离散后bode图