【含源代码】北邮dsp-MATLAB试验一重叠相加和重叠保留

重叠相加法与重叠保存法的原理实现侯凯（吉林大学 通信工程学院 吉林 长春 130012）0概述线性卷积是求离散系统响应的主要方法之一，许多重要应用都建立在这一理论基础上，如卷积滤波等。

用圆周卷积计算线性卷积的方法归纳如下:将长为N 2的序列x(n)延长到L,补L -N 2个零，将长为N 1的序列h(n)延长到L,补L -N 1个零。

如果L ≥N1+N2-1,则圆周卷积与线性卷积相等，此时,可有FFT 计算线性卷积，方法如下：a.计算X(k)=FFT[x(n)]b.求H(k)=FFT[h(n)]c.求Y(k)=H(k)Y(k) k=0～L -1d.求y(n)=IFFT[Y(k)] n=0～L -1可见，只要进行二次FFT,一次IFFT 就可完成线性卷积计算。

上述结论适用于x(n)、h(n)两序列长度比较接近或相等的情况，如果x(n)、h(n)长度相差较多。

例如，h(n)为某滤波器的单位脉冲响应,长度有限，用来处理一个很长的输入信号x(n)，或者处理一个连续不断的信号，按上述方法，h(n)要补许多零再进行计算，计算量有很大的浪费，或者根本不能实现。

为了保持快速卷积法的优越性，可将x(n)分为许多段后处理，每小段的长与h(n)接近，其处理方法有两种：重叠相加法和重叠保留法。

1重叠相加法——由分段卷积的各段相加构成总的卷积输出假定x i (n)表示图中第i 段x(n)序列如下图：22()(1)1()0i x n iN n i N x n ≤≤+-⎧=⎨⎩则输入序列可表为：()()i i x n x n ∞=-∞=∑图1 长序列分段滤波于是输出可分解为： ()()*()()*()()i i i i i y n x n h n x n h n y n ∞∞=-∞=-∞===∑∑其中 ()()*()i i y n x n h n =由此表明，只要将x(n)的每一段分别与h(n)卷积，然后再将这些卷积结果相加起来就可得到输出序列，这样，每一段的卷积都可用上面讨论的快速卷积来计算。

matlab中加法运算Matlab是一种强大的数值计算和数据可视化工具，它提供了丰富的函数和工具箱，方便用户对数据进行加法运算和相关的数学操作。

本文将介绍在Matlab中进行加法运算的方法和技巧。

我们需要了解Matlab中的基本数据类型。

Matlab支持多种数据类型，包括数字、字符串、矩阵等。

对于数字类型，可以直接使用加号（+）进行加法运算。

例如，我们可以计算两个数字的和：a = 3;b = 4;c = a + b;在上述代码中，我们定义了两个变量a和b，分别赋值为3和4。

然后使用加号将它们相加，并将结果赋值给变量c。

通过输出变量c 的值，我们可以得到两个数字的和。

在Matlab中，使用分号（;）可以抑制输出结果，如果不希望输出结果，可以在代码末尾加上分号。

除了数字类型，Matlab还支持字符串类型。

对于字符串类型，加法运算表示字符串的拼接。

例如，我们可以将两个字符串拼接在一起：str1 = 'Hello';str2 = 'World';str3 = str1 + ' ' + str2;在上述代码中，我们定义了两个字符串变量str1和str2，分别赋值为'Hello'和'World'。

然后使用加号将它们拼接在一起，并添加一个空格，将结果赋值给变量str3。

通过输出变量str3的值，我们可以得到拼接后的字符串。

除了基本数据类型，Matlab还支持矩阵和向量的加法运算。

在Matlab中，可以使用矩阵和向量进行加法运算，得到对应位置上的元素相加的结果。

例如，我们可以计算两个矩阵的和：A = [1 2; 3 4];B = [5 6; 7 8];C = A + B;在上述代码中，我们定义了两个2×2的矩阵A和B，然后使用加号将它们相加，并将结果赋值给矩阵C。

通过输出矩阵C的值，我们可以得到两个矩阵对应位置上的元素相加的结果。

实验一MATLAB基本操作及运算MATLAB是一种强大的数值计算和数据可视化工具，广泛应用于科学研究、工程设计、数据分析等领域。

