能帮我检查一下我的程序吗

yeshuitang

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)

sogooda

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

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

yeshuitang

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


g=0.98

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

yeshuitang

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

sigma665

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

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

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

eight

原帖由 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 ...

试试用相反的角度来思考问题——如果别人给你这样的东西，你能看懂吗？