大多数学生可能对FFT变换后的共轭对称特性有所了解,但怎么利用共轭对称特性经IFFT转换为实数序列会些困难,因为在书中有几处会用到这样的特性,所以在这里先介绍共轭对称特性,再说明IFFT得实数序列。

FFT后的共轭对称特性为X(k)=X*(N-k)(其中k=0,1,…,N-1),严格地说,共轭对称轴在k=0处(但在MATLAB中下标或索引从1开始,此时共轭对称轴在k=1处)。以下用一个8点随机数的FFT来观察。程序清单为

x=randn(1,8);

X=fft(x);

X1=fftshift(X);

subplot 221; stem(real(X),'k');

subplot 222; stem(imag(X),'k');

subplot 223; stem(real(X1),'k');

subplot 224; stem(imag(X1),'k');

序列x经FFT后为X,又经fftshift后改为X1,其图如下。图中标出了Xr和Xi,标明X的实部和虚部,对X1也是一样标出实部和虚部。

为了说明方便起见,图中标出了Xr(1)~Xr(8)和 Xi(1)~Xi(8),在X1图中也把X的Xr(1)~Xr(8)和 Xi(1)~Xi(8)也标注出来。在X1的实部图中可以看出以Xr(1)相当于是一条对称轴的轴线,Xr(2)=Xr(8)、Xr(3)=Xr(7)、Xr(4)=Xr(6);在X1的虚部图中也可以看出Xr(1)相当于反对称的轴线,Xi(2)=-Xi(8)、Xi(3)=-Xi(7)、Xi(4)=-Xi(6)。在这8条谱线中X (1)和X(5)都不和其它谱线具有对称性的,X (1)永远可以做对称轴。

按K=0,1,…,N-1,对称关系式为X(k)=X*(N-k),若按K=1,…,N,对称关系式为X(k)=X*(N-k+2)。

在N是偶数时可以延伸把X(5),相当于X(N/2+1)也可以作为对称轴。从X的实部图中可看到Xr(4)=Xr(6) 、Xr(3)=Xr(7)、Xr(2)=Xr(8);从X的虚部图中可看到Xi(4)=-Xi(6) 、Xi(3)=-Xi(7)、Xi(2)=-Xi(8)。同上以X (1)为对称轴时得的结果一样。

对果我们要取局部的谱线进行反变换。为了保证反变换后还是实数序列,在取有限谱线时还要保留谱线的共轭对称的特性(这样选取IFFT后理论上为实数序列,但由于在MATLAB中计算机计算误差的存在,常常虚部还有数值,但都很小,所以在反变换后还用real函数取实部)。例如我们在正频率区间取X(2:3)两条谱线,为了保留谱线共轭对称性,按上的讨论在负频应取X(7:8)两条谱线。如果还要保留直流分量,则在正频率区间取X(1:3)叁条谱线。

当序列长N,取m条谱线(不包括直流分量),则取谱线有

X(1:m+1) ~ X(N-m+1:N),或X(end-m+1:end)

这是在pr3_1_1中用的关系式。若在m条谱线中已包含了直流分量,可把上关系变换一下:

X(1:m+1) ~ X(N-m+1:N),或X(end-m+1:end)

设j=m+1,则有m=j-1,代入上式得

X(1:j) ~ X(N-j+2:N),或X(end-j+2:end)

这就是在第9章中程序pr9_2_1用到的关系式。

mg01as.jpg

(34.87 KB, 下载次数: 17)

2013-11-27 10:10 上传

