利用N点复数序列求2个N点实数序列的快速傅里叶变换

伟大的马师兄

已于 2023-10-17 23:13:28 修改

阅读量705

点赞数 1

分类专栏：傅里叶分析文章标签：傅里叶分析快速傅里叶变换离散傅里叶变换

于 2023-10-14 23:54:15 首次发布

本文链接：https://blog.csdn.net/qq_38967414/article/details/133833483

版权

傅里叶分析专栏收录该内容

3 篇文章 4 订阅

订阅专栏

一、方法简介

通过一个 $N$ 点复数序列求出两个 $N$ 点实数序列的离散傅里叶变换，进一步提升快速傅里叶变换的效率。

二、方法详解

$x(n)$ 和 $y(n)$ 是实数序列，且长度都为 $N$ ，定义复数序列：

$z(n)=x(n)+jy(n)$ ， $0\leq n\leq N-1$

则序列 $x(n)$ 和 $y(n)$ 可表示为：

$x(n)=\frac{1}{2}(z(n)+z^{*}(n))$

$y(n)=\frac{1}{j2}({z(n)-z^{*}(n)})$

$z(n)$ 的离散傅里叶变换可表示为：

$Z(k)=X(k)+jY(k)$

序列 $x(n)$ 和 $y(n)$ 的离散傅里叶变换可表示为：

$X(k)=\frac{1}{2}({DFT[z(n)]+DFT[z^{*}(n)]})$

$Y(k)=\frac{1}{j2}(DFT[z(n)]-DFT[z^{*}(n)])$

因为 $DFT[x^{*}(n)]=X^{*}(N-k)$ ，所以：

$X(k)=\frac{1}{2}({Z(k)+Z^{*}(N-k)})$

$Y(k)=\frac{1}{j2}({Z(k)-Z^{*}(N-k)})$

三、MatLab代码示例

x = [8 2 3 4 5 8];
y = [3 9 2 7 6 9];

x_fft = fft(x);
y_fft = fft(y);

z = x+1i*y;

Z = fft(z);
N = length(Z);

X = zeros(1,N);
Y = zeros(1,N);

R = real(Z);
I = imag(Z);

X(1) = R(1);
Y(1) = I(1);

for j = 2 : N
    X(j) = 1/2*(Z(j)+conj(Z(N+2-j)));
    Y(j) = 1/2i*(Z(j)-conj(Z(N+2-j)));
end

运行结果：