MATLP运用.doc

下载文档

6
0
约3.89千字
约 7页
2017-02-09 发布于重庆
举报
版权申诉
保障服务

MATLP运用.doc

1、本文档共7页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

MATLP运用

卢彬彬工程力学一阶摄动公式：二阶摄动公式：新系统的值：取质量阵，刚度阵计算原系统特征值和特征向量 m0=[1 0 0;0 1 0;0 0 2]; k0=[2 -1 0;-1 3 -2;0 -2 2]; [x0,a0]=eig(k0,m0) x0 = 0.2691 -0.8782 -0.3954 0.5008 -0.2231 0.8363 0.5817 0.2992 -0.2685 a0 = 0.1392 0 0 0 1.7459 0 0 0 4.1149 原系统的特征值为孤立特征值在这种情况下，当、很小时，新系统的特征值及相应的均只有小变化。根据摄动理论，可将、按小参数展开成幂级数现在来确定，，，取根据一阶摄动公式求， m1=[1 -1 0;-1 2 1;0 1 2]; k1=[1 0 2;0 1 -1;2 -1 1]; for i=1:3 a1(i,i)=x0(:,i)*(k1-a0(i,i)*m1)*x0(:,i); end a1 a1 = 0.4875 0 0 0 -0.9223 0 0 0 -6.0652 for i=1:3 for b=1:3 if i~=b s1=0; y1=(x0(:,b)*(k1-a0(i,i)*m1)*x0(:,i))/(a0(i,i)-a0(b,b))*x0(:,b); s1=s1+y1; end end end for i=1:3 x1(:,i)=s1-0.5*(x0(:,i)*m1*x0(:,i))*x0(:,i); end x1 = 0.9617 1.4024 1.5502 -0.0937 0.3564 -0.5016 -0.8542 -0.4777 -0.1427 for i=1:3 a2(i,i)=x0(:,i)*(k1*x1(:,i)-a0(i,i)*m1*x1(:,i)-a1(i,i)*m0*x1(:,i)-a1(i,i)*m1*x0(:,i)); end a2 = 0.7213 0 0 0 0.4836 0 0 0 15.8159 根据二阶摄动公式求， for i=1:3 a2(i,i)=x0(:,i)*(k1*x1(:,i)-a0(i,i)*m1*x1(:,i)-a1(i,i)*m0*x1(:,i)-a1(i,i)*m1*x0(:,i)); end a2 = 0.7213 0 0 0 0.4836 0 0 0 15.8159 for i=1:3 for b=1:3 if i~=b s2=0; y2=(x0(:,b)*(k1*x1(:,i)-a0(i,i)*m1*x1(:,i)-a1(i,i)*m0*x1(:,i)-a1(i,i)*m1*x0(:,i))/(a0(i,i)-a0(b,b)))*x0(:,b); s2=s2+y2; end end end for i=1:3 x2(:,i)=s2-0.5*(x0(:,i)*m1*x1(:,i)+x1(:,i)*m0*x1(:,i)+x1(:,i)*m1*x0(:,i))*x0(:,i); end x2 = -0.7673 -0.5366 -1.5170 0.0519 -0.1363 1.0265 0.6459 0.1828 -0.0909 1.当=0.01时，一阶摄动得到的特征值和特征向量 e=0.01 y1=a0+e*a1 y1 = 0.1441 0 0 0 1.7367 0 0 0 4.0543 t1=x0+e*x1 t 1= 0.2787 -0.8642 -0.3799 0.4998 -0.2196 0.8313 0.5732 0.2944 -0.2699 二阶摄动得到的特征值和特征向量 2=a0+e*a1+e^2*a2 y2 = 0.1441 0