conv1.m

来自「有关dsp的,有好多子程序,还是不错的,欢迎使用」· M 代码 · 共 41 行

function y=conv1(x,h,N,M)
if M>=N
   for n=2:M+N
      y(n)=0;
      if n<=N+1
         for k=1:n-1
            y(n)=y(n)+x(k)*h(n-k);
         end
      end
      if n>N+1&&n<=M+1
          for k=1:N
              y(n)=y(n)+x(k)*h(n-k);
          end
      end
      if n>M+1
        for k=n-M:N
          y(n)=y(n)+x(k)*h(n-k);
        end
      end
   end
else
    for n=2:M+N
      y(n)=0;
      if n<=M+1
         for k=1:n-1
            y(n)=y(n)+x(k)*h(n-k);
         end
      end
      if n>M+1&&n<=N+1
          for k=n-M:n-N
             y(n)=y(n)+x(k)*h(n-k);
          end 
      end
      if n>N+1
        for k=n-M:N
          y(n)=y(n)+x(k)*h(n-k);
        end
      end
    end
end