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

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

24小時(shí)熱門(mén)版塊排行榜

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

返回列表

lavender_12

新蟲(chóng) (正式寫(xiě)手)

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

[求助] 求拐點(diǎn)

我在做高斯擬合時(shí)，手動(dòng)擬合的，求拐點(diǎn)。程序如下：
a= xlsread('E:\UV-vis\2011414\dong 20100831 klto','dong 20100831 klto','f3:g563');%0831klto 以abs 中hv（1240/波長(zhǎng)）為橫坐標(biāo)，以（2.303*abs*hv）^2 為縱坐標(biāo)計(jì)算
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] %
那么這時(shí)候會(huì)求出來(lái)一個(gè)c = [ 3.6925373244157819483595872056264, 103.34849107643072995968936356082]，但是有可能有2個(gè)拐點(diǎn)，那么第二個(gè)求出來(lái)的c其實(shí)掩蓋了第一個(gè)求出來(lái)的值，這樣，我就沒(méi)法知道第一個(gè)c值了，我該怎么做？

回復(fù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有247人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收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

榮譽(yù)版主 (著名寫(xiě)手)

彩色掛圖

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

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

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

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

贊一下

回復(fù)此樓

showmethemoney

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

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 lavender_12 的主題更新