给出一组数据_如何同matlab中提供的工具箱拟合一条曲线,[MATLAB工具箱] 怎样用一组参数同时拟合两个曲线...

icon1.gif 怎样用一组参数同时拟合两个曲线

如题,小弟想拟合曲线r1和r2,这两条曲线都用的同一组参数A1 A2 A3 A4 A5 A6 A7 Aw E1 E2 E3 E4 E5 E6 E7 Ew进行非线性最小二乘拟合,请问怎样用matlab实现一组参数同时拟合两个方程?(如果分开拟合的话每个参数必然有两个值)。下面是相关的数据及方程,其中y1和y2分别是r1和r2的实验值,r1和r2的自变量值相同,均为x1.请各位大侠能鼎力相助,最好能给出代码或者计算方法,谢谢

x1=[0.51333,0.36301,0.34194,0.28502,0.26485,0.17771];

y1=[0.50536,0.61684,0.63566,0.68357,0.70947,0.72182];

y2=[2.5363,1.16376,0.75531,0.4164,0.27812,0.1259];

syms A1 A2 A3 A4 A5 A6 A7 Aw E1 E2 E3 E4 E5 E6 E7 Ew x;

k1=A1*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E1)/4098.802);

k2=A2*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E2)/4098.802);

k3=A3*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E3)/4098.802);

k4=A4*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E4)/4098.802);

k5=A5*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E5)/4098.802);

k6=A6*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E6)/4098.802);

k7=A7*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E7)/4098.802);

r1=-(x*k1*k2*k3*(-k5-k7))/(-(x*k1*k2-k3*(-x*k1-k2-x*k4-k6))*(-k5-k7)-x*k3*k4*(-x*k1+k7));

r2=-(x^2)*k1*k3*k4*k5/((x*k1*k2-k3*(-x*k1-k2-x*k4-k6))*(-k5-k7)-x*k3*k4*(-x*k1+k7));

100 <= A1 <= 10000;

10^11 <= A2 <= 10^14;

10^12 <= A3 <= 10^14;

100 <= A4 <= 10000;

10^12 <= A5 <= 10^14;

10^12 <= A6 <= 10^14;

10^12 <= A7 <= 10^14;

100 <= Aw <= 10000;

0 <= E1 <= 100;

50 <= E2 <= 140;

50 <= E3 <= 140;

0 <= E4 <= 100;

50 <= E5 <= 150;

100 <= E6 <= 150;

0 <= E7 <= 100;

0 <= Ew <= 100;

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值