感谢@小凤同学在网上找的这段代码，我觉得写得太好了，和大家分享一下。
如果侵权，请和我联系。

function r= my_conv(a, b)
m=length(a);
n=length(b);
r=zeros(1, m+n-1);
for k = 1:m
    c = a(k)*b;
    d = r(1, k:k+n-1);
    d = d+c;
    r(1, k:k+n-1) = d;
end

运行测试：

运行环境：Win10, MATLAB 2010b

在matlab中的命令窗口运行结果如下图所示。