数字信号处理课程设计说明书dft对称性的验证及以应用.docx
- 文档编号:993800
- 上传时间:2022-10-15
- 格式:DOCX
- 页数:16
- 大小:166KB
数字信号处理课程设计说明书dft对称性的验证及以应用.docx
《数字信号处理课程设计说明书dft对称性的验证及以应用.docx》由会员分享,可在线阅读,更多相关《数字信号处理课程设计说明书dft对称性的验证及以应用.docx(16页珍藏版)》请在冰豆网上搜索。
数字信号处理课程设计说明书dft对称性的验证及以应用
课程设计任务书
学生姓名:
专业班级:
电信0801
指导教师:
工作单位:
信息工程学院
题目:
DFT对称性的验证及以应用
初始条件:
具备数字信号处理的理论知识;
具备Matlab编程能力;
了解DFT的对称原理及应用;
提供编程所需要的计算机一台
要求完成的主要任务:
(包括课程设计工作量及其技术要求,以及说明书撰写等具体要求)
1、独立编写程序验证DFT的对称性
2、用DFT的对称性用一次FFT实现两个序列的FFT变换
3、完成符合学校要求的设计说明书
时间安排:
一周,其中3天程序设计,2天程序调试
指导教师签名:
年月日
系主任(或责任教师)签名:
年月日
摘要
有限长序列在数字信号处理是很重要的一种序列,反映它的"有限长"特点的一种有用工具是离散傅里叶变换(DFT)。
离散傅里叶变换除了作为有限长序列的一种傅里叶表示法在理论上相当重要之外,而且由于存在着计算离散傅里叶变换的有效快速算法,因而离散傅里叶变换在各种数字信号处理的算法中起着核心的作用。
而离散傅立叶变换的对称性,在求实序列的离散傅立叶变换中有重要作用。
可以实现一次DFT的计算得到两个序列DFT的高效算法,而DFT可以通过一次快速FFT变换来实现。
关键词:
DFT共轭对称性matlab
1DFT基础知识
1.1离散傅立叶变换(DFT)定义
有限长序列在数字信号处理是很重要的一种序列,当然可以用Z变换和傅里叶变换来研究它,但是,可以导出反映它的"有限长"特点的一种有用工具是离散傅里叶变换(DFT)。
离散傅里叶变换除了作为有限长序列的一种傅里叶表示法在理论上相当重要之外,而且由于存在着计算离散傅里叶变换的有效快速算法,因而离散傅里叶变换在各种数字信号处理的算法中起着核心的作用。
设x(n)是一个长度为M的有限长序列,则定义x(n)的N点离散傅里叶变换为:
正变换:
=DFT[]==
反变换:
=IDFT[]==
或
=RN(k)=RN(k)
x(n)=RN(n)=RN(n)
式中,N称为DFT变换区间长度,N≥M。
DFT隐含有周期性。
1.2复共轭序列的DFT
设是的复共轭序列,长度为N,则
(1)已知
=DFT[]
则
DFT[]=
且
(2)已知
=DFT[]
则
DFT[]=
1.3DFT的共轭对称性
DFT有对称性,但由于DFT中讨论的序列及其离散傅立叶变换均为有限长序列,且定义区间为0到N-1,所以这里的对称性是指关于N/2点的对称性。
下面讨论DFT的共轭对称性质。
1.3.1有限长共轭对称序列和共轭反对称序列
长度为的有限长序列,若满足
(1.1)
称序列为共轭对称序列,一般用来表示。
若满足
(1.2)
称序列为共轭反对称序列,一般用来表示
即
=,0≤n≤N-1
=-,0≤n≤N-1
当N为偶数时,把代入式(1.1)与式(1.2),得
(1.3)
(1.4)
式(1.3)与式(1.4)说明共轭对称序列与其共轭序列以成偶对称,共轭反对称序列与其共轭序列以成奇对称。
当N为奇数时,把代入式(1.1)与式(1.2),得
(1.6),(1.6)
式(1.5)与式(1.6)说明共轭对称序列与其共轭序列以成偶对称,共轭反对称序列与其共轭序列以成奇对称。
设一长度为的有限长序列,令
则有
(1.7)
这说明任一有限长序列,都表示成一个共轭对称序列与共轭反对称序列的和,在频域下同样有类似结论
(1.8)
式中(1.9)
(1.10)
1.3.2共轭对称性分析
(1)当x(n)为长度N的复数序列时,有
=]
=(1.11)
同理可得
(1.12)
即
式(1.11)和(1.12)说明复数序列实数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭对称分量;复书序列虚数部分的离散傅立叶变换是原来序列离散傅里叶变换的共轭反对称分量。
另一方面,由式(1.7)知有限长序列可分解为共轭对称分量与共轭反对称分量,即
=+
可得其离散傅立叶变换
=(1.13)
同理可得
=(1.14)
即
上面两式说明复序列共轭对称分量序列的离散傅立叶变换是原来序列离散傅立叶变换的实数部分;复序列共轭对称分量的离散傅立叶变换是原来序列离散傅立叶变换的虚数部分。
综上可得到有限长复序列的DFT的共轭对称性质如下
①将有限长序列x(n)分成实部与虚部,即:
则:
②将有限长序列x(n)分成共轭对称部分和共轭反对称部分,即
=+,
则:
(2)当x(n)为长度N的实数序列或纯虚数序列时,有
当x(n)为实序列时,则
又据)的对称性:
有
当x(n)为纯虚序列时,则
又据)的对称性:
有
离散傅立叶变换的对称性,在求实序列的离散傅立叶变换中有重要作用。
例如,有两个实数序列和,为求其离散傅立叶变换,可以分别用和作为虚部和实部构造一个复数序列x(n),求出x(n)的离散傅立叶变换,然后根据式(1.9)和(1.10)得到的共轭对称分量和,分别对应和,从而实现一次DFT的计算可得到两个序列DFT的高效算法。
而DFT可以通过一次快速FFT变换来实现。
2程序设计与分析
本次课设计分两个部分,一个是要验证N点的DFT的对称性,另一个是要用一次快速傅立叶变换FFT实现两个序列的DFT
2.1N点DFT对称性的验证
2.1.1程序流程图
由于函数ezplot只能画出既存在SymbolicMathToolbox中又存在于总matlab工具箱中的函数,而gedc(实信号分解为循环偶分量和循环奇分量)和dft(计算离散付利叶变换)仅存在SymbolicMathToolbox中,因此需要在自己的工作目录work下创建。
此后可以直接调用这些函数。
N点的DFT的对称性验证流程图如图2-1所示
开始
求x序列的共轭对称与反对称分量
画出共轭对称与反对称分量图形
求出X(K),Xep,Xop
画出real(X(K)),imag(X(K)),Xep,Xop的图形
Xep
结束
图2-1验证对称性流程图
输入x序列
n=0:
N-1
2.1.2程序编写与结果分析
首先在目录work下创建gedc的M文件,gedc的M文件是用来生成共轭对称分量与共轭反对称分量的,程序如下:
function[xec,xoc]=gedc(x);
N=length(x);
n=0:
(N-1);
xec=0.5*(x+x(mod(-n,N)+1));
xoc=0.5*(x-x(mod(-n,N)+1));
再是在目录work下创建dft的M文件,dft为离散傅立叶变换,程序如下:
function[Xk]=dft(xn,N);
n=[0:
1:
N-1];
k=[0:
1:
N-1];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=xn*WNnk;
主程序:
(1)N=12,序列为x=[2.501.6-3-221.6-3-144.5-2]的程序设计与结果分析
程序:
figure
(1)
n=0:
11;x=input('请输入序列x=');
[xep,xop]=gedc(x);
subplot(2,1,1);
stem(n,xep);
title('共轭对称分量')
xlabel('n');
ylabel('xep');
axis([-0.5,12.5,-3,4]);
subplot(2,1,2);
stem(n,xop);
title('共轭反对称分量');
xlabel('n');
ylabel('xop');
axis([-0.5,12.5,-4,4]);
figure
(2)
X=dft(x,12);
Xep=dft(xep,12);
Xop=dft(xop,12);
subplot(2,2,1);
stem(n,real(X));
axis([-0.5,12.5,-10,10]);
title('real(X)');
xlabel('k');
subplot(2,2,2);
stem(n,imag(X));
axis([-0.5,12.5,-17,17]);
title('imag(X)');
xlabel('k');
subplot(2,2,3);
stem(n,Xep);
axis([-0.5,12.5,-10,10]);
title('DFT[xep(n)]');
xlabel('k');
subplot(2,2,4);
stem(n,imag(Xop));
axis([-0.5,12.5,-17,17]);
title('DFT[xop(n)]');
xlabel('k');
结果:
图2-2共轭对称分量与共轭反对称分量
图2-3对称性的验证图形
分析:
从图2-3可以看出复数序列实数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭对称分量;复数序列虚数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭反对称分量。
复序列共轭对称分量序列的离散傅立叶变换是原来序列离散傅立叶变换的实数部分;复序列共轭反对称分量的离散傅立叶变换是原来序列离散傅立叶变换的虚数部分。
从而验证了DFT的对称性。
(2)N=14,序列为x=(1.2).^n的程序设计与结果分析
程序:
figure
(1)
n=0:
13;x=input('请输入序列x=');
[xep,xop]=circevod(x);
subplot(2,1,1);
stem(n,xep);
title('共轭对称分量')
xlabel('n');
ylabel('xep');
axis([-0.5,14.5,0,7]);
subplot(2,1,2);
stem(n,xop);
title('共轭反对称分量');
xlabel('n');
ylabel('xop');
axis([-0.5,14.5,-6,6]);
figure
(2)
X=dft(x,14);
Xep=dft(xep,14);
Xop=dft(xop,14);
subplot(2,2,1);
stem(n,real(X));
axis([-0.5,14.5,-13,6]);
title('real(X)');
xlabel('k');
subplot(2,2,2);
stem(n,imag(X));
axis([-0.5,14.5,-25,25]);
title('imag(X)');
xlabel('k');
subplot(2,2,3);
stem(n,Xep);
axis([-0.5,14.5,-13,6]);
title('DFT[xep(n)]');
xlabel('k');
subplot(2,2,4);
stem(n,imag(Xop));
axis([-0.5,14.5,-25,25]);
title('DFT[xop(n)]');
xlabel('k');
结果:
图2-4共轭对称分量与共轭反对称分量
图2-5对称性的验证图形
分析:
从图2-5可以看出复数序列实数部分的离散傅立叶变换是原来序列离散
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 课程设计 说明书 dft 对称性 验证 应用