| 查看: 981 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
求助各位大神，参数拟合结果不对，金币全送。谢谢
|
|||
 求助各位大神，用下面的程序做的参数拟合，拟合结果和图形都不对，希望大神们帮帮忙，金币全送。谢谢function fit_nonl12 clear all;clc format long data=[ 2561.0000 10.0000 2560.0000 15.0000 2560.0000 25.0000 2558.0000 40.0000 2557.0000 60.0000 2556.0000 100.0000 2555.0000 160.0000 2554.0000 250.0000 2552.0000 400.0000 2551.0000 630.0000 2550.0000 1000.0000 ]; xspan=data(:,2); Texp=data(:,1); k0=[2.0 30.0 6.0 0.5]; lb=[0.1 0.1 0.1 0.0]; ub=[10.0 100.0 20.0 1.0]; % OPTIONS=optimset('MaxFunEvals',1000); [k,resnorm,residual,exitflag,output,lambda,jacobian]=... lsqnonlin(@ObjFunc,k0,lb,ub,OPTIONS,xspan,Texp); ci=nlparci(k,residual,jacobian); %residual; fprintf('\n\n拟合结果;\n') fprintf('\n\t参数 E1 = %.16f',k(1)) fprintf('\n\t参数 E2 = %.16f',k(2)) fprintf('\n\t参数 E3 = %.16f',k(3)) fprintf('\n\t参数 E4 = %.16f',k(4)) y=kineticsEqs(xspan,k); R2=1-sum((Texp-y).^2)./sum((Texp-mean(y)).^2); fprintf('\n\tR^2 = %.16f',R2); figure plot(xspan,kineticsEqs(xspan,k),'b',xspan,Texp,'or'),legend('计算值','实验值','Location','Best') % function f =ObjFunc(k,xspan,Texp) f=kineticsEqs(xspan,k)-Texp; function x1 =kineticsEqs(x,k) E1=k(1);E2=k(2);E3=k(3);E4=k(4); x1=E1+(((2.*pi.*E3).^(-E4)+x.^(E4).*cos(pi.*E4./2)).*E2.*(2.*pi.*E3).^(-E4))./(((2.*pi.*E3).^(-E4)+x.^(E4).*cos(pi.*E4./2)).^2+(x.*(E4).*sin(pi.*E4./2)).^2); |
回复此楼» 本帖已获得的红花（最新10朵）
why小海» 猜你喜欢
投稿chemical physical letters不送审？
已经有2人回复
-
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有226人回复
-
杂志要求提供official language-editing certificate
已经有4人回复
-
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有16人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
» 抢金币啦！回帖就可以得到:
-
南方医科大学中药学院 申请考核博士一名 （天然药化方向，天然产物分离经验优先）
+1/278
-
电力全国重点实验室双一流A类长江学者团队招2026年全日制博士1-2名/博士后
+2/84
-
博后平台选择
+1/70
-
中国石油大学（北京）国家级大人才团队博士招生2名：化学、材料、石油工程：油田化学
+2/67
-
香港科技大学计算物理及流体力学课题组招收全奖博士后及博士生（2026年9月入学）
+1/45
-
南方医科大学生物医学工程学院招收申请考核制博士生2名、博士后2名（2026）
+1/32
-
上海市“光探测材料与器件”工程技术研究中心（上海应用技术大学）招聘优秀研究人员
+1/28
-
教育部重点实验室和清华大学某国家重点实验室，联合培养硕生、博生，并长期招博士后
+1/27
-
武汉工程大学绿碳技术与智能材料课题组诚招2026年博士研究生
+2/26
-
意大利米兰理工大学急聘CSC公派留学博士生（物理或无机材料科学方向）
+2/18
-
太原理工大学集成电路学院院长团队招收2026年博士研究生
+1/9
-
【博士招生】哈工大（深圳）智能学部机器人与先进制造学院 陆文杰老师课题组
+1/7
-
《中文期刊点评》这个模块，怎么不能点评期刊了呢自动跳转到主页了
+1/5
-
海南大学化学院—功能分子器件团队2026博士/研究助理招生+博士后招聘
+1/5
-
南京林业大学国家级青年人才团队招收2026年生物质转化/炭材料/储能等方向博士生
+1/4
-
邵阳学院食品与化学工程学院硕士调剂
+1/4
-
上海理工大学顾敏院士、张轶楠教授团队 招聘 2026级 光学工程 博士生
+2/2
-
澳科大诚招2026年秋季生物材料全奖博士研究生（今日16:30线上宣讲会）
+1/1
-
苏州大学医学院纳米生物医学方向招收申请考核制博士生1名
+1/1
-
海南大学国家优青团队招聘“AI/大数据+材料”方向专任教师（事业编制）
+1/1
6楼2019-03-28 11:22:35
2楼2019-03-26 13:50:09
3楼2019-03-26 19:51:16
4楼2019-03-26 20:19:44