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(),ff,x);
dy=;
>>=ode15s('c7impode',,);
plot(t,x)

sogooda 发表于 2008-4-30 09:06

把你要求解的方程也贴出来。
初到本版，建议先看看置顶帖：进版必读、发帖经验，常见问题和聚宝盆。

[ 本帖最后由 sogooda 于 2008-4-30 09:07 编辑 ]

yeshuitang 发表于 2008-4-30 09:25

回复 楼主 的帖子

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


g=0.98

[ 本帖最后由 sogooda 于 2008-4-30 10:16 编辑 ]

yeshuitang 发表于 2008-4-30 09:31

t=0 时的初始条件是　x=1 x'=0 y=0 y'=1

sigma665 发表于 2008-4-30 10:01

回复 4楼 的帖子

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

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

如果还有背景的话，最好一并介绍下

eight 发表于 2008-4-30 11:09

原帖由 yeshuitang 于 2008-4-30 09:01 发表 http://www.chinavib.com/forum/images/common/back.gif
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(),ff,x);
dy=试试用相反的角度来思考问题——如果别人给你这样的东西，你能看懂吗？
页: [1]
查看完整版本: 能帮我检查一下我的程序吗