Matlab or Octave Script for this problem.: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
(New page: % This octave/Matlab script demonstrates the state space model of the J=0.01; b=0.1; K=0.01; R=1; L=0.5; A1=[[0,1,0];[0,-b/J,K/J];[0,-K/L,-R/L]]; B1=[0,0,1/L]'; C1=[0,1,0]; % The outp...)
 
No edit summary
Line 1: Line 1:
% This octave/Matlab script demonstrates the state space model of the
% This octave/Matlab script demonstrates the state space model of the

J=0.01;
J=0.01;
b=0.1;
b=0.1;
K=0.01;
K=0.01;
R=1;
R=1;
L=0.5;
L=0.5;
A1=[[0,1,0];[0,-b/J,K/J];[0,-K/L,-R/L]];
A1=[[0,1,0];[0,-b/J,K/J];[0,-K/L,-R/L]];
B1=[0,0,1/L]';
B1=[0,0,1/L]';
C1=[0,1,0]; % The output is the angular velocity omega.
C1=[0,1,0]; % The output is the angular velocity omega.
sys1=ss(A1,B1,C1); % Note that theta cannot be controlled.
sys1=ss(A1,B1,C1); % Note that theta cannot be controlled.
% There is a warning in octave from this.
% So we use only the theta dot, and i_a as state variables.
% So we use only the theta dot, and i_a as state variables.
A2=[[-b/J, K/J];[-K/L,-R/L]];
B2=[0,1/L]';
A2=[[-b/J, K/J];[-K/L,-R/L]];
C2=[1,0];
B2=[0,1/L]';
C2=[1,0];
sys2=ss(A2,B2,C2);
sys2=ss(A2,B2,C2);
figure(1)
figure(1)
clf;
clf;

Revision as of 12:21, 10 March 2010

% This octave/Matlab script demonstrates the state space model of the 
J=0.01;
b=0.1;
K=0.01;
R=1;
L=0.5;
A1=[[0,1,0];[0,-b/J,K/J];[0,-K/L,-R/L]];
B1=[0,0,1/L]';
C1=[0,1,0]; % The output is the angular velocity omega.
sys1=ss(A1,B1,C1);  % Note that theta cannot be controlled.
% There is a warning in octave from this.
% So we use only the theta dot, and i_a as state variables.
A2=[[-b/J, K/J];[-K/L,-R/L]];
B2=[0,1/L]';
C2=[1,0];
sys2=ss(A2,B2,C2);
figure(1)
clf;
Retrieved from "http:///class-wiki/index.php?title=Matlab_or_Octave_Script_for_this_problem.&oldid=9218"