matlab画幅值谱 功率谱

在Matlab中，我们可以使用fft函数来计算信号的幅值谱和功率谱。幅值谱表示信号在不同频率下的幅值大小，而功率谱表示信号在不同频率下的功率大小。

要计算幅值谱，首先需要输入信号。可以通过以下代码创建一个包含信号的向量：

t = 0:0.1:10; % 时间向量
x = sin(2*pi*2*t) + 0.5*sin(2*pi*5*t); % 信号向量

接下来，我们可以使用fft函数对信号进行傅里叶变换来计算幅值谱：

N = length(x); % 信号长度
X = fft(x); % 傅里叶变换
amplitude_spectrum = abs(X(1:N/2+1)); % 幅值谱

在上面的代码中，我们使用abs函数来计算傅里叶变换的幅值，并且只保留了一半的结果，因为幅值谱是对称的。

要计算功率谱，可以使用下面的代码：

power_spectrum = (1/(N^2)) * abs(X(1:N/2+1)).^2; % 功率谱

在上面的代码中，我们使用了信号的幅值谱，并且对每个频率的幅值平方进行了计算。

最后，我们可以使用plot函数将幅值谱和功率谱绘制出来：

fs = 1/(t(2)-t(1)); % 采样率
frequencies = (0:(N/2))*fs/N; % 频率向量

figure;
subplot(2,1,1);
plot(frequencies, amplitude_spectrum);
xlabel('频率 (Hz)');
ylabel('幅值');
title('幅值谱');

subplot(2,1,2);
plot(frequencies, power_spectrum);
xlabel('频率 (Hz)');
ylabel('功率');
title('功率谱');

上面的代码中，我们使用subplot函数将幅值谱和功率谱绘制在同一个图像中的不同子图中，并通过xlabel、ylabel和title函数来添加轴标签和标题。

运行以上代码，将会生成一个包含信号的幅值谱和功率谱的图像。