数字信号处理实验报告2.doc
上传人:sy****28 上传时间:2024-09-13 格式:DOC 页数:8 大小:383KB 金币:14 举报 版权申诉
预览加载中,请您耐心等待几秒...

数字信号处理实验报告2.doc

数字信号处理实验报告2.doc

预览

在线预览结束,喜欢就下载吧,查找使用更方便

14 金币

下载此文档

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

信息与通信工程学院数字信号处理实验报告实验题目:线性卷积的分段计算班级:2009211114姓名:廖宁学号:09210411日期:2011.6实验要求:实现重叠相加和重叠保留算法,完成线性卷积的分段计算(可任意指定及);试验报告包括:对算法的概括性说明;源代码及流程图;实验结果;结合教材3.5.1节作运算量分析;总结实验过程中出现的问题以及解决问题的具体措施。重叠相加法:基本步骤:对x(n)修饰,不够点数的在后面补零。分段,x(n)按点数的N分成几段x1(n),x2(n)……xk(n)等。X1(n),x2(n)……xk(n)等与h(n)进行循环卷积。得到的各段卷积结果重叠部分相加,重新连接成一个新的序列即为线性卷积结果。开始流程图:输入xn,hn和长度NL点循环卷积重叠相加输出前N点否是否完成全部数据处理是连接输出y(n)源代码:functiony=circular_conv(x1,x2,L)%循环卷积X1k=fft(x1,L);X2k=fft(x2,L);Yk=X1k.*X2k;y=ifft(Yk);function[y]=overlaplus(x,h,N)%重叠相加法lenx=length(x);M=length(h);L=N+(M-1);x=[x,zeros(1,L-1)];K=floor(lenx/N);fork=0:Kxk=x(k*N+1:k*N+N);Y(k+1,:)=circular_conv(xk,h,L);endtn=zeros(1,M-1);fork=0:KY(k+1,1:M-1)=Y(k+1,1:M-1)+tn;tn=Y(k+1,N+1:L);endY=Y(:,2:N+1)';y=(Y(:))';实验结果:X1(n)=0:10;X2(n)=[10-1];N=6;Y(n)=[1222222222-9-10]结果正确。重叠保留法:基本描述:重叠保留法和重叠相加法都用循环卷积求线性卷积,重叠保留法中,相邻的两段循环卷积结果含有M-1点的重叠,只要将结果中发生重叠的前M-1点舍去,保留未重叠的部分并输出,这就是重叠保留法的核心。开始流程图:是否处理完全部数据填入保留值,进行分段卷积分段成xk(n)输入xn,hn和N否输出y(n)是源代码:functiony=overlapsave(x,h,N)%重叠保留法Lx=length(x);M=length(h);ifN<MN=M+1;endL=N+M-1;t=zeros(1,M-1);T=ceil(Lx/N);x=[x,zeros(1,(T+1)*N-Lx)];y=zeros(1,(T+1)*N);fori=0:1:Txi=i*N+1;x_seg=[t,x(xi:xi+N-1)];t=x_seg(N+1:N+M-1);y_seg=circular_conv(x_seg,h,L);y(xi:xi+N-1)=y_seg(M:N+M-1);end实验结果:X1(n)=0:10;X2(n)=[10-1];N=6;Y(n)=[1222222222-9-10]结果正确。