基于FFT的频谱信号分析Matlab程序

基于FFT的频谱信号分析Matlab示例程序

clc
clear
close all
data=load（‘test1.txt‘）;
figure（1）
plot（data（:1）data（:2）‘r‘data（:1）data（:3）‘b‘）;
legend（‘信号1‘‘信号2‘）;
grid

%求付利叶变换
frequency=1000*（0:512）/1024;  % 实际频率点
fftdata=fft（data）;   %付利叶变换
sinfft=fftdata（1:5132）;
cosfft=fftdata（1:5133）;

%求能量最大点
[fz1f1]=max（abs（sinfft））;
[fz2f2]=max（abs（cosfft））;

% 求两信号能量最大点的频率
sin_frequency = frequency（f1）;   %得到实际频率
cos_frequency = frequency（f2）;

% 求两信号能量最大点的相位
phase1=phase（sinfft（f1））*57.3;
phase2=phase（cosfft（f1））*57.3;

fprintf（‘信号1频率=%fHz\n信号2频率=%fHz\n‘sin_frequencycos_frequency）;
fprintf（‘信号2比信号1幅值衰减=%f db\n‘20*log10（fz2/fz1））;
fprintf（‘信号2比信号1相位差=%f度\n‘phase1-phase2）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        760  2008-05-15 12:07  signal_analyszer.m

     文件      34817  2008-05-14 13:22  test1.txt

     文件        310  2009-08-30 21:03  说明.txt

----------- ---------  ---------- -----  ----

                35887                    3