地震数据处理时，功率谱和功率谱密度分别的Matlab代码

以下是地震数据处理时计算功率谱和功率谱密度的Matlab代码示例：

假设有一个地震数据的时间序列为x，采样频率为Fs。

计算功率谱：

% 计算功率谱
N = length(x); % 数据长度
X = fft(x); % 对数据进行快速傅里叶变换
Pxx = (1/(Fs*N)) * abs(X).^2; % 计算功率谱
f = linspace(0,Fs/2,N/2+1); % 频率向量
Pxx = Pxx(1:N/2+1); % 取一半的功率谱，因为对称性
% 画功率谱图
plot(f,Pxx);
xlabel('频率（Hz）');
ylabel('功率谱密度（m^2/Hz）');

计算功率谱密度：

% 计算功率谱密度
N = length(x); % 数据长度
window = hamming(N); % 窗口函数
noverlap = N/2; % 重叠长度
nfft = 2^nextpow2(N); % 傅里叶变换点数
[S,f] = pwelch(x,window,noverlap,nfft,Fs); % 计算功率谱密度
% 画功率谱密度图
plot(f,S);
xlabel('频率（Hz）');
ylabel('功率谱密度（m^2/Hz）');