% 以下是基于Matlab实现的同震反演模型
% 代码实现了使用psokinv进行同震反演
% 希望对同震研究领域的研究者有所帮助

% 加载数据
load data.mat

% 初始化参数
param.dt = 0.01; % 采样间隔
param.f0 = 10; % 激发频率
param.t0 = 5; % 激发时间
param.vel = 1500; % 介质速度
param.num_iters = 100; % 迭代次数

% 调用psokinv进行同震反演
result = psokinv(data, param);

% 输出结果
disp(result)