关于.二阶系统的幅频与相频特性曲线的实验

《关于.二阶系统的幅频与相频特性曲线的实验》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、WORD格式可编辑二阶系统的幅频、相频特性曲线一．实验目的　　在理论学习的基础上，通过本实验熟悉matlab编程，了解二阶系统的频率特性，加深对二阶系统的幅频和相频特性曲线的理解。二．实验原理　　按二阶系统的微分方程（1-1）作拉普拉斯变换，并整理后得到幅频和相频特性的表达式（1-2）和（1-3）。三．实验内容实验内容为选取为0.1、0.2、0.3、0.5、0.7、1.0时对应的取值范围在0.1~10的幅频和相频特性曲线。四．仿真实验1.幅频特性曲线代码。其中z表示，w1、w2、w3、w4、w5、w6表示换算过

2、程的变量，A1、A2、A3、A4、A5、A6表示。开始：clearallclcx=0.1:0.1:0.9;y=1.0:1.0:10;z=[xy];%w/wnw1=1-z.^2;i=0.1;w2=2*i*z;w3=w1.^2;w5=w2.^2;w6=w3+w5;w4=sqrt(w6);A1=1./w4;%专业知识整理分享WORD格式可编辑i=i+0.2;w21=2*i*z;w31=w1.^2;w51=w21.^2;w61=w31+w51;w41=sqrt(w61);A2=1./w41;%i=i+0.2;w22=2

3、*i*z;w32=w1.^2;w52=w22.^2;w62=w32+w52;w42=sqrt(w62);A3=1./w42;%i=i+0.2;w23=2*i*z;w33=w1.^2;w53=w23.^2;w63=w33+w53;w43=sqrt(w63);A4=1./w43;%i=i+0.2;w24=2*i*z;w34=w1.^2;w54=w24.^2;w64=w34+w54;w44=sqrt(w64);A5=1./w44;%i=i+0.2;w25=2*i*z;w35=w1.^2;w55=w25.^2;w65=

4、w35+w55;w45=sqrt(w65);A6=1./w45;N=plot(z,A1,z,A2,z,A3,z,A4,z,A5,z,A6);M=[0.1;0.2;0.5;1.0;2;5;10];set(gca,'xtick',M);专业知识整理分享WORD格式可编辑xlabel('w/wn');%x轴注解ylabel('A(w)');%y轴注解title('fupintexingquxian');gridon2.图形3.相频特性曲线代码。其中z表示，w3、w4、w5、w6和m3、m4、m5、m6表示换算过程的变

5、量，B1、B11、B12、B13、B14、B15表示0.1~0.9时的值，C、C11、C12、C13、C14、C15表示0.9~10时的值。开始：clearclcx=0.1:0.1:0.9;y=1.0:1.0:10;z=[xy];%w/wn%w1=1-z.^2;i=0.1;w3=x.^-1;w4=w3-x;w6=2*i*(w4.^-1);pi=180/3.14159;w5=abs(w6);专业知识整理分享WORD格式可编辑B1=-atan(w5)*pi;m3=y.^-1;m4=m3-y;m55=2*i*(m4.

6、^-1);pi=180/3.14159;m5=abs(m55)A1=-atan(m5)*pi;A2=A1-(A1+90)*2;C=[B1A2];i=0.2;w31=x.^-1;w41=w31-x;w61=2*i*(w41.^-1);pi=180/3.14159;w51=abs(w61);B11=-atan(w51)*pi;m31=y.^-1;m41=m31-y;m61=2*i*(m41.^-1);pi=180/3.14159;m51=abs(m61)A11=-atan(m51)*pi;A21=A11-(A11+

7、90)*2;C11=[B11A21];i=0.3;w32=x.^-1;w42=w32-x;w62=2*i*(w42.^-1);pi=180/3.14159;w52=abs(w62);B12=-atan(w52)*pi;m32=y.^-1;m42=m32-y;m62=2*i*(m42.^-1);pi=180/3.14159;m52=abs(m62)A12=-atan(m52)*pi;A22=A12-(A12+90)*2;C12=[B12A22];专业知识整理分享WORD格式可编辑i=0.5;w33=x.^-1;w

8、43=w33-x;w63=2*i*(w43.^-1);pi=180/3.14159;w53=abs(w63);B13=-atan(w53)*pi;m33=y.^-1;m43=m33-y;m63=2*i*(m43.^-1);pi=180/3.14159;m53=abs(m63)A13=-atan(m53)*pi;A23=A13-(A13+90)*2;C13=[B13A23];i=0.7;