fcmhump.m

来自「波浪数值模拟」· M 代码 · 共 19 行

xeta = [0:399]';
y = [0:99];


yc = 50;
sigmay = 20


hump = 0.1*exp( -0.5*  ones(40,1)*(y-yc).^2/sigmay^2 ) );

sumhump = sum(sum(hump))/(40*100);

hump = hump - sumhump;

contour(xeta,y,hump')

save -ascii fCyhump.eta   hump