西工大电气系统仿真第三四次作业精要.docxVIP

电气系统仿真第三第四次作业专业：电气工程及其自动化学生姓名：郭振超学生学号:2014301801学生班号三次1.clc;clear;num=[1 0];den=[1 2 10];h=tf(num,den) %h(s)hz=tf(num,den,0.1) %h(z)num1=[1 -1];den1=[1 1];num2=[1 2];den2=[1 4 5];sys1=tf(num1,den1);sys2=tf(num2,den2);H1=[sys1;sys2]H2=[sys1 sys2]2.clc;clear;z1=[-12];p1=[-1 -4 -5];k1=3;z2=[0];p2=[2 -4];k2=-2;z3=[];p3=[2 -4];k3=-2;h1=zpk(z1,p1,k1)h2=zpk(z2,p2,k2)h3=zpk(z3,p3,k3)%下面转化为传递函数型tf_h1=tf(h1)tf_h2=tf(h2)tf_h3=tf(h3)%下面转化为状态空间模型ss_h1=ss(h1)ss_h2=ss(h2)ss_h3=ss(h3)3.clc;clear;%G1(s)编写G1=tf([1],[1 10]);%G2(s)编写G2=tf([1],[1 1]);%G3(s)编写G3=tf([1 0 1],[1 4 4]);%G4(s)编写G4=tf([1 1],[1 6]);%H1(s)编写H1=tf([1 1 ],[1 2]);%H2(s)编写H2=tf([2],[1]);%H3(s)编写H3=tf([1],[1]);%三种方式表达式不一样但是绘图结果是一样的这是因为一个系统的传递函数表达形式不是唯一的%比较点前移Hx1=H1/G1/G2;%此时G2 G3 H2 构成的闭环传递函数设为 G23G23=feedback(G2*G3,H2);%G1 G23 G4 H3 Hx1算出G（s）G=feedback(G1*G23*G4,Hx1-H3,+1)figure;step(G) G34=feedback(G3*G4,H1,+1); G23=feedback(G2*G34,H2/G4); G=feedback(G1*G23,H3)figure;step(G) Hx1=H1*G4/G2; G23=feedback(G2*G3,Hx1-H2,+1); G=feedback(G1*G23*G4,H3)figure;step(G)4.1clc;clear;a = [-0.5572 -0.7814;0.7814 0];b = [1 -1;0 2];c = [1.9691 6.4493];d=[0 0];sys_ss = ss(a,b,c,0)sys_tf=tf(sys_ss)sys_zpk=zpk(sys_ss)[z,p,k]=ss2zp(a,b,c,d,2)4.2clc;clear;close all;a = [-0.5572 -0.7814;0.7814 0];b = [1 -1;0 2];c = [1.9691 6.4493];d=[0 0];sys = ss(a,b,c,d)sysd=c2d(sys,0.1)figure (1) ;step(sys,r-,sysd,b--)figure (2) ;bode(sys,r-) hold;bode(sysd,b--)5.clc;clear;close all;H=tf([25],[1 6 25]);step(H,r-)[y,t]=step(H);ypeak=max(y);n=length(y);for k=1:n if(y(k)==ypeak)peaktime=t(k) temp=k;%记住此时对应数组的位置 endendystead=y(n);over=(ypeak-ystead)/ypeak%上升时间，题目要求0到100%，实际定义的10%-90%，按题目要求for k=1:temp if(y(k)=ystead)risetime=t(k) break; endend%调节时间，98%-102%ydown=0.98;yup=1.02;for k=temp:n if(y(k)=yup)settime=t(k) break; endend6.clc;clear;close all;H=tf([1],[1 1]);Hd=c2d(H,0.1);figure(1);step(H,r-,Hd,b--);figure(2);impulse(H,r-,Hd,b--);figure(3);bode(H,r-,Hd,b--);7.clc;clear;num=[1 0 3 4 0 ];d