有对matlab比较熟的么？--ode32 ode45解微分方程组

foolboy007UID: 28322

解微分方程组
2F+H'=0                     (1)
F^2+F'H-G^2-F''=0      (2)
2FG+HG'-G''=0             (3)

P'+H'H-H''=0               (4)

==========定义函数=============
function dy=od45(t,y);
% y1=F  y2=F' y3=G  y4=-G' y5=-H  
dy=zeros(5,1);
dy(1)=y(2)                                    % y1'=y2
dy(2)=y(1)^2-y(2)*y(5)-y(3)^2;      %(2)
dy(3)=-y(4);                                  % y3'=-y4      
dy(4)=2*y(1)*y(3)+y(5)+y(4);         %(3)
dy(5)=-2*y(1);                             % (1)
==========================

主程序

clear;

[T,Y]=ode45(@od45,[0:.1:10],[0 0.51 1 0.616 0]);


现在的问题是 matlab倒是能运行 可是运行一会儿就报错
再就是 方程(4) 怎么搞啊？

husterUID: 2110691

晕
我们用过这个

是Singal&System里面分析电路用的，解方程没用过，还是帮老大顶顶

[ 本帖最后由 huster 于 2006-10-3 22:59 编辑 ]

蝴蝶结UID: 195566

问学数学的，我就准备找个数学系的教我matlab，上次用excel和origin画图，要疯掉了。。。

foolboy007UID: 28322

origin 和 mathematica 功能比matlab强大
也比较复杂

我现在matlab还是靠在国内学的Tubor C的那点底子琢磨出来的

nmlyUID: 109487

dont understant why u didnt just carry on ur settings to define H'' and P

[ 本帖最后由 nmly 于 2006-10-3 23:38 编辑 ]

foolboy007UID: 28322

原帖由 nmly 于 2006-10-3 23:37 发表
dont understant why u didnt just carry on ur settings to define H'' and P

那是因为我没有H' 和P的初值
没有boundary condition 怎么算？
我只有5个初值

nmlyUID: 109487

you have the initial value of F, right? then H' is not a problem

but, the initial value of P is necessary anyway......

nmlyUID: 109487

but probably, u have a final value of P?

foolboy007UID: 28322

错了 我有P的初值 是0

foolboy007UID: 28322

哦  我明白了
尝试更新代码
谢谢哈

nmlyUID: 109487

:L
okay, happy ending, just lots of trivial coding work of rk method left le....