本文将介绍MATLAB的基本操作和运算。

首先，我们需要了解MATLAB中的基本数据类型，包括数值型、字符型和逻辑型。

数值型可以是整数、实数、复数等；字符型用单引号或双引号包围字符；逻辑型用true和false表示。

MATLAB提供了各种数学运算函数，包括四则运算、三角函数、指数函数等。

例如，加法可以使用加号（+），减法可以使用减号（-），乘法可以使用乘号（*），除法可以使用除号（/）。

三角函数可以使用sin、cos、tan等函数，指数函数可以使用exp函数。

此外，还可以使用log 函数进行对数运算。

MATLAB还可以进行矩阵运算。

矩阵可以使用方括号（[]）表示，每一行用分号（;）分隔。

可以使用矩阵乘法运算符（*）进行矩阵相乘，使用点乘运算符（.）进行矩阵对应元素的运算。

矩阵还可以进行转置、逆运算等。

除了基本运算，MATLAB还提供了各种其他功能。

例如，可以使用plot函数进行数据可视化，使用subplot函数绘制多个图形。

可以使用for循环和while循环进行循环操作，使用if语句进行条件判断。

MATLAB还可以进行文件读写操作。

可以使用load函数从文件中加载数据，使用save函数将数据保存到文件中。

可以使用fopen函数打开文件，使用fclose函数关闭文件。

可以使用fprintf函数写入文本文件，使用fscanf函数读取文本文件。

还可以使用imread函数读取图像文件，使用imwrite函数保存图像文件。

MATLAB还具备向量化的能力。

向量化是指使用矩阵代替循环进行计算，能够提高代码的执行效率。

例如，可以使用点乘运算符（.）对矩阵的每个元素进行计算，而不是使用循环逐个计算。

使用向量化的方法，可以更加简洁地编写代码。

在MATLAB中还有很多强大的功能等待探索，例如符号计算、模拟仿真、深度学习等。

maltlab光谱叠加MATLAB是一种非常强大的计算机语言和程序设计环境，广泛应用于科学、工程和工业等领域。

光谱叠加是MATLAB中经常使用的一个功能，它可以将多个光谱图像叠加在一起，以便于分析和比较。

光谱叠加的基本原理是将不同光谱的数据点对应起来，然后将它们的反射率或吸光度值相加。

在MATLAB中进行光谱叠加的方法主要有以下几种：1.使用plot函数进行光谱叠加。

这种方法需要将每个光谱的数据点先读取进MATLAB中，并将它们存储为向量或矩阵形式。

然后，使用plot函数将这些图像绘制到同一个坐标系中，并按照需要进行调整和编辑。

2.使用多条曲线的叠加命令“hold on”和“hold off”。

这种方法需要使用hold on命令来开启允许多个曲线的叠加，然后用hold off命令来结束叠加。

它与plot函数相比更加灵活和方便，因为可以随时添加或删除曲线。

3.使用matshow或imagesc函数对光谱数据进行可视化。

这种方法适用于需要对大量光谱数据进行快速分析和比较的情况。

它可以将多个光谱的反射率或吸光度数据绘制成矩阵形式的图像，并进行一定程度的色彩变换和调整。

4.使用ColorOrder和LineOrder命令来设置光谱的颜色和线条样式。

这种方法可以根据需要设置不同颜色和线条样式的光谱图像，以便于更直观和清晰地比较和分析它们之间的不同。

在进行光谱叠加时，需要注意一些问题和技巧。

比如，要先对所读取的光谱数据进行必要的清理和预处理，以确保数据质量和准确性；同时要合理选择坐标系和比例尺，以避免数据重合和混淆；此外，还要进行一定的数据平滑、峰值提取和分析等操作，以更好地理解和掌握光谱数据的含义和特征。

总之，光谱叠加是MATLAB中非常重要和实用的一个功能，可以帮助研究人员更好地进行光谱分析和比较。

掌握好这种方法，可以大大提高工作效率和分析能力，进一步推动科学研究和工程应用的进步和发展。

