| 查看: 1245 | 回复: 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 |
回复此楼» 猜你喜欢
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
关于一篇论文里的用四阶龙格库塔法求微分方程的图，急。给充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