fftanalyse.m

计算脑电信号非线性参数

%本程序心电信号频谱分析%
%林宛华  中南大学生物医学工程医学图像与信号处理实验室   %
%E-mail: [e-mail]%
%2008年3月10日%
function fftanalyse()
load ECGfs1000;
t=0:0.001:2.499;
subplot(2,1,1),plot(t,ECGfs1000);Ylabel = ('幅值');Xlabel = ('时间');
y = fft(ECGfs1000,2500);
p = y.*conj(y)/2500;  %计算功率谱密度
ff = 1000*(0:1249)/2500;    %计算变换后不同点所对应的频率值
subplot(2,1,2),plot(ff,p(1:1250));