MATLAB 计算序列线性卷积的算法

在MATLAB平台上有多种方法可以实现线性卷积。本文主要提供方法三的实现。

一、实验方法

方法一：使用MATLAB内置的conv函数：这个函数直接用于计算两个向量的线性卷积。例如，y = conv(u, v)将返回向量u和v的线性卷积结果。如果u的长度为N，v的长度为M，则u和v卷积的结果的长度为L=N+M-1。此外，conv函数还支持指定形状的卷积结果输出，例如y = conv(u, v, shape)将返回由shape指定的u和v卷积结果的一部分。

方法二：使用FFT计算，流程图如下：

 方法三：通过傅里叶变换的定义，通过循环语句与移位补零实现线性卷积。

 二、实验代码（基于MATLAB R2023a）

clear all;
x=[1,2,3,4,5,3,2];
y=[1,2,3,4,7,6,4];

yf=flip(y);%翻转
L1=length(x);
L2=length(y);
total_length=L1+L2-1;%卷积和的长度

z=zeros(1,total_len