土壤水的基质势或土壤水吸力是随土壤含水率而变化的，其关系曲线称为土壤水分特征曲线或土壤持水曲线，土壤水分特征曲线表示土壤水的能量和数量之间的关系，是研究土壤水分的保持和运动所用到的反映土壤水分基本特性的曲线。
matlab教程
比如我们用van Gnuchuten方程，打开matlab软件首先建立函数，file--new-script建立新的函数文件，
function y=ojefun(x,xdata)
y=x(1)+(x(2)-x(1))./(1+(x(3)*xdata).^x(4)).^(1-1./x(4));
函数为ojefun，那么这个函数文件名必须为ojefun.m，这个很重要。
然后在command window里输入以下命令
xdata=[0 50.65 293.77 790.14 992.74 5065 10130 15195];
ydata=[0.5650 0.4013 0.2502 0.2324 0.2307 0.1926 0.1812 0.1730];
x0=[0.1,0.1,0.01,1];
[x,resnorm]=lsqcurvefit(@ojefun,x0,xdata,ydata)
然后回车便得到
x =
0.1676 0.5653 0.0475 1.5291
resnorm =
4.7675e-004
对应x1,x2,x3,x4和残差平方和。（xdata为吸力值，ydata为对应含水量，x0分别x1,x2,x3,x4的初始值）
|
发表于 2012-1-19 10:14:14
楼主
发表于 2012-3-12 17:19:26
