matalb编写一个地震波形求出傅里叶幅值谱的程序

以下是一个基于Matlab的地震波形傅里叶幅值谱的程序示例：

% 定义参数 dt = 0.01; % 采样间隔时间 T = 10; % 采样时间 f0 = 5; % 基频

% 生成地震波形 t = 0:dt:T; y = sin(2pif0*t);

% 计算傅里叶幅值谱 N = length(y); f = (0:N-1)/(dt*N); Y = fft(y); amp = abs(Y)/N;

% 绘制波形和幅值谱图 subplot(2,1,1); plot(t,y); xlabel('Time(s)'); ylabel('Amplitude'); title('Earthquake Waveform');

subplot(2,1,2); plot(f,amp); xlim([0 50]); xlabel('Frequency(Hz)'); ylabel('Amplitude'); title('Fourier Amplitude Spectrum');