數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告-文庫(kù)吧
2025-07-19 23:50 本頁面
【正文】 越大，計(jì)算量、存儲(chǔ)量也隨之增大。一般取N為2的整次冪，以便用FFT計(jì)算，若已給定N，可用補(bǔ)零方法便N為2的整次冪。C 和N確定后，即可確定所需相應(yīng)模擬信號(hào)的長(zhǎng)度 (3)分辨率反比于T，而不是N，在給定的T的情況下，靠減小來增加N是不能提高分辨率的，因?yàn)闉槌?shù) 2．譜分析步驟； A 數(shù)據(jù)準(zhǔn)備 (4) B 使用FFT計(jì)算信號(hào)的頻譜 (5) (6) C 由頻譜計(jì)算幅度譜、相位譜和功率譜 (7) (8) (9)3．實(shí)驗(yàn)中用到的一些基本函數(shù)簡(jiǎn)介y=fft(x,n) ； 計(jì)算n點(diǎn)的FFT。abs(x) ； 取絕對(duì)值。angle(z) ； 取相角。[Pxx, f]= periodogram (xn, nfft, fs, window) ；%周期圖譜估計(jì)[Pxx, f]=pwelch (xn, nfft, fs, window, noverlap)；%平均周期圖法 Pxx=psd (xn) ；功率譜密度 實(shí)驗(yàn)內(nèi)容 1. 已知序列x(n)=2sin()+cos() 0≤n100,試?yán)L制x(n)及它的頻譜圖。若x(n)=sin()+2cos(),結(jié)果又如何？【程序代碼】clear all。close all。clcN=100。n=0:N1。xn=2*sin(*pi*n)+cos(*pi*n)。XK=fft(xn,N)。magXK=abs(XK)。phaXK=angle(XK)。subplot(1,2,1)plot(n,xn)xlabel(39。n39。)。ylabel(39。x(n)39。)。title(39。x(n)N=10039。)。subplot(1,2,2)k=0:length(magXK)1。stem(k,magXK,39。.39。)。xlabel(39。k39。)。ylabel(39。|X(K)|39。)。title(39。X(K)N=10039。)。 2. 對(duì)下面信號(hào)進(jìn)行頻譜分析，求幅度譜和相位譜。 (1)，， (2)，【程序代碼】clear all。close all。clcfs1=5000。dt1=1/fs1。N1=。n1=0:N11。xn1=.^(n1*dt1)。Xk1=fft(xn1,N1)。mag1=2*abs(Xk1)/N1。pha1=angle(Xk1)。f1=n1*fs1/N1。subplot(221)stem(f1,mag1,39。fill39。)。title(39。Magnitude 139。)xlabel(39。f39。)。ylabel(39。|Xk(1)|39。)。subplot(222)stem(f1,pha1,39。fill39。)。title(39。Phase 139。)xlabel(39。f39。)。ylabel(39。|Phk(1)|39。)。fs2=128。dt2=1/fs2。N2=18。n2=0:N21。xn2=sin(n2*dt2)./(n2*dt2+eps)。Xk2=fft(xn2,N2)。mag2=2*abs(Xk2)/N2。pha2=angle(Xk2)。f2=n2*fs2/N2。subplot(223)stem(f2,mag2,39。fill39。)。title(39。Magnitude 239。)xlabel(39。f39。)。ylabel(39。|Xk(2)|39。)。subplot(224)stem(f2,pha2,39。fill39。)。title(39。Phase 239。)xlabel(39。f39。)。ylabel(39。|Phk(2)|39。) 3. 給定信號(hào)，，現(xiàn)在對(duì)采樣，采樣點(diǎn)數(shù)，采樣頻率=50Hz，設(shè)采樣序列為，編寫程序計(jì)算的頻譜，并繪圖；改變采樣頻率，得到序列，計(jì)算的頻譜，并繪圖；增大采樣點(diǎn)數(shù)，得到序列，計(jì)算的頻譜，并繪圖；采樣點(diǎn)數(shù)N=64，采樣頻率=300Hz，在采樣點(diǎn)后補(bǔ)零得到新序列，計(jì)算的頻譜，并繪圖。(Signal )【程序代碼】clear all。close all。 clcN=16。n=0:N1。f1=15。f2=18。fs=50。dt=1/fs。xn=sin(2*pi*f1*n*dt)+2*sin(2*pi*f2*n*dt)。Xk=fft(xn,N)。mag=2*abs(Xk)/N。pha=angle(Xk)。f=fs/N*n。figure(1)。subplot(121)plot(f,mag)。gridtitle(39。Magnitude39。)。xlabel(39。f39。)。ylabel(39。|Xk|39。)。subplot(122)plot(f,pha)。gridtitle(39。Phase39。)。xlabel(39。f39。)。ylabel(39。???161。167。Xk??39。)。N=16。n=0:N1。f1=15。f2=18。fs=36。 dt=1/fs。xn=sin(2*pi*f1*n*dt)+2*sin(2*pi*f2*n*dt)。Xk=fft(xn,N)。mag=2*abs(Xk)/N。pha=angle(Xk)。f=fs/N*n。figure(2)。subplot(221)plot(f,mag,39。r39。)。gridtitle(39。Magnitude(36Hz Sampling requency)39。)。xlabel(39。f39。)。ylabel(39。|Xk|39。)。subplot(222)plot(f,pha,39。r39。)。gridtitle(39。Phase(36Hz Sampling requency)39。)。xlabel(39。f39。)。ylabel(39。???161。167。Xk??39。)。fs=100。 dt=1/fs。xn=sin(2*pi*f1*n*dt)+2*sin(2*pi*f2*n*dt)。Xk=fft(xn,N)。mag=2*abs(Xk)/N。pha=angle(Xk)。f=fs/N*n。subplot(223)plot(f,mag,39。r39。)。gridtitle(39。Magnitude(100Hz Sampling Frequency)39。)。xlabel(39。f39。)。ylabel(39。|Xk|39。)。subplot(224)plot(f,pha,39。r39。)。gridtitle(39。Phase(100Hz Sampling Frequency)39。)。xlabel(39。f39。)。ylabel(39。???161。167。Xk??39。)。N=64。n=0:N1。f1=15。f2=18。fs=50。dt=1/fs。xn=sin(2*pi*f1*n*dt)+2*sin(2*pi*f2*n*dt)。Xk=fft(xn,N)。mag=2*abs(Xk)/N。pha=angle(Xk)。f=fs/N*n。figure(3)。subplot(221)pl
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1