《數(shù)字信號(hào)處理(第三版)》第10章實(shí)驗(yàn)程序(驗(yàn)證版)(西安電子科技大學(xué)出版_主編：高西全_丁玉美)

1、數(shù)字信號(hào)處理(第三版）第10章實(shí)驗(yàn)程序注意：每次實(shí)驗(yàn)前，請(qǐng)新建以下8個(gè)m文件并依實(shí)驗(yàn)需要（已注明）移動(dòng)至matlab當(dāng)前工作目錄下，便于主程序調(diào)用。各程序均已驗(yàn)證，請(qǐng)直接復(fù)制即可。文件名：tstem.m（實(shí)驗(yàn)一、二需要）程序：function tstem(xn,yn)%時(shí)域序列繪圖函數(shù)%xn:被繪圖的信號(hào)數(shù)據(jù)序列，yn:繪圖信號(hào)的縱坐標(biāo)名稱（字符串）n=0:length(xn)-1;stem(n,xn,'.');xlabel('n');ylabel('yn');axis(0,n(end),min(xn),1.2*max(xn);文件名：tplot

2、.m（實(shí)驗(yàn)一、四需要）程序：function tplot(xn,T,yn)%時(shí)域序列連續(xù)曲線繪圖函數(shù)%xn:信號(hào)數(shù)據(jù)序列,yn:繪圖信號(hào)的縱坐標(biāo)名稱（字符串）%T為采樣間隔n=0;length(xn)-1;t=n*T;plot(t,xn);xlabel('t/s');ylabel(yn);axis(0,t(end),min(xn),1.2*max(xn);文件名：myplot.m（實(shí)驗(yàn)一、四需要）程序：%(1)myplot;計(jì)算時(shí)域離散系統(tǒng)損耗函數(shù)并繪制曲線圖。 function myplot(B,A)%B為系統(tǒng)函數(shù)分子多項(xiàng)式系數(shù)向量%A為系統(tǒng)函數(shù)分母多項(xiàng)式系數(shù)向量H,W=fr

3、eqz(B,A,1000)m=abs(H);plot(W/pi,20*log10(m/max(m);grid on;xlabel('omega/pi');ylabel('幅度(dB)')axis(0,1,-80,5);title('損耗函數(shù)曲線');文件名：mstem.m（實(shí)驗(yàn)一、三需要）程序：function mstem(Xk)%mstem(Xk)繪制頻域采樣序列向量Xk的幅頻特性圖M=length(Xk);k=0:M-1;wk=2*k/M;%產(chǎn)生M點(diǎn)DFT對(duì)應(yīng)的采樣點(diǎn)頻率（關(guān)于pi歸一化值）stem(wk,abs(Xk),'.'

4、;);box on;%繪制M點(diǎn)DFT的幅頻特性圖xlabel('w/pi');ylabel('幅度');axis(0,2,0,1.2*max(abs(Xk);文件名：mpplot.m（實(shí)驗(yàn)一需要）程序：%(2)mpplot;計(jì)算時(shí)域離散系統(tǒng)損耗函數(shù)和相頻特性函數(shù)，并繪制曲線圖。function mpplot(B,A,Rs)%mpplot(B,A,Rs)%時(shí)域離散系統(tǒng)損耗函數(shù)和相頻特性繪圖%B為系統(tǒng)函數(shù)分子多項(xiàng)式系數(shù)向量%A為系統(tǒng)函數(shù)分母多項(xiàng)式系數(shù)向量%Rs為濾波器阻帶最小衰減，省略則幅頻曲線最小值取80dBif nargin<3 ymin=-80;else

5、 ymin=-Rs-20;end;%確定幅頻曲線縱坐標(biāo)最小值H,W=H,W=freqz(B,A,1000)m=abs(H);subplot(2,2,1);plot(W/pi,20*log10(m/max(m);grid on;xlabel('omega/pi');ylabel('幅度(dB)')axis(0,1,ymin,5);title('損耗函數(shù)曲線');subplot(2,2,3);plot(W/pi,p/pi);xlabel('omega/pi');ylabel('相位/pi');grid on;title

6、('(b)相頻特性曲線');文件名：mfftplot.m（實(shí)驗(yàn)一需要）程序：function mfftplot(xn,N)%mfftplot(xn,N)計(jì)算序列向量xn的N點(diǎn)fft并繪制其幅頻特性曲線Xk=fft(xn,N);%計(jì)算信號(hào)xn的頻譜的N點(diǎn)采樣%=以下為繪圖部分=k=0:N-1;wk=2*k/N;m=abs(Xk);mm=max(m);plot(wk,m/mm);grid on;xlabel('omega/pi');ylabel('幅度(dB)');axis(0,2,0,1.2);title('幅度特性曲線');文件名

7、：mstg.m（實(shí)驗(yàn)四需要）程序：function st=mstg%產(chǎn)生信號(hào)序列向量st,并顯示st的時(shí)域波形和頻譜%st=mstg 返回三路調(diào)幅信號(hào)相加形成的混合信號(hào)，長度N=1600N=1600; %N為信號(hào)st的長度Fs=10000;T=1/Fs;Tp=N*T; %采樣頻率Fs=10kHz，Tp為采樣時(shí)間t=0:T:(N-1)*T;k=0:N-1;f=k/Tp;fc1=Fs/10;%第1路調(diào)幅信號(hào)的載波頻率fc1=1000Hzfm1=fc1/10; %第1路調(diào)幅信號(hào)的調(diào)制信號(hào)頻率fm1=100Hzfc2=Fs/20; %第2路調(diào)幅信號(hào)的載波頻率fc2=500Hzfm2=fc2/10; %

8、第2路調(diào)幅信號(hào)的調(diào)制信號(hào)頻率fm2=50Hzfc3=Fs/40;%第3路調(diào)幅信號(hào)的載波頻率fc3=250Hzfm3=fc3/10; %第3路調(diào)幅信號(hào)的調(diào)制信號(hào)頻率fm3=25Hzxt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %產(chǎn)生第1路調(diào)幅信號(hào)xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %產(chǎn)生第2路調(diào)幅信號(hào)xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %產(chǎn)生第3路調(diào)幅信號(hào)st=xt1+xt2+xt3; %三路調(diào)幅信號(hào)相加fxt=fft(st,N); %計(jì)算信號(hào)st的頻譜%=以下為繪圖部分，繪制st的時(shí)域

9、波形和幅頻特性曲線=subplot(3,1,1);plot(t,st);grid;xlabel('t/s');ylabel('s(t)');axis(0,Tp/8,min(st),max(st);title('(a) s(t)的波形');subplot(3,1,2);stem(f,abs(fxt)/max(abs(fxt),'.');grid;title('(b) s(t)的頻譜');axis(0,Fs/5,0,1.2);xlabel('f/Hz');ylabel('幅度');文件名

10、：xtg.m（實(shí)驗(yàn)五需要）程序：function xt=xtg(N)%實(shí)驗(yàn)五信號(hào)x(t)產(chǎn)生,并顯示信號(hào)的幅頻特性曲線%xt=xtg(N) 產(chǎn)生一個(gè)長度為N,有加性高頻噪聲的單頻調(diào)幅信號(hào)xt,采樣頻率Fs=1000Hz%載波頻率fc=Fs/10=100Hz,調(diào)制正弦波頻率f0=fc/10=10Hz.Fs=1000;T=1/Fs;Tp=N*T;t=0:T:(N-1)*T;fc=Fs/10;f0=fc/10; %載波頻率fc=Fs/10，單頻調(diào)制信號(hào)頻率為f0=Fc/10;mt=cos(2*pi*f0*t);  %產(chǎn)生單頻正弦波調(diào)制信號(hào)mt，頻率為f0ct=cos(2*pi*fc*t);

11、  %產(chǎn)生載波正弦波信號(hào)ct，頻率為fcxt=mt.*ct;          %相乘產(chǎn)生單頻調(diào)制信號(hào)xtnt=2*rand(1,N)-1;   %產(chǎn)生隨機(jī)噪聲nt%設(shè)計(jì)高通濾波器hn,用于濾除噪聲nt中的低頻成分,生成高通噪聲%=fp=150; fs=200;Rp=0.1;As=70;     % 濾波器指標(biāo)fb=fp,fs;m=0,1;         

12、        % 計(jì)算remezord函數(shù)所需參數(shù)f,m,devdev=10(-As/20),(10(Rp/20)-1)/(10(Rp/20)+1);n,fo,mo,W=remezord(fb,m,dev,Fs);    % 確定remez函數(shù)所需參數(shù)hn=remez(n,fo,mo,W);     % 調(diào)用remez函數(shù)進(jìn)行設(shè)計(jì),用于濾除噪聲nt中的低頻成分yt=filter(hn,1,10*nt);     

13、%濾除隨機(jī)噪聲中低頻成分，生成高通噪聲ytxt=xt+yt;           %噪聲加信號(hào)fst=fft(xt,N);k=0:N-1;f=k/Tp;subplot(3,1,1);plot(t,xt);grid;xlabel('t/s');ylabel('x(t)');axis(0,Tp/5,min(xt),max(xt);title('(a) 信號(hào)加噪聲波形')subplot(3,1,2);plot(f,abs(fst)/max(abs(fst

14、);grid;title('(b) 信號(hào)加噪聲的頻譜')axis(0,Fs/2,0,1.2);xlabel('f/Hz');ylabel('幅度')10.1 系統(tǒng)響應(yīng)及系統(tǒng)穩(wěn)定性close all;clear all;clc;%內(nèi)容1：調(diào)用filter解差分方程，由系統(tǒng)對(duì)u(n)的響應(yīng)判斷穩(wěn)定性%=A=1,-0.9;B=0.05,0.05;  %系統(tǒng)差分方程系數(shù)向量B和Ax1n=1 1 1 1 1 1 1 1 zeros(1,50);  %產(chǎn)生信號(hào)x1(n)=R8(n)x2n=ones(1,128);  &

15、#160; %產(chǎn)生信號(hào)x2(n)=u(n)hn=impz(B,A,58);    %求系統(tǒng)單位脈沖響應(yīng)h(n)subplot(2,2,1);y='h(n)'tstem(hn,y);    %調(diào)用函數(shù)tstem繪圖title('(a) 系統(tǒng)單位脈沖響應(yīng)h(n)');box ony1n=filter(B,A,x1n);    %求系統(tǒng)對(duì)x1(n)的響應(yīng)y1(n)subplot(2,2,2);y='y1(n)'tstem(y1n,y);title('(b)

16、 系統(tǒng)對(duì)R8(n)的響應(yīng)y1(n)');box ony2n=filter(B,A,x2n);    %求系統(tǒng)對(duì)x2(n)的響應(yīng)y2(n)subplot(2,2,4);y='y2(n)'tstem(y2n,y);title('(c) 系統(tǒng)對(duì)u(n)的響應(yīng)y2(n)');box on%內(nèi)容2：調(diào)用conv函數(shù)計(jì)算卷積%=x1n=1 1 1 1 1 1 1 1 ;  %產(chǎn)生信號(hào)x1(n)=R8(n)h1n=ones(1,10) zeros(1,10);h2n=1 2.5 2.5 1 zeros(1,10);y21n=co

17、nv(h1n,x1n);y22n=conv(h2n,x1n);figure(2)subplot(2,2,1);y='h1(n)'tstem(h1n,y);    %調(diào)用函數(shù)tstem繪圖title('(d) 系統(tǒng)單位脈沖響應(yīng)h1(n)');box onsubplot(2,2,2);y='y21(n)'tstem(y21n,y);title('(e) h1(n)與R8(n)的卷積y21(n)');box onsubplot(2,2,3);y='h2(n)'tstem(h2n,y);

18、60;   %調(diào)用函數(shù)tstem繪圖title('(f) 系統(tǒng)單位脈沖響應(yīng)h2(n)');box onsubplot(2,2,4);y='y22(n)'tstem(y22n,y);title('(g) h2(n)與R8(n)的卷積y22(n)');box on%內(nèi)容3：諧振器分析%=un=ones(1,256);    %產(chǎn)生信號(hào)u(n)n=0:255;xsin=sin(0.014*n)+sin(0.4*n); %產(chǎn)生正弦信號(hào)A=1,-1.8237,0.9801;B=1/100.49,0,-1/1

19、00.49;  %系統(tǒng)差分方程系數(shù)向量B和Ay31n=filter(B,A,un);    %諧振器對(duì)u(n)的響應(yīng)y31(n)y32n=filter(B,A,xsin);    %諧振器對(duì)u(n)的響應(yīng)y31(n)figure(3)subplot(2,1,1);y='y31(n)'tstem(y31n,y);title('(h) 諧振器對(duì)u(n)的響應(yīng)y31(n)');box onsubplot(2,1,2);y='y32(n)'tstem(y32n,y);title(

20、9;(i) 諧振器對(duì)正弦信號(hào)的響應(yīng)y32(n)');box on10.2 時(shí)域采樣與頻域采樣時(shí)域采樣理論驗(yàn)證程序exp2a.m：%= close all;clear all;clc;Tp=64/1000;             %觀察時(shí)間Tp=64微秒%產(chǎn)生M長采樣序列x(n)% Fs=1000;T=1/Fs;  Fs=1000;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.128;alph=pi*50*20.5;omega=pi*5

21、0*20.5;xnt=A*exp(-alph*n*T).*sin(omega*n*T);Xk=T*fft(xnt,M);  %M點(diǎn)FFTxnt)yn='xa(nT)'subplot(3,2,1);tstem(xnt,yn);           %調(diào)用自編繪圖函數(shù)tstem繪制序列圖box on;title('(a) Fs=1000Hz');k=0:M-1;fk=k/Tp;subplot(3,2,2);plot(fk,abs(Xk);title('

22、;(a) T*FTxa(nT),Fs=1000Hz');xlabel('f(Hz)');ylabel('幅度');axis(0,Fs,0,1.2*max(abs(Xk)%= % Fs=300Hz和 Fs=200Hz的程序與上面Fs=1000Hz完全相同。 頻域采樣理論驗(yàn)證程序exp2b.m：%=close all;clear all;clc;M=27;N=32;n=0:M;%產(chǎn)生M長三角波序列x(n)xa=0:floor(M/2);  xb= ceil(M/2)-1:-1:0; xn=xa,xb;Xk=fft(xn,1024); %1

23、024點(diǎn)FFTx(n), 用于近似序列x(n)的TFX32k=fft(xn,32)  ;%32點(diǎn)FFTx(n)x32n=ifft(X32k); %32點(diǎn)IFFTX32(k)得到x32(n)X16k=X32k(1:2:N); %隔點(diǎn)抽取X32k得到X16(K)x16n=ifft(X16k,N/2); %16點(diǎn)IFFTX16(k)得到x16(n)subplot(3,2,2);stem(n,xn,'.');box ontitle('(b) 三角波序列x(n)');xlabel('n');ylabel('x(n)&#

24、39;);axis(0,32,0,20)k=0:1023;wk=2*k/1024; subplot(3,2,1);plot(wk,abs(Xk);title('(a)FTx(n)');xlabel('omega/pi');ylabel('|X(ejomega)|');axis(0,1,0,200)k=0:N/2-1;subplot(3,2,3);stem(k,abs(X16k),'.');box ontitle('(c) 16點(diǎn)頻域采樣');xlabel('k');ylabel('

25、|X_1_6(k)|');axis(0,8,0,200)n1=0:N/2-1;subplot(3,2,4);stem(n1,x16n,'.');box ontitle('(d) 16點(diǎn)IDFTX_1_6(k)');xlabel('n');ylabel('x_1_6(n)');axis(0,32,0,20)k=0:N-1;subplot(3,2,5);stem(k,abs(X32k),'.');box ontitle('(e) 32點(diǎn)頻域采樣');xlabel('k');yla

26、bel('|X_3_2(k)|');axis(0,16,0,200)n1=0:N-1;subplot(3,2,6);stem(n1,x32n,'.');box ontitle('(f) 32點(diǎn)IDFTX_3_2(k)');xlabel('n');ylabel('x_3_2(n)');axis(0,32,0,20)10.3 用FFT對(duì)信號(hào)做頻譜分析exp301.m：%= clear all;close all;clc;%實(shí)驗(yàn)內(nèi)容(1)%=x1n=ones(1,4);    %產(chǎn)生序列向量

27、x1(n)=R4(n)M=8;xa=1:(M/2);  xb=(M/2):-1:1; x2n=xa,xb;    %產(chǎn)生長度為8的三角波序列x2(n)x3n=xb,xa;X1k8=fft(x1n,8);        %計(jì)算x1n的8點(diǎn)DFTX1k16=fft(x1n,16);      %計(jì)算x1n的16點(diǎn)DFTX2k8=fft(x2n,8);        %計(jì)算

28、x1n的8點(diǎn)DFTX2k16=fft(x2n,16);        %計(jì)算x1n的16點(diǎn)DFTX3k8=fft(x3n,8);        %計(jì)算x1n的8點(diǎn)DFTX3k16=fft(x3n,16);        %計(jì)算x1n的16點(diǎn)DFT%以下繪制幅頻特性曲線subplot(2,2,1);mstem(X1k8); %繪制8點(diǎn)DFT的幅頻特性圖title('(1a) 8點(diǎn)D

29、FTx_1(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X1k8)subplot(2,2,3);mstem(X1k16); %繪制16點(diǎn)DFT的幅頻特性圖title('(1b)16點(diǎn)DFTx_1(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X1k16)figure(2)subplot(2,2,1);mstem(X2k8); %繪制8點(diǎn)DFT的幅頻特性圖title('(2a) 8

30、點(diǎn)DFTx_2(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X2k8)subplot(2,2,2);mstem(X2k16); %繪制16點(diǎn)DFT的幅頻特性圖title('(2b)16點(diǎn)DFTx_2(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X2k16)subplot(2,2,3);mstem(X3k8); %繪制8點(diǎn)DFT的幅頻特性圖title('(3a) 8點(diǎn)DFTx_3

31、(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X3k8)subplot(2,2,4);mstem(X3k16);%繪制16點(diǎn)DFT的幅頻特性圖title('(3b)16點(diǎn)DFTx_3(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X3k16) exp302.m：%實(shí)驗(yàn)內(nèi)容(2) 周期序列譜分析%=close all;clear all;clc;N=8;n=0:N-1; 

32、;       %FFT的變換區(qū)間N=8x4n=cos(pi*n/4);x5n=cos(pi*n/4)+cos(pi*n/8);X4k8=fft(x4n);       %計(jì)算x4n的8點(diǎn)DFTX5k8=fft(x5n);       %計(jì)算x5n的8點(diǎn)DFTN=16;n=0:N-1;      %FFT的變換區(qū)間N=16x4n=cos(pi*n/4);x5

33、n=cos(pi*n/4)+cos(pi*n/8);X4k16=fft(x4n);      %計(jì)算x4n的16點(diǎn)DFTX5k16=fft(x5n);      %計(jì)算x5n的16點(diǎn)DFTfigure(3)subplot(2,2,1);mstem(X4k8); %繪制8點(diǎn)DFT的幅頻特性圖title('(4a) 8點(diǎn)DFTx_4(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X4k8)su

34、bplot(2,2,3);mstem(X4k16); %繪制16點(diǎn)DFT的幅頻特性圖title('(4b)16點(diǎn)DFTx_4(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X4k16)subplot(2,2,2);mstem(X5k8); %繪制8點(diǎn)DFT的幅頻特性圖title('(5a) 8點(diǎn)DFTx_5(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X5k8)subplot(2

35、,2,4);mstem(X5k16); %繪制16點(diǎn)DFT的幅頻特性圖title('(5b)16點(diǎn)DFTx_5(n)');xlabel('/');ylabel('幅度');axis(0,2,0,1.2*max(abs(X5k16) exp303.m：%實(shí)驗(yàn)內(nèi)容(3) 模擬周期信號(hào)譜分析%=close all;clear all;clc;figure(4)Fs=64;T=1/Fs;N=16;n=0:N-1;        %FFT的變換區(qū)間N=16x6nT=cos(8*

36、pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);   %對(duì)x6(t)16點(diǎn)采樣X6k16=fft(x6nT);      %計(jì)算x6nT的16點(diǎn)DFTX6k16=fftshift(X6k16);      %將零頻率移到頻譜中心Tp=N*T;F=1/Tp;   %頻率分辨率Fk=-N/2:N/2-1;fk=k*F;       %產(chǎn)生16點(diǎn)DFT對(duì)應(yīng)的采樣點(diǎn)頻率（以零頻

37、率為中心）subplot(3,1,1);stem(fk,abs(X6k16),'.');box on %繪制8點(diǎn)DFT的幅頻特性圖title('(6a) 16點(diǎn)|DFTx_6(nT)|');xlabel('f(Hz)');ylabel('幅度');axis(-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k16)N=32;n=0:N-1;        %FFT的變換區(qū)間N=16x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)

38、+cos(20*pi*n*T);   %對(duì)x6(t)32點(diǎn)采樣X6k32=fft(x6nT);      %計(jì)算x6nT的32點(diǎn)DFTX6k32=fftshift(X6k32);      %將零頻率移到頻譜中心Tp=N*T;F=1/Tp;   %頻率分辨率Fk=-N/2:N/2-1;fk=k*F;       %產(chǎn)生16點(diǎn)DFT對(duì)應(yīng)的采樣點(diǎn)頻率（以零頻率為中心）subplot(3,1,2);st

39、em(fk,abs(X6k32),'.');box on %繪制8點(diǎn)DFT的幅頻特性圖title('(6b) 32點(diǎn)|DFTx_6(nT)|');xlabel('f(Hz)');ylabel('幅度');axis(-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k32)N=64;n=0:N-1;        %FFT的變換區(qū)間N=16x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); 

40、  %對(duì)x6(t)64點(diǎn)采樣X6k64=fft(x6nT);      %計(jì)算x6nT的64點(diǎn)DFTX6k64=fftshift(X6k64);      %將零頻率移到頻譜中心Tp=N*T;F=1/Tp;   %頻率分辨率Fk=-N/2:N/2-1;fk=k*F;       %產(chǎn)生16點(diǎn)DFT對(duì)應(yīng)的采樣點(diǎn)頻率（以零頻率為中心）subplot(3,1,3);stem(fk,abs(X6k64),'

41、.'); box on%繪制8點(diǎn)DFT的幅頻特性圖title('(6a) 64點(diǎn)|DFTx_6(nT)|');xlabel('f(Hz)');ylabel('幅度');axis(-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k64)10.4 IIR數(shù)字濾波器設(shè)計(jì)及軟件實(shí)現(xiàn)exp4.m：%=clear all;close all;clc;Fs=10000;T=1/Fs;   %采樣頻率%調(diào)用信號(hào)產(chǎn)生函數(shù)mstg產(chǎn)生由三路抑制載波調(diào)幅信號(hào)相加構(gòu)成的復(fù)合信號(hào)stst=mstg;%低通濾波器設(shè)計(jì)與實(shí)現(xiàn)%=f

42、p=280;fs=450;wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60;   %DF指標(biāo)（低通濾波器的通、阻帶邊界頻）N,wp=ellipord(wp,ws,rp,rs); %調(diào)用ellipord計(jì)算橢圓DF階數(shù)N和通帶截止頻率wpB,A=ellip(N,rp,rs,wp);      %調(diào)用ellip計(jì)算橢圓帶通DF系統(tǒng)函數(shù)系數(shù)向量B和Ay1t=filter(B,A,st);     %濾波器軟件實(shí)現(xiàn)% 低通濾波器設(shè)計(jì)與實(shí)現(xiàn)繪圖部分figure(2);subp

43、lot(3,1,1);myplot(B,A);  %調(diào)用繪圖函數(shù)myplot繪制損耗函數(shù)曲線yt='y_1(t)'subplot(3,1,2);tplot(y1t,T,yt); %調(diào)用繪圖函數(shù)tplot繪制濾波器輸出波形%帶通濾波器設(shè)計(jì)與實(shí)現(xiàn)%=fpl=440;fpu=560;fsl=275;fsu=900;wp=2*fpl/Fs,2*fpu/Fs;ws=2*fsl/Fs,2*fsu/Fs;rp=0.1;rs=60;N,wp=ellipord(wp,ws,rp,rs);    %調(diào)用ellipord計(jì)算橢圓DF階數(shù)N和通帶截止頻率wpB,A=ellip(N,rp,rs,wp); %調(diào)用ellip計(jì)算橢圓帶通DF系統(tǒng)函數(shù)系數(shù)向量B和Ay2t=filter(B,A,st);     %濾波器