 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 719 | 回復(fù): 4 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||||
Charles201新蟲 (小有名氣)
|
[求助]
求助matlab高手看一下菜鳥寫的非線性擬合程序，有關(guān)苯胺絡(luò)合萃取的 已有2人參與
|
||||
|
待擬合方程 原始數(shù)據(jù)： x1 x2 y 0.10 3.27 6.49 0.20 2.99 9.10 0.30 2.86 12.65 0.40 2.77 16.96 0.50 2.73 22.88 0.60 2.70 28.12 已寫的程序，一直出錯，求各位大神幫我修改一下： m文件： function y=aniline_1fun(b,x) x1=x(0:0.01:1.00); x2=x(4.00:-0.02:2.00); b1=m; b2=KE; y=b1*b2*0.0104^(b1-1)*x1/(1+b2*0.0104^b1)/(1+10^(4.67-x2)); >> x=[0.10,0.20,0.30,0.40,0.50,0.60; 3.27,2.99,2.86,2.77,2.73,2.70]'; >> y=[6.49,9.10,12.65,16.96,22.88,28.12]; >> b0=[0,5]'; >> [b,Q]=lsqcurvefit(aniline_1fun,b0,x,y) |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有72人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有68人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
Charles201: 金幣+20, ★★★很有幫助 2015-08-05 15:14:38
感謝參與，應(yīng)助指數(shù) +1
Charles201: 金幣+20, ★★★很有幫助 2015-08-05 15:14:38
|
dingd大神這么耐心地講解，我直接給你結(jié)果吧。 我這些結(jié)果是用1stOpt軟件做的，MATLAB肯定也可以做，在這里不我討論MATLAB代碼。 dingd大神說的過擬合很有道理的。 如果按照雙參數(shù)公式： y=b1*b2*0.0104^(b1-1)*x1/(1+b2*0.0104^b1)/(1+10^(4.67-x2)) 將會得到無限組結(jié)果： b1 -0.480056251311531 b2 -0.110217887023378 ---------- ------------- b1 -1.43102339222656 b2 -0.0013970972484279 ---------- ------------- b1 -0.594292430579858 b2 -0.0652159450286256 ---------- ------------- b1 -0.438212028430683 b2 -0.13357689181283 ---------- ------------- b1 1.03961552899074 b2 -118.681614122118 …… …… 盡管b1和b2變化不定，但是兩個參數(shù)按照下式計算： b3=b1*b2*0.0104^(b1-1)/(1+b2*0.0104^b1) b3一定都是等于3425.94298397334。 也就是說，樓主最原始的公式是過擬合的。 |
5樓2015-08-04 11:17:06
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
Charles201: 金幣+30, ★★★★★最佳答案 2015-08-05 15:14:26
感謝參與，應(yīng)助指數(shù) +1
Charles201: 金幣+30, ★★★★★最佳答案 2015-08-05 15:14:26
|
公式有問題吧，過擬合： y=b1*b2*0.0104^(b1-1)*x1/(1+b2*0.0104^b1)/(1+10^(4.67-x2)); 可寫成： y=b1*b2*0.0104^(b1-1)/(1+b2*0.0104^b1)/(1+10^(4.67-x2))*x1; 上式中：b1*b2*0.0104^(b1-1)/(1+b2*0.0104^b1)可以用一個參數(shù)代替，如b3 等價公式：y=b3/(1+10^(4.67-x2))*x1; 僅有一個參數(shù) |
» 本帖已獲得的紅花（最新10朵）
Charles2012樓2015-08-03 18:22:37
Charles201
新蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 31.2
- 帖子: 74
- 在線: 24.3小時
- 蟲號: 2670335
- 注冊: 2013-09-22
- 專業(yè): 分離過程
3樓2015-08-03 18:26:06
Charles201
新蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 31.2
- 帖子: 74
- 在線: 24.3小時
- 蟲號: 2670335
- 注冊: 2013-09-22
- 專業(yè): 分離過程
|
我將m文件改成了如下形式 function y=aniline_1fun(b,x) x1=x(0:0.01:1.00); x2=x(4.00:-0.02:2.00); b3=b1*b2*0.0104^(b1-1)/(1+b2*0.0104^b1); y=b3*x1/(1+10^(4.67-x2)); 主程序： x=[0.10,0.20,0.30,0.40,0.50,0.60; 3.27,2.99,2.86,2.77,2.73,2.70]'; y=[6.49,9.10,12.65,16.96,22.88,28.12]; b3=[0]'; [b,Q]=lsqcurvefit(aniline_1fun,b0,x,y) 但是運行后提示輸入?yún)?shù)的數(shù)目不足。請問這是什么原因呢？ |
4樓2015-08-03 19:39:23
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +5 | 生物工程調(diào)劑 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
|
[論文投稿]
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|