MATLAB计算地震两点相干函数代码

以下是一个简单的 MATLAB 代码，用于计算地震两点之间的相干函数：

% 设置两个地震记录文件名
file1 = 'record1.txt';
file2 = 'record2.txt';

% 读取文件中的数据
data1 = load(file1);
data2 = load(file2);

% 获取数据的长度和采样频率
n = length(data1);
fs = 100; % 假设采样频率为100Hz

% 计算两个信号的功率谱密度
Pxx = abs(fft(data1)).^2 ./ n;
Pyy = abs(fft(data2)).^2 ./ n;

% 计算两个信号的交叉功率谱密度
Pxy = fft(data1) .* conj(fft(data2)) ./ n;

% 计算相干函数
C = abs(Pxy).^2 ./ (Pxx .* Pyy);

% 绘制相干函数曲线
f = (0:n-1) * fs / n;
plot(f, C);
xlabel('频率 (Hz)');
ylabel('相干函数');

请注意，这只是一个简单的示例代码，实际情况可能会更加复杂。例如，您可能需要对数据进行滤波、去噪等预处理操作，以获得更准确的结果。同时，您也需要确保您的数据格式正确，并且您已经正确设置了采样频率等参数。