confirm.m

用Matlab实现重叠保留法计算循环卷积的程序实例

x=[1,2,3,4,5];
n=0:4;
x1=x(mod(-n,5)+1);%x1是x的循环反对称序列
[xec,xoc]=circevod(x);%用函数circevod可以求出x的圆周偶部分和圆周奇部分
X=fft(x,5);% X是x的5点DFT
X1=fft(x1,5);
X1con=conj(X1);% X1con=X[<-k>N]的共轭,可以验证X1con=X;
Xre=real(X);
X1con_re=real(X1);%验证Xre和X1con_re是否相等
Xim=imag(X);
X1con_im=-imag(X1);%验证Xim和X1con_im是否相等
Xabs=abs(X);
X1abs=abs(X1);%验证Xabs和X1abs是否相等
theta=angle(X);
theta1=-angle(X1);%验证theta和theta1是否相等