版塊導(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í)熱門版塊排行榜

返回列表

本帖產(chǎn)生 1 個(gè) 程序強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

» 猜你喜歡

一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有5人回復(fù)
323求調(diào)劑已經(jīng)有6人回復(fù)
352求調(diào)劑已經(jīng)有3人回復(fù)
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑已經(jīng)有8人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有18人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
280分求調(diào)劑一志愿085802 已經(jīng)有7人回復(fù)
0854電子信息求調(diào)劑已經(jīng)有3人回復(fù)
263求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求matlab求解一個(gè)方程組的問題已經(jīng)有6人回復(fù)
【求助】利用MATLAB求解的外圍輪廓包絡(luò)線為什么不與圖形外圍輪廓重合？怎樣才可以精確已經(jīng)有8人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】matlab求解一個(gè)簡(jiǎn)單的積分（急�。。。� 已經(jīng)有7人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問題已經(jīng)有21人回復(fù)
【求助】Matlab非線性最小二乘擬合活度系數(shù)模型（Willson、NRTL、UNIQUAC）已經(jīng)有17人回復(fù)
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【討論】CFD軟件是求解PDE的，matlab中也可以求解PDE 已經(jīng)有6人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【求助】請(qǐng)教蟲友關(guān)于matlab差分方程的求解和畫圖程序已經(jīng)有4人回復(fù)
【求助】求MAtlab求解程序！�。� 已經(jīng)有20人回復(fù)

1樓 2011-11-27 11:23:40

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

回帖置頂 ( 共有1個(gè) )

dbb627

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

專家經(jīng)驗(yàn): +4
程序強(qiáng)帖: 2
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

★
jjdg(金幣+1, 程序強(qiáng)帖+1): 辛苦了 2011-11-28 01:15:06

CODE:

X=[-7.022 -2.391 0.7963 106.444 10e-4 0.365

-7.022 -2.391 0.7963 106.444 10e-3 0.479

-7.022 -2.391 0.7963 106.444 10e-2 0.592

-7.096 -2.524 1.9082 98.851  10e-4 0.427

-7.096 -2.524 1.9082 98.851  10e-3 0.575

-7.096 -2.524 1.9082 98.851  10e-2 0.711

-6.648 -2.392 1.9277 135.257 10e-4 0.595

-6.648 -2.392 1.9277 135.257 10e-3 0.701

-6.648 -2.392 1.9277 135.257 10e-2 0.891];

%====回歸法==============

y= X(:,6)./(1-X(:,6));

XX=[X(:,1).*X(:,5) X(:,2).*X(:,5) X(:,3).*X(:,5) X(:,4).*X(:,5) X(:,5)];

[B,BINT,R,RINT,STATS] = REGRESS(y,XX)

rcoplot(R,RINT)

%=========最小二乘法=========

F=@(a,x)(a(1)*x(:,1)+a(2)*x(:,2)+a(3)*x(:,3)+a(4)*x(:,4)+a(5)).*x(:,5)./...

    (1+(a(1)*x(:,1)+a(2)*x(:,2)+a(3)*x(:,3)+a(4)*x(:,4)+a(5)).*x(:,5));

[par,RESNORM,RESIDUAL,EXITFLAG]= lsqcurvefit(F,rand(1,5), X(:,1:5), X(:,6))

%=====結(jié)果比較====

YY=XX*B;

Te=YY./(1+YY);

f=F(par,X(:,1:5));

plot(X(:,6),'*-');hold on;plot(Te,'ro-')

hold on;plot(f,'<-g')

legend('原數(shù)據(jù)','回歸法','最小二乘')

'回歸法',
Te =

0.0151
0.1329
0.6052
0.0251
0.2049
0.7204
0.0770
0.4547
0.8929

'最小二乘'
f =

0.0726
0.4390
0.8867
0.1919
0.7037
0.9596
0.5254
0.9172
0.9910

贊一下(1人)

回復(fù)此樓

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-11-27 14:22:50

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

普通回帖

dbb627

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

專家經(jīng)驗(yàn): +4
程序強(qiáng)帖: 2
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

★ ★
aimforao(金幣+30): 不錯(cuò)，這么短的時(shí)間內(nèi)做出來，以后要向你學(xué)習(xí)了。 2011-11-27 13:50:12
jjdg(金幣+2): 辛苦了 2011-11-27 14:12:34

