| 查看: 1063 | 回复: 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小海» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有198人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有11人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有6人回复
» 抢金币啦！回帖就可以得到:
-
当新能源不再是一个选择题
+1/90
-
重庆大学刘坤教授课题组招聘等离子体农业应用相关方向科研助理（农学、等离子体、电子
+1/85
-
五邑大学 化学工程与技术调剂 4月10号九点截止！！！！！
+1/41
-
浙江大学高分子系毛峥伟教授招收科研助理
+1/38
-
山东航空学院安全工程专业085702第三轮接收研究生调剂，欢迎咨询
+1/36
-
大家好，请教个问题，连续挤出胶皮，有没有在线降低气味VOC的方式。
+1/20
-
大连海洋大学资源与环境专硕（0857）接收调剂学生
+1/19
-
福建师范大学环资学院环境微生物技术组招收2026级硕士（4月10日，线上面上）
+1/11
-
省重点实验室招收调剂生，海洋生物专业名额充足
+1/10
-
省重点实验室招收调剂生，海洋生物专业，名额充足，可快速安排复试并出结果
+1/9
-
重庆大学附属涪陵医院2026年度全职博士后招聘简章（年薪30-65万）
+1/8
-
河南农大刘立杰课题组化学学硕07调剂，线上面试
+1/8
-
赵松方教授联合山大于伟泳杰青、延世大Jong-Hyun Ahn院士招博士后 、博士研究生
+1/8
-
二次调剂名额13名，华北理工大学矿物加工工程方向接收调剂研究生
+1/7
-
Postdoc position in continuum mechanics of brain tissue at KTH
+1/5
-
鲁东大学水产学院，水产学硕接收调剂3-4人
+1/5
-
长江师范学院化学化工学院化学工程方向有大量调剂名额
+1/2
-
上海市过程安全与智能控制重点实验室招收安全工程研究生调剂-上岸既是高质量就业！
+1/2
-
山西大学招收2026级申请考核制化学材料类博士研究生一名
+1/1
-
昆明理工大学陈江照教授团队招收博士研究生3名
+1/1
4楼2019-03-26 20:19:44
2楼2019-03-26 13:50:09
3楼2019-03-26 19:51:16
5楼2019-03-27 20:41:08
