Question

求matlab程序解释，越详细越好。是关于卷积的源代码。

function a=myconv(b,c)
bs=size(b);
cs=size(c);
i=any(bs-cs);
if i
error('error')
end
i=any(~(bs-1));
if ~i
error('error')
end
ko=0;
if bs(1)>bs(2)
b=b';
c=c';
ko=1;
end
bs=size(b);
cs=size(c);
ss=2*bs(2)-1;
a=zeros(1,ss);
for i=1:cs(2)
q=zeros(1,i-1);
p=zeros(1,ss-cs(2)+1-i);
ba=[q,c,p];
ma=b(i)*ba;
a=a+ma;
end
if ko
a=a';
end
end

Answer 1

个人认为这个程序有问题
i=any(bs-cs);
if i
error('error')
end
这里明显是没道理的，如果bs和cs不一样，i就会是1，然后这里就会报错。这很明显不合理。
如果b=[1 1],c=[1 1 1]，那么bs=[1 2],cs=[1 3]，从而i=1，程序到这里就报错结束了，但是很明显b和c是可以做卷积的，结果是[1 2 2 1]
MATLAB里面的内置卷积源程序是conv.m，卷积说白了就是多项式乘法，标准程序里是用FIR滤波器写的，很精辟。你可以试试conv([1 1],[1 1 1 ])，结果就是[1 2 2 1]，但是你这个myconv就不行。我不知道"my"是谁，你还是问本人吧。

追问

不知道您能不能给个卷积的程序看看，多谢