引用回帖:

1樓: Originally posted by aimforao at 2011-11-27 11:23:40:
ie=(a*x1+b*x2+c*x3+d*x4+e)*x5/(1+(a*x1+b*x2+c*x3+d*x4+e)*x5);
x1 x2 x3 x4 x5 ie
-7.022 -2.391 0.7963 106.444 10e-4 0.365
-7.022 -2.391 0.7963 106.444 10e-3 0.479
-7.022 -2.3 ...

基于你的方程形式比較特殊，兩種方法，回歸法和最小二乘法

CODE:

X=[-7.022 -2.391 0.7963 106.444 10e-4 0.365

-7.022 -2.391 0.7963 106.444 10e-3 0.479

-7.022 -2.391 0.7963 106.444 10e-2 0.592

-7.096 -2.524 1.9082 98.851  10e-4 0.427

-7.096 -2.524 1.9082 98.851  10e-3 0.575

-7.096 -2.524 1.9082 98.851  10e-2 0.711

-6.648 -2.392 1.9277 135.257 10e-4 0.595

-6.648 -2.392 1.9277 135.257 10e-3 0.701

-6.648 -2.392 1.9277 135.257 10e-2 0.891];

%====回歸法==============

y= X(:,6)./(1-X(:,6));

XX=[X(:,1).*X(:,5) X(:,2).*X(:,5) X(:,3).*X(:,5) X(:,4).*X(:,5) X(:,5)];

[B,BINT,R,RINT,STATS] = REGRESS(y,XX)

rcoplot(R,RINT)

%=========最小二乘法=========

F=@(a,x)(a(1)*x(:,1)+a(2)*x(:,2)+a(3)*x(:,3)+a(4)*x(:,4)+a(5)).*x(:,5)./...

    (1+(a(1)*x(:,1)+a(2)*x(:,2)+a(3)*x(:,3)+a(4)*x(:,4)+a(5)).*x(:,5));

[par,RESNORM,RESIDUAL,EXITFLAG]= lsqcurvefit(F,B, X(:,1:5), X(:,6))

贊一下(1人)

回復(fù)此樓

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.

2樓2011-11-27 13:09:07

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

dbb627

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

專家經(jīng)驗(yàn): +4
程序強(qiáng)帖: 2
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

★ ★ ★
余澤成(金幣+3): 謝謝參與應(yīng)助！ 2011-12-02 17:24:57

回歸法結(jié)果
B =

20.1351
      0
19.7714
1.3244
      0

BINT =

6.9309 33.3393
      0       0
  -10.3720 49.9148
0.4520 2.1969
      0       0

R =

0.5595
0.7661
-0.0822
0.7194
1.0952
-0.1167
1.3857
1.5106
-0.1649

RINT =

-2.2418 3.3607
-1.9613 3.4934
-0.3545 0.1901
-2.0378 3.4766
-1.4941 3.6846
-0.3734 0.1400
-1.0478 3.8193
-0.8151 3.8363
-0.3877 0.0579

STATS =

0.8439 26.9497 0.0010 1.1444

最小二乘法

Local minimum possible.

lsqcurvefit stopped because the final change in the sum of squares relative to
its initial value is less than the default value of the function tolerance.

par =

  296.9674 65.5619  306.2606 19.8283  -34.1639

RESNORM =

0.3692

RESIDUAL =

-0.2924
-0.0400
0.2947
-0.2351
0.1287
0.2486
-0.0696
0.2162
0.1000

EXITFLAG =

   3

贊一下(1人)

回復(fù)此樓

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.

3樓2011-11-27 13:11:55

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

aimforao

禁蟲 (正式寫手)

本帖內(nèi)容被屏蔽

4樓2011-11-27 13:12:13

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

aimforao

禁蟲 (正式寫手)

本帖內(nèi)容被屏蔽

5樓2011-11-27 14:01:51

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

aimforao

禁蟲 (正式寫手)

本帖內(nèi)容被屏蔽

7樓2011-11-27 15:41:51

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

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

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

★
aimforao(金幣+5): 這個(gè)我用過，效果不是太好。 2011-11-27 22:28:32
jjdg(金幣+1): 感謝參與 2011-11-28 01:15:33

