24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1247  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

何叶知秋mh

新虫 (初入文坛)

[求助] 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时,ξ系数皆为常数,图像应该一直为零,不应该有波峰。
闭关学习中,谢谢好心大侠们。

matlab龙格库塔法求解初值为零高阶微分方程编程数学问题请教
`UE@$[V9OG8Y86U]H6L{DDR.png


matlab龙格库塔法求解初值为零高阶微分方程编程数学问题请教-1
IC8M`SCVF][K{0CDZ0)_K_F.png

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何叶知秋mh

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by --八云-- at 2015-07-01 22:53:06
一口气看完你的标题差点断气了。。。。

写出来标题的我可以当杀手了.......
3楼2015-07-02 10:08:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

--八云--

专家顾问 (知名作家)

博览众家之长,终集百家之短

【答案】应助回帖

一口气看完你的标题差点断气了。。。。
人原本是不死的,后来他创造了死亡,故而人们歌颂他。
2楼2015-07-01 22:53:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

--八云--

专家顾问 (知名作家)

博览众家之长,终集百家之短

引用回帖:
3楼: Originally posted by 何叶知秋mh at 2015-07-02 04:08:32
写出来标题的我可以当杀手了..........

嘿嘿
人原本是不死的,后来他创造了死亡,故而人们歌颂他。
4楼2015-07-02 14:22:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见