資源描述:
《數(shù)字信號(hào)處理西電》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、-數(shù)字信號(hào)處理上機(jī)第一次實(shí)驗(yàn)實(shí)驗(yàn)一:設(shè)給定模擬信號(hào)xate1000t,t的單位是ms。(1)利用MATLAB繪制出其時(shí)域波形和頻譜圖(傅里葉變換),估計(jì)其等效帶寬(忽略譜分量降低到峰值的3%以下的頻譜)。(2)用兩個(gè)不同的采樣頻率對(duì)給定的xat進(jìn)行采樣。15000樣本秒采樣xat得到x1n。○以fs畫出x1n及其頻譜X1ej。21000樣本秒采樣xat得到x1n?!鹨詅s畫出x1n及其頻譜X1ej。比較兩種采樣率下的信號(hào)頻譜,并解釋。實(shí)驗(yàn)一MATLAB程序:(1)○1clc;fs=5000;ts=1/fs;N=1000;t=(-N:N)*ts;s
2、=exp(-abs(t));plot(t,s,'linewidth',1.5)xlabel('時(shí)間')ylabel('幅度')set(gca,'fontweight','b','fontsize',12)SPL=N*100;figuresp=fftshift(fft(s,SPL));sp=sp/max(sp)*100;freqb=-fs/2:fs/SPL:fs/2-fs/SPL;plot(freqb,abs(sp))xlabel('頻率')ylabel('頻譜幅度')--set(gca,'fontweight','b','fontsize',12
3、)yy=abs(abs(sp)-3);[aa,freqind]=min(yy);(freqind-SPL/2)*fs/SPL--○2--clc;--fs=1000;ts=1/fs;N=1000;t=(-N:N)*ts;s=exp(-abs(t));plot(t,s,'linewidth',1.5)xlabel('時(shí)間')ylabel('幅度')set(gca,'fontweight','b','fontsize',12)SPL=N*100;figuresp=fftshift(fft(s,SPL));sp=sp/max(sp)*100;freqb=
4、-fs/2:fs/SPL:fs/2-fs/SPL;plot(freqb,abs(sp))xlabel('頻率')ylabel('頻譜幅度')set(gca,'fontweight','b','fontsize',12)yy=abs(abs(sp)-3);[aa,freqind]=min(yy);(freqind-SPL/2)*fs/SPL--實(shí)驗(yàn)三:設(shè)xn1,2,2,x2n1,2,3,4,編寫MATLAB程序,計(jì)算:1(1)5點(diǎn)圓周卷積y1n;(2)6點(diǎn)圓周卷積y2n;(3)線性卷積y3n;(4)畫出的y1n,y2n和y3n時(shí)間軸對(duì)齊。a=[1,
5、2,2];b=[1,2,3,4];y1=cconv(a,b,5)y2=cconv(a,b,6)y3=conv(a,b)figure(1);subplot(311)stem(y1);gridontitle('五點(diǎn)圓周卷積y1(n)');--xlabel('n'),ylabel('y1(n)');axis([06015])subplot(312)--stem(y2);gridontitle('六點(diǎn)圓周卷積y2(n)');xlabel('n'),ylabel('y2(n)');axis([06015])subplot(313)stem(y3);grid
6、ontitle('線性卷積y3(n)');xlabel('n'),ylabel('y3(n)');axis([06015])--實(shí)驗(yàn)四:給定因果系統(tǒng):yn0.9yn1xn(1)求系統(tǒng)函數(shù)Hz并畫出零極點(diǎn)示意圖。(2)畫出系統(tǒng)的幅頻特性Hej和相頻特性。(3)求脈沖響應(yīng)hn并畫序列圖。--提示:在MATLAB中,zplane(b,a)函數(shù)可畫零極點(diǎn)圖;Freqz(b,a,N)可給出0,范圍內(nèi)--均勻間隔的N點(diǎn)頻率響應(yīng)的復(fù)振幅;Impz(b,a,N)可求Hz的逆變換(即脈沖響應(yīng))。--clca=[1,0]b=[1,-0.9]figu
7、re(1)zplane(b,a);title('零極點(diǎn)分布圖')w=[-3*pi:0.01:3*pi];[h,phi]=freqz(b,a,w);figure(2);subplot(2,1,1);plot(w,abs(h));gridon;title('幅頻特性');xlabel('f/Hz'),ylabel('H(f)');subplot(2,1,2);plot(w,phi);gridon;title('相頻特性');xlabel('f/Hz'),ylabel('W(f)');----數(shù)字信號(hào)處理第二次實(shí)驗(yàn)1.給定模擬信號(hào)xt2sin4t5co
8、s8t,對(duì)其進(jìn)行采樣,用DFT(FFT)進(jìn)行信號(hào)頻譜分析。(1)確定最小采樣頻率和最小采樣點(diǎn)數(shù)。(2)若以t0.01nn0:N1秒進(jìn)行采