24小时热门版块排行榜    

查看: 307  |  回复: 1

nwpusolid

铁虫 (正式写手)

[求助] 大神帮忙看看我编写的这个程序对不对??

编写一个求解常微分方程的程序,函数如下所示


以下是自己编写的函数
function [dfy]=mytt(t,fy)
dfy=[fy(2);4*(-fy(1)-t)/(sqrt(1-fy(1)^2))];

在commond窗口调用ode45函数进行计算
[t,yy]=ode45('mytt',[-1,0],[1;0])
但是自己的结果确实完成不符合论文中的结果的 yy中有很多nan的数值。。
不知道自己的程序问题出在了哪里,希望大神可以帮忙看看 多谢

图2是理论结果,程序中yita自己取得是4.但是不影响

1.jpg



2.jpg
回复此楼

» 猜你喜欢

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

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

cyberdstar

新虫 (初入文坛)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与! 2013-04-21 19:21:30
你的方程输入参数有问题啊,还有你的方程在初始条件下奇异处理下就好了啊
2楼2013-04-21 10:03:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nwpusolid 的主题更新
信息提示
请填处理意见