版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

lavender_12

新蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 6.6
散金: 72
紅花: 3
帖子: 407
在線: 267.6小時
蟲號: 1054230
注冊: 2010-07-08
性別: MM
專業(yè): 環(huán)境污染化學(xué)

[求助] 求拐點

我在做高斯擬合時，手動擬合的，求拐點。程序如下：
a= xlsread('E:\UV-vis\2011414\dong 20100831 klto','dong 20100831 klto','f3:g563');%0831klto 以abs 中hv（1240/波長）為橫坐標(biāo)，以（2.303*abs*hv）^2 為縱坐標(biāo)計算
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值了，我該怎么做？

回復(fù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有53人回復(fù)
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)

1樓 2012-05-24 15:18:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

先畫圖知道x的大概范圍，然后使用
evalin(symengine, 'numeric::solve(Y), x =0...10)')
獲得該范圍內(nèi)的解。
solve好像默認(rèn)只給最先找到的解

贊一下

回復(fù)此樓

showmethemoney

2樓2012-05-24 16:32:14