处理(第三版)(高西全)第2章时域离散信号和系统.docVIP

成 绩： 《数字信号处理》 作业与上机实验 （第二章） 班 级： 电实13-1BF 学 号： 姓 名： 曾威 任课老师： 李民 完成时间： 2014/3/23 信息与通信工程学院 2013—2014学年第 2 学期 第2章 时域离散信号和系统的频率分析 1、设计两个数学信号处理系统： 系统初始状态为零。分别用这两个系统对数字信号： 进行处理。该信号为缓慢变化的指数信号（）上叠加了一个正弦干扰噪声序列，我们希望通过该系统对进行处理来消除这个正弦干扰噪声。 .应用dtft子程序分析信号的频谱，并用MATLAB工具画出 (1)matlab代码如下： %dtft函数 function [ X,w ] = dtft( x,n,dw,k ) X=x*(exp(-1j*dw)).^(n*k); w=dw*k; end %应用dtft子程序分析信号x(n)的频谱 n=0:140; x=1.02.^n+0.5*cos(2*pi*n/8+pi/4); dw=pi/500; k=-1500:1500; [ X,w ] = dtft( x,n,dw,k ); %调用dtft函数 magX=abs(X); %信号x(n)的幅度谱 angX=angle(X); %信号x(n)的相位谱 subplot(2,1,1); plot(w/pi,magX); axis([0,1,0,800]); title(信号x(n)幅频特性); xlabel(w); ylabel(幅度); subplot(2,1,2); plot(w/pi,angX); axis([0,1,-4,4]); title(信号x(n)相频特性); xlabel(w); ylabel(相位); (2)信号的频谱图见图一: 图一 信号的频谱图 . 应用Hmp子程序分析系统一与系统二的频谱特性，画出频谱图（）。 (1)matlab代码如下： %Hmp函数 function [H,w]=Hmp(b,a,dw,k) M=length(b)-1; N=length(a)-1; ib=0:M;ia=0:N; H=(b*(exp(-j*dw)).^(ib*k))./(a*(exp(-j*dw)).^(ia*k)); w=dw*k; end %系统一的频谱特性 b=[1/8,1/8,1/8,1/8,1/8,1/8,1/8,1/8]; a=1; k=-100:100; dw=(2*pi/80); [X,w]=Hmp(b,a,dw,k); magX=abs(X); figure(1) subplot(2,1,1);plot(w/pi,mag); axis([0,1,0,1.5]); title(系统一的幅频特性); xlabel(w) ylabel(magX) angX=angle(X); subplot(2,1,2);plot(w/pi,angX); axis([0,1,-4,4]); title(系统一的相频特性); xlabel(w) ylabel(angX) %系统二的频谱特性 b=[1,-1.4142,1]; a=[1,-1.3576,0.9216]; k=-100:100; dw=(2*pi/80); [X1,w1]=Hmp(b,a,dw,k); magX=abs(X1); figure(2) subplot(2,1,1);plot(w1/pi,magX); axis([0,1,0,1.5]); title(系统二的幅频特性); xlabel(w); ylabel(magX); angX=angle(X1); subplot(2,1,2);plot(w1/pi,angX); axis([0,1,-4,4]); title(系统二的相频特性); xlabel(w); ylabel(angX); (2)系统一与系统二的频谱图见图二、三: 图二 系统一的频谱图 图三 系统二的频谱图 . 分别对这两个系统在微机上基于迭代法编程实现该信号处理算法，求解处理后的信号，画出波形。应用dtft子程序分析信号的频谱，并用MATLAB工具画出画出频谱图（）。 (1)matlab代码如下： %dtft函数 function [ X,w ] = dtft( x,n,dw,k ) X=x*(exp(-1j*dw)).^(n*k); w=dw*k