版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

小木虫zb

木虫 (正式写手)

应助: 11 (小学生)
金币: 1560.9
红花: 5
帖子: 345
在线: 142.1小时
虫号: 1227944
注册: 2011-03-10
性别: GG
专业: 分离过程

[求助] 求助matlab拟合方程组

遇到个方程组，不知道怎么用matlab拟合其中的参数，求助大神。k1=0.45967,n1=1/0.33824,k2=0.40308,n2=1/0.38477，拟合求解n11,n12,k12,n22,n21,k21
方程组如下
q1=(k1*c1^(1/n11 +1/n1 ))/(c1^(1/n11 )+k12*c2^(1/n12 ) )
q2=(k2*c2^(1/n22 +1/n2 ))/(c2^(1/n22 )+k21*c1^(1/n21 ) )
源数据如下
c1 q1 c2 q2
0.04556 0.04352 0.099 0.03592
0.10253 0.08542 0.12602 0.08207
0.2056 0.17075 0.27154 0.16137
0.6889 0.30201 0.87702 0.27525
1.51941 0.38387 1.81968 0.34115
2.40038 0.45855 2.72298 0.41266
3.37753 0.51955 3.82469 0.45595
6.32386 0.60045 6.74548 0.54047
9.51436 0.64661 9.89625 0.59228
12.75747 0.68528 13.15723 0.62842
16.07689 0.7131 16.59229 0.63979

回复此楼

» 猜你喜欢

津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名已经有0人回复
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有295人回复
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）已经有6人回复
【修正版】单环用磁约束低速电子实验方案（简化版）已经有0人回复
桂林理工大学物理学专业招收调剂，还有三个名额！！！已经有22人回复
考博自荐已经有4人回复
山东大学第二批博士研究生招生已经有0人回复
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学已经有0人回复
《电磁学》教材推荐已经有1人回复
【急招】合肥工大核聚变材料计算方向2026级工程博士生已经有3人回复

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

matlab数据拟合问题求助已经有7人回复
matlab数据拟合求助！！希望高手帮忙写一下matlab的程序谢谢了！！已经有15人回复
matlab拟合求助已经有11人回复
matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较已经有4人回复
高金求助matlab解微分方程组已经有12人回复
非线性方程组matlab拟合已经有9人回复
急切求助用matlab曲线拟合色散方程已经有10人回复
求助matlab拟合参数的问题已经有2人回复
【求助】matlab拟合二元线性方程已经有3人回复
【求助】使用MATLAB拟合出按照公式的系数已经有19人回复
【求助】求助matlab数值拟合已经有8人回复
【求助】求教多元指数方程如何解有软件否？已经有6人回复
【求助】matlab函数拟合，拟合结果并保存，方便调用已经有4人回复
【求助】matlab曲面拟合表达式已经有10人回复
【求助】如何使得matlab拟合的曲线强制经过一个点已经有8人回复
【求助】指数函数方程求解已经有5人回复
【求助】求MAtlab求解程序！！！已经有20人回复
【求助】求教matlab解非线性方程组已经有9人回复
【求助】使用Matlab拟合反应动力学方程问题已经有7人回复

1楼 2013-01-09 11:07:22

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

回帖支持 ( 显示支持度最高的前 50 名 )

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-01-09 15:58:11
小木虫zb: 金币+5, ★★★很有帮助 2013-01-10 08:53:05

用1stOpt求解试试：

CODE:

Constant  k1=0.45967,n1=1/0.33824,k2=0.40308,n2=1/0.38477;

Variable c1,q1,c2,q2;

SharedModel;

Function q1=(k1*c1^(1/n11 +1/n1 ))/(c1^(1/n11 )+k12*c2^(1/n12 ));

         q2=(k2*c2^(1/n22 +1/n2 ))/(c2^(1/n22 )+k21*c1^(1/n21 ));

Data;

  c1            q1          c2            q2

0.04556        0.04352        0.099        0.03592

0.10253        0.08542        0.12602        0.08207

0.2056        0.17075        0.27154        0.16137

0.6889        0.30201        0.87702        0.27525

1.51941        0.38387        1.81968        0.34115

2.40038        0.45855        2.72298        0.41266

3.37753        0.51955        3.82469        0.45595

6.32386        0.60045        6.74548        0.54047

9.51436        0.64661        9.89625        0.59228

12.75747        0.68528        13.15723        0.62842

16.07689        0.7131        16.59229        0.63979

均方差(RMSE): 0.0335802018362698
残差平方和(SSE): 0.0248078590180215
相关系数(R): 0.997283667098149
相关系数之平方(R^2): 0.994574712660733
决定系数(DC): 0.994566132239987
F统计(F-Statistic): 36.2326353309996

参数最佳估算
-------------------- -------------
n11 0.241571501600743
k12 0.158571836379447
n12 0.218728458244284
n22 -0.168403298883327
k21 0.136718576047343
n21 -0.186531568585624

赞一下(2人)

回复此楼

3楼2013-01-09 13:56:50

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

forcal

禁虫 (初入文坛)

★ ★
csgt0: 金币+2, 应助指数+1, 谢谢 2013-01-14 16:20:27

本帖内容被屏蔽

10楼2013-01-12 16:16:50

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

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
计算强帖: 2
应助: 367 (硕士)
贵宾: 0.457
金币: 7812
散金: 8897
红花: 20
沙发: 7
帖子: 2810
在线: 515.9小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
小木虫zb: 金币+15, ★★★★★最佳答案 2013-01-10 09:52:54

CODE:

function zbzbzbzb

data=[ 0.04556        0.04352        0.099        0.03592   %  c1,q1,c2,q2

 0.10253        0.08542        0.12602        0.08207

 0.2056        0.17075        0.27154        0.16137

 0.6889        0.30201        0.87702        0.27525

 1.51941        0.38387        1.81968        0.34115

 2.40038        0.45855        2.72298        0.41266

 3.37753        0.51955        3.82469        0.45595

 6.32386        0.60045        6.74548        0.54047

 9.51436        0.64661        9.89625        0.59228

 12.75747        0.68528        13.15723        0.62842

 16.07689        0.7131        16.59229        0.63979];

xdata=[data(:,1),data(:,3)];

ydata=[data(:,2),data(:,4)];

x0=[0.2,0.2,0.2,9314,106,0.6]; %n11,n12,k12,n22,n21,k21

[x,resnorm] = lsqcurvefit(@myfun,x0,xdata,ydata)

y=myfun(x,xdata)

end

function y=myfun(x,xdata)

c1=xdata(:,1);

c2=xdata(:,2);

n11=x(1);

n12=x(2);

k12=x(3);

n22=x(4);

n21=x(5);

k21=x(6);

k1=0.45967;

n1=1/0.33824;

k2=0.40308;

n2=1/0.38477;

q1=(k1*c1.^(1/n11 +1/n1 ))./(c1.^(1/n11 )+k12*c2.^(1/n12 ) );

q2=(k2*c2.^(1/n22 +1/n2 ))./(c2.^(1/n22 )+k21*c1.^(1/n21 ) );

y=[q1,q2];

end

赞一下(1人)

回复此楼

showmethemoney

2楼2013-01-09 11:45:59

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

普通回帖

小木虫zb

木虫 (正式写手)

应助: 11 (小学生)
金币: 1560.9
红花: 5
帖子: 345
在线: 142.1小时
虫号: 1227944
注册: 2011-03-10
性别: GG
专业: 分离过程

引用回帖:

2楼: Originally posted by csgt0 at 2013-01-09 11:45:59
function zbzbzbzb
data=;
xdata=;
ydata=;
x0=; %n11,n12,k12,n22,n21,k21
= lsqcurvefit(@myfun,x0,xdata,ydata)
y=myfun(x,xdata)
end

function y=myfun(x,xdata)
c1=xdata(:,1);
c2=xdata(:,2); ...

运行不出来，出现错误

360截图20130110092051698.jpg

赞一下

回复此楼

4楼2013-01-10 09:21:45

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