function xdot = OD(t,x)
global A B C Da K
Aa=A+(B*K*C);
Ba=B*Da;
u=K*x;
xdot=(Aa*x)+(Ba*u);
end
When running this code it show error in the line u=K*x; as not enough inputs. Here A,B,C,Da,K are all 2x2 matrix.

8 commentaires

Vashist Hegde
Vashist Hegde le 30 Juin 2020
The code looks fine. The error does not seem to be for line:5
Are you sure that while running this function, you are giving t,x as arguments?
Sargondjani
Sargondjani le 30 Juin 2020
Look at the workspace when you get the error (with "dbstop if error"). Are K and x defined?
Rik
Rik le 30 Juin 2020
Did you run this function with the green button (or with the f5 hotkey)?
Reeshabh Kumar Ranjan
Reeshabh Kumar Ranjan le 30 Juin 2020
Did you save the file before running it again?
Stephen23
Stephen23 le 1 Juil 2020
Rather than using global variables you should parameterize the function:
Gopika R
Gopika R le 1 Juil 2020
@Vashist Hedge: Yes, the arguments are t and x.
Gopika R
Gopika R le 1 Juil 2020
@Rik: I ran it with the play button, the green one.
Walter Roberson
Walter Roberson le 1 Juil 2020
When you use the green run button, then where are you expecting that MATLAB will look for the values of t and x to use inside the routine?

Connectez-vous pour commenter.

 Réponse acceptée

Ameer Hamza
Ameer Hamza le 1 Juil 2020

0 votes

It seems that you are trying to solve a matrix differential equation. If you are using ode45, one possible way this error can occur if you call it like this
ode45(@(t,x) OD, [0 10], [1; -1; 0; 0])
The correct syntax is
ode45(@(t,x) OD(t,x), [0 10], [1; -1; 0; 0]) % or
ode45(@OD, [0 10], [1; -1; 0; 0])
The following shows working code with random values for matrices
global A B C Da K
A = rand(4);
B = rand(4);
C = rand(4);
Da = rand(4);
K = rand(4);
ode45(@(t,x) OD, [0 10], [1; -1; 0; 0])
function xdot = OD(t,x)
global A B C Da K
Aa=A+(B*K*C);
Ba=B*Da;
u=K*x;
xdot=(Aa*x)+(Ba*u);
end
Also, using global variables is not a good coding practice. Following shows how to pass these matrices to function OD without using global
A = rand(4);
B = rand(4);
C = rand(4);
Da = rand(4);
K = rand(4);
ode45(@(t,x) OD(t,x,A,B,C,Da,K), [0 10], [1; -1; 0; 0])
function xdot = OD(t,x,A,B,C,Da,K)
Aa=A+(B*K*C);
Ba=B*Da;
u=K*x;
xdot=(Aa*x)+(Ba*u);
end

Plus de réponses (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by