计算机仿真课程设计(3)

Gczw =

1/z*(1+363/125/z)*(1-6703/10000/z)*(1-3011/10000/z)*(1+477/2500/z)*(925332561593647/2251799813685248+2353474703694299/4503599627370496/z)

>> f7=limit(Gezw,z,inf)-1

f7 =

b1-1

>> f8=subs(Gezw,z,-2.904)-1

f8 =

238144/131769*b1-29768000/47832147*b2+3721000000/17363069361*b3-465125000000/6302794178043*b4+58140625000000/2287914286629609*b5-1

>> f9=subs(Gezw,z,0.6703)-1

f9 =

10870209/44930209*b1+108702090000/301167190927*b2+1087020900000000/2018723680783681*b3+10870209000000000000/13531504832293013743*b4+108702090000000000000000/90701676890860071119329*b5-1

>> f10=subs(Gezw,z,0.3011)-1

f10 =

48846121/9066121*b1+488461210000/27298090331*b2+4884612100000000/82194549986641*b3+48846121000000000000/247487790009776051*b4+488461210000000000000000/745185735719435689561*b5-1

>> f11=subs(Gezw,z,-0.1908)-1

f11 =

8862529/227529*b1-22156322500/108531333*b2+55390806250000/51769445841*b3-138477015625000000/24694025666157*b4+346192539062500000000/11779050242756889*b5-1

>> [b1j,b2j,b3j,b4j,b5j]=solve(f7,f8,f9,f10,f11)

b1j =1

b2j =

27824558815573855050308402/17509966474003480259529201

b3j =

1958599880988746129641979/2501423782000497179932743

b4j =

-19421775618281782311371837/93386487861351894717489072

b5j =

-7276583151737082934420395/124515317148469192956652096

>> B=double([b1j,b2j,b3j,b4j,b5j])

B =

1.0000 1.5891 0.7830 -0.2080 -0.0584

>> Gezw=subs(Gezw,[b1 b2 b3 b4 b5],B)

Gezw =

(1-1/z)^2*(1+3578267065776849/2251799813685248/z+3526291609466701/4503599627370496

/z^2-7492981345843367/36028797018963968/z^3-2105496282460953/36028797018963968/z^4)

13、求针对单位速度信号输入的最少拍无波纹控制器Dw(z)并说明Dw(z)的可实现性。

>>Gz=0.68436*z^-1*(1+2.904*z^-1)*(1-0.6703*z^-1)*(1-0.3011*z^-1)*(1+0.1908*z^-1)/(1-z^-1)^2/(1-0.8187*z^-1)/(1-0.3679*z^-1)/(1-0.2019*z^-1)

Gz =

17109/25000/z*(1+363/125/z)*(1-6703/10000/z)*(1-3011/10000/z)*(1+477/2500/z)/(1-1/z)^2/(1-8187/10000/z)/(1-3679/10000/z)/(1-2019/10000/z)

>> Dwz=Gczw/Gezw/Gz

Dwz =

25000/17109*(925332561593647/2251799813685248+2353474703694299/4503599627370496/z)/(1+3578267065776849/2251799813685248/z+3526291609466701/4503599627370496/z^2-7492981345843367/36028797018963968/z^3-2105496282460953/36028797018963968/z^4)*(1-8187/10000/z)*(1-3679/10000/z)*(1-2019/10000/z)

14、用程序仿真方法分析单位速度信号输入时闭环系统动态性能和稳态性能。

>> [N,D]=numden(simplify(Gczw));

>> numc=sym2poly(N)

numc =1.0e+029 *

0.5783 1.9635 0.2602 -1.6051 0.1282 0.0822

>> denc=sym2poly(D)

denc =

1.0e+029 *

1.4074 0 0 0 0 0 0

>> t=0:0.2:2

t =

0 0.2000 0.4000 0.6000 0.8000 1.0000

>> u=t

u =

0 0.2000 0.4000 0.6000 0.8000 1.0000

>> dlsim(numc,denc,u)