meina861010

不好意思啊，下面是代碼，能幫忙看看是哪運行出錯了嗎？
data = {{1, 18.3}, {2, 28.6}, {3, 25}, {4, 15}, {5, 7}, {6, 3.3}, {7, 1.5}, {8, 0.7}, {9, 0.4}, {10, 0.2}, {11, 0.1}, {12, 0.001}, {13, 0.002}, {14, 0.001}, {15, 0.001}, {16, 0.001}};
nlm = NonlinearModelFit[data, (1/Sqrt[x! Cosh[r]] (1/2*Exp[I*\[Theta]]*Tanh[r])^(x/2)* Exp[-(1/2)*((Abs[\[Alpha]])^2 + (Conjugate[\[Alpha]])^2*Exp[I*\[Theta]]*Tanh[r])]*(1 + x + x^2 + x^3 + x^4 + x^5 + x^6 + x^7 + x^8 + x^9 + x^10 + x^11 + x^12 + x^13 + x^14 + x^15 + x^16)*((\[Alpha] + Conjugate[\[Alpha]]*Exp[I*\[Theta]]*Tanh[r])/(2 Exp[I*\[Theta]]*Tanh[r])^(1/2))), {\[Alpha], r, \[Theta]}, x];
Normal[nlm]
Show[ListPlot[data], Plot [nlm[x], {x, 0, 20}], Frame -> True]

贊一下

回復(fù)此樓

18樓2015-09-18 13:57:01

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

xzczd

木蟲 (小有名氣)

應(yīng)助: 72 (初中生)
金幣: 1941.4
紅花: 11
帖子: 290
在線: 144.2小時
蟲號: 3305534
注冊: 2014-07-03

引用回帖:

18樓: Originally posted by meina861010 at 2015-09-18 13:57:01
不好意思啊，下面是代碼，能幫忙看看是哪運行出錯了嗎？
data = {{1, 18.3}, {2, 28.6}, {3, 25}, {4, 15}, {5, 7}, {6, 3.3}, {7, 1.5}, {8, 0.7}, {9, 0.4}, {10, 0.2}, {11, 0.1}, {12, 0.001}, {13, 0.002}, ...

阿爾法和西塔的取值可能是復(fù)數(shù)？但是你用于擬合的原始數(shù)據(jù)又全是實的？這個模型真的沒有問題嗎？

如果真的沒有問題，那么：
1. 這大概又是模型本身的強非線性導(dǎo)致的擬合失敗，可能需要精心尋找合適的初值，關(guān)于這個你仔細看看FindFit（NonlinearModelFit其實是這個函數(shù)的殼兒）的自帶幫助的“可能存在的問題“一節(jié)。
2. 如果我沒記錯的話FindFit等一系列擬合函數(shù)實際上是默認擬合參數(shù)在實數(shù)域，所以如果你的參數(shù)是復(fù)數(shù)的話就需要把實部和虛部分開設(shè)。關(guān)于這個可以參看這個帖子：http://mathematica.stackexchange.com/a/43237/1871，里面有高手已經(jīng)寫了一個針對復(fù)數(shù)域擬合的程序包，你可以試試看。

不過，我稍微試了一下，不給初值的情況下擬合效果依舊欠佳，恐怕你這模型還是要找到合理的初值才有救。至于什么初值合適，我想你這個模型創(chuàng)造者肯定比我清楚。

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

meina861010

小木蟲Mathematica版塊已毀（當然原本也不咋的），建議大家前往百度貼吧或Stackexchange。

19樓2015-09-18 14:25:37

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

meina861010

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 257
紅花: 5
帖子: 138
在線: 45.1小時
蟲號: 1075736
注冊: 2010-08-16
性別: MM
專業(yè): 生物化學(xué)

送紅花一朵

引用回帖:

19樓: Originally posted by xzczd at 2015-09-18 14:25:37
阿爾法和西塔的取值可能是復(fù)數(shù)？但是你用于擬合的原始數(shù)據(jù)又全是實的？這個模型真的沒有問題嗎？

如果真的沒有問題，那么：
1. 這大概又是模型本身的強非線性導(dǎo)致的擬合失敗，可能需要精心尋找合適的初值，關(guān)于 ...

謝謝詳細的解答，我去看看，有問題再向你請教。

贊一下

回復(fù)此樓

20樓2015-09-18 14:42:33

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

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6	SUICHILD 2026-03-12	6/300	2026-03-17 09:24 by 霧散后相遇lc
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 機械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 材料371求調(diào)劑 +9	鱷魚? 2026-03-11	11/550	2026-03-13 22:53 by JourneyLucky
[考研] 一志愿中科院，化學(xué)方向，295求調(diào)劑 +4	一氧二氮 2026-03-11	4/200	2026-03-13 22:35 by JourneyLucky
[考研] 336求調(diào)劑 +6	Iuruoh 2026-03-11	6/300	2026-03-13 22:06 by JourneyLucky
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629