MATLAB之卷积

卷积的物理意义

输入信号通过一个新兴因果移不变系统后的输出为信号与该系统的冲击响应的卷积。任何一个信号都可以表示为单位冲击信号之和。

相关卷积定理

卷积与互相关，自相关函数的关系

条件：假设系统的输入，输出，和单位脉冲响应分别为x(n), y(n), h(n).
Rxy(m) = Rxx(m)*h(m)
对上式进行傅里叶变换，得：
Sxy(w) = Sxx(w)*H(jw)

卷积与单位脉冲响应的关系

x(t)*δ(t) = x(t)
x(t)*δ(t-t0) = x(t-t0)

卷积的计算

卷积向量表示形式

一方不动，另一方前后倒置，从左向右移动，每次移动后重叠部分的乘积，再组合起来即输出信号

                  h1  h2  h3
x4  x3  x2  x1

已知：
输入信号：x=[x1,x2,x3,x4] 维度：4
系统响应：h=[h1,h2,h3] 维度：3
输出信号：y=[y1,y2,y3,y4,y5,y6] 维度：4+3-1=6

其中：

y1=x1*h1
y2=x1*h2+x2*h1
y3=x1*h3+x2*h2+x3*h1
y4=x2*h3+x3*h2+x4*h1
y5=x3*h3+x4*h2
y6=x4*h3

卷积矩阵表达形式

将上式写成矩阵形式：
法一：

法二 ：

基于法一，用MATLAB求X矩阵的三种方法

% 功能：用上述的法一用三种方法求X大矩阵

% 编辑者：lily
% 日期：2019,4,14

clear;
clc;
close all;
% ======================= input signal ==========================
x=[1,2,3,4];
h=[5,6,7];
%求矩阵的向量个数
xn=size(x,2);%求矩阵列数，即几个
hn=size(h,2)