ode45解二阶微分方程组的问题

rainbowolf

大家帮忙看看这个方程:
   %y(1)=x;
   %y(2)=x';
   %y(3)=x1;
  %y(4)=x2;
则dtdy是y的倒数
  dtdy(1)=y(1)';
  dtdy(2)=y(2)';
  dtdy(3)=y(3)';
  dtdy(4)=y(4)';
则:
   dtdy(1)=y(2);
  dtdy(2)=(ku-k*y(1)-N1*k*y(1)+N1*k*y(3))/m;
  dtdy(3)=.................;?
  dtdy(4)=.................;?
  后边两个化简不出来?  有能帮忙看看的吗?

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
看看help 文档，或找一本Matlab书看看就能解决的问题．
By xjzuo
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

[ 本帖最后由 xjzuo 于 2007-5-21 14:19 编辑 ]

zhlong

你应该是有x,x1和x2三个变量，所以应该设
y(1)=x;
y(2)=x';
y(3)=x1;
y(4)=x1';
y(5)=x2;
y(6)=x2';
我也没化简过，你试试！

rainbowolf

不是啊! 你看看
dtdy=[x';x'';x1';x2']的啊,只需要x,x',x1,和x2啊  y=[x;x';x1;x2]

[ 本帖最后由 ChaChing 于 2010-8-12 18:27 编辑 ]

zhlong

你给的图中的原方程不是又三个变量么？

rainbowolf

是的啊
但是如果按你的思路怎么编? 能不能写一点?  我刚学 从书上看到的  比较生硬
对拉  我主要是求出 x 的值  其中x1  x2 都是中间变量

[ 本帖最后由 ChaChing 于 2010-8-12 18:34 编辑 ]

zhlong

直接利用后面两个方程把x2消去，然后利用剩下的两个方程

[ 本帖最后由 ChaChing 于 2010-8-12 18:37 编辑 ]

rainbowolf

x2能消去吗?  好像不能吧

[ 本帖最后由 eight 于 2007-5-19 23:32 编辑 ]

zhlong

把刚贴的图中最后一个式子代入原第二个方程，那么原第二个方程不是没有了x2吗？
那你的 dtdy(3)=.................;?不就好写了吗？

它不好写不就是因为里面含x2'么？

rainbowolf

谢谢  我看看

yufeng

你编程好像有问题，导数怎么能这么表示呢？

quanwang2003

应该比较简单吧