《汽车系统动力学Matlab作业报告》由会员分享，可在线阅读，更多相关《汽车系统动力学Matlab作业报告（18页珍藏版）》请在金锄头文库上搜索。
1、汽车系统动力学 Matlab作业报告小组成员：组内任务分配二、 Matlab 程序与图形1、不同转向特性车辆在不同车速下的系统特征根m=1000;I=1500;a1=1.15;b1=1.35;Caf=53000;Car=53000;i=1;R=;for uc=10:5:100;D=(I*(Caf+Car)+m*(a1A2*Caf+b1A2*Car)/(m*I*uc);S=(a1+b142*Caf*Car/(m*I*ucA2)+(b1*Car-a1*Caf)/I;P=1 D S;r=roots(P);R(i,1)=r(1,1);R(i,2)=r(2,1);i=i+1;endplot(real(R
2、(:,1),imag(R(:,1),bo);holda2=1.25;b2=1.25;t=1;S=;for uc=10:5:100P=m 0;0 I;Q=(Caf+Car)/uc,m*uc+(a2*Caf-b2*Car)/uc;(a2*Caf-b2*Car)/uc,(a2A2*Caf+b 2A2*Car)/uc;R=Caf;a2*Caf;A=-PA(-1)*Q;d=eig(A);i=imag(d);r=real(d);S(t,1)=r(1);S(t,2)=i(1);t=t+1;endplot(S(:,1),S(:,2),*)a3=1.35;b3=1.15;for uc=10:5:100P=m 0
3、;0 I;Q=(Caf+Car)/uc,m*uc+(a3*Caf-b3*Car)/uc; (a3*Caf-b3*Car)/uc,(a3A2*Caf+b3A2*Car)/uc;R=Caf;a3*Caf;A=-pA(-1)*Q； d=eig(A);i=imag(d);r=real(d);S(t,1)=r(1)；S(t,2)=i(1)； t=t+1;end grid on plot(S(:,1),S(:,2),d);axis(-14 2 0 3);xlabel(实轴(Re);ylabel(虚轴(Im);text(-8,2.8,不足车专向);text(0,0.2,过多转向);text(-3,0.2,中
4、性转向)set(gca,FontName,Heketica,FontSize,10)title(不同转向特性车辆在不同车速下的系统特征根,FontSize,12);E一第唱小同转向特性车特在不同钿卜的杂设特征根2.1、具有不同转向特性车辆的横摆角速度幅频和相频响应m=1000;I=1500;a1=1.15;a2=1.25;a3=1.35;b1=1.35;b2=1.25;b3=1.15;Caf=53000;Car=53000;D=;C=;M=;uc=20;i=1;for f=0.1:0.1:10w=2*pi*f;E=1 0;0 1;P=m 0;0 I;R1=Caf;a1*Caf;R2=Caf;a
5、2*Caf;R3=Caf;a3*Caf;Q1=(Caf+Car)/uc m*uc+(a1*Caf-b1*Car)/uc;(a1*Caf-b1*Car)/uc (a1A2*Caf+b1A2*Car)/uc;Q2=(Caf+Car)/uc m*uc+(a2*Caf-b2*Car)/uc;(a2*Caf-b2*Car)/uc (a2A2*Caf+b2A2*Car)/uc;Q3=(Caf+Car)/uc m*uc+(a3*Caf-b3*Car)/uc;(a3*Caf-b3*Car)/uc (a3A2*Caf+b3A2*Car)/uc;A1=-inv(P)*Q1;A2=-inv(P)*Q2;A3=-in
6、v(P)*Q3;B1=inv(P)*R1;B2=inv(P)*R2;B3=inv(P)*R3;Hw1=-inv(A1-1i*w*E)*B1;Hw2=-inv(A2-1i*w*E)*B2;Hw3=-inv(A3-1i*w*E)*B3;D(i,1)=angle(Hw1(2)*180/pi;D(i,2)=abs(Hw1(2);C(i,1)=angle(Hw2(2)*180/pi;C(i,2)=abs(Hw2(2);M(i,1)=angle(Hw3(2)*180/pi;M(i,2)=abs(Hw3(2);i=i+1;endsubplot(2,1,1)f=0.1:0.1:10;semilogx(f,D(
7、:,2),k,f,C(:,2),r-.,f,M(:,2)grid onset(gca,Xtick,0.1,0.3,1.5,3,10)set(gca,FontName,Helvetica,FontSize,10)legend(不足车t向,中性转向,过多转向)title( 具有不同转向特性车辆的横摆角速度幅频响应,FontSize,12);xlabel(频率 /Hz)ylabel(横摆角速度增益/( /S)/( )subplot(2,1,2)f=0.1:0.1:10;semilogx(f,D(:,1),k,f,C(:,1),r-.,f,M(:,1)m=1000;I=1500;a1=1.15;a2=
8、1.25;a3=1.35;b1=1.35;b2=1.25;b3=1.15;Caf=53000;Car=53000;D=;C=;M=;uc=20;i=1;for f=0.1:0.1:10w=2*pi*f;E=1 0;0 1;P=m 0;0 I;R1=Caf;a1*Caf;R2=Caf;a2*Caf;R3=Caf;a3*Caf;Q1=(Caf+Car)/uc m*uc+(a1*Caf-b1*Car)/uc;(a1*Caf-b1*Car)/uc (a1A2*Caf+b1A2*Car)/uc;Q2=(Caf+Car)/uc m*uc+(a2*Caf-b2*Car)/uc;(a2*Caf-b2*Car)
9、/uc (a2A2*Caf+b2A2*Car)/uc;Q3=(Caf+Car)/uc m*uc+(a3*Caf-b3*Car)/uc;(a3*Caf-b3*Car)/uc (a3A2*Caf+b3A2*Car)/uc;A1=-inv(P)*Q1;A2=-inv(P)*Q2;A3=-inv(P)*Q3;B1=inv(P)*R1;B2=inv(P)*R2;B3=inv(P)*R3;Hw1=-inv(A1-1i*w*E)*B1;Hw2=-inv(A2-1i*w*E)*B2;Hw3=-inv(A3-1i*w*E)*B3;D(i,1)=angle(Hw1(2)*180/pi;D(i,2)=abs(Hw1
10、(2);C(i,1)=angle(Hw2(2)*180/pi;C(i,2)=abs(Hw2(2);M(i,1)=angle(Hw3(2)*180/pi;M(i,2)=abs(Hw3(2); i=i+1;endsubplot(2,1,1)f=0.1:0.1:10;semilogx(f,D(:,2),k,f,C(:,2),r-.,f,M(:,2)grid onset(gca,Xtick,0.1,0.3,1.5,3,10)set(gca,FontName,Helvetica,FontSize,10) legend(不足车t向,中性转向,过多转向) title( 具有不同转向特性车辆的横摆角速度幅频响
11、应xlabel(频率 /Hz)ylabel(横摆角速度增益/( /S)/( )subplot(2,1,2)f=0.1:0.1:10;semilogx(f,D(:,1),k,f,C(:,1),r-.,f,M(:,1)grid onset(gca,Xtick,0.1,0.3,1.5,3,10)set(gca,FontName,Helvetica,FontSize,10)legend(不足车t向,中性转向,过多转向)title( 具有不同转向特性车辆的横摆角速度相频响应xlabel(频率 /Hz)ylabel(横摆角速度相位()grid onset(gca,Xtick,0.1,0.3,1.5,3,1
12、0)set(gca,FontName,Helvetica,FontSize,10)legend(不足车t向,中性转向,过多转向)title( 具有不同转向特性车辆的横摆角速度相频响应xlabel(频率 /Hz),FontSize,12);,FontSize,12);,FontSize,12);ylabel(横摆角速度相位()H行不同转向特件4刊的横找条速唆喟有不附具有不同转向特件车辆的横摆俸速度相频晌应0.31.53频率fHz2.2、 具有不同转向特性车辆的侧向加速度和相频响应m=1000;I=1500;a1=1.15;a2=1.25;a3=1.35;b1=1.35;b2=1.25;b3=1.
13、15;Caf=53000;Car=53000;D=;C=;M=;uc=20;L=a1+b1;i=1;for f=0.1:0.1:10w=2*pi*f;E=1 0;0 1;Vi=w*I*Caf;Rr=L*Caf*Car/uc;Vr1=(L*b1*Caf*Car/uc-m*a1*Caf*uc);Vr2=(L*b2*Caf*Car/uc-m*a2*Caf*uc);Vr3=(L*b3*Caf*Car/uc-m*a3*Caf*uc);Ri1=w*m*a1*Caf;Ri2=w*m*a2*Caf;Ri3=w*m*a3*Caf;Dr1=-wA2*m*I+LA2*Caf*Car/(ucA2)+m*(b1*Car
14、-a1*Caf);Dr2=-wA2*m*I+LA2*Caf*Car/(ucA2)+m*(b2*Car-a2*Caf);Dr3=-wA2*m*I+LA2*Caf*Car/(ucA2)+m*(b3*Car-a3*Caf);Di1=w*(I*(Caf+Car)+m*(a1A2*Caf+b1A2*Car)/uc;Di2=w*(I*(Caf+Car)+m*(a2A2*Caf+b2A2*Car)/uc;Di3=w*(I*(Caf+Car)+m*(a3A2*Caf+b3A2*Car)/uc;Hay1=1i*w*(Vr1+1i*Vi)/(Dr1+1i*Di1)+uc*(Rr+1i*Ri1)/(Dr1+1i*Di1);Hay2=1i*w*(Vr2+1i*Vi)/(Dr2+1i*Di2)+uc*(Rr+1i*Ri2)/(Dr2+1i*Di2);Hay3=1i*w*(Vr3+1i*Vi)/(Dr3+1i*Di3)+uc*(Rr+1i*Ri3)/(Dr3+1i*Di3);D(i,1)=angle(Hay1)*180