实验要求1.每个实验进行之前须充分预习准备，实验完成后一周内提交实验报告；2.填写实验报告时，分为实验题目、实验目的、实验内容、实验结果、实验小结五项；3.实验报告要求：实验题目、实验目的、实验内容、实验结果四项都可打印；但每次实验的实验内容中的重要代码（或关键函数）后面要用手工解释其作用。

实验小结必须手写！（针对以前同学书写实验报告时候抄写代码太费时间的现象，本期实验报告进行以上改革）。

实验一信号、系统及系统响应实验目的：1. 掌握使用MATLAB进行函数、子程序、文件编辑等基本操作；2. 编写一些数字信号处理中常用序列的3. 掌握函数调用的方法。

实验内容：1.在数字信号处理的基本理论和MATLAB信号处理工具箱函数的基础上，可以自己编写一些子程序以便调用。

（1）单位抽样序列δ(n-n0)的生成函数impseq.m（2）单位阶跃序列u(n-n0)的生成函数stepseq.m（3）两个信号相加的生成函数sigadd.m（4）两个信号相乘的生成函数sigmult.m（5）序列移位y(n)=x(n-n0)的生成函数sigshift.m（6）序列翻褶y(n)=x(-n)生成函数sigfold.m（7）奇偶综合函数evenodd.m（8）求卷积和2.产生系列序列，并绘出离散图。

（1） x1(n)=3δ(n-2)-δ(n+4) -5≤n≤5（2） x3(n)=cos(0.04πn)+0.2w(n) 0≤n≤50其中：w(n)是均值为0，方差为1 的白噪声序列。

3.设线性移不变系统的抽样响应h(n)=(0.9)^n*u(n),输入序列x(n)=u(n)-u(n-10),求系统的输出y(n).实验二 系统响应及系统稳定性1.实验目的（1）掌握 求系统响应的方法。

（2）掌握时域离散系统的时域特性。

（3）分析、观察及检验系统的稳定性。

2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。

Dsp-matlab实验实验一：重叠相加法和重叠保留法的实现设计报告课题名称：学生姓名：级：班班内序号：学号：2015/06/15 日期：目录一、实验原理·········································二、Matlab源代码·································运行结果Matlab三、····························结果分析Matlab四、····································五、遇到的难题与解决方法····························参考文献·························································一、实验原理1、算法来源DFT 是连续傅里叶变换在时域和频域上都离散的形式，将时域信号的采样变换为在离散时间傅里叶变换频域的采样。

在Matlab中进行图像融合与图像叠加的方法与技巧引言：随着数字图像处理和计算机视觉领域的发展，图像融合和图像叠加变得越来越重要。

图像融合是指将多幅图像合成为一幅具有更清晰、更丰富信息的图像，而图像叠加则是在保留所叠加图像的原始信息的同时，使图像更加丰富和易于理解。

Matlab作为一种强大的科学计算工具，提供了丰富的图像处理函数和工具箱，可以很方便地进行图像融合与图像叠加。

一、图像融合的方法与技巧1. 融合算法图像融合的基本方法有加权平均法、空间域融合法、频域融合法、小波融合法等。

加权平均法是最简单的方法，通过计算图像像素的平均值来融合。

空间域融合法是通过对直接融合的图像进行空间域操作来提取融合结果。

频域融合法则是通过将图像转换到频域，然后进行频域操作来实现融合。

小波融合法是基于小波变换的方法，利用小波分析的多尺度分解能力对图像进行分析和融合。

根据具体需求和图像的特点，选择合适的融合算法是非常重要的。

2. 图像预处理在进行图像融合之前，通常需要进行图像预处理，以提高融合结果的质量。

常用的图像预处理方法包括灰度拉伸、直方图均衡化、滤波等。

灰度拉伸是通过对图像的像素值进行线性变换，将图像像素值的范围拉伸到合适的范围内，从而增加图像的对比度。

直方图均衡化则是将图像的像素值在灰度直方图上均匀分布，以增强图像的细节。

滤波是通过对图像进行滤波操作，如低通滤波、高通滤波等，以去除图像中的噪声和不需要的细节。

3. 图像融合的策略图像融合的策略可以根据具体需求来选择。

常见的策略包括全局融合和局部融合。

全局融合是将所有图像的信息进行融合，得到整体的融合结果。

