實(shí)驗(yàn)一：利用DFT研究信號的幅頻特性班級19通信一班姓名潘佳隆學(xué)號201909522114【實(shí)驗(yàn)?zāi)康摹?、了解不同點(diǎn)數(shù)DFT的區(qū)別與聯(lián)系?！緦?shí)驗(yàn)內(nèi)容及結(jié)果】《數(shù)字信號處理》教材87-88頁，運(yùn)行例3.1.2中的程序。（1）Matlab程序清單：%DFT的MATLB計(jì)算xn=[1111];%輸入時域序列向量xn=R8(n)Xk16=fft(xn,16);%計(jì)算xn的16點(diǎn)DFTXk32=fft(xn,32);%計(jì)算xn的32點(diǎn)DFT%以下為繪圖部分k=0:15;wk=2*k/16;%產(chǎn)生16點(diǎn)DFT對應(yīng)的采樣點(diǎn)頻率(關(guān)于π歸一化值)subplot(3,2,1);stem(wk,abs(Xk16),'.');%繪制16點(diǎn)DFT的幅頻特性圖title('(a)16點(diǎn)DFT的幅頻特性圖潘佳隆20190952214');xlabel('ω/π');ylabel('幅度')subplot(3,2,5);stem(wk,angle(Xk16),'.');%繪制16點(diǎn)DFT的相頻特性圖line([0,2],[0,0]);title('(b)16點(diǎn)DFT的相頻特性圖')xlabel('ω/π');ylabel('相位');axis([0,2,-3.5,3.5])k=0:31;wk=2*k/32;%產(chǎn)生32點(diǎn)DFT對應(yīng)的采樣點(diǎn)頻率(關(guān)于π歸一化值)subplot(3,2,2);stem(wk,abs(Xk32),'.');%繪制32點(diǎn)DFT的幅頻特性圖title('(c)32點(diǎn)DFT的幅頻特性圖');xlabel('ω/π');ylabel('幅度')subplot(3,2,6);stem(wk,angle(Xk32),'.');%繪制32點(diǎn)DFT的相頻特性圖line([0,2],[0,0]);title('(d)32點(diǎn)DFT的相頻特性圖');xlabel('ω/π');ylabel('相位');axis([0,2,-3.5,3.5])運(yùn)行結(jié)果（Matlab圖形）：【實(shí)驗(yàn)成績】教師批閱：年月日實(shí)驗(yàn)二：巴特沃斯低通模擬濾波器的設(shè)計(jì)班級19通信一班姓名潘佳隆學(xué)號201909522114【實(shí)驗(yàn)?zāi)康摹?、了解巴特沃斯低通模擬濾波器的設(shè)計(jì)過程?！緦?shí)驗(yàn)內(nèi)容及結(jié)果】《數(shù)字信號處理》教材176-177頁，運(yùn)行例6.2.2中的程序。（1）Matlab程序清單：wp=2*pi*5000;ws=2*pi*12000;Rp=2;As=30;%設(shè)置濾波器參數(shù)[N,wc]=buttord(wp,ws,Rp,As,'s');%計(jì)算濾波器階數(shù)N和3dB截止頻率[B,A]=butter(N,wc,'s');%計(jì)算濾波器系統(tǒng)函數(shù)分子分母多項(xiàng)式系數(shù)k=0:511;fk=0:14000/512:14000;wk=2*pi*fk;Hk=freqs(B,A,wk);subplot(2,2,1);plot(fk/1000,20*log10(abs(Hk)));gridonxlabel('頻率(kHz)');ylabel('幅度(dB)')axis([0,14,-40,5])title('潘佳隆201909522114')（2）運(yùn)行結(jié)果（Matlab圖形）：【實(shí)驗(yàn)成績】教師批閱：年月日實(shí)驗(yàn)三：IIR數(shù)字濾波器的設(shè)計(jì)班級19通信一班姓名潘佳隆學(xué)號201909522114【實(shí)驗(yàn)?zāi)康摹?、掌握脈沖響應(yīng)不變法設(shè)計(jì)IIR數(shù)字濾波器?！緦?shí)驗(yàn)內(nèi)容及結(jié)果】《數(shù)字信號處理》教材199-200頁，運(yùn)行例6.3.2中的程序。（1）Matlab程序清單：%T=1sT=1;Fs=1/T;wp=0.2*pi/T;ws=0.35*pi/T;rp=1;rs=10;[N1,wc1]=buttord(wp,ws,rp,rs,'s');[B1,A1]=butter(N1,wc1,'s');%繪制帶通濾波器的損耗函數(shù)曲線fk=0:Fs/2/512:Fs/2;wk=2*pi*fk;Hk=freqs(B1,A1,wk);subplot(2,2,1);plot(fk,20*log10(abs(Hk)));gridonxlabel('f(Hz)');ylabel('幅度(dB)')axis([0,0.5,-50,5]);title('(a)模擬濾波器（T=1s）潘佳隆201909522114');[Bz1,Az1]=impinvar(B1,A1,Fs);Hz=freqz(Bz1,Az1,wk);subplot(2,2,3);plot(wk/pi,20*log10(abs(Hz)));gridon;xlabel('\omega/\pi');ylabel('幅度(dB)')axis([0,1,-50,5]);title('(b)數(shù)字濾波器（T=1s）');%T=0.1sT=0.1;Fs=1/T;wp2=wp/T;ws2=ws/T;rp=1;rs=10;[N2,wc2]=buttord(wp2,ws2,rp,rs,'s');[B2,A2]=butter(N2,wc2,'s');%繪制帶通濾波器的損耗函數(shù)曲線fk=0:Fs/2/512:Fs/2;wk=2*pi*fk;Hk=freqs(B2,A2,wk);subplot(2,2,2);plot(fk,20*log10(abs(Hk)));gridonxlabel('f(Hz)');ylabel('幅度(dB)')axis([0,5,-50,5]);title('(c)模擬濾波器（T=0.1s）');[Bz2,Az2]=impinvar(B2,A2,Fs);Hz=freqz(Bz2,Az2,wk);subplot(2,2,4);plot(wk/pi,20*log10(abs(Hz)));gridon;xlabel('\omega/\pi');ylabel('幅度(dB)')axis([0,1,-50,5]);title('(d)數(shù)字濾波器（T=0.1s）');（2）運(yùn)行結(jié)果（Matlab圖形）：【實(shí)驗(yàn)成績】教師批閱：年月日實(shí)驗(yàn)四：FIR數(shù)字濾波器的設(shè)計(jì)班級19通信一班姓名潘佳隆學(xué)號201909522114【實(shí)驗(yàn)?zāi)康摹?、了解用Kaiser窗函數(shù)法設(shè)計(jì)FIR數(shù)字濾波器。【實(shí)驗(yàn)內(nèi)容及結(jié)果】《數(shù)字信號處理》教材233-234頁，運(yùn)行例7.2.2中的程序。（1）Matlab程序清單：%用kaiser窗函數(shù)設(shè)計(jì)線性相位低通FIRDFfp=1500;fs=2500;rs=40;Fs=10000;wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;Bt=ws-wp;%計(jì)算過渡帶寬度alph=0.5842*(rs-21)^0.4+0.07886*(rs-21);%根據(jù)(7.2.16)式計(jì)算kaiser窗的控制參數(shù)αN=ceil((rs-8)/2.285/Bt);%根據(jù)(7.2.17)式計(jì)算kaiser窗所需階數(shù)Nwc=(wp+ws)/2/pi;%計(jì)算理想高通濾波器通帶截止頻率(關(guān)于π歸一化）hn=fir1(N,wc,kaiser(N+1,alph));%調(diào)用kaiser計(jì)算低通FIRDF的h(n)%以下是繪圖部分M=1024;hk=fft(hn,M);n=0:N;subplot(2,2,1);stem(n,hn,'.');line([0,30],[0,0])xlabel('n');ylabel('h(n)');k=1:M/2;w=2*(0:M/2-1)/M;su