水波仿真源码程序,程序员在编程的过程中可以参考学习使用,希望对IT程序员有用,此源码程序简单易懂、方便阅读,有很好的学习价值!
% Examples:
clc;
clear;
close all;
% Examples:
nx = 101; xmin = 0; xmax = 100; x = linspace(xmin,xmax,nx);
ny = 51; ymin = 0; ymax = 50; y = linspace(ymin,ymax,ny);
wind_data.U = 10; wind_data.thetaU = 0; wind_data.X = 1e6;
[s,Tp,fm,B,Sk,kx,ky] = sea_surface(x,y,wind_data,'PM','none');
figure(1)
subplot(211),mesh(x,y,s),ylabel('y(m)'),xlabel('x(m)')
subplot(212),mesh(kx,ky,Sk),ylabel('ky (1/m)'),xlabel('kx (1/m)')
[s,Tp,fm,B,Sk,kx,ky] = sea_surface(x,y,wind_data,'PM','cos2');
figure(2)
subplot(211),mesh(x,y,s),ylabel('y(m)'),xlabel('x(m)')
subplot(212),mesh(kx,ky,Sk),ylabel('ky (1/m)'),xlabel('kx (1/m)')
[s,Tp,fm,B,Sk,kx,ky] = sea_surface(x,y,wind_data,'PM','mits');
figure(3)
subplot(211),mesh(x,y,s),ylabel('y(m)'),xlabel('x(m)')
subplot(212),mesh(kx,ky,Sk),ylabel('ky (1/m)'),xlabel('kx (1/m)')
[s,Tp,fm,B,Sk,kx,ky] = sea_surface(x,y,wind_data,'PM','hass');
figure(4)
subplot(211),mesh(x,y,s),ylabel('y(m)'),xlabel('x(m)')
subplot(212),mesh(kx,ky,Sk),ylabel('ky (1/m)'),xlabel('kx (1/m)')