请帮我用matlab设计一个频散曲线提取的代码，算法为f-k方法

好的，下面是用MATLAB实现f-k算法提取频散曲线的代码：

% 读入数据
data=load('data.txt');
[nz,nx]=size(data);

% 参数设置
dz=1; % 深度采样间隔
dx=1; % 横向采样间隔
dt=0.001; % 采样时间间隔
fnyquist=1/(2*dt); % 奈奎斯特频率
nf=512; % 频率采样点数
df=fnyquist/nf; % 频率采样间隔
fmax=fnyquist; % 最大频率
f=[0:df:fmax-df]'; % 频率轴

% FK变换
p=2*pi*fft2(data);

% 滤波
w=2*pi*f;
kx=[0:1/(nx*dx):(nx-1)/(2*nx*dx),-nx/(2*dx):1/(nx*dx):-1/(dx)];
kz=[0:1/(nz*dz):(nz-1)/(2*nz*dz),-nz/(2*dz):1/(nz*dz):-1/(dz)];
[KX,F]=meshgrid(kx,w);
[KZ,F]=meshgrid(kz,w);
K2=KX.^2+KZ.^2;
k2min=min(K2(:));
k2max=max(K2(:));
D=exp(-1i*K2/(4*k2min));
p_fk=p.*D;

% 反FK变换
data_fk=ifft2(p_fk);

% 提取频散曲线
fdom=zeros(nz,1);
for iz=1:nz
    [~,imax]=max(abs(data_fk(iz,:)));
    fdom(iz)=f(imax);
end

% 绘制频散曲线
figure;
plot(fdom);
xlabel('Depth');
ylabel('Frequency (Hz)');
title('Dispersion Curve');

其中，data.txt是输入数据文件，每行代表一条地震记录。代码中，我们首先读入数据，然后设置参数，包括采样间隔、频率采样点数、最大频率等。接着，我们利用FFT2函数进行FK变换，并进行滤波操作。然后，我们将反FK变换后的数据用于提取频散曲线。最后，我们绘制频散曲线，即波速与频率的关系曲线。