如何在不使用内置函数的情况下在 Matlab 上创建自己的卷积代码答案

【问题标题】：How do i create my own convolution code on Matlab without using the built-in function如何在不使用内置函数的情况下在 Matlab 上创建自己的卷积代码
【发布时间】：2019-09-03 16:52:38
【问题描述】：

我正在尝试在 Matlab 上创建一个程序来显示 2 个信号的卷积，但不使用“conv”函数。

输入是??????(????)=??????xp(-4t+4)*????(????−1)+0.5????(?? ??−5)+2????(??????−10)

我计算出 h(t) = 4*exp(-4t)*u(t)。我正在尝试在 MATLAB 上编写一个程序来计算和显示代表两个信号的两个数组的卷积，但不使用内置的“conv”函数。我有一个工作代码，但它包含“conv”，所以我试图在不使用它的情况下这样做。

t=linspace(0,11,1101);  

%create more than 500 time-series values

vs=exp(-4*(t-1)).*heaviside(t-1)+0.5*heaviside(t-5)+2*dirac(t-10);

%Input signal array

h=4*exp(-4*t).*heaviside(t);

%Input signal array

vc=conv(vs,h);

t1=linspace(0,22,2201);

%Convolution of two arrays

subplot(3,1,1)

plot(t,vs)

xlabel('t')

ylabel('v_s(t)')

title('Input signal')

subplot(3,1,2)

plot(t,h)

xlabel('t')

ylabel('h(t)')

title('System impulse response')

subplot(3,1,3)

plot(t1,vc)

xlabel('t')


ylabel('v_c(t)')

title('Output signal')

【问题讨论】：

请在此处发布之前搜索您的问题。这里有数百个问题，来自人们发布他们的卷积代码并询问它有什么问题。您可以选择其中任何一个，根据发布的答案修复它们，然后完成。

标签： matlab convolution

【解决方案1】：

尝试用这个替换 conv 行：

N = length(vs);
M = length(h);   
lout=N+M-1;
vc=zeros(1,lout); 
for i = 1:N
    for k = 1:M
      vc(i+k-1) = vc(i+k-1) + h(k)*vs(i);
    end
end

是卷积公式的基本matlab实现

【讨论】：