【正文】 )。 title(39。Y039。)。 xlabel(39。 Y0=abs(fftshift(Y))。)。 title(39。y39。)。 xlabel(39。 subplot(2,1,1)。 w=Kaiser(N,beta)。 0 0 . 5 1505隨機信號0 500 1000 5 0050利用公式0 500 1 0 0 5 00利用函數(shù) p e r i o d o g r a m0 500 5 0050利用函數(shù) p s d0 500 1 0 0 5 00利用函數(shù) p w e l c h 33 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])。 29 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) 作出各濾波器的單位脈沖響應 2) 作出各濾波器的幅頻響應并比較各濾波器的通帶紋波和阻帶紋波。)。 title(39。Y039。)。 xlabel(39。 Y0=abs(fftshift(Y))。)。 title(39。y39。)。 xlabel(39。 subplot(2,1,1)。 w=Kaiser(N,beta)。 0 5 10 15 20 25 30 35 40 45 5000 . 51wy時域波形 1 5 0 1 0 0 5 0 0 50 100 1500102030WY0頻譜波形 凱撒窗 N=51。頻譜波形 39。)。 26 ylabel(39。W39。 plot([128:127],Y0)。 subplot(2,1,2)。時域波形 39。)。 ylabel(39。w39。 stem([0:N1],w)。 Y=fft(w,256)。 0 5 10 15 20 25 30 35 40 45 5000 . 51wy時域波形 1 5 0 1 0 0 5 0 0 50 100 1500102030WY0頻譜波形三角形窗 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)))。 24 0 5 10 15 20 25 30 35 40 45 5000 . 51wy時域波形 1 5 0 1 0 0 5 0 0 50 100 1500102030WY0頻譜波形 漢寧窗 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。 23 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。 end norm_wc=1。 if(abs(norm_wr1)abs(norm_wr2))0 22 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。 0 5 10 15 20 25 30 35 40 45 50 0 . 500 . 5ny(n)零輸入響應0 5 10 15 20 25 30 35 40 45 50101ny(n)x ( n ) = c o s ( p i / 1 0 . * n )0 5 10 15 20 25 30 35 40 45 50101ny(n)x ( n ) = c o s ( p i / 5 . * n )0 5 10 15 20 25 30 35 40 45 50 0 . 500 . 5ny(n)x ( n ) = c o s ( p i / 1 0 * 7 . * n ) 0 0 . 5 1 1 . 5 2 2 . 5 300 . 51t時域波形 1 0 8 6 4 2 0 2 4 6 8 1000 . 51幅度譜w理論值計算值 21 (1) 確定 DFT 計算的各 參數(shù)（抽樣間隔，截短長度，頻譜分辨率等）； 答：抽樣間隔 T=，截斷長度 Tp=3， N=300,取 512 fsam=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)。零輸入響應 39。)。ylabel(39。n39。stem(n,y)。 y3=filter(b,a,x3,zi)。 y1=filter(b,a,x1,zi)。 x2=cos(pi/5.*n)。 y=filter(b,a,x,zi)。 x=zeros(N)。 a=[1,]。 N=50。)。title(39。w39。plot(w,abs(h))。h=freqs(b,a,w)。)。 title(39。t39。 plot(h)。 h=impulse(sys,t)。)。title(39。t39。)。 plot(real(z),imag(z),39。)。 plot(real(p),imag(p),39。 z=roots(b)。 sys=tf(b,a)。s39。b=[1 ]*k。 j5，增益 k =1