数字信号处理与解答介绍.doc

  1. 1、本文档共31页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
一、数字信号处理(确定性信号) 对于一个LTI系统,设其输入序列为矩形冲激信号x(n)=u(n)-u(n-10),而冲激相应为,用MATLAB求解输出信号。可以直接调用卷积函数来实现。 解: clear all x=[1,1,1,1,1,1,1,1,1,1]; n=[0:9]; y=0.9.^n; z=conv(x,y); N=length(z); stem(0:N-1,z); 图像如下: 编程求两个序列之间的相关系数。设序列x(k)={3,11,7,0,-1,4,2},n=[-3,-2,-1,0,1,2,3],将x进行移位再加上一个白噪声信号,即y(k)=x(k-2)+w(k),其中k属于n,需要计算x序列与y序列之间的相关系数,可以使用卷积来实现。 解: clear all x=[3,11,7,0,-1,4,2]; nx=[-3:3]; [y,ny]=sigshift(x,nx,2); w=randn(1,length(y)); nw=ny; [y,ny]=sigadd(y,ny,w,nw); [x,nx]=sigfold(x,nx); [rxy,nrxy]=conv_m(y,ny,x,nx); subplot(1,1,1); stem(nrxy,rxy) axis([-5,10,-50,250]); xlabel(延迟量1); ylabel(rxy); title(噪声序列的互相关) 图像如下: 利用filter函数计算冲激相应和单位阶跃响应。设离散系统由下列差分方程表示:。 解: 冲激响应: clear all a1=[1,-1,0.9]; b1=[1]; n=0:100; x1=[1 zeros(1,100)]; %补零 y1filter=filter(b1,a1,x1); stem(n,y1filter); title(冲激响应); xlabel(x); ylabel(y); 阶跃响应: clear all a1=[1,-1,0.9]; b1=[1]; n=0:100; x2=ones(1,101); %全一矩阵 y1filter=filter(b1,a1,x2); stem(n,y1filter); title(阶跃响应); xlabel(x); ylabel(y); 图像如下: 编程求解有限时宽复指数序列,n=[0,10]的频谱及能量谱。 解: n=0:10;x=(0.9*exp(j*pi/3)).^n; k=-200:200;w=(pi/100)*k; X=x*(exp(-j*pi/100)).^(n*k); magX=abs(X);angX=angle(X); subplot(2,1,1); plot(w/pi,magX);grid xlabel(以w/pi为单位的频率); ylabel(|X|) title(幅度部分) subplot(2,1,2); plot(w/pi,angX/pi);grid xlabel(以w/pi为单位的频率); ylabel(弧度/pi) title(相角部分) figure; Y=sum(abs(x).^2); stem(Y); title(能量); 图像如下: 已知信号,随机正态白噪声为的均值为0,标准误差为1。试比较采样点数分别为45、64点时的频谱(幅度谱)图像的差异。 解: clear all; fs=100; N=45; n=0:N-1; t=n/fs; x=2*sin(4*pi*t)+5*cos(8*pi*t); y=x+randn(size(x)); X=fft(y,N); magX=abs(X); f=n*fs/N; subplot(2,1,1); plot(f,magX); xlabel(频率Hz); ylabel(45点幅度谱); clear all; fs=100; N=64; n=0:N-1; t=n/fs; x=2*sin(2*pi*2*t)+5*cos(2*pi*4*t); y=x+randn(size(x)); X=fft(y,N); magX=abs(X); f=n*fs/N; subplot(2,1,2); plot(f,magX); xlabel(频率Hz); ylabel(64点幅度谱);

文档评论(0)

挑战不可能 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档