如何在不使用内置函数的情况下在 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实现

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

How do i create my own convolution code on Matlab without using the built-in function

matlab

convolution