一.基本内容: 基于FFT的功率谱分析程序设计与应用
1) 对一个人为产生的信号,采用FFT变换方法进行功率谱分析
已知信号x(n)=100.0*COS(2*3.14*SF*n/FS) +50.0*COS(4*3.14*SF*n/FS) ,n=0,1,2 ……N-1
式中: N---采样点数,必须为2x
SF---信号频率,SF=10 Hz
FS---采样频率
其FFT变换结果X(k)可用下面提供的FFT子程序求出,计算功率谱的公式为:
G(k)=2(XR(k)2 +XI(k)2)/N , k=0,1,2 ……N/2-1
式中: XR(k)--- X(k)的实部
XI(k)--- X(k)的虚部
请用VB,VC或C++Builder编译器编程,或采用MATLAB计算。处理结果(时域波形或频谱)采用窗口显示,并且拷贝至作业稿中。
注意:除FFT分析程序可参考附件所附的程序外,其余程序必须自己设计(小波或EMD除外),后面限选题目要求与此相同。
该信号的处理结果如下:
其matlab程序如下:
SF=10.0;%信号频率
N=1024;%采样点数
n=0:1023;%时间轴N个点
FS=300;%采样频率
t=n/FS;
x=100.0*cos(2*3.14*SF*t) +50.0*cos(4*3.14*SF*t);
figur
原创力文档

文档评论(0)