What's wrong? Would you help me revise it?
2 views (last 30 days)
Show older comments
MATLAB有错吗? (1)dx=[dx(1);dx(2);dx(3);dx(4);dx(5);dx(6)]; % Is it wrong?
(2)[t,x]=ode45('ODE45_fun',tspan,x0); %Is it wrong?
(3)What's wrong?
The element type "name" must be terminated by the matching end-tag "</name>". Could not parse the file: d:\matlab7\toolbox\ccslink\ccslink\info.xml
(4)What's wrong? Would you help me revise it?:
Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. 点击 funfun\private\odearguments at 110 后出现110-118: f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. [m,n] = size(f0); if n > 1 error('MATLABdearguments:FoMustReturnCol',... [funstring(ode) ' must return a column vector.']) elseif m ~= neq msg = sprintf('an initial condition vector of length %d.',m); error('MATLABdearguments:SizeIC', ['Solving ' funstring(ode) ' requires ' msg]); end
(5)What's wrong? Would you help me revise it?:
Error in ==> ode45 at 173 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ... 点击 funfun\private\odearguments at 173 后出现173-178: [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ... options, threshold, rtol, normcontrol, normy, hmax, htry, htspan, ... dataType] = ... odearguments(FcnHandlesUsed, solver_name, odeFcn, tspan, y0, ... options, varargin); nfevals = nfevals + 1;
programme:
ODE45_fun:
function dx=ODE45_fun(t,x)
a1=1;a2=1;e1=9;e2=7;f1=56;f2=98;g1=76;g2=665;
a=a1+a2; e=e1+e2; f=f1+f2; g=g1+g2;
dx(1)=x(3)-x(2)+x(4)*x(2)-2*x(5)*x(2)-x(2)*x(2)-x(5)-e-f-g;
dx(2)=x(1)-x(3)+x(4)*x(2)-x(2)*x(2)-x(5)*x(5)-x(5)*x(2)-e*f-g;
dx(3)=x(2)-x(1)+x(5)*x(1)-x(4)*x(3)+x(3)+x(4)+e*g+f;
x(6)=x(4)*x(5)+x(2)*x(3)+e*f+e*g; % it's x(6),not dx(6)
dx(4)=x(5)*a+x(3)*x(1)+x(3)+e*f*g*x(6)
dx(5)=x(4)*a+x(2)*x(1)-e*f*x(1)-x(2)*x(2)*x(5)-x(2)*x(5)*x(5)-e*f-f*g
dx=[dx(1);dx(2);dx(3);dx(4);dx(5);dx(6)]; % should be revised as dx=[dx(1);dx(2);dx(3);dx(4);dx(5)]; ? But,x(6)?
ODE45_main:
tspan=[0,10];
x0=[0;0.1;2;0;0.1;2];
[t,x]=ode45('ODE45_fun',tspan,x0);
data=[t,x];
save ODE45_data.txt data -ascii
subplot(2,3,1),plot(t,x(1))
subplot(2,3,2),plot(t,x(2))
subplot(2,3,3),plot(t,x(3))
subplot(2,3,4),plot(t,x(4))
subplot(2,3,5),plot(t,x(5))
subplot(2,3,6),plot(t,x(6))
0 Comments
Answers (0)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!