作者h062471033 (Branco)
看板MATLAB
标题Re: [问题] 联立微分方程式
时间Thu Oct 24 15:14:48 2013
各位版上的高手，我在解联立微分方程式，
但是解完的结果是NaN，不知道哪边出问题，
是程式码有问题?还是推导降阶有误?还是方程式本身无解? 谢谢!
% {x''+x'+x=-y+c c是常数项，x(t) y(t) 前面的系数忽略
% {y''+y'+y=-x
% y(1)=x => dydt(1)= x' = y(2)
% y(2)=x' => dydt(2)= x'' = -x'-x-y+c = -y(2)-y(1)-y(3)+c
% y(3)=y => dydt(3)= y' = y(4)
% y(4)=y' => dydt(4)= y'' = -y'-y-x = -y(4)-y(3)-y(1)
function dydt=beamdisp(t,y)
a=11988;
b=2.198;
c=115972
d=12078;
e=10.99;
f=6441;
cc=32.62;
dydt(1) = y(2);
dydt(2) = -a.*y(1)-b.*y(2)-c.*y(3)+cc;
dydt(3) = y(4);
dydt(4) = -d.*y(3)-e.*y(4)-f.*y(1);
dydt = dydt';
(command window)
[t,y]=ode45(@beamdisp,[0 1],[0 0 0 0])
y1=y(:,1);
plot(t,y1,'r')
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.115.66.117
1F：推 zupo:推倒降阶有误,我看出来是有错的,现在人在台湾休假 10/28 23:01
2F：→ zupo:礼拜五帮你解看看 10/28 23:01