科技大學數字信號處理課程設計(完整版)

　　

【正文】 l(39。（2）T0=。f=1/T/N*(N/2:(N/21))。exp(t)的幅度譜39。z=[39。w1=w(1:length(y)/2)。t=(0:N1)*Ts。convert39。(1). 三角脈沖信號x1(t)= t=2::2。x=(1+(t./2)).*(nm)。t=10:10。n=(t=2)。ylabel(39。ylabel(39。ylabel(39。u=(t=0)。subplot(3,1,3)。subplot(3,1,2)。subplot(3,1,1)。 離散時間序列的時域分析及信號的運算、作圖并與理論值進行比較：（1）x(n)=2δ(n+n0) （2） x(n)=() n [sin()+cos()] n=4:4。x=m./t。 ylabel(39。n=5*exp(t)5*exp(3*t)。t39。axis([2,6,])。（2）t=1::1。)。x=n.*m。x(t)39。plot(t,x)。x=().^n。stem([0:length(x)1],x)。stem([0:length(h)1],h)。stem([0:length(y)1],y)。m=exp(2*t)。x(n))39。y(n)39。h(n)=x(n)*y(n)39。x=(1+(t./2)).*(nm)。m=(t=2)。plot(2*t,x)。x=tripuls(t,2,0)。,Fw,39。x=exp(t)。plot(w1,mag(1:length(y)/2))。N=39。)。stem(f,abs(fftshift(X)))。 N=20。Frequency(Hz)39。 N=241。 xlabel(39。)。 stem(n,x)。stem(omega,real(fftshift(X)))。Tp=3。xlabel(39。subplot(2,1,2)。w39。3． 有限長序列， 0≤n≤31，分別用N=32，N=60，N=120點DFT計算其頻譜。 stem(omega,abs(fftshift(X)))。subplot(2,1,2)。 )。每種窗函數有其自身的特性，不同的窗函數適用于不同的應用。窗函數的頻譜由一個主瓣和幾個旁瓣組成，主瓣以時域信號的每個頻率成份為中心。k=1。subplot(1,2,1)。xlabel(39。subplot(1,2,2)。Grid(2) 1個極點s =0，增益k =1； a=[1,0]。z=roots(b)。)。h=impulse(sys,t)。)。sys=tf(b,a)。plot(real(z),imag(z),39。)。title(39。s39。)。title(39。t39。b=[1 ]*k。plot(real(p),imag(p),39。t39。plot(h)。h=freqs(b,a,w)。)。sys=tf(b,a)。plot(real(z),imag(z),39。)。title(39。w39。a=[1,]。y1=filter(b,a,x1,zi)。ylabel(39。xlabel(39。)。)。)。 plot(w,abs(H))。Magnitude response39。 X=T*fft(x,N)。)。幅度譜39。計算值39。Ar=60。norm_wr1=(((wr1^2)(t))/(B*wr1))。s39。w=linspace(0,2*pi,500)。ylabel(39。xlabel(39。)。Y039。Y=fft(w,256)。)。W39。漢寧窗N=51。)。Y0=abs(fftshift(Y))。title(39。stem([0:N1],w)。時域波形39。ylabel(39。Y=fft(w,256)。)。W39。凱撒窗N=51。)。Y0=abs(fftshift(Y))。title(39。n=0:M。 [H2,w]=freqz(b2,wc,512)。r39。stem(n1,real(h2),39。axis([0 25 ])。high39。 [H6,w]=freqz(b6,wc,512)。figure(4)subplot(3,1,1)。subplot(3,1,3)。y3=conv(x,h3)。subplot(3,2,2)。axis([0 48 1 2])。 (1) 繪出此隨機信號的時域波形； (2) 試分別用周期圖法、平均周期圖法和Welch法分析該序列的功率譜估計。隨機信號39。 grid on。 grid on。 grid on。)。xlabel(39。)。Y039。title(39。Y0=abs(fftshift(Y))。)。N=51。[Pxx3,F3] = pwelch(xn, window39。[Pxx2,F2] = psd(xn, 1024,Fs, window, noverlap)。[Pxx1,F1] = periodogram(xn,window,1024,Fs)。grid on。 % 抽樣頻率t = 0:1/Fs:1。stem(n2,y5)。axis([0 48 1])。y5=conv(x,h5)。g39。axis([0 25 ])。h6=ifft(H6,N1)。 [H5,w]=freqz(b5,wc,512)。high39。)。g39。h2=ifft(H2,N1)。n1=0:N11。)。xlabel(39。y39。w=Kaiser(N,beta)。ylabel(39。時域波形39。stem([0:N1],w)。)。subplot(2,1,2)。w39。)。xlabel(39。y39。w=1/2*(1cos(2*pi*k/(N1)))。ylabel(39。時域波形39。stem([0:N1],w)。title(39。Y0=abs(fftshift(Y))。)。)。plot(w,20*log10(abs(h)))。[b_LP,a_LP]=butter(N,norm_wc,39。norm_wc1=(((wc1^2)(t))/(B*wc1))。wc1=(2/T)*tan(omegac1/2)。 axis([10,10,0,])。xlabel(39。y=1./(j*w+2)。plot(t,x)。(1) 確定DFT計算的各參數（抽樣間隔，截短長度，頻譜分辨率等）； 答：抽樣間隔T=，截斷長度Tp=3，N=300,取512fsam=100。Frequency(rad)39。y(n)39。x(n)=cos(pi/5.*n)39。stem(n,y2)。)。)。y3=filter(b,a,x3,zi)。x=zeros(N)。title(39。)。h=impulse(sys,t)。)。z=roots(b)。j5，增益k =1；