数字信号课程设计、.doc

《数字信号处理》课程设计 学院：信息工程学院 班级：…………… 姓名：…………… 学号：…………… 第二章 信号的时域分析 任一信号可以在时域对其进行分析和描述，利用傅立叶变换理论也可以对其进行频域分析，以便更好地对信号进行存储、传输和处理，达到提取有用信号的目的。其中离散非周期信号的DTFT可表示为。 利用DFT计算序列的频谱； 解：N=51; n=-(N-1)/2:(N-1)/2; %若N为偶，n = -N / 2: (N / 2 - 1); x1=(n=0); x=((0.5).^n).*x1; X=fft(x,N); subplot(2,1,1); stem(n,fftshift(x)); ylabel(x [n]); xlabel(Time n); omega = 2 * pi / N * n; subplot (2, 1, 2 ); stem (omega, real (fftshift(X))); ylabel(X [k]); xlabel(Frequency (rad)); 2、有限长序列， 0≤n≤31，分别用N=32，N=60，N=120点DFT计算其频谱。 解： N=32; n=0:N-1; x=cos(3*pi/8*n); X=1/N*fft(x,N); omega=2*pi/N*(n-N/2); subplot(2,1,1); stem(omega,abs(fftshift(X))); axis([-pi,pi,0,1]); ylabel(Magnitude); xlabel(Frequency (rad) ); subplot(2,1,2); stem(omega,angle(fftshift(X))); axis([-pi,pi,-4,4]); ylabel(Phase); xlabel(Frequency (rad) ); N=60; n=0:N-1; x=cos(3*pi/8*n); X=1/N*fft(x,N); omega=2*pi/N*(n-N/2); subplot(2,1,1); stem(omega,abs(fftshift(X))); axis([-pi,pi,0,1]); ylabel(Magnitude); xlabel(Frequency (rad) ); subplot(2,1,2); stem(omega,angle(fftshift(X))); axis([-pi,pi,-4,4]); ylabel(Phase); xlabel(Frequency (rad) ); N=120; n=0:N-1; x=cos(3*pi/8*n); X=1/N*fft(x,N); omega=2*pi/N*(n-N/2); subplot(2,1,1); stem(omega,abs(fftshift(X))); axis([-pi,pi,0,1]); ylabel(Magnitude); xlabel(Frequency (rad) ); subplot(2,1,2); stem(omega,angle(fftshift(X))); axis([-pi,pi,-4,4]); ylabel(Phase); xlabel(Frequency (rad) ); 第三章 系统的分析与设计 信息系统就是为了有效的传输和交换信息而对载荷信息的信号进行加工、处理和变换的设备的总称。系统根据输入、输出信号的不同分为连续系统和离散系统。其分析方法有卷积法，零输入、零状态法，变换域法等。 1、 已知一个LTI系统的差分方程为： (1)初始条件y(-1)=1,y(-2)=2，输入x(n)=u(n)，计算系统的零输入响应； (2)当下面三个信号分别通过系统，分别计算系统的响应： (3)指出这是一个什么特性的系统。 解： N=50; n=0:N-1; a=[1,-1.143,0.4128]; b=[0.0675,0.1349,0.0675]; x=zeros(N); zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi); x1=cos(pi/10.*n); x2=cos(pi/5.*n); x3=cos(pi/10*7.*n); y1=filter(b,a,x1,zi); y2=filter(b,a,x2,zi); y3=filter(b,a,x3,zi); figure(1) subplot(4,1,1);stem(n,y); xlabel(n);ylabel(y(n));title(零输入响应); subplot(4,1,2);stem(n,y1); xlabel(n);ylab