符号说明 运行环境：MATLAB2016a 附录 1 程序 1：数据拟合 x1=[]; x2=[]; c=[]; %plot3(x1,x2,c) x0=[ 1 1 1 1]; d=length(x0); f=@(a,x)0.08*(ones(1,30)-(x(1,:)./7).*(ones(1,30)-a(1)*(ones(1,30)-((x(1,:)./7).^a(2))).*exp(a(3)* sqrt((x(2,:))./2982))))+ones(1,30)*a(4); a=lsqcurvefit(f,x0,[x1;x2],c) %%作图 B=0.08*(ones(1,30)-(x1./7).*(ones(1,30)-a(1)*(ones(1,30)-((x1./7).^a(2))).*exp(a(3)*sqrt((x2)./29 82))))+a(4)*ones(1,30); plot3(x1,x2,B)