MATELAB实验仿真.docx
- 文档编号:29920985
- 上传时间:2023-08-03
- 格式:DOCX
- 页数:23
- 大小:318.25KB
MATELAB实验仿真.docx
《MATELAB实验仿真.docx》由会员分享,可在线阅读,更多相关《MATELAB实验仿真.docx(23页珍藏版)》请在冰豆网上搜索。
MATELAB实验仿真
实验2
1.
num=[52530];den=[16108];sys=tf(num,den);t=0:
0.01:
12;
figure
(1)impulse(sys,t);grid
xlabel('t');ylabel('c(t)');title('impulsereponse');
figure
(2)
step(sys,t);gridxlabel('t');ylabel('c(t)');title('stepreponse');
figure(3)
u=0*t;lsim(sys,u,t,0);gridxlabel('t');ylabel('c(t)');title('zeroreponse');
2.
forwn=2:
2:
12
num=[0,wn*wn];den=[12*0.7*wnwn*wn];sys=tf(num,den);t=0:
0.01:
4;
ifwn==2
y1=(sys);
elseifwn==4
y2=(sys);
elseifwn==6
y3=(sys);
elseifwn==8
y4=(sys);
elseifwn==10
y5=(sys);
elseifwn==12
y6=(sys);
end
end
figure
(1)
step(y1,t,y2,t,y3,t,y4,t,y5,t,y6,t)
legend('wn=2','wn=4','wn=6','wn=8','wn=10','wn=12')
grid;xlabel('t');ylabel('c(t)');title('stepreponse');
forwn=2:
2:
12
num=[0,wn*wn];den=[12*0.7*wnwn*wn];sys=tf(num,den);t=0:
0.01:
4;
ifwn==2
y1=(sys);
elseifwn==4
y2=(sys);
elseifwn==6
y3=(sys);
elseifwn==8
y4=(sys);
elseifwn==10
y5=(sys);
elseifwn==12
y6=(sys);
end
end
figure
(1)
u=0*t;lsim(y1,y2,y3,y4,y5,y6,u,t,0);grid
xlabel('t');ylabel('c(t)');title('zeroreponse');
forwn=2:
2:
12
num=[0,wn*wn];den=[12*0.7*wnwn*wn];sys=tf(num,den);
t=0:
0.01:
2.5;
ifwn==2
y1=(sys);
elseifwn==4
y2=(sys);
elseifwn==6
y3=(sys);
elseifwn==8
y4=(sys);
elseifwn==10
y5=(sys);
elseifwn==12
y6=(sys);
end
end
figure
(1)
impulse(y1,y2,y3,y4,y5,y6,t);gridxlabel('t');ylabel('c(t)');title('impulsereponse');
3.
fore=0.2:
0.1:
2.0
wn=6;num=[0,wn*wn];den=[12*e*wnwn*wn];sys=tf(num,den);
t=0:
0.01:
20;
ife==0.2
y1=(sys);
elseife==0.4
y2=(sys);
elseife==0.6
y3=(sys);
elseife==0.8
y4=(sys);
elseife==1.0
y5=(sys);
elseife==1.5
y6=(sys);
elseife==2.0
y7=(sys);
end
end
figure
(1)
step(y1,t,y2,t,y3,t,y4,t,y5,t,y6,y7,t);grid
xlabel('t');ylabel('c(t)');title('stepreponse');
fore=0.2:
0.1:
2.0
wn=6;num=[0,wn*wn];den=[12*e*wnwn*wn];sys=tf(num,den);
t=0:
0.01:
20;
ife==0.2
y1=(sys);
elseife==0.4
y2=(sys);
elseife==0.6
y3=(sys);
elseife==0.8
y4=(sys);
elseife==1.0
y5=(sys);
elseife==1.5
y6=(sys);
elseife==2.0
y7=(sys);
end
end
figure
(1)
u=0*t;lsim(y1,y2,y3,y4,y5,y6,y7,u,t,0);grid
xlabel('t');ylabel('c(t)');title('zeroreponse');
fore=0.2:
0.1:
2.0
wn=6;num=[0,wn*wn];den=[12*e*wnwn*wn];sys=tf(num,den);
t=0:
0.01:
6;
ife==0.2
y1=(sys);
elseife==0.4
y2=(sys);
elseife==0.6
y3=(sys);
elseife==0.8
y4=(sys);
elseife==1.0
y5=(sys);
elseife==1.5
y6=(sys);
elseife==2.0
y7=(sys);
end
end
figure
(1)
impulse(y1,y2,y3,y4,y5,y6,y7,t);grid
xlabel('t');ylabel('c(t)');title('impulsereponse');
实验3
1.
forwn=2:
2:
12
num=[0,wn*wn];den=[12*0.7*wnwn*wn];sys=tf(num,den);
ifwn==2
y1=(sys);
elseifwn==4
y2=(sys);
elseifwn==6
y3=(sys);
elseifwn==8
y4=(sys);
elseifwn==10
y5=(sys);
elseifwn==12
y6=(sys);
end
end
bode(y1,y2,y3,y4,y5,y6)
grid
forwn=2:
2:
12
num=[0,wn*wn];den=[12*0.7*wnwn*wn];sys=tf(num,den);
ifwn==2
y1=(sys);
elseifwn==4
y2=(sys);
elseifwn==6
y3=(sys);
elseifwn==8
y4=(sys);
elseifwn==10y1-y6
y5=(sys);
elseifwn==12
y6=(sys);
end
end
nyquist(y1,'k--',y2,'k*',y3,'k:
',y4,'k+',y5,'kp',y6,'k');grid
y1
y6
2
fore=0.2:
0.1:
2.0
wn=6;num=[0,wn*wn];den=[12*e*wnwn*wn];sys=tf(num,den);
t=0:
0.01:
8;
ife==0.2
y1=(sys);elseife==0.4
y2=(sys);elseife==0.6
y3=(sys);elseife==0.8
y4=(sys);elseife==1.0
y5=(sys);elseife==1.5
y6=(sys);elseife==2.0
y7=(sys);
end
end
bode(y1,'k--',y2,'k:
',y4,'k+',y5,'kp',y6,'kh',y7,'k');grid
fore=0.2:
0.1:
2.0
wn=6;
num=[0,wn*wn];
den=[12*e*wnwn*wn];
sys=tf(num,den);
t=0:
0.01:
8;
ife==0.2
y1=(sys);
elseife==0.4
y2=(sys);
elseife==0.6
y3=(sys);
elseife==0.8
y4=(sys);
elseife==1.0
y5=(sys);
elseife==1.5
y6=(sys);
elseife==2.0
y7=(sys);
end
end
nyquist(y1,'k--',y2,'k:
',y4,'k+',y5,'kp',y6,'kh',y7,'k');grid
fore=0.2:
0.1:
2.0
wn=6;
num=[0,wn*wn];
den=[12*e*wnwn*wn];
sys=tf(num,den);
ife==0.6
y3=(sys);
end
end
nyquist(y3)
grid
实验4
1.
num=[0,1];
den=conv(conv([10],[11]),[12]);
sys=tf(num,den);
rlocus(sys)
2.
num=[0,1];
den=conv(conv([100],[12]),[15]);
sys=tf(num,den);
rlocus(sys)
3.
num=[0,1];den=conv(conv([10],[14]),[1420]);
sys=tf(num,den);
rlocus(sys)
实验5
1.
num=[5,10];
den=conv(conv([10],[11]),[0.511]);
sys=tf(num,den);
margin(sys)
2.
num=[1,2];
den=conv([149],[149]);
sys=tf(num,den);
rlocus(sys)
3.
(1)
num=[02.7];
den=[1540];
sys=tf(num,den);
margin(sys)
(2)
num=[02.7];
den=[15-40];
sys=tf(num,den);
margin(sys)
实验7
1.
k=25;
G0=zpk([],[0-1-4],25);
[h0,r,wx,wc]=margin(G0)
wm=2.5;L=bode(G0,wm);
Lwc=20*log10(L);a=10^(-0.1*Lwc);
T=1/(wm*sqrt(a));
phi=asin((a-1)/(a+1));
Gc=(1/a)*tf([a*T1],[T,1]);
Gc=a*Gc;G=Gc*G0;bode(G,'K',G0,'K--');
grid;
[h,r,wx,wc]=margin(G)
h0=0.8000
r=-5.0230
wx=2.0000
wc=2.2318
h=1.1522
r=3.3696
wx=2.6902
wc=2.4998
>>Gc
Transferfunction:
0.508s+1
------------
0.3149s+1
2.
G0=zpk([],[0-10-2],k);
[h0,r,wx,wc]=margin(G0)
wm=2;L=bode(G0,wm);Lwc=20*log10(L);
a=10^(-0.1*Lwc);
T=1/(wm*sqrt(a));
phi=asin((a-1)/(a+1));
Gc=(1/a)*tf([a*T1],[T,1]);
Gc=a*Gc;G=Gc*G0;bode(G,'K',G0,'K--');
grid;
[h,r,wx,wc]=margin(G)
h0=12
r=60.4231
wx=4.4721
wc=0.9070
h=7.6292
r=85.4485
wx=8.8283
wc=2
>>Gc
Transferfunction:
1.442s+1
------------
0.1733s+1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATELAB 实验 仿真