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 Kommentare

Vashist Hegde
Vashist Hegde am 30 Jun. 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 am 30 Jun. 2020
Look at the workspace when you get the error (with "dbstop if error"). Are K and x defined?
Rik
Rik am 30 Jun. 2020
Did you run this function with the green button (or with the f5 hotkey)?
Reeshabh Kumar Ranjan
Reeshabh Kumar Ranjan am 30 Jun. 2020
Did you save the file before running it again?
Stephen23
Stephen23 am 1 Jul. 2020
Rather than using global variables you should parameterize the function:
Gopika R
Gopika R am 1 Jul. 2020
@Vashist Hedge: Yes, the arguments are t and x.
Gopika R
Gopika R am 1 Jul. 2020
@Rik: I ran it with the play button, the green one.
Walter Roberson
Walter Roberson am 1 Jul. 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?

Melden Sie sich an, um zu kommentieren.

 Akzeptierte Antwort

Ameer Hamza
Ameer Hamza am 1 Jul. 2020

0 Stimmen

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

Weitere Antworten (0)

Community Treasure Hunt

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

Start Hunting!

Translated by