matlab写一段程序，模拟水陆检鬼波压制

在MATLAB中，我们可以编写一个简单的程序来模拟水陆交替检测（也称为鬼影波抑制），通常用于信号处理领域，尤其是在地震数据分析中。这个过程通常涉及到滤波、去噪以及信号重构等步骤。以下是一个简化的示例，我们将创建一个随机信号并应用一种常见的鬼影波抑制技术——自适应Wiener滤波。

% 导入所需库
import numpy as np
from scipy.signal import wiener

% 创建模拟信号（假设是地震数据）
length_data = 1000; % 数据长度
freq_water = 50; % 水体频率
freq_terrain = 100; % 土壤频率
water_amplitude = 0.5; % 水体幅度
terrain_amplitude = 1; % 土壤幅度

% 生成水体和土壤成分的正弦信号
signal_water = water_amplitude * sin(2 * pi * freq_water * (0:length_data) / length_data);
signal_terrain = terrain_amplitude * sin(2 * pi * freq_terrain * (0:length_data) / length_data);

% 合成原始信号
original_signal = signal_water + signal_terrain;

% 应用自适应Wiener滤波
% 这里我们假设我们知道水体成分的信息，可以用其作为噪声模型
noise_floor = mean(signal_water); % 假设噪声等于水体成分的平均值
wiener_filter = wiener(original_signal, [0 0], noise_floor); % 第二个参数设置滤波窗口

% 显示原始信号和经过滤波后的信号
figure;
subplot(2,1,1);
plot(original_signal);
title('Original Signal');
xlabel('Sample Index');
ylabel('Amplitude');

subplot(2,1,2);
plot(wiener_filter);
title('Filtered Signal with Adaptive Wiener Filter');
xlabel('Sample Index');
ylabel('Amplitude');