下载本文档
版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领
文档简介
1、matlab的 fft 函数2010-04-26 22:16相关语法：y = fft(x)y = fft(x,n)y = fft(x,dim)y = fft(x,n,dim)定义如下：相关的一个例子：fs = 1000;%采样频率t = 1/fs;%采样时间l = 1000;%总的采样点数t = (0:l-1)*t;%时间序列（时间轴）%产生一个幅值为0.7 频率为50hz正弦 +另外一个信号的幅值为1 频率为120hz的正弦信号x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);y = x + 2*randn(size(t);%混入噪声信号plot(fs*t(1
2、:50),y(1:50)%画出前 50 个点title(signal corrupted with zero-mean random noise)xlabel(time (milliseconds)nfft = 2nextpow2(l); %y = fft(y,nfft)/l;%求得最接近总采样点的2n, 这里应该是 210=1024进行 fft变换（除以总采样点数，是为了后面精确看出原始信号幅值）f = fs/2*linspace(0,1,nfft/2+1);%频率轴 ( 只画到fs/2即可，由于y 为实数，后面一半是对称的)% 画出频率幅度图形，可以看出 50hz 幅值大概 0.7,120
3、hz 幅值大概为 1. plot(f,2*abs(y(1:nfft/2+1)title(single-sided amplitude spectrum of y(t) xlabel(frequency (hz)ylabel(|y(f)|)ps：前段时间，只是为了自己看明白， 没有管太多，刚才上空间发现几位的留言，感觉惭愧。其实要掌握 matlab 中 fft 的用法，主要有两点注意的地方：1、从公式上看， matlab 的 fft 序号是从 1 到 n，但是绝大多数教材上是从0 到n-1。2、y=fft （x）之后，这个 y 是一个复数，它的模值应该除以（ length(x)2），才能得到各个频率信号实际幅值。（注：专业文档是经验性极
温馨提示
- 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
评论
0/150
提交评论