声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 870|回复: 5

[综合讨论] 能帮我检查一下我的程序吗

[复制链接]
发表于 2008-4-30 09:01 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
function dy=c7impode(t,x)
dx=inline(['[0.98*cosx(1)-2*p(2)+x(3)*x(2)^2,'...
      ' 0.98*sinx(1)+x(3)*p(1)+2*x(4)*x(2)]'],
...,'p','x');
    ff=optimset;dx1=fsolve(dx,x([1,3]),ff,x);
dy=[x(2);dx1(1);x(4);dx1(2)];
>>[t,x]=ode15s('c7impode',[0,2],[1,0,0,1]);
plot(t,x)
回复
分享到:

使用道具 举报

发表于 2008-4-30 09:06 | 显示全部楼层
把你要求解的方程也贴出来。
初到本版,建议先看看置顶帖:进版必读、发帖经验,常见问题和聚宝盆。

[ 本帖最后由 sogooda 于 2008-4-30 09:07 编辑 ]
 楼主| 发表于 2008-4-30 09:25 | 显示全部楼层

回复 楼主 的帖子

g*cosx=2y''-yx' ^2
  gsinx=-yx''-2y'x'


g=0.98

[ 本帖最后由 sogooda 于 2008-4-30 10:16 编辑 ]
 楼主| 发表于 2008-4-30 09:31 | 显示全部楼层
t=0 时的初始条件是 x=1 x'=0 y=0 y'=1
发表于 2008-4-30 10:01 | 显示全部楼层

回复 4楼 的帖子

请详细说明白你的问题,不要别人问一点,你就说一点

包括你的方程,最好以图片的形式,附件上来
你的程序,尽量完整
你的出错提示

如果还有背景的话,最好一并介绍下
发表于 2008-4-30 11:09 | 显示全部楼层
原帖由 yeshuitang 于 2008-4-30 09:01 发表
function dy=c7impode(t,x)
dx=inline(['[0.98*cosx(1)-2*p(2)+x(3)*x(2)^2,'...
      ' 0.98*sinx(1)+x(3)*p(1)+2*x(4)*x(2)]'],
...,'p','x');
    ff=optimset;dx1=fsolve(dx,x([1,3]),ff,x);
dy=[x(2);dx ...
试试用相反的角度来思考问题——如果别人给你这样的东西,你能看懂吗?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-11 15:52 , Processed in 0.066460 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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