24小时热门版块排行榜    

查看: 2242  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

qs秦秦

新虫 (小有名气)

[求助] MATLAB拟合土壤水分特征曲线 已有1人参与

想实现这篇论文的程序,但是老是出现错误,刚接触这个东西,不知道怎么去改,求会matlab的虫友们帮帮忙,谢谢了!

MATLAB拟合土壤水分特征曲线
AFEHPDQU}Q[FF~5F7$TC@{Q.jpg


MATLAB拟合土壤水分特征曲线-1
A9X4]%_4K~6{NCQF)4U5UGJ.jpg
回复此楼

» 收录本帖的淘帖专辑推荐

程序

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qs秦秦

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 月只蓝 at 2014-07-28 20:03:15
需要指出的是,截图中的计算步骤,在输入y1=... 时,”1./1.5003“这一项的点除是多余的。

多谢啦
4楼2014-08-09 15:36:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
qs秦秦: 金币+5, ★★★★★最佳答案 2014-08-09 15:37:14
一下所有程序复制进一个m文件,点击运行即可。
CODE:
function feixianxingnihe_3
clear all;clc
format long


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];
xspan=xdata;    %x的数据,在此输入
Texp=ydata;    %T的数据,在此输入

x0=[0.1 0.1 0.1 1];
k0=x0;
lb=-[1 1 1 1]*1e9;
ub=[1 1 1 1]*1e9;


%-------------------------------------------------------------------------

% 使用函数lsqnonlin()进行参数估计

OPTIONS=optimset('MaxFunEvals',1000);
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
    lsqnonlin(@ObjFunc,k0,lb,ub,OPTIONS,xspan,Texp);

ci = nlparci(k,residual,jacobian);
%residual;
fprintf('\n\n拟合结果:\n')
fprintf('\n\t参数 c1 = %.16f',k(1))
fprintf('\n\t参数 c2 = %.16f',k(2))
fprintf('\n\t参数 c3 = %.16f',k(3))
fprintf('\n\t参数 c4 = %.16f',k(4))
fprintf('\n\t残差平方和 = %.6e',resnorm)
y=KineticsEqs(xspan,k);
R2=1-sum((Texp-y).^2)./sum((Texp-mean(y)).^2);
%fprintf('\n\t相关系数之平方R^2 = %.6e',R2);

mm=max(xspan)-min(xspan);
xspan1=min(xspan):0.01:max(xspan)+0.05*mm;
figure(3)
plot(xspan1,KineticsEqs(xspan1,k),'b',xspan,Texp,'or'),legend('拟合曲线','实验数据','Location','Best'),...
    axis([-10 16000 0.1 0.6])


%-------------------------------------------------------------------------

function f = ObjFunc(k,xspan,Texp)
f=KineticsEqs(xspan,k)-Texp;

%------------------------------------------------------------------------
function xt = KineticsEqs(x,n)
k1=n(1);
k2=n(2);
k3=n(3);
k4=n(4);
xt=k1+(k2-k1)./(1+(k3*x).^k4).^(1-1/k4);

计算结果:
拟合结果:

        参数 c1 = 0.1674744224584986
        参数 c2 = 0.5653038321277274
        参数 c3 = 0.0475760372012390
        参数 c4 = 1.5279673014244484
        残差平方和 = 4.767805e-004>>
MATLAB拟合土壤水分特征曲线-2
附图1.png

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2014-07-28 20:00:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

需要指出的是,截图中的计算步骤,在输入y1=... 时,”1./1.5003“这一项的点除是多余的。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
3楼2014-07-28 20:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuanbeichuan

新虫 (初入文坛)

大神你好。
请问 请问怎样用excel 拟合。
用matlab水平不高
5楼2016-09-07 11:39:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 一志愿华中农业071010,总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 8/400 2026-03-20 15:58 by babero
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
信息提示
请填处理意见