| 查看: 822 | 回复: 1 | ||
lavender_12新虫 (正式写手)
|
[求助]
求拐点
|
|
我在做高斯拟合时，手动拟合的，求拐点。程序如下： a= xlsread('E:\UV-vis\2011414\dong 20100831 klto','dong 20100831 klto','f3:g563');%0831klto 以abs 中hv（1240/波长）为横坐标，以（2.303*abs*hv）^2 为纵坐标计算 plot(a(:,1),a(:,2)) cftool(a(:,1),a(:,2)) syms x a1 = -32.75 ; b1 = 5.169 ; c1 = 0.001421 ; a2 = 215.7 ; b2 = 5.317 ; c2 = 0.1951 ; a3 = 227 ; b3 = 5.057 ; c3 = 0.2606 ; a4 = 168.8 ; b4 = 4.735 ; c4 = 0.2747 ; a5 = 156.4 ; b5 = 4.391 ; c5 = 0.3278 ; a6 = 0 ; b6 = 4.29 ; c6 = 6.22e-005 ; a7 = 151.7 ; b7 = 3.947 ; c7 = 0.4023 ; y=a1*exp(-((x-b1)/c1)^2) + a2*exp(-((x-b2)/c2)^2) + ... a3*exp(-((x-b3)/c3)^2) + a4*exp(-((x-b4)/c4)^2) + ... a5*exp(-((x-b5)/c5)^2) + a6*exp(-((x-b6)/c6)^2) + ... a7*exp(-((x-b7)/c7)^2) Y=diff(y,x,2); x=solve(Y,x); y=eval(y); c=[x,y] % 那么这时候会求出来一个c = [ 3.6925373244157819483595872056264, 103.34849107643072995968936356082]，但是有可能有2个拐点，那么第二个求出来的c其实掩盖了第一个求出来的值，这样，我就没法知道第一个c值了，我该怎么做？ |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有130人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
2楼2012-05-24 16:32:14
