科技大學數(shù)字信號處理課程設(shè)計(更新版)

　　

【正文】 a=conv([1 +5j],[1 ])。plot(w,abs(h))。t39。title(39。)。s39。title(39。)。plot(real(z),imag(z),39。sys=tf(b,a)。)。h=impulse(sys,t)。)。z=roots(b)。j5，增益 k =1 ；a=conv([1 5j],[1 5j])。plot(h)。t39。plot(real(p),imag(p),39。b=[1]*k。xlabel(39。)。*39。sys=tf(39。FFT 產(chǎn)生離散的頻譜，出現(xiàn)在FFT 每個譜線的是在每個譜線上的連續(xù)卷積頻譜。如若信號中有許多遠離被測頻率的強干擾頻率分量，應(yīng)選擇旁瓣衰減速度較快的窗函數(shù)；如果強干擾頻率分量緊鄰被測頻率時，應(yīng)選擇旁瓣峰值較小的窗函數(shù)；如果被測信號含有兩個或兩個以上的頻率成份，應(yīng)選用主瓣很窄的窗函數(shù)；如果是單一頻率信號，且要求幅度精度較高，則推薦用寬主瓣的窗函數(shù)。2. 若信號持續(xù)時間無限，且無解析表達式，如何利用DFT分析其頻譜？答：先用窗函數(shù)截取，將信號變?yōu)橛邢揲L的信號，然后利用DFT分析其頻譜。 axis([pi,pi,4,4])。ylabel(39。 N=32。legend(39。r.39。)。 T=1/fsam。X[k]39。x [n]39。 n=(N1)/2:(N1)/2。)。 t=0:T:T0。ylabel(39。t=0:T:T0。Frequency(Hz)39。 N=32。fs=39。頻率（弧度/秒）39。mag=Ts*abs(y)。)。syms t wxt=sym(39。m=(t=2)。x=(1+(t./2)).*(nm)。t=10:10。xlabel(39。xlabel(39。xlabel(39。x=m.*u。y=x*h39。h39。x39。m=x.*y。x(t)39。xlabel(39。（4）w0=2*pi。x=t(t=0)。 plot(t,x)。x=(t=0)。)。t=0::10。 x=cos(w0*t).*cos(w1*t)。（5）t=10::10。t39。h=[0,1,2,3]。Time index n39。Time index n39。Time index n39。subplot(3,1,1)。subplot(3,1,2)。subplot(3,1,3)。（5）已知信號x(t)=(1+t/2)[u(t+2)u(t2)], 求x(t+2),x(t2),x(t),x(2t),x(t)t=10:10。x=(1+(t./2)).*(nm)。m=(t=2)。plot(t,x)。Fw=fourier(xt,t,w)。fs=6。*Ws/length(y)。幅度譜39。]。x=cos((pi/8)*t+pi/4)。magnitude39。f=1/T/N*(N/2:(N/21))。 4．產(chǎn)生一個淹沒在噪聲中的信號x(t)，例如由50Hz和120Hz的正弦信號以及一個零均值的隨機噪聲疊加而成。Xm=fft(x,N)/N。)。 X=fft(x,N)。Time n39。Frequency (rad)39。 X=T*fft(x,N)。)。幅度譜39。計算值39。X=1/N*fft(x,N)。 xlabel(39。)。 柵欄現(xiàn)象：補零增加序列長度4. 在利用DFT分析連續(xù)信號頻譜時，如何選擇窗函數(shù)？ 答：加窗處理對頻譜分析主要有兩個方面的影響：頻譜中出現(xiàn)多余的高頻分量；譜線變成了具有一定寬度的譜峰。5. 討論序列后補零對頻譜分析結(jié)果的影響。如果時間序列的長度不是周期的整數(shù)倍，窗函數(shù)的連續(xù)頻譜將偏離主瓣的中心，頻率偏移量對應(yīng)著信號頻率和FFT 頻率分辨率的差異，這個偏移導致了頻譜中出現(xiàn)旁瓣，所以，窗函數(shù)的旁瓣特性直接影響著各頻譜分量向相鄰頻譜的泄漏寬度。sys=tf(b,a)。plot(real(z),imag(z),39。)。title(39。)。hold on。零極點圖39。)。sys=tf(39。*39。)。xlabel(39。k=1。subplot(1,2,1)。xlabel(39。subplot(1,2,2)。grid(5) 零點在s =0. 5，極點在s=0 .1177。p=roots(a)。o39。grid %運行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng)%沖激響應(yīng)t=0::3。沖激響應(yīng)h(t)39。)。sys=tf(39。*39。)。xlabel(39。subplot(2,2,3)。grid 離散系統(tǒng)分析1. 已知一個LTI系統(tǒng)的差分方程為： (1)初始條件y(1)=1,y(2)=2，輸入x(n)=u(n)，計算系統(tǒng)的零輸入響應(yīng)； (2)當下面三個信號分別通過系統(tǒng)，分別計算系統(tǒng)的響應(yīng)： (3)指出這是一個什么特性的系統(tǒng)。x1=cos(pi/10.*n)。xlabel(39。)。)。)。stem(n,y3)。x(n)=cos(pi/10*7.*n)39。Magnitude39。 T=1/fsam。)。r.39。legend(39。omegar1=*pi。wr2=(2/T)*tan(omegar2/2)。else norm_wr=abs(norm_wr1)。[b_BP,a_BP]=lp2bp(b_LP,a_LP,sqrt(t),B)。xlabel(39。Y=fft(w,256)。)。W39。漢明窗N=51。)。Y0=abs(fftshift(Y))。title(39。stem([0:N1],w)。時域波形39。ylabel(39。w=*cos(2*pi*k/(N1))+*cos(4*pi*k/(N1))。y39。xlabel(39。)。)。Y0=abs(fftshift(Y))。title(39。stem([0:N1],w)。時域波形39。ylabel(39。 3）若當輸入為，計算各濾波器的輸出并作出響應(yīng)波形N=10。 H1_db=20*log10(abs(H1))。 H3_db=20*log10(abs(H3))。stem(n1,real(h1))。stem(n1,real(h3),39。 H4_db=20*log10(abs(H4))。% hanning window b6=fir1(M,wc/pi,39。,w,H6_db,39。r39。x=1+2*cos(pi/4*n1)+cos(pi/2*n1)。figure(5)subplot(3,2,1)。axis([0 48 1])。stem(n2,y6)。subplot(3,3,1)。xlabel(39。 xlabel (39。 xlabel (39。plot(F3, 10*log10(Pxx3))。Y=fft(w,256)。)。W39。33