實驗二連續(xù)時間系統(tǒng)的頻率響應(yīng)

1、實驗二連續(xù)時間系統(tǒng)的頻率響應(yīng)39022622 龔小川.實驗?zāi)康模?.進一步加深對連續(xù)時間系統(tǒng)頻率響應(yīng)理解；2 .掌握借助計算機計算任意連續(xù)時間系統(tǒng)頻率響應(yīng)的方法。二. 實驗原理1.本實驗的基本內(nèi)容就是將系統(tǒng)函數(shù)的幅頻特性曲線以及相頻特性曲線給畫出來。而系統(tǒng)函數(shù)H(s)mK (s引j 1n(si 1jw，則 H (jw)mK (jw Zj)j 1nPi)(jw Pi)jwZjN jejjwPiiej(jw )(jw )mj.e j 1jjie(w)即(1)計算所有零點模之積及極點模之積，兩者之商即為(2)計算所有零點相角之和及極點相角之和，兩者之差即為j ( w )(jw ) eH (s)的幅度

2、;H (s)的相角。2.通過零極點圖通過幾何的方法來計算，而且通過零極點圖可以迅速地判斷系統(tǒng)的濾波特性。通過零極點圖進行計算的方法是:(1 )在S平面上標(biāo)出系統(tǒng)的零極點位置；(2) 選擇S平面的坐標(biāo)原點為起始點，沿虛軸向上移動，計算此時各極點和零點與該點的膜和夾角；(3) 將所有零點的模相乘，再除以各極點的模，得到對應(yīng)頻率處的幅頻特性的值;(4 )將所有零點的幅角相加，減去各極點的幅角，得到對應(yīng)頻率處的相角。三. 實驗流程圖jm ?輸出輸出tempi,temp2,即即為零點為極點模之積模之積令 i=0,temp2=1in ?ternj - ternji *y+輸出 temp1/temp2,極

3、為H(s)之幅度輸出 temp1-temp2,極 為H(s)之相角四. 實驗代碼#i nclude stdio.h#in elude math.h#in elude graphics.hfloat ata nn ew(float t1,float t2);int mai n()float z102,p102,out1003,f50;int a,iout;int i,itemp,k,ktemp,j,m, n;float w,temp1,temp2,prew,pretemp1,pretemp2,ptemp,h,fout;int gdriver,gmode=0;char s10;gdriver=O;

4、/*/for(i=0;i=10;i+)a=sca nf(%f %fj, &ziO, &zi1);if(a!=0);elsefflush(stdin);/*清空輸入緩沖區(qū)，使得第二個scanf不被忽略*/for(k=0;k=10;k+)a=sca nf(%f %fj,&pkO, &pk1);if(a!=0);elsebreak;break;/*/initgraph(&gdriver,&gmode, e:tcbgi);setbkcolor(IO);setcolor(4);setli nestyle(1,0,1);recta ngle(50,20,600,420);for(ptemp=20;ptem

5、p=420;ptemp=ptemp+50) li ne(50,ptemp,600,ptemp);for(ptemp=50;ptemp=600;ptemp=ptemp+50) lin e(ptemp,20,ptemp,420);setl in estyle(0,0,1);setcolor(1);lin e(200,20,200,420);lin e(50,420,600,420);lin e(200,20,190,30);lin e(200,20,210,30);outtextxy(210,30,A);lin e(600,420,590,410);lin e(600,420,590,430);o

6、uttextxy(600,430,w);for(iout=0;iout=5;iout+)spri ntf(s,%d,iout);outtextxy(iout*50+200,430,s);for(fout=1.0;fout=8;fout=fout+1)spri ntf(s,%.1f,fout/10);outtextxy(170,420-fout*50,s);setcolor(4);/*prin tf(w幅度相角 n);*/for(w=0,j=0;w=5.0;j+)temp1=1;for(itemp=0;itemp=i-1;itemp+)temp仁temp1*sqrt(zitemp0*zitemp

7、0+(zitemp1-w)*(zitemp1-w);for(ktemp=0;ktemp0) lin e(prew*40+200,420-pretemp1*500,w*40+200,420-temp1*500);prew=w;pretemp仁tempi;w=w+0.1;for(m=0;moutm+11)break;prin tf(%.4f,%f,outm1,outm0);lin e(100,420-outm1/sqrt(2)*500,500,420-outm1/sqrt(2)*500);getch();closegraph();initgraph(&gdriver,&gmode, e:tcbgi

8、);setbkcolor(10);setcolor(4);setli nestyle(1,0,1);recta ngle(50,20,600,420);for(ptemp=20;ptemp=420;ptemp=ptemp+50)li ne(50,ptemp,600,ptemp);for(ptemp=50;ptemp=600;ptemp=ptemp+50) lin e(ptemp,20,ptemp,420);for(w=0,j=0;w=5;j+)temp2=0;for(ktemp=0;ktemp=k-1;ktemp+)temp2=temp2-ata nn ew(-pktemp0,w-pktemp

9、1);for(itemp=0;itemp=180)temp2=temp2-360;else if(temp20)lin e(prew*50+200,180+pretemp2,w*50+200,180+temp2);prew=w;pretemp2=temp2;w=w+0.1;getch();closegraph();for( j=0;j=50;j+)prin tf(%.2f,%.4f,%.2f*,outj0,outj1,outj2);prin tf(nn);for(n=0,j=0;n=50 ;n+)fn=fabs(outm1/sqrt(2)-outn 1);for(n=0,h=f0,j=0;nf

10、 n)h=f n;j=n;prin tf(Fl=%.1fn,outj0);for(n=m,h=f0,j=0;nf n)h=f n;j=n;prin tf(Fh=%.1f,outj0);getch();return 0; float ata nn ew(float t1,float t2)if(t1 0&t20)return ata n( t2/t1)/3.14*180;else if(t1 0&t20)return 180+ata n(t2/t1)/3.14*180;else if(t10)return 180+ata n(t2/t1)/3.14*180;else if(t1 0&t20)re

11、turn 90;else if(t1=0&t2 0&t2=0)return 180;else if(t1 0&t2=0)return -180;五. 實驗數(shù)據(jù)及所繪圖形零點 z1= 0;極點 p1=-1-j,p2=-1+jw|H(jw)|0(jw)w|H( jw)| .(i Iii|Ii ii si1Iia iiii i| P/5)ii；tD.i DJ. _ _- .-1111111111i11 1I1i11i:!-(0 亦 IT)1 11 1111iI1ii1,1! 0i1i11 ii!11i丨ii11iI1J1零點 z仁2+1.5j,z2=2-1.5j; 極點p仁-0.5-0.7j,p2=

12、-0.5+0.7j,p3=-2.5-2j,p4=-2.5+2j0.00.823995197.6069762.50.125927-30.1683910.10.827748188.7530952.60.117968-34.2282300.20.838357179.0547982.70.110810-38.1362930.30.853583168.2866332.80.104336-41.8993550.40.868832156.2324432.90.098449-45.5231510.50.876442142.77567830.093071-49.0127090.60.866588128.0477

13、023.10.088137-52.3725880.70.830997112.5538083.20.083592-55.6070560.80.76958397.1076083.30.079392-58.7202040.90.69100282.5233353.40.075499-61.7160311.00.60750369.3081333.50.071880-64.5984821.10.52856857.5963953.60.068509-67.3714831.20.45889247.2757383.70.065362-70.0389431.30.39964338.1349873.80.062418-72.6047581.40.35012029.9547293.90.059661-75.0727981.50.30896222.54471840.057073-77.4468961.60.27472715.7524114.10.054642-79.7308371.70.2461209.4598394.20.052356-81.9283401.80.2220603.5772464.30.050202-84.0430451.90.201673-1.9632154.40.048172-