科技大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)(存儲(chǔ)版)
【正文】 ar2=*pi。B=wc2wc1。endnorm_wc=1。Fs=1/T。frequency(rad)39。subplot(2,1,1)。title(39。)。k=0:N1。ylabel(39。plot([128:127],Y0)。頻譜波形39。xlabel(39。)。Y039。Y=fft(w,256)。)。W39。三角形窗N=51。ylabel(39。plot([128:127],Y0)。頻譜波形39。xlabel(39。)。Y039。M=N1。 h1=ifft(H1,N1)。h3=ifft(H3,N1)。axis([0 25 ])。g39。 h4=ifft(H4,N1)。high39。g39。)。y1=conv(x,h1)。stem(n2,y1)。subplot(3,2,4)。axis([0 48 1 2])。plot(t,xn)。利用公式39。利用函數(shù)periodogram39。利用函數(shù)psd39。 xlabel (39。subplot(2,1,1)。title(39。)。)。xlabel(39。y39。w=Kaiser(N,beta)。subplot(3,3,5)。plot(F2, 10*log10(Pxx2))。plot(F1, 10*log10(Pxx1))。plot([0:1023]*Fs/1024,10*log10(Pxx))。 %粗略地估計(jì)xn的功率譜，做N=1024點(diǎn)FFT： Pxx = abs(fft(xn,1024)).^2/1001。subplot(3,2,6)。stem(n2,y3)。n2=0:2*N12。axis([0 25 ])。stem(n1,real(h5),39。r39。h5=ifft(H5,N1)。 [H4,w]=freqz(b4,wc,512)。subplot(3,1,3)。figure(2)subplot(3,1,1)。 [H3,w]=freqz(b3,wc,512)。 [H1,w]=freqz(b1,wc,512)。 1) 作出各濾波器的單位脈沖響應(yīng) 2) 作出各濾波器的幅頻響應(yīng)并比較各濾波器的通帶紋波和阻帶紋波。)。title(39。subplot(2,1,1)。)。subplot(2,1,2)。w39。頻譜波形39。plot([128:127],Y0)。ylabel(39。k=0:N1。)。title(39。subplot(2,1,1)。)。subplot(2,1,2)。w39。)。xlabel(39。y39。w=boxcar(N)。grid on。)。if(abs(norm_wr1)abs(norm_wr2))0 norm_wr=abs(norm_wr2)。wr1=(2/T)*tan(omegar1/2)。omegac2=*pi。)。plot(w,abs(fftshift(X)),w,abs(y),39。t39。 N=512。ylabel(39。title(39。subplot(4,1,4)。n39。y(n)39。零輸入響應(yīng)39。stem(n,y)。y=filter(b,a,x,zi)。)。h=freqs(b,a,w)。plot(h)。t39。plot(real(p),imag(p),39。b=[1 ]*k。w39。title(39。)。plot(real(z),imag(z),39。sys=tf(b,a)。)。h=impulse(sys,t)。)。z=roots(b)。 j5，增益k =1； a=conv([1 +5j],[1 ])。plot(h)。t39。plot(real(p),imag(p),39。b=[1]*k。t39。title(39。)。s39。)。零極點(diǎn)圖39。hold on。)。沒(méi)有出現(xiàn)旁瓣的原因是旁瓣正處在窗函數(shù)主瓣兩側(cè)采樣頻率間隔處的零分量點(diǎn)。 絕大多數(shù)應(yīng)用采用漢寧（Hanning）窗即可得到滿意的結(jié)果，因?yàn)樗哂休^好的頻率分辨率和抑制頻譜泄漏的能力。 頻率泄漏：選擇適合的窗函數(shù)。Phase39。)。x=cos(3*pi/8*n)。,39。title(39。時(shí)域波形39。 x=exp(2*t)。xlabel(39。 xlabel(39。 x=().^n.*(n=0)。magnitude39。 。)。X=1/N*fft(x,N)。ylabel(39。t=0:T:T0。的結(jié)果39。ylabel(39。w=(0:length(y)1)39。ezplot(FFP,[10*pi 10*pi])axis([10*pi 10*pi 0 1])(2). 單邊指數(shù)信號(hào)x2(t)=eu(t)N=256。)。x=(1+(t./2)).*(nm)。t=10:10。n=(t=2)。)。)。)。h=conv(x,y)。xlabel(39。xlabel(39。xlabel(39。（3）已知LTI離散系統(tǒng)，x(n)=[1 1 1],h(n)=[0 1 2 3],求y(n)x=[1,1,1]。xlabel(39。)。 t=0::5。axis([2,6,])。x(t)39。數(shù)字信號(hào)處理課程設(shè)計(jì)第1章 信號(hào)的時(shí)域分析 連續(xù)信號(hào)的時(shí)域分析用Matlab產(chǎn)生下列信號(hào)并與人工分析結(jié)果進(jìn)行比較：(1) r (t)= t u (t) 1t10(2) x (t)=1+cos10 t 1t1(3) x (t)=(5e t 5e 3 t ) u (t) 1t5(4) x (t)=cos(2πt) cos(20πt) 0t5(5) x (t)=sin(t)/t 10t10（1）t=1::5。 ylabel(39。plot(t,x)。 w1=20*pi。t39。)。stem(n,m)。)。)。)。y=n.*u。Time index n39。Time index n39。Time index n39。m=(t=2)。plot(t,x)。n=(t=2)。(t+1)*Heaviside(t+1)2*t*Heaviside(t)+(t1)*Heaviside(t1)39。FFP=abs(FFw)。Ws=2*pi*fs。)。 num2str(fs) 39。T=T0/N。)。x=cos(*pi)*t+cos(*pi*t)。magnitude39。x=sin(100*pi*t)+sin(240*pi*t)+ randn(size(t))。ylabel(39。 %若N為偶，n = N / 2: (N / 2 1)。)。)。t=0:T:Tp。title(39。)。理論值39。 n=0:N1。Magnitude39。ylabel(39。3. 在利用DFT分析連續(xù)信號(hào)頻譜時(shí)，會(huì)出現(xiàn)哪些誤差？如何克服或改善這些誤差？ 答：混疊現(xiàn)象：減少時(shí)域抽樣間隔可以克服這種
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1