| 查看: 737 | 回复: 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小海» 猜你喜欢
郑州大学光电器件方向招收一名博士
已经有20人回复
-
氧化硅薄膜玻璃 求购
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有70人回复
-
vasp极化强度
已经有0人回复
-
vasp计算的一些问题
已经有0人回复
-
国科温州研究院 项目聘用有意“二战”的优秀本科毕业生
已经有21人回复
-
SCAPS-1D学习交流
已经有0人回复
» 抢金币啦！回帖就可以得到:
-
浙江工业大学生物工程学院招聘
+5/220
-
南方科技大学招聘行政助理
+1/82
-
恭喜老弟以第一作者发表三大顶刊文章！！！！
+1/74
-
博士论文盲审，接心软审稿人
+1/73
-
211博后出站求职，第一性原理+分子动力学+机器学习力场开发
+1/32
-
北航工程热物理、航空工程与环境团队招聘博后
+1/31
-
专业提供SEM-EDS/NMR/FTIR/DSC/TGA/DMA/TMA/GCMS 加扣870 559 66
+1/23
-
专业提供SEM-EDS/NMR/FTIR/DSC/TGA/DMA/TMA/GCMS 加扣870 559 66
+1/22
-
专业提供SEM-EDS/NMR/FTIR/DSC/TGA/DMA/TMA/GCMS 加扣870 559 66
+1/21
-
诚心寻找另一半，期待相亲缘分
+1/14
-
国家杰青、长江学者张卫东教授课题组招收科研助理启事
+1/14
-
山东第一医科大学纳米剂型与药物递送团队诚聘优秀博士后
+1/13
-
海归博士高层次人才协助提升科技论文质量，加快发表流程
+1/11
-
中西部211还是沿海二本求建议
+1/8
-
海归博士高层次人才协助提升科技论文质量，加快发表流程
+1/8
-
华南师范大学华南先进光电子研究院
+1/7
-
清华大学储凌课题组招聘有机化学及化学生物学方向博士后
+1/7
-
某知名农化企业招聘一位年薪30+的工艺组长（青岛）
+1/4
-
基因敲除
+1/3
-
中国科学院杭州医学所金程课题组招聘博士后
+1/2
+5/220
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
