推 zupo:推倒降階有誤,我看出來是有錯的,現在人在台灣休假 10/28 23:01
→ zupo:禮拜五幫你解看看 10/28 23:01
各位版上的高手，我在解聯立微分方程式，
但是解完的結果是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