- 1、本文档共13页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《数字信号处理》Matlab实验
学院:信息与通信工程学院
姓名: 兰梅珍
班级:2012211114
学号:2012210421
班内序号:27
实验一:实现重叠相加和重叠保留算法,完成线性卷积的分段计算(可任意指定及);
1.对算法的概括性说明;
重叠相加法和重叠保留法都是为了实现卷积的快速计算,两种方法都是以追逐段地方式通过循环卷积来完成线性卷积的计算。将输入序列x(n)进行分段,每段长为N,且N=M(M为有限长因果序列h(n)的长度),x(n)逐段与h(n)进行循环卷积,在重叠保留法中需在x(n)序列首部加入长度为M-1的0序列。在算法中,在获得N个点的输入后,进行N+M-1点循环卷积计算,之后输出N个点。通过for循环逐段进行循环卷积,使用fft和ifft计算两个有限长的N点循环卷积结果。
2.源代码及流程图;
重叠相加法
function [Y]=overl(x,h,N)
Lx=length(x);
M=length(h);
x=[x,zeros(1,N-1)];
t=zeros(1,M-1);
Y=zeros(1,Lx+M-1);
a=floor(Lx/N);
for k=0:a
A=x(k*N+1:k*N+N);
y1=fft(A,Lx+M-1);
y2=fft(h,Lx+M-1);
y3=y1.*y2;
q=ifft(y3,Lx+M-1);
Y(k*N+1:k*N+M-1)=q(1:M-1)+t(1:M-1);
Y(k*N+M:k*N+N)=q(M:N);
t(1:M-1)=q(N+1:N+M-1);
M-1点缓存end
M-1点缓存
Y(1:Lx+M-1);
+L点循环卷积 X1(n) 后M-1点 y(n)
+
L点循环卷积
L=N+M-1 h(n)
重叠保留法
function [Y]=overlsav(x,h,N)
Lx=length(x);
M=length(h);
M1=M -1;
L=N-M1;
h=[h,zeros(1,N-M)];
x=[zeros(1,M1),x,zeros(1,N-1)];
a=floor ((Lx+M1-1)/(L))+1;
Y=zeros(1,N);
for k=0:a-1
xk=x(k*L+1:k*L+N);
b=fft(xk,N);
C=fft(h,N);
Z=b.*C;
Y(k+1,:)=ifft(Z,N);
end
Y=Y(:,M:N);
Y=(Y(:));
序列衔接N+M-1点循环卷积M-1点缓存
序列衔接
N+M-1点循环卷积
M-1点缓存
Y(n)
每N点输入 h(n)
3.实验结果;
重叠相加法
x=[1,2,3,4,5,6,7,8,9,10]
h=[1,0,-1]
N=6
重叠保留法
x=[1,2,3,4,5,6,7,8,9,10]
h=[1,0,-1] N=4
4.结合教材3.5.1节作运算量分析;
有限长因果序列x(n)h(n)的长度分别为N和M,直接计算线性卷积y(n),y(n)可视为N个序列的叠加结果,序列长度为M,所以每成一个序列需完成M次乘法,则共需要完成MN次乘法运算。这N个序列依次向右移动一位故需(N-1)(M-1)次加法运算。
5.总结实验过程中出现的问题以及解决问题的具体措施。
N越大时,算法的运算量越小,即实现重叠相加法和重叠保留法的运算量越小。
实验二:周期序列的谱分析
利用DFT分析模拟信号之频谱;试验要求:
1.设定采样周期并说明原因;
由奈奎斯特抽样条件,fs=2f=2*8=16Hz,所以取fs=20Hz,T=1/f=0.05s。
2.若令,确定该序列之周期并说明原因;
因为x(n)=cos(4/5πn),2π/4=5/2,所以周期N=5.
3.绘制10个周期内的取值情况;
T=1
您可能关注的文档
最近下载
- PICC护士进修汇报心得ppt.pptx
- 无处不在的人格.docx
- 安全生产法律法规培训.ppt
- 园林植物识别与应用高职完整全套教学课件.pptx
- 糖尿病护理查房PPT课件.pptx VIP
- 2024~2025学年八年级数学上册第13章轴对称(4)—考点考题点点通「含答案」.pdf
- 部编版四年级语文下册四下语文《语文园地一》优秀教学设计精选合集-微课教案.doc VIP
- 新教材高中地理第六章自然灾害第四节地理信息技术在防灾减灾中的应用教案含解析新人教版必修第一册.docx
- Unit 3 How many(复习课件)四年级英语上册期末核心考点(译林版三起).pptx
- 华能青海分公司招聘笔试题库2024.pdf
文档评论(0)