【正文】 )。title(39。Y039。)。xlabel(39。Y0=abs(fftshift(Y))。)。title(39。y39。)。xlabel(39。subplot(2,1,1)。w=Kaiser(N,beta)。N=51。)。 xlabel (39。subplot(3,3,5)。[Pxx3,F3] = pwelch(xn, window39。 grid on。利用函數(shù)psd39。plot(F2, 10*log10(Pxx2))。[Pxx2,F2] = psd(xn, 1024,Fs, window, noverlap)。 grid on。利用函數(shù)periodogram39。plot(F1, 10*log10(Pxx1))。[Pxx1,F1] = periodogram(xn,window,1024,Fs)。 grid on。利用公式39。plot([0:1023]*Fs/1024,10*log10(Pxx))。grid on。隨機信號39。plot(t,xn)。 %粗略地估計xn的功率譜，做N=1024點FFT： Pxx = abs(fft(xn,1024)).^2/1001。 % 抽樣頻率t = 0:1/Fs:1。 (1) 繪出此隨機信號的時域波形； (2) 試分別用周期圖法、平均周期圖法和Welch法分析該序列的功率譜估計。axis([0 48 1 2])。subplot(3,2,6)。stem(n2,y5)。axis([0 48 1 2])。subplot(3,2,4)。stem(n2,y3)。axis([0 48 1])。subplot(3,2,2)。stem(n2,y1)。n2=0:2*N12。y5=conv(x,h5)。y3=conv(x,h3)。y1=conv(x,h1)。axis([0 25 ])。g39。subplot(3,1,3)。)。stem(n1,real(h5),39。axis([0 25 ])。figure(4)subplot(3,1,1)。g39。r39。h6=ifft(H6,N1)。 [H6,w]=freqz(b6,wc,512)。high39。h5=ifft(H5,N1)。 [H5,w]=freqz(b5,wc,512)。high39。 h4=ifft(H4,N1)。 [H4,w]=freqz(b4,wc,512)。high39。axis([0 25 ])。g39。subplot(3,1,3)。)。stem(n1,real(h2),39。axis([0 25 ])。figure(2)subplot(3,1,1)。g39。r39。h3=ifft(H3,N1)。 [H3,w]=freqz(b3,wc,512)。h2=ifft(H2,N1)。 [H2,w]=freqz(b2,wc,512)。 h1=ifft(H1,N1)。 [H1,w]=freqz(b1,wc,512)。n1=0:N11。n=0:M。M=N1。 1) 作出各濾波器的單位脈沖響應(yīng) 2) 作出各濾波器的幅頻響應(yīng)并比較各濾波器的通帶紋波和阻帶紋波。)。title(39。Y039。)。xlabel(39。Y0=abs(fftshift(Y))。)。title(39。y39。)。xlabel(39。subplot(2,1,1)。w=Kaiser(N,beta)。凱撒窗N=51。頻譜波形39。)。ylabel(39。W39。plot([128:127],Y0)。subplot(2,1,2)。時域波形39。)。ylabel(39。w39。stem([0:N1],w)。Y=fft(w,256)。三角形窗N=51。頻譜波形39。)。ylabel(39。W39。plot([128:127],Y0)。subplot(2,1,2)。時域波形39。)。ylabel(39。w39。stem([0:N1],w)。Y=fft(w,256)。k=0:N1。)。title(39。Y039。)。xlabel(39。Y0=abs(fftshift(Y))。)。title(39。y39。)。xlabel(39。subplot(2,1,1)。w=1/2*(1cos(2*pi*k/(N1)))。漢寧窗N=51。頻譜波形39。)。ylabel(39。W39。plot([128:127],Y0)。subplot(2,1,2)。時域波形39。)。ylabel(39。w39。stem([0:N1],w)。Y=fft(w,256)。k=0:N1。)。title(39。Y039。)。xlabel(39。Y0=abs(fftshift(Y))。)。title(39。y39。)。xlabel(39。subplot(2,1,1)。w=boxcar(N)。)。ylabel(39。frequency(rad)39。grid on。plot(w,20*log10(abs(h)))。w=linspace(0,2*pi,500)。Fs=1/T。)。[b_LP,a_LP]=butter(N,norm_wc,39。s39。endnorm_wc=1。if(abs(norm_wr1)abs(norm_wr2))0 norm_wr=abs(norm_wr2)。norm_wc1=(((wc1^2)(t))/(B*wc1))。norm_wr1=(((wr1^2)(t))/(B*wr1))。B=wc2wc1。wr1=(2/T)*tan(omegar1/2)。wc1=(2/T)*tan(omegac1/2)。Ar=60。omegar2=*pi。omegac2=*pi。 axis([10,10,0,])。計算值39。理論值39。)。xlabel(39。幅度譜39。)。plot(w,abs(fftshift(X)),w,abs(y),39。y=1./(j*w+2)。)。title(39。t39。plot(t,x)。 X=T*fft(x,N)。t=0:T:Tp。 N=512。(1) 確定DFT計算的各參數(shù)（抽樣間隔，截短長度，頻譜分辨率等）； 答：抽樣間隔T=，截斷長度Tp=3，N=300,取512fsam=100。Magnitude response39。)。ylabel(39。Frequency(rad)39。 plot(w,abs(H))。)。title(39。y(n)39。)。xlabel(39。subplot(4,1,4)。x(n)=cos(pi/5.*n)39。)。ylabel(39。n39。stem(n,y2)。)。title(39。y(n)39。)。xlabel(39。subplot(4,1,2)。零輸入響應(yīng)39。)。ylabel(39。n39。stem(n,y)。y3=filter(b,a,x3,zi)。