liming1985 发表于 2007-12-6 08:43

请教高手下面的程序错误是啥意思

opions=odeset('Reltol',,'Abstol',)
T0=903;
P0=0.1e+6;
y0=;
=ode45('yifan',,y0,opions,N1,Neb01,rou_b,L1,qn0);
h=size(y);
e=h(1);
f=y(e,:);
X1=f(1)+f(2)+f(3);


执行结果是：

opions =
            AbsTol: 1.0000e-006
               BDF: []
            Events: []
         InitialStep: []
            Jacobian: []
         JConstant: []
            JPattern: []
                Mass: []
      MassConstant: []
      MassSingular: []
            MaxOrder: []
             MaxStep: []
         NonNegative: []
         NormControl: []
         OutputFcn: []
         OutputSel: []
            Refine: []
            RelTol: 1.0000e-004
               Stats: []
          Vectorized: []
    MStateDependence: []
         MvPattern: []
      InitialSlope: []
??? Error using ==> feval
Undefined command/function 'yifan'.
Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...
Error in ==> ZUOYE22 at 30
=ode45('yifan',,y0,opions,N1,Neb01,rou_b,L1,qn0);

[ 本帖最后由 eight 于 2007-12-6 18:39 编辑 ]

eight 发表于 2007-12-6 09:18

原帖由 liming1985 于 2007-12-6 08:43 发表 http://www.chinavib.com/forum/images/common/back.gif
下面是部分程序：
function rt=yifan(t,y,flag,n,Neb0,rou_b,L,qn0)
Kp=exp(19.76-1.537e+4/T-0.5223*log(T))
Pe=y(5)*(1-y(1)-y(2)-y(3))/(n+1+y(1)+y(2));
Ps=y(5)*(y(1))/(n+1+y(1)+y(2));
Pb=y(5)*y(2 ...
请搜索版面

yufeng 发表于 2007-12-6 09:19

如果你自己编龙格库塔没有必要嵌入函数吧 我的qq604443022 交流一下
页: [1]
查看完整版本: 请教高手下面的程序错误是啥意思