fcmyhump.m

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

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


yc = 50;
sigmay = 10


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