Matlab中FFT函数的源代码示例
来源:互联网 发布:记账软件onlinedown 编辑:程序博客网 时间:2024/05/29 09:24
clear
fs=100;N=128; %采样频率和数据点数
n=0:N-1;t=n/fs; %时间序列
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); %信号
y=fft(x,N); %对信号进行快速Fourier变换
mag=abs(y); %求得Fourier变换后的振幅
f=n*fs/N; %频率序列
subplot(2,2,1),plot(f,mag); %绘出随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('N=128');grid on;
subplot(2,2,2),plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('N=128');grid on;
%对信号采样数据为1024点的处理
fs=100;N=1024;n=0:N-1;t=n/fs;
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); %信号
y=fft(x,N); %对信号进行快速Fourier变换
mag=abs(y); %求取Fourier变换的振幅
f=n*fs/N;
subplot(2,2,3),plot(f,mag); %绘出随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('N=1024');grid on;
subplot(2,2,4)
plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('N=1024');grid on;
运行结果为:
fs=100Hz,Nyquist频率为fs/2=50Hz。整个频谱图是以Nyquist频率为对称轴的。并且可以明显识别出信号中含有两种频率成分:15Hz和40Hz。由此可以知道FFT变换数据的对称性。因此用FFT对信号做谱分析,只需考察0~Nyquist频率范为内的福频特性。若没有给出采样频率和采样间隔,则分析通常对归一化频率0~1进行。另外,振幅的大小与所用采样点数有关,采用128点和1024点的相同频率的振幅是有不同的表现值,但在同一幅图中,40Hz与15Hz振动幅值之比均为4:1,与真实振幅0.5:2是一致的。为了与真实振幅对应,需要将变换后结果乘以2除以N。
- Matlab中FFT函数的源代码示例
- MATLAB实现的FFT函数
- Matlab中fft函数的用法及关键问题详解
- Matlab中fft函数的用法及关键问题详解
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的使用方法
- MATLAB中FFT的应用
- 查看matlab中函数源代码的方法
- fedora 14 下的网络相关设置方法
- first blog
- 现在为了让程序兼容vista,我们需要给程序加上一个manifest文件,可是有时候vc6它就会弹出Resource Compiler Error RC2170 的错误,可能是vc6一个bug,我个人对vc6还是很喜欢的,不过经常会有编译死掉或出错的问题,
- servlet 、GenericServlet 、HttpServlet 详解
- DLL中键盘钩子失效的问题解决
- Matlab中FFT函数的源代码示例
- Freemarker原理解析
- 几种排序方案
- JavaBean_重构_内省操作
- STL容器类简介
- 知识漏洞
- zoj 1392
- 《Java SE 6.0 编程指南》学习笔记
- Altium Designer多图纸原理图设计方法探讨