二维PM海洋谱模拟matlab代码

以下是一个简单的二维PM海洋谱模拟的Matlab代码示例：

clc; clear; close all;

Lx = 5000; % 模拟区域x轴长度（m）
Ly = 4000; % 模拟区域y轴长度（m）
Nx = 256; % 离散点个数x轴方向
Ny = 256; % 离散点个数y轴方向
dx = Lx/Nx; % x方向离散间隔（m）
dy = Ly/Ny; % y方向离散间隔（m）

kx = (-Nx/2:Nx/2-1)*2*pi/Lx; % x方向波数
ky = (-Ny/2:Ny/2-1)*2*pi/Ly; % y方向波数
[KX, KY] = meshgrid(kx, ky); % 波数网格

% 二维PM海洋谱
g = @(k) 5*exp(-1.25*(k/0.1).^(-4)); % 能量谱密度函数
S = g(sqrt(KX.^2 + KY.^2)); % 能量谱密度

% 生成随机相位谱
theta = rand(Nx, Ny)*2*pi;

% 计算波高
H = sqrt(2*S.*dx*dy);
H = ifft2(H.*exp(1j*theta));

% 绘制波面图
x = 0:dx:Lx-dx;
y = 0:dy:Ly-dy;
[X, Y] = meshgrid(x, y);
surf(X, Y, real(H));
xlabel('x (m)');
ylabel('y (m)');
zlabel('Height (m)');

这段代码生成了一个256x256的二维PM海洋波高场，并绘制了波面图。你可以根据需要调整模拟区域大小、离散点个数和能量谱密度函数等参数。