cal_psd.m

用 CC法计算延迟时间和潜入位数的matlab程序

% MATLAB 7 辅助信号处理技术与应用
% 经典功率谱估计 293
% 6.2.1  相关函数估计 293
% 6.2.2  基于直接法的功率谱估计 299
% 6.2.3  基于间接法的功率谱估计 302
% 6.2.4  基于改进直接法的功率谱估计 303
% 6.2.5  基于多窗口法的功率谱估计 309
% 6.3  基于参数建模的功率谱估计 311
% 6.3.1  基于AR模型的功率谱估计 312
% 6.3.2  基于MA模型的功率谱估计 329
% 6.3.3  基于ARMA模型的功率谱估计 331
% 6.3.4  基于最小方差的功率谱估计 334
% 6.4  基于非参数建模的功率谱估计 337
% 6.4.1  相关矩阵的特征分解 337
% 6.4.2  基于MUSIC算法的功率谱估计 339
% 6.4.3  基于特征向量的功率谱估计 342
% 6.4.4  信号与噪声子空间维数估计 344
Fs = 1000;
t = 0:1/Fs:1;
nfft = 1024;
index=1:round(nfft/2)+1;
hz=length(t)/(nfft*1);
frequency=index*hz;
xn = sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(size(t));
Pxx = abs(fft(xn,nfft).^2/1001);
plot(frequency, Pxx);