精编完整版科技大学数字信号处理毕业论文.doc

(此文档为word格式，下载后您可任意编辑修改！) 数字信号处理课程设计 第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); （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); 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],)=[0 1 2 3],求y(n) x=[1,1,1]; =exp(-1*t); x=m.*u; y=n.*u; )));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()=x(n)*y(n));xlabel(Time index n); （5）已知信号x(t)=(1+t2)[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); syms t w xt=sym((t+1)*Heaviside(t+1)-2*t*Heaviside(t)+(t-1)*Heaviside(t-1)); Fw=fourier(xt,t,w); FFw=maple(convert,Fw,piecewise); FFP=abs(FFw); ezplot(FFP,[-10*pi 10*pi]) axis([-10*pi 10*pi 0 1]) (2). 单边指数信号x2(t)=eu(t) N=256; fs=6; Ts=1fs; t=(0:N-1)*Ts; x=exp(-t); y=fft(x); mag=Ts*abs(y); Ws=2*pi*fs; w=(0:length(y)-1)*Wslength(y); X=1.sqrt(w.^2+1); w1=