实验二 连续时间系统的频率响应.docx

1、实验二 连续时间系统的频率响应实验二 连续时间系统得频率响应 龚小川一.实验目得： 1、 进一步加深对连续时间系统频率响应理解;2.掌握借助计算机计算任意连续时间系统频率响应得方法。二.实验原理、本实验得基本内容就就是将系统函数得幅频特性曲线以及相频特性曲线给画出来。而系统函数，令,则即（)计算所有零点模之积及极点模之积,两者之商即为得幅度； (2)计算所有零点相角之与及极点相角之与,两者之差即为得相角。2、通过零极点图通过几何得方法来计算,而且通过零极点图可以迅速地判断系统得滤波特性。通过零极点图进行计算得方法就是:(1)在 平面上标出系统得零极点位置；(2)选择 平面得坐标原点为起始点,沿

2、虚轴向上移动,计算此时各极点与零点与该点得膜与夹角;()将所有零点得模相乘,再除以各极点得模,得到对应频率处得幅频特性得值；（4)将所有零点得幅角相加，减去各极点得幅角,得到对应频率处得相角。三.实验流程图jm? 否就是输出temp1,即为零点相角之与令j=0,temp1=0in?否就是输出temp2,即为极点相角之与令i=0,temp2=0输出temp1-temp2,极为H(s)之相角四.实验代码include stdi、h#includmath、h#nleahi、hfloattannw（ot 1，oa t2)；t mi() flot z10,p102,out003,f5; t ，out;

3、inti,iemp，,ktem,j，m,n;float w,temp，tmp,prew,petm1,prtep2,ptem,h,fou;int gdriver,goe=0;chr 1； drive=0; /* for(=;i=10；+） a=scf(%fj,i0,&i1)； if(a!=0）； ele flsh（std）;/*清空输入缓冲区,使得第二个scanf不被忽略/ for（k=0;=10;k+) =ca(%f fj,&p0,pk)； （a！=0）; els bek; rea; */itgrp(grive，&god, e:tgi)；stbkcolor(10）； stcolor（); se

4、tinestyle(1,0,)； rectanle(50,0，0,20); for（pte=2；ptep420;ptem=pmp+50) lie(50，ptep,600,pemp); or(t=;temp=600;ptmp=ptemp+50) li（ptemp,20，ptem,420)； setlinstye(0,0,1)； setcor(1); line(20,20,200，420); lne（50,420，600,420); lie(200,20,190，3); line(0，2,210,0); outxy(10，30,A); lin（600，4,90，410）; i(00,42,59，3

5、0); outety(6,43,）; fr(iou=0；ot=5;iu+) sprint(,d,iu)； outxtxy(iot*500,s); or（fout=1、0；out8；ut=fou+1) itf(,%1f,fu/10); outxxy(1,0fou5,s）; setcolor(4); /*prinf( 幅度 相角n）；/for(=0,j;=5、0;j+） emp1=1; for（itep0;item0) lin（prew*4200，420-preemp1*500,*40+20，40-*5）； pew； preemp1=tem1; w=+0、1; for(=0;mut+) brek;

6、 prinf(%4f,f，t1,out0); ine(0，42utm1/qrt(2）*50,00,2-otm1/sqrt(）*00); tch(); loga(); inigh(&gdivr,gd, :cbgi); etbkclor(10)； stcoor(4); selietyle(1,，1)； rectnge(50，20,600,42)； o(pemp=2;ptemp=;pmp=ptep+50） le(0,ptm,60，tm); fr(mp=0;temp=600;pm=tem+50) lne(pm,20，ptemp，4); for（w=0，j=0;w=5;j+) temp=0； for(k

7、ep=;ktempk;ktemp+） temptemp-atannew（-ptep0，pktmp)；o(itep=;itemp=8) m=tem2-30； ls if(t) ln（prw*50+00,8+pretmp2,*50+20，10p);rew=w; pretemp2=emp2;w=w+0、1； getch（); losgrp（); for(j=0；j0;j+) ritf(%.2，%.f,%.2f*,tj0,outj1,otj2); prif(n）; fo（n0,j=0;=50；n+) n=fabs(outm1/sqrt(2)-ut1); for(=，=f0，j=；fn） h=f; =n

8、; printf(Fl%1n,ut0); or(n=m,h0,;nn) h=f; jn; rnt(Fh=%.1f,ut); ec（); rturn ；float aannew(lot t1,ltt) i(t&t20) retun an(2/t)3、4*180;ls if(t10&20） return 0+tan(21)/、4*8;ele if(t10&t20) retun ata(t2/t1)3、1*80; els if(t=0&t20） return 90；lse f（t1=0&20) retur -9;elsef(t=0&t2=） retrn 0; elif(0&t2=0) return

9、180; se if(t0&=) etur -180; 五实验数据及所绘图形零点1=;极点 =-1-j,p21+jw|H()|(j）wH(w)|()0、0、0999984、26712、60、381-4、4760、2、099978、71152、7、577-44、437、30、18482、5592、0、06-4、20520、40、19936366、0633、90、335472-47、8660、5、28960、55、32539-4、40、60、2525553、103、0、315814-50、83340、0、394647、16263、20、36705-、1670、80、80740、36753、30、28046-5、4350、90、170923、471、0、289812-54、5790、447226、5673、0、2819855、7391、10、470519、7523、60、274528-56、7019、20、489113、4983、0、26731-57、669