| 查看: 1589 | 回复: 3 | |||
[求助]
matlab龙格库塔法求解初值为零高阶微分方程编程数学问题请教
|
|
一： %自由落体问题的高阶偏微分方程 s=0.5*(s'')*t^2 function dy=myfun_1(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=2*y(1)/(t*t); [x,y]=ode45(@myfun_1,[0,20],[0,0]); plot(t,y(:,1),'-',t,y(:,2),'.') 问题1:初值t=0，而t在分母上； 问题2:或者什么其他算法可以实现。 难道只能把初值t=0改成一个很小的数值吗？ 二： 类似问题dy(1)=y(2); dy(2)=（右边是关于y(1)和y(2)的函数且初值为零不含非齐次项）; 用龙格库塔法y(1)始终为零怎么解决，例如我看的是黄怀纬博士论文第六章 确定m，n时，ξ系数皆为常数，图像应该一直为零，不应该有波峰。 闭关学习中，谢谢好心大侠们。  `UE@$[V9OG8Y86U]H6L{DDR.png  IC8M`SCVF][K{0CDZ0)_K_F.png |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于一篇论文里的用四阶龙格库塔法求微分方程的图，急。给充50话费。求大神
已经有4人回复
- 用四阶龙格库塔法求解常微分方程组
已经有15人回复
- 请问什么是分段龙格库塔的方法？需要用matlab软件采用这种方法求二阶微分方程的解
已经有3人回复
- mathcad怎样用四阶龙格库塔法求解二阶微分方程？
已经有8人回复
- 分数阶微分方程的求解（用差分法，matlab编程求解）
已经有9人回复
- 求matlab 用四阶龙格-库塔法求解微分方程
已经有4人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 【求助】Matlab中利用四阶龙格-库塔法求解微分方程！！！！
已经有9人回复
--八云--
专家顾问 (知名作家)
博览众家之长，终集百家之短
-
专家经验: +143 - 博学EPI: 4
- 应助: 216 (大学生)
- 金币: 4097.5
- 散金: 1408
- 红花: 70
- 帖子: 9216
- 在线: 1033.8小时
- 虫号: 1368331
- 注册: 2011-08-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 物理
2楼2015-07-01 22:53:06
3楼2015-07-02 10:08:32
--八云--
专家顾问 (知名作家)
博览众家之长，终集百家之短
-
专家经验: +143 - 博学EPI: 4
- 应助: 216 (大学生)
- 金币: 4097.5
- 散金: 1408
- 红花: 70
- 帖子: 9216
- 在线: 1033.8小时
- 虫号: 1368331
- 注册: 2011-08-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 物理
4楼2015-07-02 14:22:31
