 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1887 | 回復(fù): 1 | ||
ponder_zsm新蟲 (初入文壇)
|
[求助]
如何用給定的帶有參數(shù)的函數(shù)擬合已知數(shù)據(jù)？ 已有1人參與
|
|
已知實(shí)驗(yàn)數(shù)據(jù)和函數(shù)原型如下： data = {{8.5, 0.828}, {12.5, 0.78}, {16.5, 0.741}, {20.4, 0.714}, {23.1, 0.66}, {29.4, 0.558}, {39.4, 0.491}, {49.4, 0.454}}; model = Integrate[(1 - x/r)^2*Exp[-(Log[r/b*Sqrt[1 + (a*a)/(b*b)]])^2/(2 Log[1 + (a*a)/(b*b)])]/(r*Sqrt[2*Pi*Log[1 + (a*a)/(b*b)]]), {r,0, Infinity}]; 這是在Mathematica中輸入的格式，函數(shù)原型可以參見上傳的附件。其中a, b是參數(shù)，x是自變量。用Findfit求a和b時，Mathematica給出的結(jié)果不對： {a -> 160.339, b -> 186.173} 合理值應(yīng)該在a=299, b=247左右，在用所給值驗(yàn)證時，所得到的函數(shù)值為虛數(shù)： f[a_, b_, x_, r_] :=Integrate[(1 - x/r)^2*Exp[-(Log[r/b*Sqrt[1 + (a*a)/(b*b)]])^2/(2 Log[1 + (a*a)/(b*b)])]/(r*Sqrt[2*Pi*Log[1 + (a*a)/(b*b)]]), {r,0, Infinity}]; f[160.339, 186.173, 8.5, r] 輸出為：0.851972 - 4.76418*10^-15 I 不知道哪里可以改進(jìn)，求高人指點(diǎn)！ 另外我嘗試用1st opt，但是也不成功： Constant h=10000; Parameters a,b; Variables x,y; Function y=int(sqr(1-x/r)*exp(-sqr(ln(r*sqrt(1+(a*a)/(b*b))/b))/(2*ln(1+(a*a)/(b*b))))/(r*sqrt(2*Pi*ln(1+(a*a)/(b*b)))),r=0,h); Data; 8.5 0.828 12.5 0.78 16.5 0.741 20.4 0.714 23.1 0.66 29.4 0.558 39.4 0.491 49.4 0.454 提示函數(shù)表達(dá)式錯誤！求指點(diǎn)！！ 函數(shù).png |
回復(fù)此樓» 猜你喜歡
本人考085602 化學(xué)工程 專碩
已經(jīng)有23人回復(fù)
-
求調(diào)劑院校信息
已經(jīng)有4人回復(fù)
-
085600材料與化工306
已經(jīng)有4人回復(fù)
-
286求調(diào)劑
已經(jīng)有10人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有4人回復(fù)
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有8人回復(fù)
-
資源與環(huán)境 調(diào)劑申請(333分)
已經(jīng)有5人回復(fù)
-
280求調(diào)劑
已經(jīng)有12人回復(fù)
-
269專碩求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- origin8.0 函數(shù)關(guān)系中參數(shù)擬合求助？
已經(jīng)有4人回復(fù)
- 急！在用origin進(jìn)行非線性擬合怎樣創(chuàng)建新的擬合函數(shù)？
已經(jīng)有14人回復(fù)
- origin 里 “三次樣條函數(shù)擬合” 與“三次多項式擬合”是一回事嗎？如圖
已經(jīng)有3人回復(fù)
- 想用GULP擬合純鉬的morse勢函數(shù)參數(shù)
已經(jīng)有13人回復(fù)
- 求助該圖形用什么函數(shù)曲線擬合，
已經(jīng)有3人回復(fù)
- Aguado-Paniagua函數(shù)擬合
已經(jīng)有6人回復(fù)
- 1stopt 擬合隱函數(shù)的四參數(shù)
已經(jīng)有7人回復(fù)
- 如何用已知函數(shù)去擬合實(shí)驗(yàn)數(shù)據(jù)，得出參數(shù)
已經(jīng)有10人回復(fù)
- 關(guān)于matlab擬合
已經(jīng)有4人回復(fù)
- 在用origin進(jìn)行非線性擬合怎樣去選擇擬合函數(shù)？
已經(jīng)有6人回復(fù)
- 通過實(shí)驗(yàn)數(shù)據(jù)擬合，求解公式中的參數(shù)
已經(jīng)有14人回復(fù)
- 求助一個邏輯斯蒂函數(shù)的擬合
已經(jīng)有8人回復(fù)
- 自定義函數(shù)擬合 請假大俠幫忙！非常非常感激不盡！�。�
已經(jīng)有9人回復(fù)
- matlab高手幫忙做個自定義復(fù)雜函數(shù)擬合，含有復(fù)數(shù)，
已經(jīng)有25人回復(fù)
- 已知函數(shù)形式，如何擬合其中的多個系數(shù)。
已經(jīng)有3人回復(fù)
- matlab的殘差大小跟什么因素有關(guān)？擬合的函數(shù)選擇？還是給定值的范圍？
已經(jīng)有4人回復(fù)
- 跪求1stOpt2.5以上版本，或者請高人幫忙擬合包含復(fù)數(shù)的函數(shù)的參數(shù)
已經(jīng)有8人回復(fù)
- origin進(jìn)行非線性擬合如何選擇函數(shù)或者自定義合適的函數(shù)
已經(jīng)有19人回復(fù)
- 請問一個有理函數(shù)的曲線擬合過程如何做，給出數(shù)據(jù)，請求擬合，
已經(jīng)有14人回復(fù)
- 【公式擬合問題】如何擬合下面這樣的一個公式，多謝大牛。
已經(jīng)有6人回復(fù)
- 【求助】急求關(guān)于matlab做非線性擬合的函數(shù)
已經(jīng)有8人回復(fù)
- 【求助】matlab函數(shù)擬合，擬合結(jié)果并保存，方便調(diào)用
已經(jīng)有4人回復(fù)
【答案】應(yīng)助回帖
|
樓主你確定真值在{a->299, b->247}附近？ 下圖是{a -> 160.339, b -> 186.173}時的擬合效果圖： 這擬合效果不算差。 下圖是你給的參數(shù)對應(yīng)的圖： 囧。 然后，關(guān)于直接反代導(dǎo)致虛數(shù)的問題：這只是一點(diǎn)數(shù)值誤差罷了，如你所見，虛部的值非常微小。想要避免這個問題，你可以把全體參數(shù)都改成準(zhǔn)確值（也就是別用小數(shù)），也可以把Integrate改成NIntegrate然后再加個高一點(diǎn)的WorkingPrecision。 |
2樓2015-01-03 21:41:26
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ponder_zsm 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 本人考085602 化學(xué)工程 專碩 +21 | 不知道叫什么！ 2026-03-15 | 23/1150 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +3 | vv迷 2026-03-21 | 4/200 |
|
|
|
[考研] 269專碩求調(diào)劑 +5 | 金恩貝 2026-03-21 | 5/250 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +6 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|