时间域上实数序列位移后DFT是否可以还是实数序列？

如何使实数序列在时间域上位移后还是实数序列？

---

---

前言

过完年放松几天，今天开始充电，买了本《MATLAB信号处理》看看，学到了一个特别好的知识点，特给大家分享一下！

---

提示：以下是本篇文章正文内容

一、DFT移位性质

一个序列在时间域上位移后进行DFT运算，根据性质得是一个复数序列。但是我们现在希望通过我们的方式最后得到一个实数序列。

其中m必须为整数，或正数或负数。

因为WN破坏了原本X(K)的共轭对称，使得IFFT后得到的x(n-m)不是一个实数序列。

二、解决方案

1.m为整数的情况

m为整数时，按照常规的解决方法就就可以。
（1）计算出位移量，再算出旋转因子。
（2）原序列FFT后得到X(k)后乘旋转因子。
（3）进行IFFT变换，然后取逆变换后的实部
这样就得到了原序列位移后的实数序列。

例：设x=cos(2pi f0* n/fs+ph1),其中fs采样频率为2000 Hz, .信号频率f0为100 Hz,ph1 是信号的初始相角一pi,n是样点值,信号长度N为40。希望通过位移使该信号的初始相角为0。
MATLAB代码如下：

fs = 2000;                      % 采样频率
N = 40;                         % 信号长度
n = 0:N-1;                      % 样点序列
f0 = 100; ph1=-pi;            % 初始频率和初始相角
x=cos(2*pi*f0*n/fs+ph1);        % 余弦信号序列
X=fft(x);

d=ph1*fs/f0/(2*pi);             %计算位移量      ph1=-pi时  d=-10为整数
Ex=exp(-1j*2*pi*n*d/N);         % 计算旋转因子W^(dk) 
Y=X.*Ex;                        %FFT结果后乘旋转因子
y=ifft(Y);                      %逆变换

%作图
subplot(211);
plot(n,x,'k');
axis([0 N -1.1 1.1]); 
title('原始信号'); ylabel('幅值'); xlabel('样点'); 

subplot(212);
plot(n,y,'k');
axis([0 N -1.1 1.1]); 
title('位移后信号的实部'); ylabel('幅值'); xlabel('样点'); 

得到的结果：

结果显然初相角为0

若上题中的ph1值为-pi/3,这时位移量不是整数，用上述方法得到的结果就不准确了。
MATLAB代码如下:

fs = 2000;                      % 采样频率
N = 40;                         % 信号长度
n = 0:N-1;                      % 样点序列
f0 = 100; ph1=-pi/3;            % 初始频率和初始相角
x=cos(2*pi*f0*n/fs+ph1);        % 余弦信号序列
X=fft(x);

%使用ceil（）函数进行取整。
d=ceil(ph1*fs/f0/(2*pi));             %计算位移量      ph1=-pi/3时  d不为整数
Ex=exp(-1j*2*pi*n*d/N);         % 计算旋转因子W^(dk) 
Y=X.*Ex;                        %FFT结果后乘旋转因子
y=ifft(Y);                      %逆变换

%作图
subplot(211);
plot(n,x,'k');
axis([0 N -1.1 1.1]); 
title('原始信号'); ylabel('幅值'); xlabel('样点'); 

subplot(212);
plot(n,real(y),'k');
axis([0 N -1.1 1.1]); 
title('位移后信号的实部'); ylabel('幅值'); xlabel('样点'); 

结果如下：

可以看出距离初相角为0有一定的误差，原因是位移量d为小数，在进行计算时，只可以使用整数计算，所以用了ceil（）的取整函数，数值处理产生误差。所以我们引入了下面的方法来解决m不为整数

2.m不为整数

方法：将实数序列转换为复数序列
具体步骤：
（1）先使序列x（n）通过希尔伯特变换得到x1（n）
（2）对x1（n）进行FFT得到X1（K）
（3）X1（K）后乘旋转因子得到Y（K）
（4）Y(K)进行IFFT取实部，得到y（n），即最后位移后的序列

例：上题ph1=-pi/3；其他条件不变。

% 例：设x=cos(2*pi* f0* n/fs+ph1),其中fs采样频率为2000 Hz, 
% .信号频率f0为100 Hz,ph1 是信号的初始相角一pi/3,n是样点值,
% 信号长度N为40。希望通过位移使该信号的初始相角为0。
fs=2000;
f0=100;
N=40;
n=0:N-1;
ph1=-pi/3;
x=cos(2*pi*f0*n/fs+ph1);
x1=hilbert(x);                  %希尔伯特变换
X=fft(x1);
d=ph1*fs/f0/(2*pi);
Ex=exp(-1j*2*pi*n*d/N);
Y=X.*Ex;
y=ifft(Y);
y1=real(y);                     % 取的实部
Y1=fft(y1);                     % FFT
df = fs/N;                      % 计算频率间隔
nk=f0/df+1;                     % 信号在nk谱线上
A=abs(real(Y1(nk)))*2/N;        % 计算幅值
Theta=angle(Y1(nk));            % 计算初始相角
fprintf('A=%5.2f   Theta=%5.6f\n',A,Theta)

subplot(211);
plot(n,x,'k');
axis([0 N -1.1 1.1]); 
title('原始信号'); ylabel('幅值'); xlabel('样点'); 

subplot(212);
plot(n,y1,'k');
axis([0 N -1.1 1.1]); 
title('位移后信号的实部'); ylabel('幅值'); xlabel('样点'); 
% set(gcf,'color','w')

结果如下


初相角从x（n）的-pi/3位移后为0，并且没有误差。

---

总结

本篇文章主要围绕一个实数序列位移后进行DFT，然后IFFT后如何得到一个实数序列展开的，通过对于位移量是否为整数分为两种方法，对于d为整数则使用常规的方法，对于d不为整数则先进行希尔伯特变换然后，最后取IFFT结果的实部即可。