（6）验证基层控制方程是否正确0 1 图7 简支梁边界条件分别如下：以如图7的简支梁为例进行验证 ， 由于 ， 故 ， 则式可以化简成如下的齐次方程：（其中为传递矩阵）令 ，。
则式可以写成将式和式代入有：取式 。

33、的2、4、6列 ， 简化得：最后可以求出 。
指数矩阵的精细算法在结构动力、优化控制等问题中 ， 通过变换都可以将运动（控制）微分方程写成状态向量形式(2-20)式中 ， 为阶状态向量；是阶常数矩阵；为阶载荷向量（或控制微量） 。
当时 ， 一阶线性常系数齐次微分方程组的解可写成（2-21）当积分步长时 ， 指数矩阵（或传递矩阵）为(2-22)因此 ， 如何精确地求得指数矩阵（或传递矩阵）的值 ， 就成为这类方法的核心 。
钟氏精细算法中其要点是利用加法定理 ， 取 ， 将矩阵缩小后 ， 保证用泰勒级数展开计算的可靠性 。
（2-23）式中 ， ；表示单位矩阵 ， 同时将式（2-23）作如下形式的分解（2-24）由于(2-25)因此式（2-24）和（2-2 。

34、5）相当于循环语句for (;
)(2-26)当N次循环结束后有(2-27)T=exp(A)能获得高精度计算结果的根本原因是： 数值计算的相对误差不随递推过程的进行而扩散 。
指数矩阵T=exp(A)的计算精度取决于的计算精度以及矩阵的谱半径和积分步长的大小 ， 所以通过选取适当的N(m=2)和积分步长的大小 ，能够使计算结果达到很高的精度 ， 甚至达到计算机所能表达的满精度 。
%基梁静力学验证clear allclcL=0.1;
b=0.008;
h1=0.008;
E1=20.6e+10;
rou1=7800;
omg2=0;
%=G=Matrix_G(h1,L,E1,rou1,omg2);
%=% %简支梁% 。

35、 T1=exppm(G*0.5,20);
% T2=exppm(G*0.5,20);
% T12=zeros(6,1);
% T12(5,1)=-1000;
% % T_tol=T2*T1;
% TF=T2*T12;
% % T_select=T_tol(2,6,3,5);
% TF_select=TF(2,6,1);
% % Y_ini=zeros(6,1);
% Y_ini(3,5,1)=-inv(T_select)*TF_select;
% Y_part_left=T1*Y_ini;
% Y_part_right=Y_part_left+T12;
% Y_end=T2*Y_part_right;
%=% %简支梁 。

36、T1=exppm(G*0.5,20);
T2=exppm(G*0.5,20);
T12=zeros(6,1);
F=1e6*L/E1/h1T12(5,1)=-F;
T_tol=T2*T1;
TF=T2*T12;
T_select=T_tol(2,6,3,5);
TF_select=TF(2,6,1);
Y_ini=zeros(6,1);
Y_ini(3,5,1)=-inv(T_select)*TF_select;
Y_part_left=T1*Y_ini;
Y_part_right=Y_part_left+T12;
Y_end=T2*Y_part_right;
%=% % 悬臂梁% T=exppm(G,20);
% Y_ini=zeros(6,1);
% Y_end=zeros(6,1);
% F=1e6*L/E1/h1% Y_end(5,1)=-F;
% % Y_ini(4:6,1)=inv(T(4:6,4:6)*0;
F;
0;
% Y_end=T*Y_ini 。

来源：(未知)

【学习资料】网址：/a/2021/0413/0021924279.html

标题：约束|约束层阻尼梁动力学特性研究——毕业论文( 六 )