帮我写一个完整的matlab代码以完成二维PM海浪谱产生的海浪

以下是一个简单的MATLAB代码，用于生成二维PM海浪谱：

%设置参数
Lx = 5000; % x方向长度
Ly = 5000; % y方向长度
Nx = 256; % x方向离散点数
Ny = 256; % y方向离散点数
dx = Lx/Nx; % x方向离散间距
dy = Ly/Ny; % y方向离散间距
kx = (2*pi/Lx)*[0:(Nx/2-1) (-Nx/2):-1]; % x方向波数
ky = (2*pi/Ly)*[0:(Ny/2-1) (-Ny/2):-1]; % y方向波数
[KX,KY] = meshgrid(kx,ky); % 生成波数矩阵
g = 9.81; % 重力加速度
h = 100; % 水深
wp = sqrt(g*KX.*tanh(KX*h)); % 角频率
A = randn(Nx,Ny); % 产生随机数矩阵
B = randn(Nx,Ny); % 产生随机数矩阵

% 计算海面高度
eta = zeros(Nx,Ny);
for i=1:Nx
    for j=1:Ny
        eta(i,j) = sqrt((A(i,j)^2 + B(i,j)^2)/2)*exp(1i*2*pi*rand); % 计算海面高度
    end
end

% Fourier反变换
eta = real(ifft2(eta));

% 绘制海浪图像
x = 0:dx:Lx-dx;
y = 0:dy:Ly-dy;
[X,Y] = meshgrid(x,y);
surf(X,Y,eta);

这段代码首先设置了一些参数，比如海域的长度和宽度、离散点数、离散间距等等。然后，通过随机数生成器生成了两个随机数矩阵，用于计算海面高度。接着，计算了海面高度，并进行Fourier反变换。最后，使用surf函数绘制了海浪图像。

请注意，这只是一个简单的实现，可能需要根据你的具体需求进行调整。