[問題] matlab級數解已回收
請問大家一下
因為最近享用MATLAB解級數解的問題
現在問題如下
有一個ODE為y''+y=0 y(0)=1 y'(0)=1
求對0展開的級數解
我的作法是
syms x y a0 a1 a2 a3 a4 a5 a6 a7 a8 a9 a10
a = [a0 a1 a2 a3 a4 a5 a6 a7 a8 a9 a10];
y = sum(a.*(x).^[0:10])
y =
a0+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10
>> dy=diff(y);
>> d2y=diff(dy);
>> ode=collect(d2y+y,x)
ode =
a10*x^10+a9*x^9+(90*a10+a8)*x^8+(a7+72*a9)*x^7+(a6+56*a8)*x^6+(a5+42*a7)*x^5
+(30*a6+a4)*x^4+(a3+20*a5)*x^3+(a2+12*a4)*x^2+(a1+6*a3)*x+2*a2+a0
(很長.......)
>> initcond1=strcat(char(subs(y,x,0)),'=1')
initcond1 =
a0 = 1
and
>> initcond2=strcat(char(subs(dy,x,0)),'=1')
initcond2 =
a1 = 1
到了這步(我的問題所在)(這裡的方程式是要每個都打還是有辦法用matlab算出來
不必一個一個輸入進去)
>> odecoeff1eq = '2*a2+a0=0';
>> odecoeff2eq = '6*a3+a1=0';
還有一個問題想請問大家
如果他沒給初始條件matlab還可以解嗎??
希望大家可以幫幫我!!!謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 122.147.16.139
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
