信号处理之圆周卷积与DFT频率分辨率分析

求4点圆周卷积与线性卷积

一、 已知x1(n)=[2，4，3，1]，x2(n)=[2，1，3]，采用时域和频域两种方法分别求它俩的4点圆周卷积与线性卷积。

1. 时域方法

function y=cirshift(x,m,N)
 
if length(x)>N
    error('N length should >=x length')
end
x=[x zeros(1,N-length(x))];
n=[0:1:N-1];
n=mod(n-m,N);
y=x(n+1);

function y=circonvt(x1,x2,N)
if length(x1)>N
    error('N length should >=x1 length')
end
if length(x2)>N
    error('N length should >=x2 length')
end
x1=[x1 zeros(1,N-length(x1))];
x2=[x2 zeros(1,N-length(x2))];
m=[0:1:N-1];
x2=