数字信号处理课程设计说明书-dft对称性的验证及以应用Word下载.doc
《数字信号处理课程设计说明书-dft对称性的验证及以应用Word下载.doc》由会员分享,可在线阅读,更多相关《数字信号处理课程设计说明书-dft对称性的验证及以应用Word下载.doc(20页珍藏版)》请在冰点文库上搜索。
1DFT基础知识 1
1.1离散傅立叶变换(DFT)定义 1
1.2复共轭序列的DFT 1
1.3DFT的共轭对称性 2
1.3.1有限长共轭对称序列和共轭反对称序列 2
1.3.2共轭对称性分析 3
2程序设计与分析 6
2.1N点DFT对称性的验证 6
2.1.1程序流程图 6
2.1.2程序编写与结果分析 7
2.2用一次FFT实现两个序列的DFT 13
2.2.1程序流程图. 13
2.2.2程序编写与结果分析 13
3课程设计心得体会 16
参考文献 17
摘要
有限长序列在数字信号处理是很重要的一种序列,反映它的"
有限长"
特点的一种有用工具是离散傅里叶变换(DFT)。
离散傅里叶变换除了作为有限长序列的一种傅里叶表示法在理论上相当重要之外,而且由于存在着计算离散傅里叶变换的有效快速算法,因而离散傅里叶变换在各种数字信号处理的算法中起着核心的作用。
而离散傅立叶变换的对称性,在求实序列的离散傅立叶变换中有重要作用。
可以实现一次DFT的计算得到两个序列DFT的高效算法,而DFT可以通过一次快速FFT变换来实现。
关键词:
DFT共轭对称性matlab
1DFT基础知识
1.1离散傅立叶变换(DFT)定义
有限长序列在数字信号处理是很重要的一种序列,当然可以用Z变换和傅里叶变换来研究它,但是,可以导出反映它的"
设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[]=
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:
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)
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);
共轭反对称分量'
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]);
real(X)'
k'
subplot(2,2,2);
stem(n,imag(X));
axis([-0.5,12.5,-17,17]);
imag(X)'
subplot(2,2,3);
stem(n,Xep);
axis([-0.5,12.5,-10,10]);
DFT[xep(n)]'
subplot(2,2,4);
stem(n,imag(Xop));
DFT[xop(n)]'
结果:
图2-2共轭对称分量与共轭反对称分量
图2-3对称性的验证图形
分析:
从图2-3可以看出复数序列实数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭对称分量;
复数序列虚数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭反对称分量。
复序列共轭对称分量序列的离散傅立叶变换是原来序列离散傅立叶变换的实数部分;
复序列共轭反对称分量的离散傅立叶变换是原来序列离散傅立叶变换的虚数部分。
从而验证了DFT的对称性。
(2)N=14,序列为x=(1.2).^n的程序设计与结果分析
13;
[xep,xop]=circevod(x);
axis([-0.5,14.5,0,7]);
axis([-0.5,14.5,-6,6]);
X=dft(x,14);
Xep=dft(xep,14);
Xop=dft(xop,14);
axis([-0.5,14.5,-13,6]);
axis([-0.5,14.5,-25,25]);
imag(X)'
图2-4共轭对称分量与共轭反对称分量
图2-5对称性的验证图形
从图2-5可以看出复数序列实数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭对称分量;
2.2用一次FFT实现两个序列的DFT
2.2.1程序流程图
一次快速傅立叶变换FFT实现两个序列的DFT流程图如图2-5所示。
输入x=+j
调用fft函数
得到和
图2-6一次FFT变换实现两序列的DFT
2.2.2程序编写与结果分析
x1=input('
请输入序列x1='
x2=input('
请输入序列x2='
N=input('
请输入N='
x=x1+j*x2;
X=fft(x,N);
k=0:
N-1;
c=conj(X);
Xep=0.5*(X+c(mod(-k,N)+1));
Xop=-j*0.5*(X-c(mod(-k,N)+1));
X1=Xep
X2=Xop
stem(k,X1);
X1'
axis([-0.5,7.5,-10,40]);
stem(k,X2);
X2'
当运行程序时,会出现提示,按提示输入x1=[1352463526],x2=[12345-5-4-3-2-1],N=10,程序运行结果如下:
X1和X2分别为x1,x2的离散傅立叶变换,X1和X2的图形如图2-7所示
图2-7x1和x2的离散傅立叶变换
当直接调用dft时,程序运行结果和上面的是相同的,从而实现了用一次FFT实现了两个序列的DFT。
3课程设计心得体会
本次课程设计主要是运用本学期所学到的数字信号处理的基础知识来设计一个符合要求的matlab程序来进行DFT对称性的验证以及应用,本次设计不仅要求我们要掌握数字信号处理课程的基础知识,还要求我们对matlab编程有深刻的理解和掌握。
用新的语言去解决工程问题根本不需要先掌握某一门语言,有效的方法是先了解那门语言的一些基本函数,然后熟悉界面,就可以开始编了。
拿到一个课题,不要急于坐在电脑前开始编程,因为当你坐在电脑前都不知道该干什么时,你就是对课题了解得不够。
你首先需要的是透彻分析课题,把你要解决的问题写下来和列出各种可能情况。
接下来,就考虑看用什么样的算法去解决,等到这一切都定下来后就可以开始着手编程了,如果你不熟悉语言的话,过程中会碰到很多问题,例如,不知道用什么样的函数去实现,这时你就可以根据实际情况去找资料,看帮助文件。
等到程序完成后,调试是非常关键的一步,看到出错没有关系,利用debug去分析,相信一定可以找到问题的所在,然后逐个改正。
经过这么一次以后,你会发现,你已经对本来不熟悉的语言已经有了很大的了解,而且可以解决实际问题,这样你就不会觉得学语言很枯燥,而是觉得很有趣。
相信经过多次反复,一定可以很快掌握一门语言的基本技巧。
虽然我现在已经初步学会了如何设计符合要求的matlab程序,但是离真正能够利用已学的知识自由设计matlab程序的还有很长的一段的距离。
课设的这段时间我确实受益匪浅,不仅是因为它发生在特别的实践,更重要的是我的专业知识又有了很大的进步,进步总是让人快乐的。
参考文献
[1]刘泉等主编.信号与系统.高等教育出版社,2006年
[2]赵红怡等主编.数字信号处理及其MATLAB实现.化学工业出版社,2007年
[3]薛定宇主编.基于MATLAB的系统仿真技术与应用.清华大学出版社,2002年
[4]程佩青主编.数字信号处理教程(第三版).清华大学出版社,2007年
[5]张德丰主编.MATLAB数字信号处理与应用.清华大学出版社,2010年
[6]张德丰主编.详解MATLAB数字信号处理.电子工业出版社,2010年
[7]李正周主编.Matlab数字信号处理与应用.清华大学出版社,2008年
20