汽车系统动力学Matlab.docVIP

下载本文档

11
0
约1.03万字
约 18页
2019-01-05 发布于浙江
举报
版权申诉

汽车系统动力学Matlab.doc

1、本文档共18页，可阅读全部内容。
2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

汽车系统动力学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*(a1^2*Caf+b1^2*Car))/(m*I*uc); S=(a1+b1)^2*Caf*Car/(m*I*uc^2)+(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; end plot(real(R(:,1)),imag(R(:,1)),bo); hold a2=1.25; b2=1.25; t=1; S=[]; for uc=10:5:100 P=[m 0;0 I]; Q=[(Caf+Car)/uc,m*uc+(a2*Caf-b2*Car)/uc;(a2*Caf-b2*Car)/uc,(a2^2*Caf+b2^2*Car)/uc]; R=[Caf;a2*Caf]; A=-P^(-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 plot(S(:,1),S(:,2),*) a3=1.35; b3=1.15; for uc=10:5:100 P=[m 0;0 I]; Q=[(Caf+Car)/uc,m*uc+(a3*Caf-b3*Car)/uc; (a3*Caf-b3*Car)/uc,(a3^2*Caf+b3^2*Car)/uc]; R=[Caf;a3*Caf]; A=-P^(-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,中性转向) set(gca,FontName,Helvetica,FontSize,10) title([不同转向特性车辆在不同车速下的系统特征根],FontSize,12); 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:10 w=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 (a1^2*Caf+b1^2*Car)/uc]; Q2=[(Caf+Car)/uc m*uc+(a2*Caf-b2*Car)/uc;(a2*Caf-b2*Car)/uc (a2^2*Caf+b2^2*Car)/uc]; Q3=[(Caf+Car)/uc m*uc+(a3*Caf-b3*Car)/uc;(a3*Caf-b3*Car)/uc (a3^2*Caf+b3^2*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(2)); C(i,1)=angle(Hw2(2))*180/