而局部融合则是将不同图像的不同区域进行融合，以保留更多的细节和纹理。

根据具体应用和需求，选择合适的融合策略可以使融合结果更加符合实际需求。

4. 参数设置与调整在进行图像融合过程中，不同的算法和方法有各自的参数，根据不同的图像和具体应用，需要适时地进行参数的设置和调整。

不同相位脉冲叠加matlab在MATLAB中，可以使用不同相位的脉冲叠加来模拟复杂的波形。

首先，我们可以创建一个基本的脉冲，然后通过改变相位来叠加多个脉冲。

下面我将从不同角度来介绍如何在MATLAB中实现不同相位脉冲的叠加。

1. 创建基本脉冲。

首先，我们可以创建一个基本的脉冲，比如高斯脉冲。

可以使用MATLAB内置的函数，如`gauspuls`来生成一个高斯脉冲。

例如：matlab.t = -10:0.1:10;y = gauspuls(t,1,0.5);plot(t,y);这将生成一个以1为中心频率、0.5为带宽的高斯脉冲，并将其显示出来。

2. 创建不同相位的脉冲。

接下来，我们可以创建不同相位的脉冲。

可以通过改变脉冲的时间延迟来实现不同相位的叠加。

例如，可以使用`exp`函数来实现时间延迟，然后将不同相位的脉冲相加。

示例代码如下：matlab.t = -10:0.1:10;y1 = gauspuls(t,1,0.5); % 基本脉冲。

y2 = gauspuls(t,1,0.5).exp(1i0.5); % 相位为0.5的脉冲。

y3 = gauspuls(t,1,0.5).exp(1i1); % 相位为1的脉冲。

y_total = y1 + y2 + y3; % 叠加不同相位的脉冲。

plot(t, abs(y_total)); % 显示叠加后的波形。

在这个例子中，我们创建了两个具有不同相位的脉冲，并将它们与基本脉冲相加，最后显示出叠加后的波形。

3. 分析叠加波形。

除了显示叠加后的波形外，我们还可以对叠加波形进行进一步的分析。

比如，可以计算叠加波形的频谱、脉冲宽度等特性，并进行相应的可视化展示。

总之，在MATLAB中实现不同相位脉冲的叠加，可以通过创建基本脉冲，改变相位并相加的方式来实现。

同时也可以对叠加后的波形进行进一步的分析和处理，以满足具体的需求。

Matlab实验报告学院：信息与通信工程学院班级：201321113学号：2013210381班内序号：18姓名：石雪原实验题目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点循环卷积结果。

二．源代码和流程图重叠相加法代码function [Y]=overpl(x,h,N)Lx=length(x); % 序列长度M=length(h); %h (n )长度x=[x,zeros(1,N-1)];t=zeros(1,M-1);Y=zeros(1,Lx+M-1);a=floor(Lx/N);for k=0:aA=x(k*N+1:k*N+N);y1=fft(A,Lx+M-1); % 利用fft 进行运算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);endY(1:Lx+M-1);对应流程图图一：每段利用fft和ifft实现循环卷积图二：对X (n)逐段进行循环卷积然后相加得输出线性卷积结果y (n)重叠保留法代码fun cti on [Y]二overlpsav(x,h,N)Lx=le ngth(x);M=le ngth(h);M仁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(：))'对应流程图装成单列向量再转置成行向量 输出输入数据 x I 各段搭接长 长度及脉冲 —\度M1，有效 响应长度k数据长度L将h 延长 至循环 长度N㈡把x 前面 加上（M- 1）个零各段进 行卷积把 K+1X N 阶输出矩 阵Y 初始 化Y 中各行均去 掉前M-1个样 Z 本，转置后构V成新的Y三.实验结果重叠相加法x=[1,2,3] h二[1,2,3] N=41 5]? 3X] 1 K 12 9 0緒=1 1 10 12 9 0重叠保留法x=[1,2,3,4,5,6,7,8,9,10] h=[1,0,-1] N=4四.结果分析（有关运算量的定量分析结果）有限长因果序列x （n）h （n）的长度分别为N和M直接计算线性卷积y （n），y（n）可视为N个序列的叠加结果，序列长度为M 所以每生成一个序列需完成M次乘法，共需完成MN次乘法运算。