请选择 进入手机版 | 继续访问电脑版

土壤之家-土壤学论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

土壤考研资料sigmaplot12.0破解版下载sigmaplot教程土壤污染
土壤修复案例土壤招聘信息土壤修复技术有哪些 
查看: 8760|回复: 4

[Matlab] 土壤水分特征曲线拟合(matlab教程)

[复制链接]
发表于 2012-1-19 10:14:14 | 显示全部楼层 |阅读模式
土壤水的基质势或土壤水吸力是随土壤含水率而变化的,其关系曲线称为土壤水分特征曲线或土壤持水曲线,土壤水分特征曲线表示土壤水的能量和数量之间的关系,是研究土壤水分的保持和运动所用到的反映土壤水分基本特性的曲线。

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-2-8 22:55:03 | 显示全部楼层
欢迎大家来交流!
回复 支持 反对

使用道具 举报

发表于 2012-3-12 17:19:26 | 显示全部楼层
提示: 该帖被管理员或版主屏蔽
回复 支持 反对

使用道具 举报

发表于 2015-7-17 11:04:41 | 显示全部楼层
太好了,感谢分享
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|小黑屋|手机版|Archiver|网站地图|土壤之家

GMT+8, 2024-4-18 11:45 , Processed in 0.214919 second(s), 18 queries .

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表