matlab频谱分析功率代码

在MATLAB中，进行功率谱分析通常涉及到快速傅立叶变换（FFT）函数，如`fft()`或`pwelch()`。下面是一个基本的示例代码，用于计算信号的功率谱密度：

% 假设有一个时间序列数据 'signal'
signal = randn(1, 1000); % 生成随机模拟信号

% 对信号进行离散傅立叶变换（假设采样频率fs=1）
fs = 1;            % 采样频率
time_vec = (0:length(signal)-1)'/fs;   % 时间向量
fft_result = fft(signal);        % 进行FFT

% 计算幅度谱，并取其平方得到功率谱
power_spectrum = abs(fft_result).^2;  
power_density = power_spectrum ./ length(signal);     % 单位：瓦特/赫兹 (W/Hz)

% 获取中心频率对应的谱值，通常只关心正半轴
frequencies = (0:(length(power_density)-1)) * fs / length(signal);

% 可能需要对结果进行归一化、滤波等进一步处理
% plot(frequencies, power_density, 'b') % 绘制功率谱图