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_length);
x0=[zeros(1,L2-1),x];
y0=[yf,zeros(1,L1-1)];
z(1)=x0*y0';

for i=1:total_length-1
x0=x0(2:end);
x0=[x0,zeros(1,1)];%每次尾部补零
z(i+1)=x0*y0';
end

figure;
subplot(2,2,1);
stem(x);
title('x');
subplot(2,2,2);
stem(y);
title('y');
subplot(2,2,3);
stem(z);
title('算法实现');
w=conv(x,y);
subplot(2,2,4);
stem(w);
title('conv函数实现');
disp(z);

三、实验结果

通过对比conv函数实现的结果，测试算法的正确性。