科技大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)(已修改)
2025-06-29 21:42 本頁(yè)面
【正文】 數(shù)字信號(hào)處理課程設(shè)計(jì)第1章 信號(hào)的時(shí)域分析 連續(xù)信號(hào)的時(shí)域分析用Matlab產(chǎn)生下列信號(hào)并與人工分析結(jié)果進(jìn)行比較：(1) r (t)= t u (t) 1t10(2) x (t)=1+cos10 t 1t1(3) x (t)=(5e t 5e 3 t ) u (t) 1t5(4) x (t)=cos(2πt) cos(20πt) 0t5(5) x (t)=sin(t)/t 10t10（1）t=1::5。x=(t=0)。plot(t,x)。axis([2,6,])。（2）t=1::1。 x=1+cos(10*t)。 plot(t,x)。 ylabel(39。x(t)39。)。xlabel(39。t39。)。（3）t=0::5。x=t(t=0)。plot(t,x)。axis([2,6,])。t=0::10。m=(t=0)。n=5*exp(t)5*exp(3*t)。x=n.*m。plot(t,x)。（4）w0=2*pi。 w1=20*pi。 t=0::5。 x=cos(w0*t).*cos(w1*t)。 plot(t,x)。 ylabel(39。x(t)39。)。xlabel(39。t39。)。（5）t=10::10。m=sin(t)。x=m./t。plot(t,x)。ylabel(39。x(t)39。)。xlabel(39。t39。)。 離散時(shí)間序列的時(shí)域分析及信號(hào)的運(yùn)算、作圖并與理論值進(jìn)行比較：（1）x(n)=2δ(n+n0) （2） x(n)=() n [sin()+cos()] n=4:4。x=().^n。y=[sin(*pi*n)+cos(*pi*n)]。m=x.*y。stem(n,m)。（3）已知LTI離散系統(tǒng)，x(n)=[1 1 1],h(n)=[0 1 2 3],求y(n)x=[1,1,1]。h=[0,1,2,3]。y=conv(x,h)。subplot(3,1,1)。stem([0:length(x)1],x)。ylabel(39。x39。)。xlabel(39。Time index n39。)。subplot(3,1,2)。stem([0:length(h)1],h)。ylabel(39。h39。)。xlabel(39。Time index n39。)。subplot(3,1,3)。stem([0:length(y)1],y)。ylabel(39。y=x*h39。)。xlabel(39。Time index n39。)；（4）已知x(t)=e –2 tu (t), y(t)=e tu (t) , 求：x(t) * y(t)t=10:10。u=(t=0)。m=exp(2*t)。n=exp(1*t)。x=m.*u。y=n.*u。h=conv(x,y)。subplot(3,1,1)。stem([0:length(x)1],x)。ylabel(39。x(n))39。)。xlabel(39。Time index n39。)。subplot(3,1,2)。stem([0:length(y)1],y)。ylabel(39。y(n)39。)。xlabel(39。Time index n39。)。subplot(3,1,3)。stem([0:length(h)1],h)。ylabel(39。h(n)=x(n)*y(n)39。)。xlabel(39。Time index n39。)。（5）已知信號(hào)x(t)=(1+t/2)[u(t+2)u(t2)], 求x(t+2),x(t2),x(t),x(2t),x(t)t=10:10。m=(t=2)。n=(t=2)。x=(1+(t./2)).*(nm)。plot(t+2,x)。t=10:10。m=(t=2)。n=(t=2)。x=(1+(t./2)).*(nm)。plot(t2,x)。t=10:10。m=(t=2)。n=(t=2)。x=(1+(t./2)).*(nm)。plot(t,x)。t=10:10。m=(t=2)。n=(t=2)。x=(1+(t./2)).*(nm)。plot(2*t,x)。t=10:10。m=(t=2)。n=(t=2)。x=(1+(t./2)).*(nm)。plot(t,x)。第2章 信號(hào)的頻域分析 利用DFT分析連續(xù)信號(hào)頻譜 1．用fourier函數(shù),理論上求下列連續(xù)時(shí)間信號(hào)的頻譜。(1). 三角脈沖信號(hào)x1(t)= t=2::2。x=tripuls(t,2,0)。plot(t,x)。syms t wxt=sym(39。(t+1)*Heaviside(t+1)2*t*Heaviside(t)+(t1)*Heaviside(t1)39。)。Fw=fourier(xt,t,w)。FFw=maple(39。convert39。,Fw,39。piecewise39。)。FFP=abs(FFw)。ezplot(FFP,[10*pi 10*pi])axis([10*pi 10*pi 0 1])(2). 單邊指數(shù)信號(hào)x2(t)=eu(t)N=256。fs=6。Ts=1/fs。t=(0:N1)*Ts。x=exp(t)。y=fft(x)。mag=Ts*abs(y)。Ws=2*pi*fs。w=(0:length(y)1)39。*Ws/length(y)。X=1./sqrt(w.^2+1)。w1=w(1:length(y)/2)。plot(w1,mag(1:length(y)/2))。xlabel(39。頻率（弧度/秒）39。)。ylabel(39。幅度譜39。)。z=[39。N=39。 num2str(N) 39。fs=39。 num2str(fs) 39。的結(jié)果39。]。title(39。exp(t)的幅度譜39。)。2．用DFT計(jì)算下列信號(hào)的頻譜：（1） T0=16。 N=32。T=T0/N。t=0:T:T0。x=cos((pi/8)*t+pi/4)。X=1/N*fft(x,N)。f=1/T/N*(N/2:(N/21))。stem(f,abs(fftshift(X)))。xlabel(39。Frequency(Hz)39。)。ylabel(39。magnitude39。)。（2）T0=。 N=20。T=T0/N。t=0:T:T0。x=cos(*pi)*t+cos(*pi*t)。X=1/N
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1