1.
>>[x,y]=dsolve('D2x+2*Dx=x+2*y-exp(-t)','Dy=4*x+3*y+4*exp(-t)')
或者
>>a='D2x+2*Dx=x+2*y-exp(-t)';
>>b='Dy=4*x+3*y+4*exp(-t)';
>>[x,y]=dsolve(a,b)
自变量用t就可以,不用换成x.因为Dx,Dy都是对t求导.
2.
这个是对x的导数,所以这个dsolve里面最后一项加'x'就对了.
不过用了ezplot,就没必要定义x=[0:.01:1]了,
如果用plot,那才需要你那样定义呢.
所以这样就可以
>> y=dsolve('Dy=y+2*x','y(0)=1','x')
y =
-2*x-2+3*exp(x)
>> ezplot(y,[0,1])
另外,题目要求你数值解也画在图上,
如果你没写错的话,x=[0:1]只是0和1两个点.
那你这样画,接上面的画图
>> x=0:1
x =
0 1
>> y1=subs(y,x)
y1 =
1 4.15484548537714
>> hold on
>> plot(x,y1,'o')
>> hold off