数字信号处理实验二:DFT的共轭对称性及应用
实验目的
1、掌握实序列的DFT共轭对称性的特点;
2、学习应用实序列DFT的共轭对称性构建频域序列以保证时域序列为实数的方法;
实验原理
1.DFT的共轭对称性
其中:
2.有限长实序列的DFT的共轭对称性
设为长度为N的有限长实序列,则是圆周共轭对称序列:。此对称性又可表示为: 。
实验内容
1、试利用 DFT 的共轭对称性,设计两种高效算法,使得计算一个N 点DFT 来得到两个实序列的 N 点 DFT。设:
(1)算法一:令,计算其16点离散傅里叶变换(提示:取 N=16时,。
由分别求出和,并通过和 求IDFT的方法验证得出的结果是否正确。
clc,clear,close all
N=16;
n=0:1:N-1;
x1=cos(pi/4*n);
x2=sin(pi/8*n);
xn=x1+x2;
Xk=fft(xn,16);
X1=real(Xk);
X2=imag(Xk);
x11=ifft(X1,16);
x22=ifft(1i*X2,16);subplot(611)
stem(n,X1)
xlabel('k');
ylabel('real(X1)');subplot(612)
stem(n,X2)
xlabel('k');
ylabel('real(X2)');subplot(613)
stem(n,x1);
xlabel('n');
ylabel('x1');subplot(614)
stem(n,x11);
xlabel('n');
ylabel('x1恢复');subplot(615)
stem(n,x2);
xlabel('n');
ylabel('x2');subplot(616)
stem(n,x22);
xlabel('n');
ylabel('x2恢复');
实验结论 1-1:、与的关系?
答:
(2)算法二:令,重复(1)。
clc,clear,close all
N=16;
n=0:1:N-1;
k=0:1:N-1;
x1=cos(pi/4*n);
x2=sin(pi/8*n);
xn=x1+1i*x2;
Xk=fft(xn,16);
Xkx=conj(Xk);%取共轭
X1=1/2*(Xk+[Xkx(1),fliplr(Xkx(2:16))]);
X2=1/2*(Xk-[Xkx(1),fliplr(Xkx(2:16))]);
x11=ifft(X1,16);
x22=ifft(X2,16)*(-1i);subplot(611)
stem(n,real(X1))
xlabel('k');
ylabel('real(X1)');subplot(612)
stem(n,real(X2))
xlabel('k');
ylabel('real(X2)');subplot(613)
stem(n,x1);
xlabel('n');
ylabel('x1');subplot(614)
stem(n,x11);
xlabel('n');
ylabel('x1恢复');subplot(615)
stem(n,x2);
xlabel('n');
ylabel('x2');subplot(616)
stem(n,x22);
xlabel('n');
ylabel('x2恢复');
实验结论 1-2:、与的关系?
答:
2.有限长实序列的DFT的共轭对称性
由有限长实序列的DFT的共轭对称性可知,频域成共轭对称的序列作IDFT
后为实序列,而实数的发送可以大大简化发送设备。OFDM 正是利用这一特性来保证发往信道的序列为实数序列的。
按要求编程完成以下内容:
设 XK_in 为一频域复数序列,XK_in=[1+j,-3-j,-3+3*j,-1-3*j];
试利用实序列的 DFT 的共轭对称性公式,将频域序列
XK_in 扩展成共轭对称形式Xk,以保证其对应的时域序列xn =ifft(Xk,16)为实数序列。
(1)求频域序列Xk;并给出Xk的实部与虚部图;
clc,clear,close all
format compact
N=16;
n=0:1:N-1;
k=0:1:7;
XK_in=[1+1i,-3-1i,-3+3*1i,-1-3*1i];
XKf=conj(fliplr(XK_in));
Xk=[0,XK_in,0,0,0,0,0,0,0,XKf];subplot(211)
stem(n,real(Xk));
xlabel('k');
ylabel('real(Xk)');
subplot(212)
stem(n,imag(Xk));
xlabel('k');
ylabel('imag(Xk)');
实验结论 2-1:说明Xk的实部与虚部各有何特点;
答:实部关于N/2偶对称,虚部关于N/2奇对称。
(2)求xn =ifft(Xk,16);
clc,clear,close all
format compact
N=16;
n=0:1:N-1;
k=0:1:7;
XK_in=[1+1i,-3-1i,-3+3*1i,-1-3*1i];
XKf=conj(fliplr(XK_in));
Xk=[0,XK_in,0,0,0,0,0,0,0,XKf];
xn =ifft(Xk,16)subplot(211)
stem(n,real(xn));
xlabel('n');
ylabel('real(xn)');
subplot(212)
stem(n,imag(xn));
xlabel('n');
ylabel('imag(xn)');
实验结论 2-2:说明xn是否为实数序列,可以用xn的实部与虚部图来说明。
答:由上图可知,xn虚部恒为0,为实数序列。
实验思考
1.对序列x(n),如何通过计算N/2点DFT而得到N点DFT?
答:
(1)对于序列x(n),可由基2FFT算法,对x(n)进行奇偶序列划分来求N/2点DFT,从而求出N点DFT。
(2)特别的,若x(n)为实序列,则由,当N=偶数时,只需计算X(k)的前面N/2+1点,而N=奇数时,只需计算X(k)的前面(N+1)/2点,从而计算出N点DFT。
数字信号处理实验二:DFT的共轭对称性及应用相关推荐
- 数字信号处理实验(二)
实验目的 音频信号实验 产生两段音频信号,将声音信号进行傅立叶变换,将时域信号转为频域信号,并分别查看两段信号变化后的幅频特性,然后将两段频域信号线性叠加,合成新的频域信号,最后将此频域信号进行反傅立 ...
- 数字信号处理实验二:数字语音信号与频谱分析
文章目录 一.实验目的 二.实验过程与结果 1.数字音频信号的频谱分析 2. 数字音频信号的频谱分析 3. 数字和弦音频信号的生成 三.结果分析与实验结论 1.读数字音频信号的频谱分析 2.数字和弦音 ...
- 数字信号处理 | 实验二 MATLAB z换和z逆变换分析+求解差分方程+求解单位冲击响应+求解幅频相频特性曲线+求解零极点
1.实验目的 (1)掌握离散时间信号的z变换和z逆变换分析 (2)掌握MATLAB中利用filter函数求解差分方程: (3)掌握MATLAB中利用impz函数求解单位冲击响应h(n); (4)掌握M ...
- Matlab验证dtft共轭性质,数字信号处理实验4重点.docx
深 圳 大 学 实 验 报 告 课程名称: 数字信号处理实验 实验项目名称: 实验4 学院: 信息工程学院 专业: 电子信息工程 指导教师: 陈佳义 报告人: 学号: 班级: 实验时间: 11.12 ...
- 奇数点偶数点fft的matlab,电子科大 数字信号处理实验2_FFT的实现
电 子 科 技 大 学 实 验 报 告 学生姓名:Shrimp 学 号: 指导教师: 一.实验室名称:数字信号处理实验室 二.实验项目名称:FFT 的实现 三.实验原理: 一.FFT 算法思想: 1. ...
- matlab fft谱分析实验报告,数字信号处理实验报告-FFT算法的MATLAB实现.doc
数字信号处理实验报告-FFT算法的MATLAB实现.doc 数字信号处理 实验报告实验二FFT算法的MATLAB实现一.实验目的通过本实验的学习,掌握离散傅立叶变换的理论,特别是FFT的基本算法以及其 ...
- 数字信号处理实验matlab版答案刘舒帆,数字信号处理实验(MATLAB版) 刘舒帆,费诺,陆辉 西安电子科技大学出版社 9787560620060...
商品描述: 基本信息 书名:数字信号处理实验(MATLAB版) 原价:31.00元 作者:刘舒帆,费诺,陆辉 著 出版社:西安电子科技大学出版社 出版日期:2013-7-1 ISBN:97875606 ...
- 离散卷积的c语言编程实验,数字信号处理实验一离散卷积c语言编程.ppt
数字信号处理实验一离散卷积c语言编程 实验一 离散卷积的C语言编程实验 DSP实验室 2005 实验性质 综合设计性实验 实验目的 1 了解和认识常用的各种信号: 2 掌握卷积的定义和计算方法: 3 ...
- matlab数字信号处理实验报告,数字信号处理实验报告一 离散信号及其MATLAB实验...
数字信号处理 离散信号及其MATLAB实验 南昌航空大学实验报告 2012 年 04 月 06 日 课程名称: 数字信号处理 实验名称: 离散信号及其MATLAB实现 班级: 090423班 学号: ...
最新文章
- 腾讯云年度最强技术大会召开在即,这次只谈技术和代码
- C++ 编译 找不到标识符 问题
- Exception in thread “main“ java.lang.NoClassDefFoundError: org/apache/ibatis/io/Resources
- 原始Junit测试Spring的问题
- 没有垃圾回收的JVM
- ethtool用法 linux_Linux命令之Ethtool用法详解
- ReentrantLock实现原理
- iPhone系统更新提醒怎么关闭?
- [oracle]数据库字段类型char(n)的问题
- 物理、化学实验仪器与设备
- C# 使用RabbitMQ(二)安装
- Ubuntu22.04安装网易云音乐实操笔记
- 心脏除颤器的工作原理
- 下软件,就靠这几个网站
- 高级会计职称计算机考什么,会计高级职称考哪些科目
- 正常图片转灰色图片——灰度图
- 如何取消笔记本Fn键
- [转载]document.readyState
- 计算机里的文档怎么设置密码,文件夹怎么设置密码,教您如何给电脑上文件夹设置密码...
- 游戏联网必备: 国内外实时对战服务详细对比
热门文章
- ubuntu文件系统字体底纹含义
- 【Filecoin周报-62】倒计时1天,太空竞赛准备怎么样了?
- Flutter-防京东商城项目-提交订单、去支付页面制作-44
- Mac 使用 之 Permission denied XXX
- 中止执行后超过2年_在法院中止执行后怎么办
- Unity3d基础教程7-灯光(灯光阴影参数设置及lightmapping制作)
- ruby自动注册163邮箱
- 华为云迁移工具推荐最佳实践:Xen虚拟化迁移到华为云
- 李德毅院士:迭代的智能——从薛定谔、图灵和维纳谈开去
- python123程序设计题说句心里话a_C程序设计基础(2019年春)-中国大学mooc-试题题目及答案...