DFT.zip_DFT 性质_DFT性质仿真_matlab对称性_傅里叶 时移_频移

clear all; x1=[1 3 5 3 6 8 3 9]; N=length(x1); %求x1[n]的长度 X1=fft(x1); X1c=circshift(X1,[1,-2]);% X1c[k]= X1[k+2]; %将x1[k]左移2位 x1c=ifft(X1c); %对X1c作傅里叶反变换 n=0:length(x1)-1; R1=real(X1); I1=imag(X1); M1=abs(X1); phase=atan2(I1,R1);%angle(X); R1c=real(X1c); I1c=imag(X1c); M1c=abs(X1c); phasec=atan2(I1c,R1c);%angle(X); figure(1); subplot(3,2,1); stem( n, x1, '.' ); xlabel( 'n' ); title('原时间序列'); subplot(3,2,2); stem( n,x1c, '.' ); xlabel( 'n' ); title('频率移位+2后的时间序列'); subplot(3,2,3); stem( n, M1, '*' ); xlabel( '频率 k' ); ylabel( '幅度' ); title( '原序列幅度谱' ); subplot(3,2,4); stem( n, M1c, '*' ); xlabel( '频率 k' ); ylabel( '幅度' ); title('频率移位后+2后的幅度谱' ); subplot(3,2,5); stem( n, phase, '*' ); xlabel( '频率 k' ); ylabel( '相位' ); title( '原序列的相位谱' ); subplot(3,2,6); stem( n, phasec, '*' ); xlabel( '频率 k' ); ylabel( '相位' ); title( '频率移位+2后的相位谱' ); X1c2=circshift(X1,[1,3]);% X1c2[k]= X1[k-3]; %将x1[k]右移3位 x1c2=ifft(X1c2); %对X1c作傅里叶反变换 n=0:length(x1)-1; R1=real(X1); I1=imag(X1); M1=abs(X1); phase=atan2(I1,R1);%angle(X); R1c2=real(X1c2); I1c2=imag(X1c2); M1c2=abs(X1c2); phasec2=atan2(I1c2,R1c2);%angle(X); figure(2); subplot(3,2,1); stem( n, x1, '.' ); xlabel( 'n' ); title('原时间序列'); subplot(3,2,2); stem( n,x1c2, '.' ); xlabel( 'n' ); title('频率移位-3后的时间序列'); subplot(3,2,3); stem( n, M1, '*' ); xlabel( '频率 k' ); ylabel( '幅度' ); title( '原序列幅度谱' ); subplot(3,2,4); stem( n, M1c2, '*' ); xlabel( '频率 k' ); ylabel( '幅度' ); title('频率移位后-3后的幅度谱' ); subplot(3,2,5); stem( n, phase, '*' ); xlabel( '频率 k' ); ylabel( '相位' ); title( '原序列的相位谱' ); subplot(3,2,6); stem( n, phasec2, '*' ); xlabel( '频率 k' ); ylabel( '相位' ); title( '频率移位-3后的相位谱' );