 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1672 | 回復(fù): 3 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
wsbl實(shí)習(xí)版主
|
[交流]
Levenberg-Marquardt 算法的例子代碼 已有3人參與
|
|||
|
#include "stdio.h" #include #include void out(double*p,int n,int m)//輸出矩陣 {int i,j; for (i=0;i if ((m*i+j)%m==0) printf("\n"  ;printf("%16.8f",*(p+m*i+j)); } printf("\n\n" ;} void LM(double *c, double *sigma, int N, double sigma_H2O,double r[3]) { int i,j,maxtime=100; //如果不收斂，允許的最大迭代次數(shù) double s,snext,det,temp0; double x[]={.5,.5},xnext[2],ep[2]; double f[6]; double f_x1[6]; double f_x2[6]; double m[4]; //矩陣 double temp1[4];//逆矩陣 double temp2[2]; double temp3[2]; double alph,beita; j=0; for (i=0;i f_x1=sigma_H2O*log(1+c/x[1]); f_x2=-sigma_H2O*x[0]*c/(x[1]*x[1])/(1+c/x[1]); } s=.0; for (i=0;i temp0=.0; for (i=0;i m[0]=temp0; temp0=.0; for (i=0;i m[1]=temp0; m[2]=temp0; temp0=.0; for (i=0;i m[3]=temp0; //out(m,2,2);//到此步正確。 alph=(fabs(*m)+fabs(*(m+3)))/2; //printf("%8.5e\n",alph); beita=2.0; m[0]=m[0]+alph; m[3]=m[3]+alph; //out(m,2,2);//到此步正確。 det=1.0/(m[0]*m[3]-m[1]*m[2]); //printf("迭%15.10f\n",det); temp1[0]=det*m[3]; temp1[1]=-det*m[1]; temp1[2]=-det*m[1]; temp1[3]=det*m[0]; //out(temp1,2,2);//到此步正確。 //temp2=mult(temp,f,2,N,1); temp0=.0; for (i=0;i temp2[0]=temp0; temp0=.0; for (i=0;i temp2[1]=temp0; //out(temp2,2,1);//到此步正確。 //temp3=mult(temp1,temp2,2,2,1); temp0=.0; for (i=0;i<2;i++) { temp0=temp0+temp1*temp2; } temp3[0]=temp0; temp0=.0; for (i=0;i<2;i++) { temp0=temp0+temp1[2+i]*temp2; } temp3[1]=temp0; //out(temp3,2,1);//到此步正確。 xnext[0]=x[0]-temp3[0]; xnext[1]=x[1]-temp3[1]; //out(xnext,2,1);//到此步正確。 snext=0; for (i=0;i snext=snext+f*f; } //printf("%8.5e",snext);//到此步正確。 ep[0]=xnext[0]-x[0]; ep[1]=xnext[1]-x[1]; //out(ep,2,1);//到此步正確。 while (((fabs(ep[0])+fabs(ep[1]))/2)>1.0e-6 && j if (snext>=s) alph=alph*beita; else alph=alph/beita; //printf("%8.3f\n",(float)j); x[0]=xnext[0]; x[1]=xnext[1]; s=snext; for (i=0;i f_x1=sigma_H2O*log(1+c/x[1]); f_x2=-sigma_H2O*x[0]*c/(x[1]*x[1])/(1+c/x[1]); } s=.0; for (i=0;i temp0=.0; for (i=0;i m[0]=temp0; temp0=.0; for (i=0;i m[1]=temp0; m[2]=temp0; temp0=.0; for (i=0;i m[3]=temp0; m[0]=m[0]+alph; m[3]=m[3]+alph; //out(m,2,2);//到此步正確。 det=1.0/(m[0]*m[3]-m[1]*m[2]); //printf("迭%15.10f\n",det); temp1[0]=det*m[3]; temp1[1]=-det*m[1]; temp1[2]=-det*m[1]; temp1[3]=det*m[0]; //out(temp1,2,2);//到此步正確。 //temp2=mult(temp,f,2,N,1); temp0=.0; for (i=0;i temp2[0]=temp0; temp0=.0; for (i=0;i temp2[1]=temp0; //out(temp2,2,1);//到此步正確。 //temp3=mult(temp1,temp2,2,2,1); temp0=.0; for (i=0;i<2;i++) { temp0=temp0+temp1*temp2; } temp3[0]=temp0; temp0=.0; for (i=0;i<2;i++) { temp0=temp0+temp1[2+i]*temp2; } temp3[1]=temp0; //out(temp3,2,1);//到此步正確。 xnext[0]=x[0]-temp3[0]; xnext[1]=x[1]-temp3[1]; //out(xnext,2,1);//到此步正確。 snext=0; for (i=0;i snext=snext+f*f; } ep[0]=xnext[0]-x[0]; ep[1]=xnext[1]-x[1]; } r[0]=xnext[0]; r[1]=xnext[1]; r[2]=(double)j; } void main() { //double p[] = {0.02,0.04,0.08,0.12,0.16,0.2,0.24,0.28};//mol/L //double q[] = {0.06810,0.06442,0.0601,0.05678,0.05467,0.05237,0.05108,0.05025};//N/m double p[] = {0.02,0.04,0.08,0.12,0.16,0.2}; double q[] = {0.06810,0.0644,0.0601,0.0568,0.0547,0.0524}; double r[3]; LM(p,q,6,0.07197,r); printf("迭代次數(shù)：%5.0f\n",r[2]); if (r[2]<100) { printf("\n希什科夫模型sigma_H2O-sigma_H2O*a*ln(1+x/b)中\(zhòng)n\n參數(shù)a的最優(yōu)值為:%9.9f\n",r[0]); printf("\n參數(shù)b的最優(yōu)值為:%9.9f\n\n",r[1]); } else printf("\n用Levenberg法失敗\n\n" ;for(int i=0;i<1;i++)//延時(shí)幾秒。 { printf("按回車鍵結(jié)束程序:" ;getchar(); } } |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 軟件學(xué)習(xí) |
» 猜你喜歡
285求調(diào)劑
已經(jīng)有5人回復(fù)
-
309求調(diào)劑
已經(jīng)有19人回復(fù)
-
08工科275分求調(diào)劑
已經(jīng)有12人回復(fù)
-
esi高被引論文是不是能對中標(biāo)有所加分和幫助呢
已經(jīng)有3人回復(fù)
-
求調(diào)劑推薦 材料 304
已經(jīng)有19人回復(fù)
-
265求調(diào)劑
已經(jīng)有11人回復(fù)
-
292分，材料與化工，申請調(diào)劑
已經(jīng)有9人回復(fù)
-
調(diào)劑
已經(jīng)有5人回復(fù)
-
英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑
已經(jīng)有12人回復(fù)
-
295材料工程專碩求調(diào)劑
已經(jīng)有13人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- Visual C++ CAD應(yīng)用程序開發(fā)技術(shù)源代碼（基于OPenGL技術(shù)開發(fā)的）
已經(jīng)有76人回復(fù)
- Levenberg-Marquardt算法能否用于計(jì)算1范數(shù)定義下的最小值問題？
已經(jīng)有5人回復(fù)
- 求解一元三次方程
已經(jīng)有8人回復(fù)
- 求助，自己編代碼實(shí)現(xiàn)別人的算法，可以申請軟件著作權(quán)嗎
已經(jīng)有4人回復(fù)
- 求助神經(jīng)網(wǎng)絡(luò) ---Levenberg-Marquardt算法
已經(jīng)有6人回復(fù)
- 基于L-M算法的BP神經(jīng)網(wǎng)絡(luò)
已經(jīng)有3人回復(fù)
- 【求助】Levenberg-Marquardt法和LBFGS擬牛頓法，求C++版源碼及使用方法
已經(jīng)有3人回復(fù)
- 【求助】MATLAB中BP神經(jīng)網(wǎng)絡(luò)的訓(xùn)練算法具體是怎么樣的？
已經(jīng)有5人回復(fù)
lingshuning
實(shí)習(xí)版主
- 應(yīng)助: 89 (初中生)
- 金幣: 1320.9
- 散金: 3426
- 紅花: 63
- 帖子: 3433
- 在線: 287.6小時(shí)
- 蟲號: 1515504
- 注冊: 2011-11-29
- 性別: GG
- 專業(yè): 工業(yè)工程與管理
4樓2011-12-03 19:18:50
lrs2008
管理員
觀天
- 應(yīng)助: 7 (幼兒園)
- 金幣: 3457.1
- 散金: 5505
- 紅花: 13
- 帖子: 424
- 在線: 251.9小時(shí)
- 蟲號: 746075
- 注冊: 2009-04-11
- 性別: GG
- 專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)
3樓2011-10-16 09:21:30
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑推薦 材料 304 +19 | 荷包蛋hyj 2026-03-26 | 19/950 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +10 | yelck 2026-04-01 | 11/550 |
|
|
|
[考研] 309求調(diào)劑 +8 | gajsj 2026-03-25 | 10/500 |
|
|
|
[考研] 266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑 +4 | 哇呼哼呼哼 2026-04-01 | 5/250 |
|
|
|
[考研] 288求調(diào)劑 一志愿哈工大 材料與化工 +26 | 洛神哥哥 2026-03-31 | 26/1300 |
|
|
|
[考研] 一志愿北交材料工程總分358 +5 | cs0106 2026-04-01 | 7/350 |
|
|
|
[考研] 311（085601）求調(diào)劑 +12 | liziyeyeye 2026-03-28 | 13/650 |
|
|
|
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +4 | 1段時(shí)間 2026-03-31 | 5/250 |
|
|
|
[考研] 0856求調(diào)劑 +9 | 楒桉 2026-03-28 | 9/450 |
|
|
|
[考研] 一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑 +9 | s1145 2026-03-31 | 9/450 |
|
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +11 | dongfang59 2026-03-27 | 11/550 |
|
|
|
[考研] 0856 335分 | 封裝or激光加工 老師看看我！！ +9 | cccchenso 2026-03-29 | 9/450 |
|
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 13/650 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
|
[考研] 085601一志愿西北工業(yè)大學(xué)初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|