利用蒙特卡洛算法求路径的综合满意度的程序源代码.docx

??

?

??

利用蒙特卡洛算法求路径的综合满意度的程序源代码

?

??

?

?

?

?

?

?

?

???

?

?

?

?

?

#include“Stdio.h”

#include“Conio.h”

#include＜stdlib.h＞

#include＜math.h＞

#defineBETA0.90/*可能性置信水平*/

#defineGOAL11/*目标节点编号*/

#defineN_TEST5000/*testnumberforMonteCarlo*/

#defineVEX_NUM12

#defineROADSECTION_NUM16

#defineMAX_VEX50

#defineSR_MAX5/*maxscenarionumberofanUncertPnt*/

#defineSTKSIZ（8*sizeof（void*）-2）

#defineCUTOFF8/*testingshowsthatthisisgoodvalue*/

struct_Scenario/*Scenario（SR），ascenarioofanUncertPnt*/

{

floatPr；/*SR’spossibilitytohappen*/

floatIs1；/*SR’simpact1*/

floatIs2；/*SR’simpact2*/

floatIs3；/*SR’simpact3*/

floatIs4；/*SR’simpact4*/

}；

struct_UncertPnt/*UncertPnt（UP），anUncertPntofaScheme*/

{

intSR_No；/*selectedscenarioNO.ofaUncertPnt*/

intSR_Num；/*scenarionumberofaUncertPnt*/

struct_ScenarioScenario［SR_MAX］；/*scenariosofaUncertPnt*/

floatIs1max；/*SR’smaximpact1*/

floatIs2max；/*SR’smaximpact2*/

floatIs3max；/*SR’smaximpact3*/

floatIs4max；/*SR’smaximpact4*/

floatIs1min；/*SR’sminimpact1*/

floatIs2min；/*SR’sminimpact2*/

floatIs3min；/*SR’sminimpact3*/

floatIs4min；/*SR’sminimpact4*/

floatDmin_max；

}；

struct_UncertPntRoadSection［VEX_NUM］［VEX_NUM］；

floattmpFitness；

structGraph{

intvexnum；

intGarr［MAX_VEX］［MAX_VEX］；

}；

structGraphG；

intvisited［MAX_VEX］；

floatGlobalIs1max=0.0，GlobalIs2max=0.0，GlobalIs3max=0.0，GlobalIs4max=0.0；

floatGlobalIs1min=1.0，GlobalIs2min=1.0，GlobalIs3min=1.0，GlobalIs4min=1.0；

floatDmin_max=0.0；

structPath{

floatfitness；

intlength；

intpath［MAX_VEX］；

intSR_No［MAX_VEX］；

}；

structPathpathes；

intpathNum=0；

structCurrentPath{

intpath；

intlength；

}currentPath；

/*Declarationofproceduresusedbythisalgorithm*/

intcomp（void*a，void*b）

{

if（（*（float*）a）＞（*（floa