科技大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)-免費(fèi)閱讀
【正文】 )。subplot(2,1,2)。w39。 grid on。noverlap=500。noverlap=500。window=boxcar(1001)。)。 Fs = 1000。subplot(3,2,5)。stem(n2,y2)。y4=conv(x,h4)。stem(n1,real(h6),39。stem(n1,real(h4))。 H6_db=20*log10(abs(H6))。,hamming(N+1))。% HP% rectangular window b4=fir1(M,wc/pi,39。r39。,w,H3_db,39。 H2_db=20*log10(abs(H2))。N1=2048。頻譜波形39。plot([128:127],Y0)。ylabel(39。beta=4。)。title(39。subplot(2,1,1)。Y039。)。xlabel(39。頻譜波形39。plot([128:127],Y0)。ylabel(39。k=0:N1。)。title(39。subplot(2,1,1)。)。subplot(2,1,2)。w39。gain(dB)39。h=freqz(b,a,w)。)。norm_wr2=(((wr2^2)(t))/(B*wr2))。T=5*pi*10^(7)。,0)。)。w=(N/2:N/21)*(2*pi/N)*fsam。subplot(2,1,1)。)。xlabel(39。ylabel(39。title(39。subplot(4,1,3)。n39。y(n)39。y2=filter(b,a,x2,zi)。b=[,]。)。沖激響應(yīng)h(t)39。grid %運(yùn)行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng)%沖激響應(yīng)t=0::3。o39。p=roots(a)。grid(6) 零點(diǎn)在s =0. 5，極點(diǎn)在s=0 .1177。subplot(2,2,3)。xlabel(39。)。*39。sys=tf(39。)。零極點(diǎn)圖39。hold on。)。沖激響應(yīng)h(t)39。grid %運(yùn)行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng) %沖激響應(yīng)t=0::3。o39。p=roots(a)。grid(3) 2個(gè)共軛極點(diǎn)s=177。subplot(1,2,2)。xlabel(39。subplot(1,2,1)。k=1。plot(h)。t39。plot(real(p),imag(p),39。b=[1]*k。旁瓣在主瓣的兩側(cè)以一定的間隔衰減至零。要選擇正確的窗函數(shù)，必須先估計(jì)信號(hào)的頻譜成份。1. 既然可直接由Fourier變換的定義計(jì)算連續(xù)信號(hào)的傅里葉變換，為何利用DFT分析連續(xù)信號(hào)的頻譜？ 答：因?yàn)橛邢揲L(zhǎng)序列的DFT分析的是有限長(zhǎng)序列，特別適合數(shù)字系統(tǒng)，且存在著快速算法，所以仍用DFT分析連續(xù)信號(hào)的頻譜。 stem(omega,angle(fftshift(X)))。 axis([pi,pi,0,1])。要求：(1) 確定DFT計(jì)算的各參數(shù)； (2) 進(jìn)行理論值與計(jì)算值比較，分析各信號(hào)頻譜分析的計(jì)算精度； (3) 詳細(xì)列出利用DFT分析離散信號(hào)頻譜的步驟； (4) 寫出實(shí)驗(yàn)原理。)。plot(w,abs(fftshift(X)),w,abs(y),39。t39。 N=512。ylabel(39。 ylabel(39。 利用DFT分析離散序列頻譜1． DFT計(jì)算序列的頻譜；N=51。f (Hz)39。 T=T0/N。)。T=T0/N。xlabel(39。2．用DFT計(jì)算下列信號(hào)的頻譜：（1） T0=16。 num2str(N) 39。xlabel(39。y=fft(x)。piecewise39。plot(t,x)。t=10:10。n=(t=2)。plot(t+2,x)。)。)。)。n=exp(1*t)。ylabel(39。ylabel(39。ylabel(39。y=[sin(*pi*n)+cos(*pi*n)]。ylabel(39。)。plot(t,x)。（3）t=0::5。 x=1+cos(10*t)。plot(t,x)。xlabel(39。m=(t=0)。 plot(t,x)。m=sin(t)。)。y=conv(x,h)。)。)。)；（4）已知x(t)=e –2 tu (t), y(t)=e tu (t) , 求：x(t) * y(t)t=10:10。stem([0:length(x)1],x)。stem([0:length(y)1],y)。stem([0:length(h)1],h)。m=(t=2)。plot(t2,x)。n=(t=2)。第2章 信號(hào)的頻域分析 利用DFT分析連續(xù)信號(hào)頻譜 1．用fourier函數(shù),理論上求下列連續(xù)時(shí)間信號(hào)的頻譜。FFw=maple(39。Ts=1/fs。X=1./sqrt(w.^2+1)。)。title(39。X=1/N*fft(x,N)。)。stem(f,abs(fftshift(X)))。確定分析長(zhǎng)度和取樣速度，計(jì)算信號(hào)的頻譜；計(jì)算其功率譜密度并作圖，指出50Hz和120Hz的正弦成分以及噪聲；詳細(xì)列出檢測(cè)信號(hào)的步驟和原理。 f=((N1)/2:(N1)/2)/N/T。 title(39。 omega = 2 * pi / N * n。)。)。subplot(2,1,1)。w=(N/2:N/21)*(2*pi/N)*fsam。)。,0)。omega=2*pi/N*(nN/2)。Frequency (rad)39。 xlabel(39。對(duì)于同一個(gè)窗函數(shù)，增加長(zhǎng)度N雖然可以減少主瓣寬度，但是不減少旁瓣泄露。 答：通過(guò)補(bǔ)零增加序列長(zhǎng)度時(shí)，同時(shí)也減少了頻譜分析的譜線間隔，從而計(jì)算出的頻譜將會(huì)顯示更多的細(xì)節(jié)。選取適當(dāng)?shù)拇昂瘮?shù)能夠提高頻譜分辨率。p=roots(a)。o39。grid %運(yùn)行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng) %沖激響應(yīng)t=0::3。沖激響應(yīng)h(t)39。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)。j5，增益k =1； a=conv([1 +5j],[1 ])。z=roots(b)。)。h=impulse(sys,t)。)。title(39。s39。)。title(39。t39。plot(w,abs(h))。 N=50。x2=cos(pi/5.*n)。n39。subplot(4,1,2)。title(39。ylabel(39。xlabel(39。)。)。t=0:T:Tp。title(39。)。理論值39。omeg
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1