振型矩阵与正则振型矩阵

最新推荐文章于 2023-04-06 15:23:44 发布

lijil168

最新推荐文章于 2023-04-06 15:23:44 发布

阅读量6.9k

点赞数 1

分类专栏： matlab

本文链接：https://blog.csdn.net/lijil168/article/details/67644068

版权

matlab 专栏收录该内容

56 篇文章 20 订阅

订阅专栏

%由振型矩阵mode_vec（特征列向量组成）、质量矩阵M求正则振型矩阵mode_vec_n
function mode_vec_n= mode_n(mode_vec,M)
M_p=mode_vec'*M*mode_vec; %求主质量阵
n=length(mode_vec);
mode_vec_n=zeros(n);
for i=1:n
mode_vec_n(:,i)=mode_vec(:,i)/sqrt(M_p(i,i));
end
%mode_vec_n'*M*mode_vec_n=I
%mode_vec_n^(-1)=mode_vec_n'*M
%I*q_n"+[w_n^2]*q_n=[0]
%x=mode_vec_n*q_n
%正则坐标 q_n=mode_vec_n^(-1) *x

%%%%%%%%---------------------------------------------------------------------------------------------------

%SH760小轿车空载主要参数
clear;
m=1340;
a=1.54;
b=1.29;
l=a+b;
Ic=2395; %绕质心的转动惯量
rou=sqrt(Ic/m);
k1=40*1000;
k2=44*1000;
M=[m*(b^2+rou^2)/l^2,m*(a*b-rou^2)/l^2;m*(a*b-rou^2)/l^2,m*(a^2+rou^2)/l^2];
K=[k1,0;0,k2];

%用matlab特征值分解法求主振型------------------------------------------------

[eig_vec,eig_val] = eig(inv(M)*K);