matlab互功率谱

matlab互功率谱
互功率谱（Cross-Power Spectrum）是信号处理中的一种方法，用于分析两个信号之间的相互关系。

在MATLAB中，可以使用xcorr 函数来计算两个信号的互功率谱。

以下是一个简单的示例：
1. 首先，生成两个正弦波信号x和y：
matlab
Fs = 1000; % 采样频率
T = 1/Fs; % 采样周期
L = 1000; % 信号长度
t = (0:L-1)*T; % 时间向量
x = sin(2*pi*50*t) + sin(2*pi*120*t); % 信号x
y = sin(2*pi*60*t) + sin(2*pi*150*t); % 信号y
2. 使用xcorr函数计算互功率谱：
matlab
Rxy = xcorr(x, y); % 计算互相关函数
Pxy = abs(fft(Rxy)).^2; % 计算互功率谱
f = Fs*(0:(L/2))/L; % 计算频率向量
figure;
plot(f, Pxy); % 绘制互功率谱图
xlabel('Frequency (Hz)');
ylabel('Power Spectrum');
title('Cross-Power Spectrum of x and y');
在这个示例中，我们首先生成了两个正弦波信号x和y，然后使用xcorr函数计算它们的互相关函数Rxy。

接下来，我们对Rxy进行傅里叶变换，得到频域表示，然后计算其平方模得到互功率谱Pxy。

最后，我们绘制了互功率谱图。