数字信号处理 时域离散随机信号处理 实验报告及答案 (高西全 丁玉美)

实验报告配套教材：

书名：数字信号处理—时域离散随机信号处理
作者：高西全 丁玉美
出版社：西安电子科技大学出版社

实验报告概述：

DSP实验报告书 姓名 学号 实验时间 课题名称 离散信号的分析 设计原理图及源程序 1.新建一个m文件stepseq.m，保存到C：/matlab7/work路径下。 具体函数段如下： function [x,n]=stepseq(n0,n1,n2) %stepseq 阶跃函数 n=[n1:n2]; %序列号 x=[(n-n0)>=0];。 %n-n0>=0时值为1 2利用上述函数完成实验内容1，具体程序如下： n=[0:20]; %n的取值为0-20，步数为1 x1=n.*(stepseq(0,0,20)-stepseq(10,0,20));%x1=n*[u(n)-u(n-10)] x2=10*exp(-0.3*(n-10)).*(stepseq(10,0,20)-stepseq(20,0,20));%x2=10*exp[-0.3*(n-10)]*[u(n- 10)-u(n-20) x=x1+x2; subplot(2,1,1); %将图形窗口分为2x1个子图，在第一个子图处绘制图形 stem(n,x); %画出图形 xlable('n');ylabel('x(n)'); axis([0,20,-1,11]); %为了使波形顶部避开图框，改变图框坐标 3.新建两个函数文件sigshift.m和sigadd.m，保存到C：/matlab7/work路径下。 具体函数段如下：sigshift.m程序—— function [y,n]=sigshift(x,m,n0); %sigshift序列移位y(n)=x(n-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;%把x1补0，扩展为具有y长度的y1 y2(find(n>=min(n2)&(n<=max(n2))==1))=x2;%把x2补0，扩展为具有y长度的y2 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)'); 仿真波形图 实验结果 如上图 问题讨论 利用第二问的两个m文件编写程序如下： n=-4:2;x=[1,-2,4,6,-5,8,10]; [x21,n21]=sigshift(2.*exp(0.5*n).*x,n,0); [x22,n22]=sigshift(cos(0.1*pi*n).*x,n,-2); [x2,n2]=sigadd(x21,n21,x22,n22); subplot(1,1,1); stem(n2,x2); xlabel('n');ylabel('x_1(n)'); 教师评分 操作成绩 报告成绩 教师签名 日 期