有他们相关的Matlab代码嘛

当然，以下是它们相关的Matlab代码：

1. 时频谱

% 定义信号参数
fs = 1000;      % 采样率
t = 0:1/fs:1;   % 时间向量
x = sin(2*pi*50*t) + sin(2*pi*150*t);  % 信号

% 计算短时傅立叶变换
N = length(x);
w = hann(N);    % 汉宁窗
noverlap = N/2; % 重叠长度
nfft = 1024;    % FFT长度
[S,F,T] = spectrogram(x,w,noverlap,nfft,fs,'yaxis');

% 绘制时频谱图
surf(T,F,10*log10(abs(S)),'edgecolor','none');
axis tight;
view(0,90);
xlabel('Time (s)');
ylabel('Frequency (Hz)');
title('Spectrogram');

2. 功率谱

% 定义信号参数
fs = 1000;      % 采样率
t = 0:1/fs:1;   % 时间向量
x = sin(2*pi*50*t) + sin(2*pi*150*t);  % 信号

% 计算功率谱
N = length(x);
X = fft(x,N);
P = abs(X).^2/N/fs;

% 绘制功率谱图
f = (0:N-1)*fs/N;
plot(f,P);
xlabel('Frequency (Hz)');
ylabel('Power');
title('Power Spectrum');

3. 振幅谱

% 定义信号参数
fs = 1000;      % 采样率
t = 0:1/fs:1;   % 时间向量
x = sin(2*pi*50*t) + sin(2*pi*150*t);  % 信号

% 计算振幅谱
N = length(x);
X = fft(x,N);
A = abs(X)/N/fs*2;

% 绘制振幅谱图
f = (0:N-1)*fs/N;
plot(f,A);
xlabel('Frequency (Hz)');
ylabel('Amplitude');
title('Amplitude Spectrum');

这些Matlab代码可以帮助你更好地理解时频谱、功率谱和振幅谱的计算和绘制方法，也可以作为你处理地震数据的参考。