DSP实验一报告.docVIP

DSP实验报告书 姓名 吴晓虹 学号 200907211106 实验时间 2012.09.17 实验名称 实验一 离散信号的分析 实验目的 1.观察常用离散时间信号的图形，掌握离散时间信号的基本序列运算。 2.理解离散时间系统的时域特性，加深对离散时间系统差分方程的理解。 3.熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。 实验内容 1.在给出的区间上产生并画出下面序列： 2.设，产生并画出下列序列的样本。 实验步骤及源程序 1.新建一个m文件stepseq.m，保存到D：/matlab6p5/work路径下。 具体函数段如下： function [x,n]=stepseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)=0]; 2利用上述函数完成实验内容1。 新建一个m文件lin1.m，保存到D：/matlab6p5/work路径下。具体函数段如下： n=[0:20]; x1=n.*(stepseq(0,0,20)-stepseq(10,0,20)); x2=10*exp(-0.3*(n-10)).*(stepseq(10,0,20)-stepseq(20,0,20)); x=x1+x2; subplot(2,1,1);stem(n,x); xlable(n);ylabel(x(n));axis([0,20,-1,11]); 3.新建两个函数文件sigshift.m和sigadd.m，保存到D：/matlab6p5/work路径下。 具体函数段如下：sigshift.m程序—— function [y,n]=sigshift(x,m,n0); n=m-n0; y=x; sigadd.m程序—— function [y,n]=sigadd(x1,n1,x2,n2); m=[min(n1),min(n2)]; p=[max(n1),max(n2)]; n=min(m):1:max(p); y1=zeros(1,length(n)); y2=y1; y1(find(n=min(n1)(n=max(n1))==1))=x1; y2(find(n=min(n2)(n=max(n2))==1))=x2; y=y1+y2; 4．运用上述两个函数完成实验内容2。具体程序如下： n=-4:2;x=[1,-2,4,6,-5,8,10]; [x21,n21]=sigshift(5*x,n,-5); [x22,n22]=sigshift(4*x,n,-4); [x23,n23]=sigadd(x21,n21,x22,n22); [x2,n2]=sigadd(x23,n23,3*x,n); subplot(2,1,1);stem(n2,x2); xlabel(n);ylabel(x_1(n)); 实验结果 分析 Lin1.m的文件运行结果 Lin2.m的文件运行结果 问题讨论 设，产生并画出下列序列的样本。 sigshift.m程序—— function [y,n]=sigshift(x,m,n0); n=m-n0; y=x; sigadd.m程序—— function [y,n]=sigadd(x1,n1,x2,n2); m=[min(n1),min(n2)]; p=[max(n1),max(n2)]; n=min(m):1:max(p); y1=zeros(1,length(n)); y2=y1; y1(find(n=min(n1)(n=max(n1))==1))=x1; y2(find(n=min(n2)(n=max(n2))==1))=x2; y=y1+y2; sigmult.m程序—— function [y,n]=sigmult(x1,n1,x2,n2); m=[min(n1),min(n2)]; p=[max(n1),max(n2)]; n=min(m):1:max(p); y1=zeros(1,length(n)); y2=y1; y1(find(n=min(n1)(n=max(n1))==1))=x1; y2(find(n=min(n2)(n=max(n2))==1))=x2; y=y1.*y2; 函数程序： n=-4:2;x=[1,-2,4,6,-5,8,10];n4=-10:10;x3=2*exp(0.5*n4);x4=cos(0.1*pi); [x31,n31]=sigmult(x3,n4,x,n); [x32,n32]=sigshift(x,n,-2); [x3