科技大学字信号处理课程设计.docVIP

PAGE 33数字信号处理课程设计第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．用fourier函数,理论上求下列连续时间信号的频谱。(1). 三角脉冲信号x1(t)= t=-2:0.1:2;x=tripuls(t,2,0);plot(t,x);sym