python中对序列FFT后结果不是共轭对称的问题

我对雷克子波进行了fft后，发现结果并不是共轭对称的。

这导致我把fft后结果取前一半区间，然后共轭对称到后一半，再反变换回去的结果就有点异常了。如果直接把fft后的全序列做ifft，结果就是对的。

第一张图是fft后，把前一半区间取共轭对称，再和后一半区间相减，可以看到误差还是挺大的，而且越到后面越大。

我在别的论坛看到有人用matlab也遇到过这个问题。

请问这是什么原因？

附上代码：

import numpy as np
from scipy import fft, ifft
import matplotlib.pyplot as plt
import math

f = 50
t = np.arange(-0.512, 0.512, 0.001)
n = int(len(t)/2)

rick = (1 - 2 * (math.pi * f * t)**2) * math.e ** (-1 * (math.pi * f * t)**2)

yy = fft(rick)

y1 = yy[: n]
y2 = np.conj(y1)[:: -1]
y3 = yy[n:]
yyy = np.append(y1, y2)

plt.plot(ifft(yy))
plt.plot(ifft(yyy))
plt.show()

plt.plot(y2-y3)
plt.show()

误差。
直接反变换和共轭对称后再反变换的结果。

python中对序列FFT后结果不是共轭对称的问题相关推荐

python中的序列总结：列表，元组，字符串
首先python中的序列有列表,元组,字符串. 列表List是python里面最基本的数据结构.序列中每个元素都是从索引(下标)从0开始,依次叠加. List操作的方法很多,只能熟悉基本常用的这个方法 ...
python 中的序列
python 中的序列序列概述 Python 内置了多种序列,比如列表.元组和字符串. 列表和元组的主要不同在于:列表是可以修改的,而元组不可以. 在需要处理一系列值时,序列很有用.在数据库中,你可 ...
Python中的序列操作
Python中的序列操作分类: python undefined 官方手册:https://docs.python.org/3.7/library/stdtypes.html#sequence-ty ...
python中的序列类型数据结构元素的切片操作_浅析python中的分片与截断序列
浅析python中的分片与截断序列序列概念在分片规则里list.tuple.str(字符串)都可以称为序列,都可以按规则进行切片操作切片操作注意切片的下标0代表顺序的第一个元素,-1代表倒序的 ...
python中的序列类型数据结构元素的切片操作_PythonI/O进阶学习笔记_4.自定义序列类(序列基类继承关系/可切片对象/推导式)...
前言: 本文代码基于python3 Content: 1.python中的序列类分类 2. python序列中abc基类继承关系 3. 由list的extend等方法来看序列类的一些特定方法 4. l ...
已知一棵二叉树的中序序列和后序序列，写一个建立该二叉树的二叉链表存储结构的算法...
已知一棵二叉树的中序序列和后序序列,写一个建立该二叉树的二叉链表存储结构的算法 #define N 10 //二叉树节点的个数 char postorderstr[]={};//后序序列 char i ...
PAT甲级1020 Tree Traversals：[C++题解]树的遍历、由中序序列和后序序列递归建树
文章目录题目分析题目链接题目分析题意重述:给定一棵二叉树的后序遍历序列和中序遍历序列,让求层次遍历的序列. 分析: 后序遍历:先左子树.右子树 ,最后再遍历根结点. 中序遍历:先左子树,再根 ...
Python 中的序列类型支持哪些公共操作
一.序列类型支持哪些公共操作 Python 中的序列类型分为:可变序列和不可变序列,都支持我们序列类型的所有操作可变序列:列表不可变序列:字符串和元祖 1.支持通过数字索引取值 a = (&quo ...
7-10 先序序列创建二叉树，输出先序序列、中序序列、后序序列并输出叶子结点数 (10 分)
7-10 先序序列创建二叉树,输出先序序列.中序序列.后序序列并输出叶子结点数 (10 分) 对于给定的二叉树,输出其先序序列.中序序列.后序序列并输出叶子结点数. 输入格式: 二叉树的先序遍历序列. ...

python中对序列FFT后结果不是共轭对称的问题

python中对序列FFT后结果不是共轭对称的问题相关推荐

最新文章

热门文章