zhangxiaodao

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1929.4
散金: 31
紅花: 1
帖子: 228
在線: 286小時
蟲號: 871428
注冊: 2009-10-13

[求助] 根據(jù)數(shù)據(jù)和方程的形式，求方程中的參數(shù)

RT，數(shù)據(jù)如下列，然后方程的形式為y=a*exp[-p/(x+q)],求方程中的p和q，謝謝各位大俠，急求啊~~
x=1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
2
2.1
2.2
2.3
2.4
y=10363972.36
5697599.795
3269017.372
1948866.775
1202604.284
765578.5464
501320.0508
336809.4864
231640.3371
162754.7914
116618.904
85080.3473
63109.81473
47535.43773

回復此樓

» 猜你喜歡

287求調劑已經有5人回復
22408 344分求調劑一志愿華電計算機技術已經有3人回復
289求調劑已經有3人回復
317求調劑已經有10人回復
廣西大學家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調劑）已經有3人回復
學校已經提交到NSFC，還能修改嗎？已經有8人回復
08工學調劑已經有5人回復
281求調劑（0805）已經有25人回復
085600材料與化工已經有6人回復
265求調劑已經有8人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

【500金幣本版活動】“經驗力場分子模擬+第一性原理計算問題討論互助及資源共享” 已經有82人回復
【求助】matlab 求解微分方程中的未知參數(shù) 已經有20人回復
求非線性方程中的參數(shù)問題已經有13人回復

1樓 2011-12-19 22:50:41

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

715211229

木蟲 (正式寫手)

應助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時
蟲號: 1357002
注冊: 2011-07-27
性別: GG
專業(yè): 傳動機械學

★ ★
臭水溝(金幣+2): 謝謝交流~~ 2011-12-20 11:03:15

matlab ：ctool工具箱
%%%擬合效果很差

Fit computation did not converge:
Maximum number of function evaluations exceeded. Increasing
MaxFunEvals (in fit options) may allow for a better fit, or
the current equation may not be a good model for the data.

Fit found when optimization terminated:

General model:
   f(x) = a*exp(-p/(x+q))
Coefficients (with 95% confidence bounds):
   a =    11.28  (-2.505e+004, 2.507e+004)
   p =    -7.497  (-3191, 3176)
   q =    -0.4499  (-151.6, 150.7)

Goodness of fit:
  SSE: 1.309e+014
  R-square: -0.159
  Adjusted R-square: -0.3697
  RMSE: 3.45e+006

贊一下(1人)

回復此樓

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

zhangxiaodao

我是蝸牛

2樓2011-12-20 05:28:41

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

715211229

木蟲 (正式寫手)

應助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時
蟲號: 1357002
注冊: 2011-07-27
性別: GG
專業(yè): 傳動機械學

cftool工具箱

回復此樓

我是蝸牛

3樓2011-12-20 05:29:19

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

zhangxiaodao

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1929.4
散金: 31
紅花: 1
帖子: 228
在線: 286小時
蟲號: 871428
注冊: 2009-10-13

送鮮花一朵

引用回帖:

樓: Originally posted by 715211229 at 2011-12-20 05:28:41:
matlab ：ctool工具箱
%%%擬合效果很差

Fit computation did not converge:
Maximum number of function evaluations exceeded. Increasing
MaxFunEvals (in fit options) may allow ...

謝謝

回復此樓

4樓2011-12-20 09:00:53

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

715211229

木蟲 (正式寫手)

應助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時
蟲號: 1357002
注冊: 2011-07-27
性別: GG
專業(yè): 傳動機械學

★
臭水溝(金幣+1): 謝謝交流~~ 2011-12-20 11:03:26

你看下這個帖子，我這個可能是錯誤的
http://www.gaoyang168.com/bbs/viewthread.php?tid=3866180

贊一下(1人)

回復此樓

我是蝸牛

5樓2011-12-20 10:48:24

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

dbb627

榮譽版主 (著名寫手)

