convolution.m

来自「convolution with matlab」· M 代码 · 共 20 行

u=[1 0 0 0 1 1 0 1 1 0];
lu=length(u);
v=[1 1 0 1];
lv=length(v);
k=lu+lv-1;
m=0;n=0;
for  i = 1:1:k
    w(i) = 0;                      
    for j = 1:1:lv
        if(((i-j)<=0) | ((i-j)>lu))
            w(i) =w(i)+ 0 * v(j);
        else
            w(i) =w(i)+ u(i - j) * v(j);
        end
    end
end

[w]