用matlab拟合多元函,matlab 多元函数拟合 - 计算模拟 - 小木虫 - 学术 科研 互动社区...

于 2021-03-16 12:06:25 发布
阅读量1.7k 收藏 3
点赞数
文章标签： 用matlab拟合多元函

cfdc6361174d9cfc8169daed13483a0c.gif

月只蓝

一、MATLAB代码和结果如下，图形结果见附图1。CODE:

function feixianxingnihe_3

clear all;clc

format long

data=[298,1.51941000000000;473,1.47000000000000;673,1.38000000000000;873,1.32000000000000;1073,1.28000000000000];

xspan=data(:,1);    %x的数据，在此输入

Texp=data(:,2);    %T的数据，在此输入

k0=[10 10];

lb=-[1 1]*1e9;

ub=[1 1]*1e9;

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

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

OPTIONS=optimset('MaxFunEvals',1000);

[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...

lsqnonlin(@ObjFunc,k0,lb,ub,OPTIONS,xspan,Texp);

fprintf('\n\n拟合结果:\n')

fprintf('\n\t参数 A = %.16f',k(1))

fprintf('\n\t参数 a = %.16f',k(2))

y=KineticsEqs(xspan,k);

R2=1-sum((Texp-y).^2)./sum((Texp-mean(y)).^2);

fprintf('\n\t相关系数之平方R^2 = %.16f',R2);

figure(1)

mm=max(xspan)-min(xspan);

xspan1=min(xspan)-0.05*mm:0.01:max(xspan)+0.05*mm;

plot(xspan1,KineticsEqs(xspan1,k),'b',xspan,Texp,'or'),legend('计算值','实验值','Location','Best')

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

function f = ObjFunc(k,xspan,Texp)

f=KineticsEqs(xspan,k)-Texp;

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

function xt = KineticsEqs(x,k)

global b

A=k(1);a=k(2);

xt=A/(3*137.8)+A./x.^0.5.*2/a^0.5.*atan((-0.5+0.5*(1+a*137.8./x.^2).^0.5).^0.5)-2*2^0.5/3*A*137.8^0.5./x.^1.5./(1+(1+a*137.8./x.^2).^0.5).^1.5;

CODE:

拟合结果:

参数 A = 502.1696785220468100

参数 a = 23729.5421378699030000

相关系数之平方R^2 = 0.9274532902868037>>

二、软件1stOpt计算代码和结果CODE:

均方差(RMSE): 0.0241425766168092

残差平方和(SSE): 0.0029143200284925

相关系数(R): 0.963175718167507

相关系数之平方(R^2): 0.927707464067493

决定系数(DC): 0.927453222207639

卡方系数(Chi-Square): 0.00102747119728203

F统计(F-Statistic): 38.4980600874317

参数        最佳估算

----------        -------------

a                  502.170022602279

a1                 23730.0107275906

====== 结果输出 =====

No        实测值y        计算值y

1        1.51941        1.5410049

2        1.47        1.4301342

3        1.38        1.3659321

4        1.32        1.3285206

5        1.28        1.3042527

w134h1122189_1405857421_285.png

附图1.png

0014f08102b5cebd97bfef69c59ca4cf.png

qj2006

引用回帖:

月只蓝 at 2014-07-20 19:57:04

一、MATLAB代码和结果如下，图形结果见附图1。

function feixianxingnihe_3

clear all;clc

format long

data=;

xspan=data(:,1);    %x的数据，在此输入

Texp=data(:,2);    %T的数据，在此输入

k0= ...

谢谢你，我这还有一组数据，麻烦你再给我拟合一次。

采用关系式：

y=A/(3*137.8)+A/x^0.5*2/a^0.5*atan((-0.5+0.5*(1+a*137.8/x^2)^0.5)^0.5)-2*2^0.5/3*A*137.8^0.5/x^1.5/(1+(1+a*137.8/x^2)^0.5)^1.5

x=[298,473,673,873,1073];

y=[1.74322,1.6,1.47,1.4,1.37];

cfdc6361174d9cfc8169daed13483a0c.gif

月只蓝

引用回帖:

qj2006 at 2014-07-20 20:04:45

谢谢你，我这还有一组数据，麻烦你再给我拟合一次。

采用关系式：

y=A/(3*137.8)+A/x^0.5*2/a^0.5*atan((-0.5+0.5*(1+a*137.8/x^2)^0.5)^0.5)-2*2^0.5/3*A*137.8^0.5/x^1.5/(1+(1+a*137.8/x^2)^0.5)^1.5

x=;

y ...

均方差(RMSE): 0.0166330356441843

残差平方和(SSE): 0.00138328937370353

相关系数(R): 0.992803768743302

相关系数之平方(R^2): 0.985659323230904

决定系数(DC): 0.985526033247728

卡方系数(Chi-Square): 0.000442346340380584

F统计(F-Statistic): 206.195148060588

参数        最佳估算

----------        -------------

A                  525.805717177879

a1                 7213.56939658239

====== 结果输出 =====

No        实测值y        计算值y

1        1.74322        1.7557642

2        1.6        1.5687255

3        1.47        1.4677259

4        1.4        1.4127136

5        1.37        1.3790017

0014f08102b5cebd97bfef69c59ca4cf.png

qj2006

引用回帖:

月只蓝 at 2014-07-20 20:07:47

均方差(RMSE): 0.0166330356441843

残差平方和(SSE): 0.00138328937370353

相关系数(R): 0.992803768743302

相关系数之平方(R^2): 0.985659323230904

决定系数(DC): 0.985526033247728

卡方系数(Chi-Square): ...

非常感谢

0014f08102b5cebd97bfef69c59ca4cf.png

qj2006

引用回帖:

月只蓝 at 2014-07-20 20:07:47

均方差(RMSE): 0.0166330356441843

残差平方和(SSE): 0.00138328937370353

相关系数(R): 0.992803768743302

相关系数之平方(R^2): 0.985659323230904

决定系数(DC): 0.985526033247728

卡方系数(Chi-Square): ...

现有以下4组数据要分别进行拟合,采用关系式：

y=A/(3*137.8)+A/x^0.5*2/a^0.5*atan((-0.5+0.5*(1+a*137.8/x^2)^0.5)^0.5)-2*2^0.5/3*A*137.8^0.5/x^1.5/(1+(1+a*137.8/x^2)^0.5)^1.5

x=[298,473,673,873,1073];

1y=[1.74024,1.63,1.49,1.4,1.36];

2y=[1.73395,1.58,1.45,1.37,1.34];

3y=[1.55275,1.51,1.41,1.34,1.3];

4y=[1.55959,1.52,1.42,1.35,1.33];

拟合之后把代码，参数，图发过来谢谢！万分感谢送你金币20个！

cfdc6361174d9cfc8169daed13483a0c.gif

月只蓝

引用回帖:

qj2006 at 2014-08-30 21:53:54

现有以下4组数据要分别进行拟合,采用关系式：

y=A/(3*137.8)+A/x^0.5*2/a^0.5*atan((-0.5+0.5*(1+a*137.8/x^2)^0.5)^0.5)-2*2^0.5/3*A*137.8^0.5/x^1.5/(1+(1+a*137.8/x^2)^0.5)^1.5

x=;

1y=;

2y=;

3y=;

4y= ...

第一组：

Variable x,y;

Function y=A1/(3*137.8)+A1/x^0.5*2/a^0.5*atan((-0.5+0.5*(1+a*137.8/x^2)^0.5)^0.5)-2*2^0.5/3*A1*137.8^0.5/x^1.5/(1+(1+a*137.8/x^2)^0.5)^1.5;

Data;

298        1.74024000000000

473        1.63000000000000

673        1.49000000000000

873        1.40000000000000

1073        1.36000000000000

均方差(RMSE): 0.0307537469550395

残差平方和(SSE): 0.00472896475887301

相关系数(R): 0.976679183783147

相关系数之平方(R^2): 0.953902228035314

决定系数(DC): 0.953374176421571

卡方系数(Chi-Square): 0.00151939891296933

F统计(F-Statistic): 62.079067211712

参数        最佳估算

----------        -------------

a1                 528.632286321485

a                  7309.23339319526

====== 结果输出 =====

No        实测值y        计算值y

1        1.74024        1.7635719

2        1.63        1.5764231

3        1.49        1.4752668

4        1.4        1.4201226

5        1.36        1.3863080

第二组：

Variable x,y;

Function y=A1/(3*137.8)+A1/x^0.5*2/a^0.5*atan((-0.5+0.5*(1+a*137.8/x^2)^0.5)^0.5)-2*2^0.5/3*A1*137.8^0.5/x^1.5/(1+(1+a*137.8/x^2)^0.5)^1.5;

Data;

298        1.73395000000000

473        1.58000000000000

673        1.45000000000000

873        1.37000000000000

1073        1.34000000000000

均方差(RMSE): 0.0180817726996743

残差平方和(SSE): 0.00163475251981344

相关系数(R): 0.992357331948071

相关系数之平方(R^2): 0.984773074271095

决定系数(DC): 0.984577352758231

卡方系数(Chi-Square): 0.000536567599397937

F统计(F-Statistic): 194.019415042202

参数        最佳估算

----------        -------------

a1                 514.710671370125

a                  5704.97616409959

====== 结果输出 =====

No        实测值y        计算值y

1        1.73395        1.7468852

2        1.58        1.5480091

3        1.45        1.4425041

4        1.37        1.3859055

5        1.34        1.3516123

第三组：

Variable x,y;

Function y=A1/(3*137.8)+A1/x^0.5*2/a^0.5*atan((-0.5+0.5*(1+a*137.8/x^2)^0.5)^0.5)-2*2^0.5/3*A1*137.8^0.5/x^1.5/(1+(1+a*137.8/x^2)^0.5)^1.5;

Data;

298        1.55275000000000

473        1.51000000000000

673        1.41000000000000

873        1.34000000000000

1073        1.30000000000000

均方差(RMSE): 0.0295054872943402

残差平方和(SSE): 0.00435286890238237

相关系数(R): 0.952330923708382

相关系数之平方(R^2): 0.906934188251261

决定系数(DC): 0.906570847157658

卡方系数(Chi-Square): 0.00149940028805953

F统计(F-Statistic): 29.2352531356999

参数        最佳估算

----------        -------------

a1                 510.648956991685

a                  21291.2744324343

====== 结果输出 =====

No        实测值y        计算值y

1        1.55275        1.5788868

2        1.51        1.4607601

3        1.41        1.3926929

4        1.34        1.3532335

5        1.3        1.3277581

第四组：

Variable x,y;

Function y=A1/(3*137.8)+A1/x^0.5*2/a^0.5*atan((-0.5+0.5*(1+a*137.8/x^2)^0.5)^0.5)-2*2^0.5/3*A1*137.8^0.5/x^1.5/(1+(1+a*137.8/x^2)^0.5)^1.5;

Data;

298        1.55959000000000

473        1.52000000000000

673        1.42000000000000

873        1.35000000000000

1073        1.33000000000000

均方差(RMSE): 0.0269603623925428

残差平方和(SSE): 0.00363430570168618

相关系数(R): 0.955014904538155

相关系数之平方(R^2): 0.912053467890022

决定系数(DC): 0.911828190854424

卡方系数(Chi-Square): 0.0012313365926436

F统计(F-Statistic): 31.1116349675785

参数        最佳估算

----------        -------------

a1                 518.480586764386

a                  25383.1634997523

====== 结果输出 =====

No        实测值y        计算值y

1        1.55959        1.5836851

2        1.52        1.4725106

3        1.42        1.4079438

4        1.35        1.3702032

5        1.33        1.3456515

w173h1122189_1409448993_123.jpg

附图1.jpg

w170h1122189_1409448993_995.jpg

附图2.jpg

w172h1122189_1409448994_214.jpg

附图3.jpg

w171h1122189_1409448996_851.jpg

附图4.jpg

，

0014f08102b5cebd97bfef69c59ca4cf.png

qj2006

引用回帖:

月只蓝 at 2014-08-31 09:36:39

第一组：

Variable x,y;

Function y=A1/(3*137.8)+A1/x^0.5*2/a^0.5*atan((-0.5+0.5*(1+a*137.8/x^2)^0.5)^0.5)-2*2^0.5/3*A1*137.8^0.5/x^1.5/(1+(1+a*137.8/x^2)^0.5)^1.5;

Data;

298        1.74024000000000

473 ...

万分感谢!

weixin_39933414
关注
  • 0
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值