📄 linear system ii.m
字号:
%-- 11/11/08 5:13 PM --%
a=[1 0 0;1 1 1;20 2 0]
a=[1 0 0;1 1 1;20 2 0];
a=[1 0 0;1 1 1;20 2 0]
a=[1 0 0;1 1 1;20 2 0];
b=[]
simulink
%-- 11/11/08 9:09 PM --%
T=[72 4 -174 24;-6 79 121 12;-6 -11 20 1;0 -12 1 1];
A=inv[T]
T=[72 4 -174 24;-6 79 121 12;-6 -11 20 1;0 -12 1 1];
A=inv(T)
A =
-0.0332 0.0463 -0.6118 0.8543
0.0043 -0.0014 0.0536 -0.1406
-0.0107 0.0147 -0.1434 0.2245
0.0629 -0.0321 0.7869 -0.9115
A =
-0.0332 0.0463 -0.6118 0.8543
0.0043 -0.0014 0.0536 -0.1406
-0.0107 0.0147 -0.1434 0.2245
0.0629 -0.0321 0.7869 -0.9115
A =
-0.0332 0.0463 -0.6118 0.8543
0.0043 -0.0014 0.0536 -0.1406
-0.0107 0.0147 -0.1434 0.2245
0.0629 -0.0321 0.7869 -0.9115
A =
-0.0332 0.0463 -0.6118 0.8543
0.0043 -0.0014 0.0536 -0.1406
-0.0107 0.0147 -0.1434 0.2245
0.0629 -0.0321 0.7869 -0.9115
T=[72 4 -174 24;-6 79 121 12;-6 -11 20 1;0 -12 1 1];
>> A=inv(T);
T=[72 4 -174 24;-6 79 121 12;-6 -11 20 1;0 -12 1 1];
>> A=inv(T);
T=[72 4 -174 24;-6 79 121 12;-6 -11 20 1;0 -12 1 1];
>> A=inv(T)
%-- 11/12/08 1:15 PM --%
control
\control
\tool\box
\tool\box\control
%-- 11/12/08 2:45 PM --%
a=[0 1 0 0 ;3 0 0 2;0 0 0 1;0 -2 0 0];
b=a*a
a=[0 1 0 0 ;3 0 0 2;0 0 0 1;0 -2 0 0];
b=a*a*a
a=[0 1 0 0 ;3 0 0 2;0 0 0 1;0 -2 0 0];
b=[0;0;0;1];
p=[-2;-1+i;-1-i];
K=place(a,b,p)
a=[0 1 0 0 ;3 0 0 2;0 0 0 1;0 -2 0 0];
b=[0;0;0;1];
p=[-2;-1+i;-1-i;-3];
K=place(a,b,p)
eig(a-b*K)
a=[0 1 0 0 ;3 0 0 2;0 0 0 1;0 -2 0 0];
c=[0 0 1 0];
p=[-2;-2;-1;-1];
K=place(a,b,p)
a=[0 1 0 0 ;3 0 0 2;0 0 0 1;0 -2 0 0];
c=[0 0 1 0];
p=[-2 -2 -1 -1];
K=place(a,b,p)
a=[0 3 0 0;1 0 0 -2;0 0 0 0;0 2 1 0];
c=[0;0;1;0];
p=[-2 -2 -1 -1];
K=place(a,c,p)
a=[0 3 0 0;1 0 0 -2;0 0 0 0;0 2 1 0];
c=[0;0;1;0];
p=[-2;-2;-1;-1];
K=place(a,c,p)
%-- 11/12/08 3:32 PM --%
a=[0 3 0 0;1 0 0 -2;0 0 0 0;0 2 1 0];
>> c=[0;0;1;0];
>> p=[-2;-2;-1;-1];
>> K=place(a,c,p)
a=[0 3 0 0;1 0 0 -2;0 0 0 0;0 2 1 0];
b=[0;0;1;0];
p=[-2;-2;-1;-1];
K=place(a,b,p)
a=[0 3 0 0;1 0 0 -2;0 0 0 0;0 2 1 0];
b=[0;0;1;0];
p=[-2;-2;-1;-1];
K=place(a,b,p)
a=[0 3 0 0;1 0 0 -2;0 0 0 0;0 2 1 0];
>> b=[0;0;1;0];
>> p=[-2 -2 -1 -1];
>> K=place(a,b,p)
%-- 11/12/08 3:44 PM --%
a=[0 3 0 0;1 0 0 -2;0 0 0 0;1 0 -4 0];
b=[0;0;1;0];
p=[-2+i;-2-i;-1;-1];
K=place(a,b,p)
a=[0 3 0 0;1 0 0 -2;0 0 0 0;1 0 -4 0];
b=[0;0;1;0];
p=[-2+i;-2-i;-1;-1]
K=place(a,b,p)
%-- 11/12/08 6:05 PM --%
a=[0 1 0 0;3 0 0 2;0 0 0 1;-21.5 -12.5 4 -7];
b=[0;0;0;1];
c=[-21.5 -10.5 4 -7];
d=[1];
[N1,d1]=ss2tf(a,b,c,d)
a=[0 1 0 0;3 0 0 2;0 0 0 1;-21.5 -12.5 4 -7];
b=[0;0;0;1];
c=[0 0 1 0];
[N1,d1]=ss2tf(a,b,c)
a=[0 1 0 0;3 0 0 2;0 0 0 1;-21.5 -12.5 4 -7];b=[0;0;0;1];
b=[0;0;0;1];
c=[0 0 1 0];
d=[0];
[N1,d1]=ss2tf(a,b,c,d)
a=[0 1 0 0;3 0 0 2;0 0 0 1;-21.5 -12.5 4 -7];b=[0;0;0;1];
b=[10/3;1;-4;-6];
%-- 11/12/08 6:34 PM --%
a=[0 1 0 0;3 0 0 2;0 0 0 1;-21.5 -12.5 4 -7];b=[0;0;0;1];
b=[-10/3;-1;4;6];
c=[0 0 1 0];
d=[1];
[N1,d1]=ss2tf(a,b,c,d)
a=[0 1 0 0;3 0 0 2;0 0 0 1;-21.5 -12.5 4 -7];b=[0;0;0;1];
b=[-10/3;-1;4;6];
c=[-21.5 -10.5 4 -7];
d=[0];
[N1,d1]=ss2tf(a,b,c,d)
%-- 11/13/08 8:40 PM --%
function xdot=zhao_f(time,x)
x11=x(1);
x12=x(2);
x21=x(3);
x22=x(4);
x31=x(5);
x32=x(6);
x41=x(7);
x42=x(8);
x51=x(9);
x52=x(10);
x61=x(11)
x62=x(12)
x71=x(13);
x72=x(14);
x81=x(15);
x82=x(16);
x91=x(17);
x92=x(18);
dx11=(x21-x11)+(x31-x11)+(x41-x11);
dx12=(x22-x12)+(x32-x12)+(x42-x12);
dx21=(x11-x21)+(x31-x21)+(x41-x21);
dx22=(x12-x22)+(x32-x22)+(x42-x22);
dx31=(x11-x31)+(x21-x31)+(x41-x31)+(x51-x31);
dx32=(x12-x32)+(x22-x32)+(x42-x32)+(x52-x32);
dx41=(x11-x41)+(x21-x41)+(x31-x41);
dx42=(x12-x42)+(x22-x42)+(x32-x42);
dx51=(x31-x51)+(x61-x51);
dx52=(x32-x52)+(x62-x52);
dx61=(x51-x61)+(x71-x61)+(x81-x61)+(x91-x61);
dx62=(x52-x62)+(x72-x62)+(x82-x62)+(x92-x62);
dx71=(x61-x71)+(x81-x71)+(x91-x71);
dx72=(x62-x72)+(x82-x72)+(x92-x72);
dx81=(x61-x81)+(x71-x81)+(x91-x81);
dx82=(x62-x82)+(x72-x82)+(x92-x82);
dx91=(x61-x91)+(x71-x91)+(x81-x91);
dx92=(x62-x92)+(x72-x92)+(x82-x92);
xdot=[dx11 dx12 dx21 dx22 dx31 dx32 dx41 dx42 dx51 dx52 dx61 dx62 dx71 dx72 dx81 dx82 dx91 dx92]';
clear;
clc;
x0=[-8 0 -6 2 -4 0 -6 -2 0 0 4 0 6 2 8 0 6 -2]';
dt=0.01;
tfinal=10;
[tt,yy]=ode45('zhao_f',[0:dt:tfinal],x0);
%
% figure;
% subplot(1,2,1);
% axis([-8 8 -8 8]);
% axis square;
% hold on;
% plot(-8,0,'r*',-6,-2,'r*',-6,2,'r*',-4,0,'r*');
% plot(0,0,'r*',4,0,'r*',6,2,'r*',6,-2,'r*',8,0,'r*');
% plot([-8,-6,-4,-6,-8],[0,2,0,-2,0],'r:');
% plot([-8,8],[0,0],'r:');
% plot([4,6,8,6,4],[0,2,0,-2,0],'r:');
% plot([-6,-6],[2,-2],'r:');
% plot([6,6],[2,-2],'r:');
subplot(1,2,1);
axis square;
hold on;
for i=0:8
plot(tt,yy(:,2*i+1),'r');
plot(tt,yy(:,2*i+2),'b');
hold on;
end
clear;
clc;
x0=[-8 0 -6 2 -4 0 -6 -2 0 0 4 0 6 2 8 0 6 -2]';
dt=0.01;
tfinal=10;
[tt,yy]=ode45('zhao_f',[0:dt:tfinal],x0);
%
% figure;
% subplot(1,2,1);
% axis([-8 8 -8 8]);
% axis square;
% hold on;
% plot(-8,0,'r*',-6,-2,'r*',-6,2,'r*',-4,0,'r*');
% plot(0,0,'r*',4,0,'r*',6,2,'r*',6,-2,'r*',8,0,'r*');
% plot([-8,-6,-4,-6,-8],[0,2,0,-2,0],'r:');
% plot([-8,8],[0,0],'r:');
% plot([4,6,8,6,4],[0,2,0,-2,0],'r:');
% plot([-6,-6],[2,-2],'r:');
% plot([6,6],[2,-2],'r:');
subplot(1,2,1);
axis square;
hold on;
for i=0:8
plot(tt,yy(:,2*i+1),'r');
plot(tt,yy(:,2*i+2),'b');
hold on;
end
clear;
clc;
x0=[-8 0 -6 2 -4 0 -6 -2 0 0 4 0 6 2 8 0 6 -2]';
dt=0.01;
tfinal=10;
[tt,yy]=ode45('zhao_f',[0:dt:tfinal],x0);
%
% figure;
% subplot(1,2,1);
% axis([-8 8 -8 8]);
% axis square;
% hold on;
% plot(-8,0,'r*',-6,-2,'r*',-6,2,'r*',-4,0,'r*');
% plot(0,0,'r*',4,0,'r*',6,2,'r*',6,-2,'r*',8,0,'r*');
% plot([-8,-6,-4,-6,-8],[0,2,0,-2,0],'r:');
% plot([-8,8],[0,0],'r:');
% plot([4,6,8,6,4],[0,2,0,-2,0],'r:');
% plot([-6,-6],[2,-2],'r:');
% plot([6,6],[2,-2],'r:');
subplot(1,2,1);
axis square;
hold on;
for i=0:8
plot(tt,yy(:,2*i+1),'r');
plot(tt,yy(:,2*i+2),'b');
hold on;
end
kk
1+2
clear
clear;
clear; clr
clear;
clc;
x0=[-8 0 -6 2 -4 0 -6 -2 0 0 4 0 6 2 8 0 6 -2]';
dt=0.01;
tfinal=10;
[tt,yy]=ode45('zhao_f',[0:dt:tfinal],x0);
%
% figure;
% subplot(1,2,1);
% axis([-8 8 -8 8]);
% axis square;
% hold on;
% plot(-8,0,'r*',-6,-2,'r*',-6,2,'r*',-4,0,'r*');
% plot(0,0,'r*',4,0,'r*',6,2,'r*',6,-2,'r*',8,0,'r*');
% plot([-8,-6,-4,-6,-8],[0,2,0,-2,0],'r:');
% plot([-8,8],[0,0],'r:');
% plot([4,6,8,6,4],[0,2,0,-2,0],'r:');
% plot([-6,-6],[2,-2],'r:');
% plot([6,6],[2,-2],'r:');
subplot(1,2,1);
axis square;
hold on;
for i=0:8
plot(tt,yy(:,2*i+1),'r');
plot(tt,yy(:,2*i+2),'b');
hold on;
end
clc
2+5
2+4;
a
a=2-=3
clc
%-- 11/18/08 10:11 AM --%
R1=A(u+w)[b1*k*c*(b+x)+b2*c(k*m*v+k1*p*a)+b3*(kmvx+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))]+uw(r1*c*(k*m*v+k1*p*a)+r2(kmvx+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x)));
R1=A(u+w)[b1*k*c*(b+x)+b2*c(k*m*v+k1*p*a)+b3*(kmvx+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))]+uw(r1*c*(k*m*v+k1*p*a)+r2(kmvx+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x)))
R1=A(u+w)[b1*k*c*(b+x)+b2*c(k*m*v+k1*p*a)+b3*(kmvx+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))]+uw(r1*c*(k*m*v+k1*p*a)+r2(kmvx+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x)));
R2=u*(u+w)*a*(b+x)c*(u+k1+k);
R1=A*(u+w)*[b1*k*c*(b+x)+b2*c(k*m*v+k1*p*a)+b3*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))]+u*w[r1*c*(k*m*v+k1*p*a)+r2(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))];
u=0.02;w=0.001;A=10000;b1=0.00001;b2=0.000005;b3=0.00001;k=0.09;k1=0.01;
p=0.85;v=0.3;m=0.85;x=0.15;r1=0.75;r2=0.1;a=0.47;b=0.78;c=0.17;
R1=A(u+w)[b1*k*c*(b+x)+b2*c(k*m*v+k1*p*a)+b3*(kmvx+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))]+uw(r1*c*(k*m*v+k1*p*a)+r2(kmvx+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x)));R2=u*(u+w)*a*(b+x)c*(u+k1+k);
R1=A*(u+w)*[b1*k*c*(b+x)+b2*c(k*m*v+k1*p*a)+b3*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))]+u*w[r1*c*(k*m*v+k1*p*a)+r2(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))];R2=u*(u+w)*a*(b+x)c*(u+k1+k);
u=0.02;w=0.001;A=10000;b1=0.00001;b2=0.000005;b3=0.00001;k=0.09;k1=0.01;
p=0.85;v=0.3;m=0.85;x=0.15;r1=0.75;r2=0.1;a=0.47;b=0.78;c=0.17;
R2=u*(u+w)*a*(b+x)c*(u+k1+k)
u*(u+w)*a*(b+x)c*(u+k1+k)
u*(u+w)*a*(b+x)c*(u+k1+k);
%-- 11/18/08 11:02 AM --%
u=0.02;w=0.001;A=10000;b1=0.00001;b2=0.000005;b3=0.00001;k=0.09;k1=0.01;
p=0.85;v=0.3;m=0.85;x=0.15;r1=0.75;r2=0.1;a=0.47;b=0.78;c=0.17;
R2=u*(u+w)*a*(b+x)c*(u+k1+k);
R1=A*(u+w)*[b1*k*c*(b+x)+b2*c(k*m*v+k1*p*a)+b3*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))]+u*w[r1*c*(k*m*v+k1*p*a)+r2(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))];R2=u*(u+w)*a*(b+x)c*(u+k1+k);u=0.02;w=0.001;A=10000;b1=0.00001;b2=0.000005;b3=0.00001;k=0.09;k1=0.01;
u=0.02;w=0.001;A=10000;b1=0.00001;b2=0.000005;b3=0.00001;k=0.09;k1=0.01;
u
R2=u*(u+w)*a*(b+x)c*(u+k1+k);
R2=u*(u+w)*a*(b+x)*c*(u+k1+k);
R2
R1=A*(u+w)*[b1*k*c*(b+x)+b2*c(k*m*v+k1*p*a)+b3*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))]+u*w[r1*c*(k*m*v+k1*p*a)+r2(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))];
R1=A*(u+w)*[b1*k*c*(b+x)+b2*c*(k*m*v+k1*p*a)+b3*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))]+u*w[r1*c*(k*m*v+k1*p*a)+r2(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))];
R1=A*(u+w)*[b1*k*c*(b+x)+b2*c*(k*m*v+k1*p*a)+b3*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))]+u*w*[r1*c*(k*m*v+k1*p*a)+r2*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))];
R1=A*(u+w)*[b1*k*c*(b+x)+b2*c*(k*m*v+k1*p*a)+b3*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))]+u*w*[r1*c*(k*m*v+k1*p*a)+r2*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))]
R1/R2
A=1000
R1=A*(u+w)*[b1*k*c*(b+x)+b2*c(k*m*v+k1*p*a)+b3*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v(1-m)*(b+x))]+u*w[r1*c*(k*m*v+k1*p*a)+r2(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))];
R1=A*(u+w)*[b1*k*c*(b+x)+b2*c*(k*m*v+k1*p*a)+b3*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))]+u*w*[r1*c*(k*m*v+k1*p*a)+r2*(k*m*v*x+k1*(1-p)*a*(b+x)+k1*p*a*x+k*v*(1-m)*(b+x))]
R2
R2=u*(u+w)*a*(b+x)*c*(u+k1+k)
%-- 11/18/08 11:50 AM --%
clear
text('ll')
text('13')
gtext('13')
%-- 11/19/08 10:12 AM --%
simulink
%-- 11/20/08 6:22 PM --%
log(8/1920)/log(2)
b1=10-log(8/1920)/log(2)
b2=10-log(12/1920)/log(2)
b3=10-log(20/1920)/log(2)
1920^(1/3)
t=1920^(1/3)
b1=10-log(8/t)/log(2)
b2=10-log(12/t)/log(2)
b3=10-log(20/t)/log(2)
8*12
96*20
196*12
%-- 11/20/08 9:09 PM --%
hold on
fig_count
gtext('p=0.1'),
gtext('p=0.8'),
figure(2)
hold off
gtext('p=0.1'),
gtext('p=0.8'),
gtext('S'),
hold on
gtext('S'),
%-- 11/20/08 10:20 PM --%
gtext('p=0.1'),
gtext('p=0.8'),
hold on
gtext('p=0.1'),gtext('p=0.2'),
gtext('sigama=0.1'),gtext('sigama=0.8'),
gtext('sigama=0.6'),
%-- 11/21/08 12:43 PM --%
gtext('v=0.1'),gtext('v=0.8')
%-- 12/16/08 11:02 AM --%
A=
A=[-0.08 -0.03 0 -0.157;-0.73 -0.377 1 0;0 -8.65 -0.5 0;0 0 1 0 ]
eig(A)
B=1 0 0 0; 0 1 0 0 ]'
B=1 0 0 0; 0 1 0 0 ]`
B=1 0;0 1; 0 0; 0 0 ]
B=[1 0;0 1; 0 0; 0 0 ]
rank(ctrl (A,B)\)
rank(ctrl (A,B))
rank(ctrl(A,B))
rank(ctrb(A,B))
C=[1 0 0 0;0 0 1 2]
ran(obsr(A,C))
ran(obsv(A,C))
rank(obsv(A,C))
M=ctrb(A,B)
G90
[NUM,DEN] = SS2TF(A,B,C,D,2
D=[0 0; 0 0]
[NUM,DEN] = SS2TF(A,B,C,D,2)
A=[-0.08 -0.03 0 -0.157;-0.73 -0.377 1 0;0 -8.65 -0.5 0;0 0 1 0 ]
%-- 12/17/08 6:06 PM --%
A=[-0.08 -0.03 0 -0.157;-0.73 -0.377 1 0;0 -8.65 -0.5 0;0 0 1 0 ]
[NUM,DEN] = SS2TF(A,B,C,D,2
[NUM,DEN] = SS2TF(A,B,C,D,2)
A=[-0.08 -0.03 0 -0.157;-0.73 -0.377 1 0;0 -8.65 -0.5 0;0 0 1 0 ]
[NUM,DEN] = SS2TF(A,B,C,D,2
[NUM,DEN] = SS2TF(A,B,C,D,2)
[NUM,DEN] = SS2TF(A,B,C,D,2
[NUM,DEN] = SS2TF(A,B,C,D,2)
[NUM,DEN] = SS2TF(A,B,C,D,2
[NUM,DEN] = SS2TF(A,B,C,D,2)
%-- 12/19/08 11:05 AM --%
[NUM,DEN] = SS2TF(A,B,C,D,2)
A =
A=[-0.08 -0.03 0 -0.157;-0.73 -0.377 1 0;0 -8.65 -0.5 0;0 0 1 0 ]
B=[1 0;0 1; 0 0; 0 0 ]
C=[1 0 0 0;0 0 1 2]
M=ctrb(A,B)
rank(ctrl (A,B)\)
rank(ctrb(A,B))
rank(obsv(A,C))
D=[0 0; 0 0]
[NUM,DEN] = SS2TF(A,B,C,D,2)
eig(A)
p11=[ 1.0000 0.9570 8.8868 0.6961 0.9914
]
roots(p11)
p123=[ -8.6500 -17.9920 -1.3840]
roots(p123)
C2-[0 0 1 2]
C2=[0 0 1 2]
B2=[ 0;1; 0; 0 ]
G22=
[NUM,DEN] = SS2TF(A,B,C,D,2)
D22=0
[NUM,DEN] = SS2TF(A,B2,C2,D,2)
C1
C2
[NUM,DEN] = SS2TF(A,B2,C2,D)
[NUM,DEN] = SS2TF(A,B2,C2,D,1)
[NUM,DEN] = SS2TF(A,B2,C2,D22,1)
aig
eig(A)
[NUM,DEN] = SS2TF(A,B2,C2,I,1)
I=0
[NUM,DEN] = SS2TF(A,B2,C2,I,1)
[NUM,DEN] = SS2TF(A,B2,C2,D22,1)
I=1
[NUM,DEN] = SS2TF(A,B2,C2,I,1)
H=[0;0;0;0]
[NUM,DEN] = SS2TF(A,H,C2,I,1)
[NUM,DEN] = SS2TF(A,-B2,F,I,1)
pp=[-0.4444 + 2.9180i -0.4444 - 2.9180i -1 -1]
KK=place(A,B,pp)
KK=place(A,B2,pp)
KK=place(A,B,pp)
B2
B4=[0;0;0;1]
KK=place(A,B4,pp)
pp=[-0.4444 + 2.9180i -0.4444 -1 -1]
KK=place(A,B4,pp)
pp=[-0.4444 + 2.9180i -0.4444 - 2.9180i -0.0341 + 0.3356i -1]
KK=place(A,B2,pp)
-0.0341 + 0.3356i-0.0341 + 0.3356i
pp=[-0.4444 + 2.9180i -0.4444 - 2.9180i -0.0341 + 0.3356i -0.0341 + 0.3356i]
KK=place(A,B2,pp)
A
A1=[]
A1=[1 1 2 0;2 0 2 1;1;2 0 0; 1 0 0 2]
A1=[1 1 2 0;2 0 2 1;1 2 0 0; 1 0 0 2]
eig(A1)
pp=[-1,-2,-1,-1]
KK=place(A1,B2,pp)
pp=[ -1.7835
-0.7025
-1,-1]
-1.7835
-0.7025
pp=[ -1.7835
-0.7025
,-1,-1]
pp=[ -1.7835 -0.7025
- 1 -1]
pp=[ -1.7835 -0.7025 -1 -1]
KK=place(A1,B2,pp)
A
-0.0800 -0.0300 0 -0.1570
-0.7300 -0.3770 1.0000 0
0 -8.6500 -0.5000 0
0 0 1.0000 0
B2
rank(ctrb(A,B2))
p5=[-1 -1 -1+i -1-i]
KK=place(A,B2,p5)
p5=[-1 ;-1; -1+i; -1-i]
KK=place(A,B2,p5)
B2
A
KK=place(A,B2,p5)
pp=[-0.4444 + 2.9180i; -0.4444 - 2.9180i ; -0.0341 + 0.3356i -0.0341 -0.3356i]
pp=[-0.4444 + 2.9180i; -0.4444 - 2.9180i ; -0.0341 + 0.3356i ;-0.0341 -0.3356i]
0
B
A=[-0.08 -0.03 0 -0.157;-0.73 -0.377 1 0;0 -8.65 -0.5 0;0 0 1 0 ];
p22=[-2;-3;-1+i;-1-i]
k=place(A,B2,p22)
p22=[-1,-1,-2,-3]
p22=[-1;-1;-2;-3]
k=place(A,B2,p22)
p22=[-1;-2;-1+i;-1-i];
k=place(A,B2,p22)
p22=[-1;-2;-1+i;-1-i];
p22=[-1;-2;-1+i;-1-i]
k=place(A,B2,p22)
p22=[-1;-2;-8;-5]
k=place(A,B2,p22)
p22=[-1;-2;-1+i;-1-i]
k=place(A,B2,p22)
B2*k
A+B_{2}*F
A+B2*F
A-B2*k
[NUM,DEN] = SS2TF(A,B2,,I,1)
f=place(A,B2,p22)
[NUM,DEN] = SS2TF(A,B2,f,I,1)
AF=A-B2*k
[NUM,DEN] = SS2TF(AF,B2,f,I,1)
I
C2
C2=[0 0 1 2]
[NUM,DEN] = SS2TF(AF,B2,c2,D,1)
D
[NUM,DEN] = SS2TF(AF,B2,C2,D,1)
[NUM,DEN] = SS2TF(AF,B2,C2,0,1)
[NUM,DEN] = SS2TF(A,B2,C2,0,1)
H=[0;0;0;0]
[NUM,DEN] = SS2TF(A,H,C2,1)
[NUM,DEN] = SS2TF(A,-B2,f,1)
G22
[NUM,DEN] = SS2TF(A,B2,C2,1)
[NUM,DEN] = SS2TF(A,B2,C2,0)
[NUM,DEN] = SS2TF(A,-H,f,0)
[NUM,DEN] = SS2TF(A,-B2,f,1)
[NUM,DEN] = SS2TF(AF,B2,f,1)
%-- 12/20/08 8:24 PM --%
A
A=[-0.08 -0.03 0 -0.157;-0.73 -0.377 1 0;0 -8.65 -0.5 0;0 0 1 0 ]
B=[1 0;0 1; 0 0; 0 0 ]
C=[1 0 0 0;0 0 1 2]
[NUM,DEN] = SS2TF(A,B,C,D,2)
D=[0 0; 0 0]
[NUM,DEN] = SS2TF(A,B,C,D,2)
[NUM,DEN] = SS2TF(A,B,C,D,1)⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -