科技大學(xué)數(shù)字信號處理課程設(shè)計-文庫吧在線文庫

　　

【正文】 誤差。對頻帶較寬或含有多個頻率成份的信號則采用連續(xù)采樣。如果原始信號的頻譜成份與FFT 中的譜線完全一致，這種情況下采樣數(shù)據(jù)的長度為信號周期的整數(shù)倍，頻譜中只有主瓣。s39。)。title(39。t39。sys=tf(39。*39。)。xlabel(39。k=1。subplot(1,2,1)。xlabel(39。subplot(1,2,2)。grid(4) 2個共軛極點s=0 .1177。p=roots(a)。o39。grid %運行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng) %沖激響應(yīng)t=0::3。沖激響應(yīng)h(t)39。)。hold on。零極點圖39。)。xlabel(39。k=1。subplot(2,2,1)。xlabel(39。subplot(2,2,2)。grid%幅頻特性w=0::5。幅頻特性39。 zi=filtic(b,a,[1,2])。figure(1)subplot(4,1,1)。title(39。ylabel(39。xlabel(39。)。)。)。Tp=3。xlabel(39。subplot(2,1,2)。w39。第4章 數(shù)字濾波器設(shè)計 IIR數(shù)字濾波器設(shè)計2．設(shè)計一個數(shù)字帶通濾波器，使之滿足以下指標：（1） 通帶（2） =, =（3） =, 60dBomegac1=*pi。wc2=(2/T)*tan(omegac2/2)。norm_wc2=(((wc2^2)(t))/(B*wc2))。s39。axis([0,2*pi,300,50])。 窗函數(shù)的特性分析1．分析并繪出各種窗函數(shù)的時域特性；矩形窗N=51。ylabel(39。plot([128:127],Y0)。頻譜波形39。xlabel(39。)。Y039。Y=fft(w,256)。)。W39。布拉克窗N=51。)。Y0=abs(fftshift(Y))。title(39。xlabel(39。)。Y039。Y=fft(w,256)。)。W39。 FIR數(shù)字濾波器設(shè)計2. 分別使用矩形窗、哈明窗、海寧窗設(shè)計一個N=10的FIR低通和高通濾波器， 截頻為。% LP% rectangular window b1=fir1(M,wc/pi,boxcar(N))。% hanning window b3=fir1(M,wc/pi,hanning(N))。)。axis([0 25 ])。,boxcar(N+1))。 H5_db=20*log10(abs(H5))。figure(3)c=plot(w,H4_db,w,H5_db,39。subplot(3,1,2)。)。y6=conv(x,h6)。subplot(3,2,3)。axis([0 48 1 2])。 % 抽樣時間xn = cos(20*pi*t) + cos(40*pi*t) + randn(size(t))。subplot(3,3,2)。subplot(3,3,3)。subplot(3,3,4)。, noverlap, 1024,Fs)。beta=4。ylabel(39。plot([128:127],Y0)。頻譜波形39。ylabel(39。時域波形39。stem([0:N1],w)。利用函數(shù)pwelch39。)。)。)。 xlabel(39。 第5章 隨機信號功率譜估計1．某隨機信號由兩余弦信號與噪聲構(gòu)成：x[k]= cos (20πk) +cos (40πk) + s [k]，s[k]為均值為0，方差為1的高斯白噪聲。stem(n2,y4)。axis([0 48 1])。y2=conv(x,h2)。axis([0 25 ])。)。,hanning(N+1))。% hamming window b5=fir1(M,wc/pi,39。)。subplot(3,1,2)。figure(1)c=plot(w,H1_db,w,H2_db,39。% hamming window b2=fir1(M,wc/pi,hamming(N))。wc=pi/3。)。subplot(2,1,2)。w39。)。xlabel(39。y39。w=bartlett(N)。)。title(39。subplot(2,1,1)。)。subplot(2,1,2)。w39。)。xlabel(39。y39。w=*cos(2*pi*k/(N1))。ylabel(39。時域波形39。stem([0:N1],w)。)。[b,a]=bilinear(b_BP,a_BP,Fs)。N=buttord(norm_wc,norm_wr,Ap,Ar,39。t=wc1*wc2。Ap=。,39。title(39。時域波形39。 x=exp(2*t)。 title(39。figure(2)[H,w]=freqz(b,a)。n39。y(n)39。x(n)=cos(pi/10.*n)39。stem(n,y1)。)。x3=cos(pi/10*7.*n)。n=0:N1。xlabel(39。)。零極點圖39。hold on。)。幅頻特性39。grid%幅頻特性w=0::5。subplot(2,2,2)。xlabel(39。subplot(2,2,1)。k=1。xlabel(39。)。*39。sys=tf(39。)。零極點圖39。hold on。)。沖激響應(yīng)h(t)39。grid %運行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng) %沖激響應(yīng)t=0::3。o39。p=roots(a)。)。h=impulse(sys,t)。)。z=roots(b)。第3章 系統(tǒng)分析與設(shè)計 連續(xù)系統(tǒng)分析2． 研究具有以下零極點連續(xù)系統(tǒng)： (1) 1個極點s=，增益 k =1；a=[1,]。6. 窗函數(shù)對頻譜分辨率有何影響？如何提高頻譜分辨率？ 答：輸入數(shù)據(jù)通過一個窗函數(shù)相當(dāng)于原始數(shù)據(jù)的頻譜與窗函數(shù)頻譜的卷積。為了在兩者之間取得平衡，應(yīng)根據(jù)實際信號特性采用合適的窗函數(shù)。Frequency (rad) 39。 )。subplot(2,1,1)。 axis([10,10,0,])。xlabel(39。y=1./(j*w+2)。plot(t,x)。；確定DFT計算的各參數(shù)（抽樣間隔，截短長度，頻譜分辨率等）； 答：抽樣間隔T=，截斷長度Tp=3，N=300,取512fsam=100。 subplot(2,1,2)。 subplot(2,1,1)。幅度譜39。 stem(f,abs(fftshift(Xm)))。 50赫茲的頻率成分對應(yīng)50和50兩個坐標點，120赫茲的頻率成分對應(yīng)120和120兩個坐標點的頻率T0=1。xlabe