作业帮 > 数学 > 作业

用matlab完成:从x轴上不同点出发,绘制出微分方程dy/dx=x^2+y^2的解的曲线

来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/05/18 23:13:55
用matlab完成:从x轴上不同点出发,绘制出微分方程dy/dx=x^2+y^2的解的曲线
请给出matlab语句
syms y(x);
dy=diff(y);
solution=dsolve(dy==x^2+y^2,y(1)==0);%假设从x=1出发,则y(1)=0,从x=k,出发,则将y(1)==0改为y(k)==0.
xi=[1:10];
yi=subs(solution,xi);
h=plot(xi,yi);
再问: 不是每次画出一条曲线,而是通过x轴上不同点的解的曲线画在同一张图里,比如过y=0、x=0的解、过y=0、x=1的解等等
再答: 这微分方程忘了怎么解了,取值1到10感觉有点怪,不过语句应该是这么写没错,你可以吧语句里面的改为Dy==x^2,这个方程比较容易理解,然后运行下感受下结果。y1=real(y1)是取复数的实部。 clear all;clc; for i=1:10; y(i)=dsolve('Dy==x^2+y^2','y(i)==0','x'); y1(i)=subs(y(i),i); end; x=[1:10]; y1=real(y1); plot(x,y1,'ro');
再问: 我这里有理论的结果图但没有程序,试了你的代码,显示的不是正确结果……⊙﹏⊙b汗……第一次的是折线,第二次的是点图,但都与结果图不同…………不知道为神马…………
附:这是mathematica中得到的方程在边值条件y[0]=0时的图形,在其他边值条件y[k]=0条件下应该也有类似图形……原题就是要求将之放在一张图中观察……
再答: 这个难度超出我能力了,不好意思,这个方程比较特殊吧,加上那些条件,这个边值条件好像是极限的性质,我matlab也是做绘图比较多,这些东西比较少弄,你另谋高就吧,祝你好运。
再问: 好吧…………anyway,谢谢解答