科技大学数号处理课程设计.doc

数字信号处理课程设计 第1章 信号的时域分析 1.1 连续信号的时域分析 用Matlab产生下列信号并与人工分析结果进行比较： (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:0.01:5; x=(t=0); plot(t,x); axis([-2,6,-0.1,1.1]); （2）t=-1:0.001:1; x=1+cos(10*t); plot(t,x); ylabel(x(t));xlabel(t); （3） t=0:0.001:5; x=t(t=0); plot(t,x); axis([-2,6,-0.1,1.1]); t=0:0.1: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:0.001:5; x=cos(w0*t).*cos(w1*t); plot(t,x); ylabel(x(t));xlabel(t); （5） t=-10:0.1:10; m=sin(t); x=m./t; plot(t,x); ylabel(x(t));xlabel(t); 1.2 离散时间序列的时域分析及信号的运算 1.使用Matlab产生下列序列、作图并与理论值进行比较： （1）x(n)=2δ(n+n0) （2） x(n)=(0.9) n [sin(0.25πn)+cos(0.25πn)] n=-4:4; x=(0.9).^n; y=[sin(0.25*pi*n)+cos(0.25*pi*n)]; m=x.*y; stem(n,m); （3）已知LTI离散系统，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(x);xlabel(Time index n); subplot(3,1,2);stem([0:length(h)-1],h); ylabel(h);xlabel(Time index n); subplot(3,1,3);stem([0:length(y)-1],y); ylabel(y=x*h);xlabel(Time index n)； （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(x(n)));xlabel(Time index n); subplot(3,1,2);stem([0:length(y)-1],y); ylabel(y(n));xlabel(Time index n); subplot(3,1,3);stem([0:length(h)-1],h); ylabel(h(n)=x(n)*y(n));xlabel(Time index n); （5）已知信号x(t)=(1+t/2)[u(t+2)-u(t-2)], 求x(t+2),x(t-2),x(-t),x(2t),-x(t) t=-10:10; m=(t=2); n=(t=-2); x=(1+(t./2)).*(n-m); plot(t+2,x); t=-10:10; m=(t=2); n=(t=-2); x=(1+(t./2)).*(n-m); plot(t-2,x); t=-10:10; m=(t=2); n=(t=-2); x=(1+(t./2)).*(n-m); plot(-t,x); t=-10:10; m=(t=2); n=(t=-2); x=(1+(t./2)).*(n-m); plot(2*t,x); t=-10:10; m=(t=2); n=(t=-2); x=(1+(t./2)).*(n-m); plot(t,-x); 信号的频域分析 2.1 利用DFT分析连续信号频谱 1．