一樓的問題本身就是多解的，1stOpt求解：

CODE:

Variable x1,x2,x3,x4,x5,ie;

Function  ie=(a*x1+b*x2+c*x3+d*x4+e)*x5/(1+(a*x1+b*x2+c*x3+d*x4+e)*x5);

Data;

x1     x2     x3     x4      x5    ie

-7.022 -2.391 0.7963 106.444 10e-4 0.365

-7.022 -2.391 0.7963 106.444 10e-3 0.479

-7.022 -2.391 0.7963 106.444 10e-2 0.592

-7.096 -2.524 1.9082 98.851  10e-4 0.427

-7.096 -2.524 1.9082 98.851  10e-3 0.575

-7.096 -2.524 1.9082 98.851  10e-2 0.711

-6.648 -2.392 1.9277 135.257 10e-4 0.595

-6.648 -2.392 1.9277 135.257 10e-3 0.701

-6.648 -2.392 1.9277 135.257 10e-2 0.891

結(jié)果1：
均方差(RMSE): 0.202527890633088
殘差平方和(SSE): 0.369157918358593
相關(guān)系數(shù)(R): 0.898894702339435
相關(guān)系數(shù)之平方(R^2): 0.808011685893902
決定系數(shù)(DC): -0.78116114903423

參數(shù) 最佳估算
---------- -------------
a 171.747254365496
b 476.920810961258
c 347.506046605747
d 19.8838381108433
e 31.1764490928192

結(jié)果2：
均方差(RMSE): 0.202527890633088
殘差平方和(SSE): 0.369157918358592
相關(guān)系數(shù)(R): 0.898894697992795
相關(guān)系數(shù)之平方(R^2): 0.808011678079559
決定系數(shù)(DC): -0.781161149034225

參數(shù) 最佳估算
---------- -------------
a 473.628802180385
b -461.290822287828
c 253.57821566365
d 19.6210369282456
e 10.4930047968312

結(jié)果3：
均方差(RMSE): 0.202527890633088
殘差平方和(SSE): 0.369157918358593
相關(guān)系數(shù)(R): 0.898894694897679
相關(guān)系數(shù)之平方(R^2): 0.808011672515192
決定系數(shù)(DC): -0.781161149034227

參數(shù) 最佳估算
---------- -------------
a 260.476868413158
b 743.522098490703
c 371.156955344628
d 17.8126496010584
e 1493.31184889042

結(jié)果4：
均方差(RMSE): 0.202527890633089
殘差平方和(SSE): 0.369157918358595
相關(guān)系數(shù)(R): 0.898894706202555
相關(guān)系數(shù)之平方(R^2): 0.808011692838977
決定系數(shù)(DC): -0.781161149034242

參數(shù) 最佳估算
---------- -------------
a 12.9772017982061
b 963.05053780153
c 396.215255571043
d 20.0489254291642
e 22.2696191795072

贊一下(1人)

回復(fù)此樓

8樓2011-11-27 19:51:16

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

aimforao

禁蟲 (正式寫手)

本帖內(nèi)容被屏蔽

9樓2011-11-29 17:07:30

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

dbb627

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

專家經(jīng)驗(yàn): +4
程序強(qiáng)帖: 2
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

引用回帖:

9樓: Originally posted by aimforao at 2011-11-29 17:07:30:
[par,RESNORM,RESIDUAL,EXITFLAG]= lsqcurvefit(F,rand(1,5), X(:,1:5), X(:,6))

這段代碼是什么意思啊，我看不懂，可以告訴一下我嗎？

你不會(huì)看幫助文檔嗎 help lsqcurvefit 就有他的用法介紹，解釋每項(xiàng)的含義太麻煩了吧

贊一下

回復(fù)此樓

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.

10樓2011-11-29 17:30:01

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

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

返回列表

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

[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11	楓橋ZL 2026-03-18	13/650	2026-03-22 20:26 by edmund7
[考研] 求調(diào)劑 +6	十三加油 2026-03-21	6/300	2026-03-22 17:00 by i_cooler
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +6	皮皮皮123456 2026-03-22	6/300	2026-03-22 16:52 by i_cooler
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 求助 +5	夢(mèng)里的無言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-21 00:39 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機(jī)用戶 2026-03-17	9/450	2026-03-20 14:24 by 無懈可擊111
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao