Erreur lors du passage des arguments à des poignées fonction

Question

Supposons que nous ayons la fonction suivante:

function f=lorenz(t,x,a,b,c) 
    % solve differential equation like this 
    %dx/dt=a*(y-x) 
    %dy/dt=-x*z+b*x-y 
    %dz/dt=xy-c*z/3 
    f=zeros(3,1);% preallocate result 
    f(1)=a*(x(2)-x(1)); 
    f(2)=-x(1)*x(3)+b*x(1)-x(2); 
    f(3)=x(1)*x(2)-c*x(3)/3; 
    end 

Pour l'exécution de ce programme, laissez-nous utiliser le fichier test suivant:

% test program 
x0=[-2 -3.5 21];% initial point 
a=input(' enter first coefficient : '); 
b=input(' enter second coefficient: '); 
c=input(' enter third coefficient : '); 
[t,x]=ode45(@(x) lorenz(x,a,b,c),[0 10],x0); 
plot(t,x(:,1),'r'); 
title(' solution of x part'); 
grid on 

J'ai essayé passer des arguments à la poignée de fonction,

test_program 
enter first coefficient : 10 
enter second coefficient: 28 
enter third coefficient : -8 

mais il me donne l'erreur suivante:

Error using @(x)lorenz(x,a,b,c) 
Too many input arguments. 

Error in odearguments (line 87) 
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. 

Error in ode45 (line 113) 
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... 

Error in test_program (line 6) 
[t,x]=ode45(@(x) lorenz(x,a,b,c),[0 10],x0); 

une solution consiste à utiliser des variables globales, comme ceci:

function f=lorenz(t,x) 
    % solve differential equation like this 
    %dx/dt=a*(y-x) 
    %dy/dt=-x*z+b*x-y 
    %dz/dt=xy-c*z/3 
    global a 
    global b 
    global c 
    f=zeros(3,1);% preallocate result 
    f(1)=a*(x(2)-x(1)); 
    f(2)=-x(1)*x(3)+b*x(1)-x(2); 
    f(3)=x(1)*x(2)-c*x(3)/3; 
    end 

, mais il prend trop de temps à courir.

Comment puis-je résoudre ce problème? Ce que je veux est de passer des arguments différents, si j'écris dans le code quelque chose comme ça

a=input('enter the coefficient : ') 

alors cela sera répété plusieurs fois.

Answer 1

N'utilisez pas de variables globales.

Le correctif est très simple, ajoutez le t en entrée trop:

[t,x] = ode45(@(t,x) lorenz(t,x,a,b,c),[0 10],x0);