二自由度悬架传递函数-多阻尼对比

clc;clear;close all; C=4700; %减振器阻尼系数，设置4组，可以对比 K= 100000; % N/m 弹簧刚度 Kt=800000; %N/m 轮胎刚度 m1=100; %kg 簧下质量 m2=1000; % kg簧上质量 num=[Kt*C,Kt*K,0]; den=[m2*m1,(m2+m1)*C,m2*K+m2*Kt+m1*K,C*Kt,K*Kt]; f=0:0.1:100; w=f*2*pi; H = freqs(num, den, w); figure semilogx(f,abs(H),'LineWidth',1); grid on; xlabel('频率 hz'); ylabel('簧上加速度/路面变化率'); %% 正弦信号输入 sys = tf(num,den); sim_f=1.5; %输入信号频率 sim_A=1; %m/s 路面高度变化率幅值 t=0:0.001:10; q_v=sim_A*sin(2*pi*sim_f*t); %输入信号 a_z2=lsim(sys,q_v,t); %运行传递函数 figure plot(t,q_v,'LineWidth',1) hold on;grid on; plot(t,a_z2,'LineWidth',1); xlabel('时间 s'); ylabel('响应值'); legend('输入','输出'); title('1.5hz输入，路面变化率到簧上加速度') %% 不同阻尼下传递曲线对比 C2=[1700,3200,4700,6200,7700]; %减振器阻尼系数，设置5组，可以对比 for i=1:5 num2=[Kt*C2(i),Kt*K,0]; den2=[m2*m1,(m2+m1)*C2(i),m2*K+m2*Kt+m1*K,C2(i)*Kt,K*Kt]; H2(:,i) = freqs(num2, den2, w); end figure for i=1:5 semilogx(f,abs(H2(:,i)),'LineWidth',1); hold on; end grid on; xlabel('频率 hz'); ylabel('簧上加速度/路面变化率'); legend('C=1700','C=3200','C=4700','C=6200','C=7700'); title('不同阻尼下，路面变化率到簧上加速度传递特性')