%《MATLAB及其在理工課程中的應(yīng)用指南》（第二版）第七章例7－3－3程序
% 二自由度無阻尼振動的傳統(tǒng)建模和解耦求解程序
% 西安電子科技大學(xué)出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
m1=1; m2=9; k1 = 4; k2=2; 			% 輸入各原始參數(shù)
x0 = [1;0]; xd0 = [0;-1]; tf= 20; 
M = [m1,0;0,m2];                    % 構(gòu)成參數(shù)矩陣
K = [k1+k2, -k2; -k2, k2]; 		
[u,L] = eig (K,M );					% 求廣義特征向量和特征值
t=linspace(0,tf,101);x=zeros(2,101);% 時間分割和輸出變量初始化
for s=1:2						    % 分別處理兩個特征值
    alfa=sqrt(u(:,s)'*M*u(:,s));	% 解耦后的質(zhì)量
    u(:,s)=u(:,s)/alfa;			    % 把特征向量[ui]歸一化
    w(s)=sqrt(L(s,s));		        % 分別求對應(yīng)于兩個特征值的分量
    xt=u(:,s)*(u(:,s)'*M*x0*cos(w(s)*t)+u(:,s)'*M*xd0/w(s)*sin(w(s)*t));
    x=x+xt;					        % 把兩個分量累加
end
for r=1:2			                % 繪圖
    subplot(2,1,r)
    plot(t,x(r,:)),grid
end
set(gcf,'color','w')                % 置圖形背景色為白色