experiment2_428done.m

student s experiment majored in automation

%experiment2_428
%=====================(1)==========
num=5*[1 2.4];den=conv(conv([1 1],[1 1]),conv([1 3 4],[1 0 1]));
[a,b,c,d]=tf2ss(num,den)
sys=ss(a,b,c,d)
%=====================(2)==========
clear;
num=[1 -0.5 2];den=[1 0.4 1];
[z,p,k]=tf2zp(num,den) %[z p k]=tf2zp(num,den)空格也可以分隔输出参数

%
%=================================(3)========
clear;clc
a=[0 1 0;0 0 1;-6 -11 -6];% 张雪平参数录入情况无保留；传递函数的最终形式无显示。
b=[1 0;2 -1;0 2];
c=[1 -1 0;2 1 -1];
d=[0 0;0 0];
[num1,den1]=ss2tf(a,b,c,d,1);
[num2,den2]=ss2tf(a,b,c,d,2);
g11=tf(num1(1,:),den1)
g12=tf(num1(2,:),den1)
g21=tf(num2(1,:),den2)
g22=tf(num2(2,:),den2)
[z11,p,k11]=tf2zp(num1(1,:),den1);
[z12,p,k12]=tf2zp(num1(2,:),den1);
[z21,p,k21]=tf2zp(num2(1,:),den2);
[z22,p,k22]=tf2zp(num2(2,:),den2);
z=[z11,z12,z21,z22]
p=p'
k=[k11,k12,k21,k22]
%=====================(4)==========
[G,H]=c2d(a,b,0.1)% 刘彦秋用c2dm；命令有遗漏；没有每道题的标号。
dsys=ss(G,H,c,d,0.1)

%=====================(5)==========
a=[0 1;-25 -4];
b=[0 0;0 0];        %b=[0 0;0 0];
c=[2 1];
d=[0 1];
[num1,den1]=ss2tf(a,b,c,d,1)
g11=tf(num1,den1)   % g12=tf(num1(2,:),den1)
[num2,den2]=ss2tf(a,b,c,d,2)
g21=tf(num2,den2)   % g22=tf(num2(2,:),den2)

%=====================(6)==========
%saturation function
x1=saturation(1,0.5)
x2=saturation(1,1.5)
%-----------------------画出饱和函数图像---------
clear;figure(1)
i=1;
for u=-3:0.01:3;
y=saturation(u,1);
yy(1,i)=y;
uu(1,i)=u;
i=i+1;
end
plot(uu,yy,'g')
grid;


%------------------------%刘维亮%--------出图速度稍慢些
figure(2)
for u=-3:0.01:3;
y=saturation(u,1);
plot(u,y,'g');hold on;
end
grid;

%=====================(7)==========
dnum=[0.5 -.01];
dden=[1 -0.15 -0.135];
G=tf(dnum,dden,0.5)
[z,p,k]=tf2zp(dnum,dden);
dsys=zpk(z,p,k,0.5)