Question

跪求！！！！急需~~~常微分方程matlab求解~~~~求指点…………

我现在需要求dy/dx=0.2231*exp(0.0613*x)-0.0287*y的特解
解法如下~~得到通解后的特解个数与定义X个数~~跪求！！！！
y=dsolve('Dy=0.2231*exp(0.0613*x)-0.0287*y','y(0)= 0 ','x');
x=0:5;
fun=inline('0.2231*exp(0.0613*x-0.0287*y');
[x,y]=ode15s(fun,[0:5], 0 );
我是想用matlab求解微分方程的通解……刚才有点晕写错啦
继续跪求

Answer 1

dy/dx=0.2231*exp(0.0613*x)-0.0287*y
dy/dx+0.0287*y=0.2231*exp(0.0613*x)
所以可以设特解为y=C*exp(0.0613*x)
带入上式后得到：0.0613C+0.0287C=0.2231
所以C=2.4788888888888888888888888888889
所以一个特解为:y=2.4788888888888888888888888888889*exp(0.0613*x)

Answer 2

大哥，你
y=dsolve('Dy=0.2231*exp(0.0613*x)-0.0287*y','y(0)= 0 ','x')

就求出来通解了，
你还要ode15s干什么。。。
如果你要求y(x=0:5)
直接
for x=0:5
f(x+1)=eval(y);
end
即可
f(1)是y(x=0),f(2)是y(x=1)...