| 查看: 1284 | 回复: 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小海» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有83人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
» 抢金币啦！回帖就可以得到:
-
安徽工程大学徐珍珍课题组招聘博士后
+1/424
-
澳门大学智慧城市物联网国重“基础设施智能监测”及“智慧城市”研究方向博士后招聘
+1/228
-
经济危机不是生产过剩
+1/137
-
郑州招聘博士
+1/86
-
四川大学先进制造技术团队招硕士博士
+1/84
-
原子层沉积（ALD）等微纳代工服务：18817872921
+1/83
-
四川大学卢伟课题组招收2027年秋季入学博士生2名
+1/79
-
南开大学环境科学与工程学院 环境与人体健康方向 招博士研究生、硕士研究生各1名
+1/78
-
北科大力学专业招生硕博研究生-能源工程偏微分方程应用及AI赋能工业软件开发方向
+2/68
-
华南理工大学黄哲昊课题组招收博士研究生
+2/58
-
中国科学院大学陈凯凯课题组招收2027级博士生
+2/50
-
香港城市大学（东莞）雷进程课题组博士研究生、研究助理招聘
+1/39
-
电子科技大学（深圳）高等研究院Beng Hau Tan 课题组现面向海内外公开招聘多名博士后
+1/33
-
澳大利亚皇家墨尔本理工招收phd
+1/32
-
华南理工大学 郎超课题组 招收2027级博士研究生
+1/29
-
安徽工程大学徐珍珍课题组招聘博士后
+1/28
-
安徽工程大学徐珍珍课题组招聘博士后
+1/28
-
北科大力学专业招生硕博研究生-能源工程偏微分方程应用及AI赋能工业软件开发方向
+2/14
-
社科后期
+1/10
-
蛋白互作实验怎么选？CoIP、Pull‑down、PLA 一文理清证据层级
+1/2
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