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

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

qj2006

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 25.6
散金: 486
红花: 1
帖子: 1111
在线: 265.4小时
虫号: 2930214
注册: 2014-01-13
性别: GG
专业: 材料物理化学

[求助] matlab 多元函数拟合已有1人参与

现有以下数据要进行拟合,采用关系式：
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.51941,1.47,1.38,1.32,1.28];
拟合之后把代码，参数，图发过来谢谢！

回复此楼

» 猜你喜欢

0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有2人回复
070200求调剂，一志愿某211，288分已经有18人回复
物理学I论文润色/翻译怎么收费? 已经有285人回复
法国博士后职位已经有0人回复
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。已经有0人回复
** 已经有1人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有3人回复
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？已经有0人回复
B区学生调剂-兰州交通大学材料科学与工程学院已经有2人回复

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

1楼 2014-07-20 17:26:57

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68129.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

引用回帖:

6楼: Originally posted by 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

附图1.jpg

附图2.jpg

附图3.jpg

附图4.jpg

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

qj2006

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

7楼2014-08-31 09:36:39

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

查看全部 10 个回答

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68129.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
qj2006: 金币+20, ★★★★★最佳答案 2014-09-05 10:36:21

一、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

附图1.png

赞一下(1人)

回复此楼

» 本帖已获得的红花（最新10朵）

山鹰HPU

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2014-07-20 19:57:04

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

qj2006

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 25.6
散金: 486
红花: 1
帖子: 1111
在线: 265.4小时
虫号: 2930214
注册: 2014-01-13
性别: GG
专业: 材料物理化学

引用回帖:

2楼: Originally posted by 月只蓝 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];

赞一下

回复此楼

3楼2014-07-20 20:04:45

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68129.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

引用回帖:

3楼: Originally posted by 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

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

4楼2014-07-20 20:07:47

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

查看全部 10 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 一志愿211，化学310分，本科重点双非，求调剂 +13	努力奋斗112 2026-04-08	13/650	2026-04-08 21:17 by 学员tURuqU
[考研] 275求调剂 +4	1624447980 2026-04-08	5/250	2026-04-08 15:29 by 哦哦123
[考研] 求调剂 +9	月@163.com 2026-04-07	11/550	2026-04-08 14:48 by qlm5820
[考研] 275 求调剂 +8	Lei812514 2026-04-07	8/400	2026-04-08 12:46 by chemisry
[考研] 材料专硕322分 +12	哈哈哈吼吼吼哈 2026-04-02	12/600	2026-04-08 11:43 by 1753564080
[考研] 388求调剂 +6	四川王涛 2026-04-07	8/400	2026-04-08 00:17 by JourneyLucky
[考研] 农学，求调剂，314分 +4	访客记录可爱 2026-04-04	4/200	2026-04-07 21:07 by 等岸
[考研] 化工调剂303分，过四级 +34	栖梧待风 2026-04-02	34/1700	2026-04-07 12:30 by 1018329917
[考研] 286求调剂 +20	Faune 2026-04-06	20/1000	2026-04-07 11:33 by 诗与自由
[考研] 08600生物与医药-327 +9	18755400796 2026-04-05	9/450	2026-04-06 22:35 by 52305043001
[考研] 化学357分，考研调剂 +11	.Starry. 2026-04-04	12/600	2026-04-06 06:28 by houyaoxu
[考研] 材料专硕322分 +10	哈哈哈吼吼吼哈 2026-04-04	10/500	2026-04-05 21:22 by 学员8dgXkO
[考研] 计算机11408，286分求调剂 +7	木子念晞 2026-04-05	7/350	2026-04-05 19:02 by chy09050039
[考研] 一志愿江南大学085501机械工程专硕326分，本科佳木斯大学 +5	顾若浮生 2026-04-03	9/450	2026-04-05 09:57 by 1753564080
[考研] 调剂 +4	是可乐不是可乐 2026-04-04	4/200	2026-04-04 19:41 by 唐沐儿
[考研] 338求调剂 +7	晟功? 2026-04-03	7/350	2026-04-03 16:46 by wxiongid
[考研] 081200-11408-276学硕求调剂 +5	崔wj 2026-04-03	5/250	2026-04-03 15:06 by arrow8852
[考研] 320求调剂 +3	农业工程与信息� 2026-04-03	3/150	2026-04-03 11:40 by 土木硕士招生
[考研] 求调剂 +7	Aniyaio 2026-04-02	7/350	2026-04-02 16:42 by zzsw+
[考研] 266求调剂 +4	学员97LZgn 2026-04-02	4/200	2026-04-02 13:03 by yulian1987

24小时热门版块排行榜

qj2006

[求助] matlab 多元函数拟合 已有1人参与

» 猜你喜欢

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

月只蓝

【答案】应助回帖

» 本帖已获得的红花（最新10朵）

月只蓝

【答案】应助回帖

» 本帖已获得的红花（最新10朵）

qj2006

月只蓝

【答案】应助回帖

[求助] matlab 多元函数拟合已有1人参与