科技大學數(shù)字信號處理課程設計(已改無錯字)
2023-07-18 21:42:48 本頁面
【正文】 l(39。y(n)39。)。title(39。x(n)=cos(pi/10*7.*n)39。)。figure(2)[H,w]=freqz(b,a)。 plot(w,abs(H))。xlabel(39。Frequency(rad)39。)。ylabel(39。Magnitude39。)。 title(39。Magnitude response39。)。(1) 確定DFT計算的各參數(shù)（抽樣間隔，截短長度，頻譜分辨率等）； 答：抽樣間隔T=，截斷長度Tp=3，N=300,取512fsam=100。Tp=3。 N=512。 T=1/fsam。t=0:T:Tp。 x=exp(2*t)。 X=T*fft(x,N)。subplot(2,1,1)。plot(t,x)。xlabel(39。t39。)。title(39。時域波形39。)。w=(N/2:N/21)*(2*pi/N)*fsam。y=1./(j*w+2)。subplot(2,1,2)。plot(w,abs(fftshift(X)),w,abs(y),39。r.39。)。title(39。幅度譜39。)。xlabel(39。w39。)。legend(39。理論值39。,39。計算值39。,0)。 axis([10,10,0,])。第4章 數(shù)字濾波器設計 IIR數(shù)字濾波器設計2．設計一個數(shù)字帶通濾波器，使之滿足以下指標：（1） 通帶（2） =, =（3） =, 60dBomegac1=*pi。omegac2=*pi。omegar1=*pi。omegar2=*pi。Ap=。Ar=60。T=5*pi*10^(7)。wc1=(2/T)*tan(omegac1/2)。wc2=(2/T)*tan(omegac2/2)。wr1=(2/T)*tan(omegar1/2)。wr2=(2/T)*tan(omegar2/2)。B=wc2wc1。t=wc1*wc2。norm_wr1=(((wr1^2)(t))/(B*wr1))。norm_wr2=(((wr2^2)(t))/(B*wr2))。norm_wc1=(((wc1^2)(t))/(B*wc1))。norm_wc2=(((wc2^2)(t))/(B*wc2))。if(abs(norm_wr1)abs(norm_wr2))0 norm_wr=abs(norm_wr2)。else norm_wr=abs(norm_wr1)。endnorm_wc=1。N=buttord(norm_wc,norm_wr,Ap,Ar,39。s39。)。[b_LP,a_LP]=butter(N,norm_wc,39。s39。)。[b_BP,a_BP]=lp2bp(b_LP,a_LP,sqrt(t),B)。Fs=1/T。[b,a]=bilinear(b_BP,a_BP,Fs)。w=linspace(0,2*pi,500)。h=freqz(b,a,w)。plot(w,20*log10(abs(h)))。axis([0,2*pi,300,50])。grid on。xlabel(39。frequency(rad)39。)。ylabel(39。gain(dB)39。)。 窗函數(shù)的特性分析1．分析并繪出各種窗函數(shù)的時域特性；矩形窗N=51。w=boxcar(N)。Y=fft(w,256)。subplot(2,1,1)。stem([0:N1],w)。xlabel(39。w39。)。ylabel(39。y39。)。title(39。時域波形39。)。subplot(2,1,2)。Y0=abs(fftshift(Y))。plot([128:127],Y0)。xlabel(39。W39。)。ylabel(39。Y039。)。title(39。頻譜波形39。)。漢明窗N=51。k=0:N1。w=*cos(2*pi*k/(N1))。Y=fft(w,256)。subplot(2,1,1)。stem([0:N1],w)。xlabel(39。w39。)。ylabel(39。y39。)。title(39。時域波形39。)。subplot(2,1,2)。Y0=abs(fftshift(Y))。plot([128:127],Y0)。xlabel(39。W39。)。ylabel(39。Y039。)。title(39。頻譜波形39。)。漢寧窗N=51。k=0:N1。w=1/2*(1cos(2*pi*k/(N1)))。Y=fft(w,256)。subplot(2,1,1)。stem([0:N1],w)。xlabel(39。w39。)。ylabel(39。y39。)。title(39。時域波形39。)。subplot(2,1,2)。Y0=abs(fftshift(Y))。plot([128:127],Y0)。xlabel(39。W39。)。ylabel(39。Y039。)。title(39。頻譜波形39。)。布拉克窗N=51。k=0:N1。w=*cos(2*pi*k/(N1))+*cos(4*pi*k/(N1))。Y=fft(w,256)。subplot(2,1,1)。stem([0:N1],w)。xlabel(39。w39。)。ylabel(39。y39。)。title(39。時域波形39。)。subplot(2,1,2)。Y0=abs(fftshift(Y))。plot([128:127],Y0)。xlabel(39。W39。)。ylabel(39。Y039。)。title(39。頻譜波形39。)。三角形窗N=51。w=bartlett(N)。Y=fft(w,256)。subplot(2,1,1)。stem([0:N1],w)。xlabel(39。w39。)。ylabel(39。y39。)。title(39。時域波形39。)。subplot(2,1,2)。Y0=abs(fftshift(Y))。plot([128:127],Y0)。xlabel(39。W39。)。ylabel(39。Y039。)。title(39。頻譜波形39。)。凱撒窗N=51。beta=4。w=Kaiser(N,beta)。Y=fft(w,256)。subplot(2,1,1)。stem([0:N1],w)。xlabel(39。w39。)。ylabel(39。y39。)。title(39。時域波形39。)。subplot(2,1,2)。Y0=abs(fftshift(Y))。plot([128:127],Y0)。xlabel(39。W39。)。ylabel(39。Y039。)。titl
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1