matlab 共轭对称分量,利用共轭对称特性经IFFT转换为实数序列相关推荐

  1. 【信号重构】经傅里叶逆变换(IFFT)后得到实数序列-含Matlab程序

    我们处理的大部分数据都是实数序列.设有实数序列x(n),把x(n)经FFT转到频域中为X(k),然后在频域中进行处理,处理完后经IFFT变成实数序列.本实例将说明如何经IFFT后得到实数序列. 一.实 ...

  2. 共轭反对称分量 matlab,共轭对称分量共轭反对称分量-Read.PPT

    共轭对称分量共轭反对称分量-Read * 序列傅里叶变换 第五讲 变换域分析 时域分析方法 变换域分析方法: 连续时间信号与系统 Laplace变换 Fourier变换 离散时间信号与系统 z变换 F ...

  3. Matlab如何进行利用离散傅里叶逆变换iDFT 从频谱恢复时域信号

    文章目录 1. 定义 2. 变换和处理 3. 函数 4. 实例演示 例1:单频正弦信号(整数周期采样) 例2:含有直流分量的单频正弦信号 例3:正弦复合信号 例4:含有随机干扰的正弦信号 例5:实际案 ...

  4. 实数序列频谱的共轭对称性(DFT与IDFT仿真实现)

    一.基础知识 1.傅里叶变换:通俗来讲,是以时间为自变量的信号与以频率为自变量的"频谱函数"之间的某种转换关系. DFT:即离散傅里叶变换,对离散序列进行傅里叶变换.设x(n)为长 ...

  5. 【数字信号处理】傅里叶变换性质 ( 序列傅里叶变换共轭对称性质 | x(n) 分解为实部序列与虚部序列 | 实部傅里叶变换 | 虚部傅里叶变换 | 共轭对称傅里叶变换 | 共轭反对称傅里叶变换 )

    文章目录 一.前置概念 1.序列对称分解定理 2.傅里叶变换 3.傅里叶变换的共轭对称分解 二.序列傅里叶变换共轭对称性质 0.序列傅里叶变换共轭对称性质 x(n) 分解为实部序列与虚部序列 x(n) ...

  6. 先验信息,先验概率,先验分布,后验概率,后验分布,共轭分布,共轭先验分布

    任何一个学科,最基本的就是概念,概念一定要清楚,清晰,否则概念都模棱两可的话,这之上的一切建筑都不牢固.很多概念可能长时间不使用就会变得模糊,所以在这里记录一下,输出是最好的记忆. 先验与后验的区别主 ...

  7. Matlab如何进行利用离散傅里叶变换DFT (快速傅里叶变换FFT)进行频谱分析

    文章目录 1. 定义 2. 变换和处理 3. 函数 4. 实例演示 例1:单频正弦信号(整数周期采样) 例2:单频正弦信号(非整数周期采样) 例3:含有直流分量的单频正弦信号 例4:正弦复合信号 例5 ...

  8. 【数字信号处理】傅里叶变换性质 ( 共轭对称与共轭反对称图像示例 | 实序列中共轭对称是偶对称 | 实序列中共轭反对称是奇对称 )

    文章目录 一.共轭对称与共轭反对称图像示例 1.共轭对称序列图示 2.共轭反对称序列图示 3.总结 一.共轭对称与共轭反对称图像示例 序列 x(n)=0.8nu(n)x(n) = 0.8^n u(n) ...

  9. 【数字信号处理】傅里叶变换性质 ( 共轭对称、共轭反对称 与 偶对称、奇对称关联 | 序列对称分解定理 )

    文章目录 一.共轭对称.共轭反对称 与 偶对称.奇对称关联 二.序列对称分解定理 证明过程 总结 一.共轭对称.共轭反对称 与 偶对称.奇对称关联 实序列 : 偶对称 : x(n)=x(−n)x(n) ...

最新文章

  1. 数字人民币实现可控匿名交易?产业升级离不开安全可信的“数字底座”
  2. Docker运行sonarqube-(代码质量检测平台)
  3. jq 实现发送验证码倒计时功能
  4. centos 6.1_64位系统 双网卡邦定
  5. python怎么把数据写入txt-python(如何将数据写入本地txt文本文件)
  6. 程序员必须尽早作打算
  7. Redis作为缓存服务器
  8. 计算机与交互式白板通过USB数据线,选购交互式电子白板注意哪些事项【详细介绍】...
  9. 三、spring中高级装配(1)
  10. IDEA导入JDK源码
  11. jquery 图片裁剪 java,jQuery实现图片上传和裁剪插件Croppie
  12. 网络地址转换--PAT配置
  13. 58欧氏空间05——对称变换和对称矩阵、实对称矩阵的标准形、正交相似
  14. Nwafu-OJ-1511 Problem 封闭的圆圈
  15. HDU 5454 Excited Database (2015年沈阳赛区网络赛E题)
  16. 移动端长按文本选择复制
  17. 详解ArcGIS中添加经纬网操作步骤
  18. oracle for update wait 解析
  19. 【mask2former】环境配置+cityscapes全景分割
  20. STM32CUBEIDE使用盗版jlink仿真器

热门文章

  1. 如何更新深度linux系统,深度操作系统 15.4 RC更新详情
  2. ajax入门学习(一)
  3. 出版社--图书管理系统
  4. 哈工大赛尔 | 事理图谱:事件演化的规律和模式
  5. marlin代码解析
  6. 试产机器发现元器件损毁 风险排查过程
  7. 用JSSE定制SSL连接
  8. 自动化测试 | 解决方案聚焦:如何使用PXI仪器进行高级数字测试
  9. HTML5的离线应用
  10. Wi-Fi6路由器芯片型号