找回密码
 注册
关于网站域名变更的通知
查看: 384|回复: 2
打印 上一主题 下一主题

求y(t)数值解并绘制其图像,与解析解对比精度。

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-5-15 13:36 | 只看该作者 回帖奖励 |正序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
% V& {8 K3 |4 n" [6 h/ U4 y9 P% r0 k
微分方程y^'' (x)-(2-1/x) y^' (x)+(1-1/x)y(x)=x^2 e^(-5x),满足边界条件y(1)=π,y(π)=1,求y(t)数值解并绘制其图像,与解析解对比精度。

该用户从未签到

3#
发表于 2020-5-15 16:04 | 只看该作者
典型的边值微分方程问题(BVP ODE),Matlab可用bvp4c命令试试,下面使用更简单的1stOpt求解:0 B) Z" F9 g0 l+ d4 g" n  {
9 b( O$ y* l  w2 @: w
Variable x=[1,pi],y=[pi,1];
' z9 H! Q! n- h$ r$ b' R- ]Plot x[x],y,y'[y2];
2 v, R, f5 y  i( U0 c. l- u" DODEFunction y''=x^2*exp(-5*x)-(1-1/x)*y+(2-1/x)*y';

该用户从未签到

2#
发表于 2020-5-15 14:46 | 只看该作者
你这是带边界条件的微分方程,采用有限差分法求解数值解,采用dsolve求解解析解。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-24 11:11 , Processed in 0.140625 second(s), 25 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表