| 查看: 1098 | 回复: 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小海» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有214人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有0人回复
» 抢金币啦！回帖就可以得到:
-
青基祈福
+1/118
-
西北农林科技大学朱杰教授团队招收 生物物理学专业2026级博士生1-2名（5.15报名截止）
+1/82
-
南方科技大学物理系 CSLM Lab 招生招聘 （软物质/生物物理，Physics of AI）
+1/42
-
谷歌学术引用量增长
+1/39
-
谷歌学术引用量增长
+1/36
-
华南师范大学（211）- 生物光子学研究院招收生物学、物理学、神经科学等专业博士
+1/36
-
北京航空航天大学拟招收2026博士研究生
+1/29
-
苏州科技大学国家高层次青年人才团队诚聘副教授/讲师（事业编）和科研助理数名
+1/27
-
双一流高校-南京林业大学-化学工程学院-国家海外优青团队招青年师资（正式事业编制）
+1/24
-
【陕西师范大学】催化化学课题组2026年招收博士后/讲师/副高
+1/17
-
南昌大学 物理与材料学院招收2026年申请考核博士1名 化学、材料、催化相关专业
+1/13
-
海南大学生物医学工程学院 光学诊疗团队 博士招生
+1/13
-
求siesta、Transiesta教程资料
+1/10
-
211高校2026年入学博士名额
+1/7
-
哈尔滨工业大学（深圳）——急招能动/材料/机械——26年9月份博士
+1/7
-
博士后招聘（高薪40万+）
+1/5
-
【招收2026年秋博士研究生】哈工大（深圳）医工学院任卫红老师（15号截止！！！）
+1/5
-
请教：国际会议报账的问题
+1/3
-
【已入职｜急需医学类（全路径）博士、】有偿访谈招募
+1/3
-
坐标济南——高校新药研发科研岗教师招聘——招聘化学合成方向硕士
+1/1
2楼2019-03-26 13:50:09
3楼2019-03-26 19:51:16
4楼2019-03-26 20:19:44
5楼2019-03-27 20:41:08
6楼2019-03-28 11:22:35