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

Question

我正在尝试在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')

Answer 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实现