- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
一、四种信号傅立叶表示:
1. 周期为的连续时间周期信号:
频谱特点: 离散非周期谱
2. 连续时间非周期信号
频谱特点: 连续非周期谱
3. 离散非周期信号
频谱特点: 周期为2(的连续谱
4. 周期为N 的离散周期信号
频谱特点:周期为N的离散谱
二、离散时间傅立叶变换定义与计算
如果绝对可加的,即,则其离散傅立叶变换定义为:
的 离散时间傅立叶逆变换可以表示为:
离散时间傅立叶变换有两个重要的特性:
周期性:离散时间傅立叶变换是的周期函数,其周期为。
三、离散时间傅立叶变换的特性:
(1)线性:
(2) 时移:
(3) 频移:
(4) 反转:
(5) 实序列的对称性:
(6) 卷积:
(7) 相乘:
(8) 能量:
详细性质参见课本
四、离散傅立叶变换
由傅立叶分析知,一个周期函数可以由其各个谐波分量的线性组合得到,这就是离散傅立叶级数(DFS),将其推广至有限持续时间序列,并产生一个新的变换,这就是傅立叶变换(DFT)
离散傅立叶变换的定义如下所示:
离散傅立叶变换的函数实现:
function [Xk]=dft(xn,N)
%计算离散傅立叶变换
%Xk为在0=nN-1间的DFT系数数组
%xn为N点有限持续时间序列
%N为DFT的长度
n=0:N-1; %n的行向量
k=0:N-1; %k的行向量
WN=exp(-j*2*pi/N); %Wn因子
nk=n*k; %产生一个含nk值的N乘N维矩阵
WNnk=WN.^nk
Xk=xn*WNnk;%DFT系数行向量
五、快速傅立叶变换
当序列长度N很大时,直接计算序列的离散傅立叶变换的计算量非常大,因此引入了快速傅立叶变换(FFT),所有这些高效算法合起来统称为快速傅立叶变换算法。
Matlab中自带有FFT函数,试用help函数察看其具体使用方法。
例题1:已知带有测量噪声的信号其中f1=50Hz,f2=120Hz,为均值为0的随机信号,数据采样频率为1000Hz,数据点数N=1024,绘制信号的频谱图和无噪声信号的频谱图。
fs=1000;
N=1024;
n=0:N-1;
t=n/fs;
f1=50;f2=120;
%信号没有噪声情况
x=sin(2*pi*f1*t)+sin(2*pi*f2*t);
y=fft(x,N);
mag=abs(y);
f=(0:length(y)-1)*fs/length(y);
subplot(212);
plot(f(1:N/2),mag(1:N/2));
xlabel(Frequence(Hz));
ylabel(Magnitude);
title(N=1024 没有噪声);
grid
%信号有噪声情况
x=x+2*randn(1,length(t));
y=fft(x,N);
mag=abs(y);
f=(0:length(y)-1)*fs/length(y);
subplot(211);
plot(f(1:N/2),mag(1:N/2));
xlabel(Frequence(Hz));
ylabel(Magnitude);
title(N=1024 有噪声);
grid
原创力文档


文档评论(0)