| 查看: 797 | 回复: 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值了，我该怎么做？ |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有92人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
2楼2012-05-24 16:32:14
