matlab中求二阶微分方程 的解析解和数值解，并在同一图形中画出解析解和数值解的图形进行比较
2个回答
展开全部
y=dsolve('(1+x^2)*D2y=2*x*Dy','y(0)=1','Dy(0)=3',x)
y =
x*(x^2 + 3) + 1
数值解法：
首先编写函数myfun.m
function dy=myfun(x,y)
dy=zeros(2,1)
dy(1)=y(2);
dy(2)=2*x*y(2)/(1+x^2)
调用ode45：
[x45,y45]=ode45('myfun',[0 10],[1 3])
plot(x45,y45,'+')
hold on
x1=0:0.5:10;
y1=subs(y,x1);
plot(x1,y1,'o')
y =
x*(x^2 + 3) + 1
数值解法：
首先编写函数myfun.m
function dy=myfun(x,y)
dy=zeros(2,1)
dy(1)=y(2);
dy(2)=2*x*y(2)/(1+x^2)
调用ode45：
[x45,y45]=ode45('myfun',[0 10],[1 3])
plot(x45,y45,'+')
hold on
x1=0:0.5:10;
y1=subs(y,x1);
plot(x1,y1,'o')
本回答被网友采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
系科仪器
2024-08-02 广告
科仪器致力于为微纳薄膜领域提供精益级测量及控制仪器，包括各种光谱椭偏、激光椭偏、反射式光谱等，从性能参数、使用体验、价格、产品可靠性及工艺拓展性等多个维度综合考量，助客户提高研发和生产效率，以及带给客户更好的使用体验。...
点击进入详情页
本回答由系科仪器提供
展开全部
我怎么求解析解它老报错啊，怪了。。。这个微分方程还算比较简单的，根本不用求数值解的方法来求解。
我再想想。
我再想想。
追问
好吧 谢啦
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
|
广告 您可能关注的内容 |
