数字信号处理
学 院 理学院 年级专业 2011级电子信息工程01班 指导老师 刘x 学 号 201101830105 姓 名 郑xx
实验一 时域离散信号的产生及时域处理
实验目的:了解Matlab软件数字信号处理工具箱的初步使用方法。掌握其简单的Matlab
语言进行简单的时域信号分析。
实验内容:[1.1]已知两序列 x1=[0,1,2,3,4,3,2,1,0];n1=[-2:6];
x2=[2,2,0,0,0,-2,-2],n2=[2:8].求他们的和ya及乘积yp.
程序如下
x1=[0,1,2,3,4,3,2,1,0];ns1=-2; x2=[2,2,0,0,0,-2,-2];ns2=2;
nf1=ns1+length(x1)-1;nf2=ns2+length(x2)-1; ny=min(ns1,ns2):max(nf1,nf2); xa1=zeros(1,length(ny));xa2=xa1; xa1(find((ny>=ns1)&(ny<=nf1)==1))=x1; xa2(find((ny>=ns2)&(ny<=nf2)==1))=x2; ya=xa1+xa2 yp=xa1.*xa2
subplot(4,4,1),stem(ny,xa1,'.') subplot(4,1,2),stem(ny,xa2,'.') line([ny(1),ny(end)],[0,0]) subplot(4,1,3),stem(ny,ya,'.') line([ny(1),ny(end)],[0,0]) subplot(4,1,4),stem(ny,yp,'.') line([ny(1),ny(end)],[0,0]) ya =
0 1 2 3 6 5 yp =
2 1 0 -2 -2
0 0 0 0 8 6 0 0 0 0 0
[1.2]编写产生矩形序列的程序。并用它截取一个复正弦序列,最后画出波形。 程序如下:clear;close all n0=input('输入序列起点:n0='); N=input('输入序列长度:N='); n1=input('输入位移:n1='); n=n0:n1+N+5; u=[(n-n1)>=0];
x1=[(n-n1)>=0]-[(n-n1-N)>=0]; x2=[(n>=n1)&(n<(N+n1))]; x3=exp(j*n*pi/8).*x2; subplot(2,2,1);stem(n,x1,'.'); xlabel('n');ylabel('x1(n)'); axis([n0,max(n),0,1]); subplot(2,2,3);stem(n,x2,'.'); xlabel('n');ylabel('x2(n)'); axis([n0,max(n),0,1]);
subplot(2,2,2);stem(n,real(x3),'.'); xlabel('n');ylabel('x3(n)的实部'); line([n0,max(n)],[0,0]);
axis([n0,max(n),-1,1]);
subplot(2,2,4);stem(n,imag(x3),'.'); xlabel('n');ylabel('x3(n)的虚部'); line([n0,max(n)],[0,0]); axis([n0,max(n),-1,1]); 输入序列起点:n0=5 输入序列长度:N=5 输入位移:n1=5
[1.3]利用已知条件,利用MATLAB生成图形。 程序如下 clear;close all N=24;M=8;
m=input('输入位移值:m='); if (m<1|m>=N-M+1)
fprintf('输入数据不在规定范围内!');break end
相关推荐: