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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (4277)
>文献求助 (466)
>虫友互识 (375)
>导师招生 (365)
>考博 (203)
>博后之家 (182)
>硕博家园 (162)
>休闲灌水 (154)
>论文投稿 (131)
>基金申请 (94)
>考研 (72)
>招聘信息布告栏 (66)
>教师之家 (60)
>绿色求助(高悬赏) (47)
>公派出国 (46)
>找工作 (38)

返回列表

小木虫zb

木虫 (正式写手)

应助: 11 (小学生)
金币: 1560.1
红花: 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

回复此楼

» 猜你喜欢

求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有105人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

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

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.5小时
虫号: 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
金币: 7807
散金: 8897
红花: 20
沙发: 7
帖子: 2809
在线: 515.8小时
虫号: 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.1
红花: 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的回帖