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

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