代码如下:
x=1:100;
y=1:100;
Q1=0.14*1000*1000*1000/24/60/60;
u1=2.637037037;
H=80+2.4*1.5*1.8/u1;
[X,Y]=meshgrid(x,y);
sy=0.16.*X.*(1.+0.0001.*X).^0.5;
sz=0.12*X;
Sc=(Q1/(2*pi*u1).*sy.*sz).*exp(-(Y.^2/(sy.^2)+2*H^2./(sz.^2)));
surf(X,Y,Sc);
shading flat
现在可以出结果 也不是NAN了,但是全是0 T-T
追答我这不全是零啊。