解矩陣微分方程組一例
瀏覽：386353 評論：1
摘要：方程組形式
有一個二階微分方程組：
[M]{D2X}+[K]{X}={0}
[M]--為對角矩陣n×n
[k]--為對稱矩陣n×n
程序：
% 有一個二階微分方程組：
% [M]{D2X}+[K]{X}={0}
% [M]--為對角矩陣n×n
% [k]--為對稱矩陣n×n
clear;clc;close all;
n=5;
rand('state',0);
% \copyright: zjliu
% Author's email: zjliu2001@163.com
M=rand(n);
K=rand(n);
Df=inline('[x(n+1:end,1);-inv(M)*K*x(1:n,1)]',...
't','x','flag','n','M','K');
[t,x]=ode45(Df,[0,10],rand(n,1),[],n,M,K);
plot(t,x(:,1:n));
for k=1:n;
eval(['Le',num2str(k),'=[''X',num2str(k),'''];']);
end
ss='Le1';
for k=2:n;
ss=[ss,',Le',num2str(k)];
end
eval(['legend(',ss,',0);']);
有一個二階微分方程組：
[M]{D2X}+[K]{X}={0}
[M]--為對角矩陣n×n
[k]--為對稱矩陣n×n
程序：
% 有一個二階微分方程組：
% [M]{D2X}+[K]{X}={0}
% [M]--為對角矩陣n×n
% [k]--為對稱矩陣n×n
clear;clc;close all;
n=5;
rand('state',0);
% \copyright: zjliu
% Author's email: zjliu2001@163.com
M=rand(n);
K=rand(n);
Df=inline('[x(n+1:end,1);-inv(M)*K*x(1:n,1)]',...
't','x','flag','n','M','K');
[t,x]=ode45(Df,[0,10],rand(n,1),[],n,M,K);
plot(t,x(:,1:n));
for k=1:n;
eval(['Le',num2str(k),'=[''X',num2str(k),'''];']);
end
ss='Le1';
for k=2:n;
ss=[ss,',Le',num2str(k)];
end
eval(['legend(',ss,',0);']);
感謝蘿卜網友
技術鄰APP
工程師必備
工程師必備
- 項目客服
- 培訓客服
- 平臺客服
TOP
1