專家經驗: +4
仿真EPI: 6
應助: 289 (大學生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學
管轄: 計算模擬

【答案】應助回帖

★
感謝參與，應助指數(shù) +1
zhangxiaodao(金幣+20): 2011-12-20 13:52:42
xiegangmai(金幣+1): 謝謝應助！ 2011-12-20 23:02:49

General model:
   f(x) = a*exp(-p/(x+q))
Coefficients (with 95% confidence bounds):
   a =          1  (1, 1)
   p =       -42  (-42, -42)
   q =       1.5  (1.5, 1.5)

Goodness of fit:
  SSE: 1.767e-007
  R-square: 1
  Adjusted R-square: 1
  RMSE: 0.0001267

初值a=10000 p=-2  q=-1

贊一下(1人)

回復此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

6樓2011-12-20 11:22:06

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

dbb627

榮譽版主 (著名寫手)

專家經驗: +4
仿真EPI: 6
應助: 289 (大學生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學
管轄: 計算模擬

★
xiegangmai(金幣+1): 辛苦了！ 2011-12-20 23:03:02

引用回帖:

5樓: Originally posted by 715211229 at 2011-12-20 10:48:24:
你看下這個帖子，我這個可能是錯誤的
http://www.gaoyang168.com/bbs/viewthread.php?tid=3866180

圖如下

贊一下(1人)

回復此樓

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

zhangxiaodao

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

7樓2011-12-20 11:23:32

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

zhangxiaodao

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1929.4
散金: 31
紅花: 1
帖子: 228
在線: 286小時
蟲號: 871428
注冊: 2009-10-13

送鮮花一朵

引用回帖:

樓: Originally posted by dbb627 at 2011-12-20 11:23:32:
圖如下

非常感謝，您擬合的太棒了~~
ps: 怎么送金幣呢？然后能講程序發(fā)給我嗎，我是MATLAB門外漢。。。。

贊一下

回復此樓

8樓2011-12-20 13:48:57

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

dbb627

榮譽版主 (著名寫手)

專家經驗: +4
仿真EPI: 6
應助: 289 (大學生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學
管轄: 計算模擬

CODE:

 x=[1.1

1.2

1.3

1.4

1.5

1.6

1.7

1.8

1.9

2

2.1

2.2

2.3

2.4];

y=[10363972.36

5697599.795

3269017.372

1948866.775

1202604.284

765578.5464

501320.0508

336809.4864

231640.3371

162754.7914

116618.904

85080.3473

63109.81473

47535.43773];

st_ = [1000 -2 -1];

ft_ = fittype('a*exp(-p/(x+q))','dependent',{'y'},'independent',{'x'},'coefficients',{'a', 'p','q'});

[cf_,good]= fit(x,y,ft_ ,'Startpoint',st_)

h_ = plot(cf_,'fit',0.95);

legend off;  % turn off legend from plot method call

set(h_(1),'Color',[1 0 0],...

     'LineStyle','-', 'LineWidth',2,...

     'Marker','none', 'MarkerSize',6);

hold on,plot(x,y,'*')

cf_ =

   General model:
   cf_(x) = a*exp(-p/(x+q))
   Coefficients (with 95% confidence bounds):
   a =          1  (1, 1)
   p =       -42  (-42, -42)
   q =       1.5  (1.5, 1.5)

good =

         sse: 1.7672e-007
   rsquare: 1
         dfe: 11
adjrsquare: 1
      rmse: 1.2675e-004

贊一下

回復此樓

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

zhangxiaodao

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

9樓2011-12-20 14:38:54

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

zhangxiaodao

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1929.4
散金: 31
紅花: 1
帖子: 228
在線: 286小時
蟲號: 871428
注冊: 2009-10-13

送鮮花一朵

引用回帖:

樓: Originally posted by dbb627 at 2011-12-20 14:38:54:
[code]
x=[1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
2
2.1
2.2
2.3
2.4];
y=[10363972.36
5697599.795
3269017.372
1948866.775
1202604.284
765578.5464
501320.0508
336809.4864
231640.3 ...

非常感謝dbb627版主，可能您也知道我問您要源程序的目的了，我手里還有一些數(shù)據(jù)，是跟上面一樣的方程，希望得到a，p，q值，請問能直接套用這個程序嗎，或者需要做哪些調整，謝謝您。

贊一下

回復此樓

10樓2011-12-20 15:41:40

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

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

[考研] 求調劑 +4	Mqqqqqq 2026-03-19	4/200	2026-03-20 14:15 by 星空星月
[考研] 295復試調劑 +7	簡木ChuFront 2026-03-19	7/350	2026-03-20 14:11 by 星空星月
[考研] 274求調劑 +8	S.H1 2026-03-18	8/400	2026-03-20 11:53 by 學員8dgXkO
[考研] 一志愿南昌大學，327分，材料與化工085600 +5	Ncdx123456 2026-03-19	5/250	2026-03-20 11:15 by wangy0907
[考研] 能源材料化學課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	14/700	2026-03-20 09:30 by kkcoco25
[考研] 081700化工學碩調劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 復試調劑 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考研] 266求調劑 +5	陽陽哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽陽哇塞
[考研] 328求調劑，英語六級551，有科研經歷 +4	生物工程調劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調劑
[考研] 一志愿中海洋材料工程專碩330分求調劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 0703化學調劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 0854可跨調劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 301求調劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 326求調劑 +4	諾貝爾化學獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學獎覬?/a>
[考研] 321求調劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 0856求調劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 289求調劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 招收0805（材料）調劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 311求調劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky