C語言實現(xiàn)運籌學中的馬氏決策算法實例

更新時間：2017年09月16日 09:01:45 作者：reghi

這篇文章主要介紹了C語言實現(xiàn)運籌學中的馬氏決策算法,簡單介紹了馬氏決策的概念,并結(jié)合實例形式分析了C語言實現(xiàn)馬氏決策算法的具體實現(xiàn)技巧,需要的朋友可以參考下

本文實例講述了C語言實現(xiàn)運籌學中的馬氏決策算法。分享給大家供大家參考，具體如下：

一、概述

馬氏決策(Markov decision)是馬爾可夫決策過程(Markov Decision Processes，簡記為MDP)的簡稱，是研究隨機序貫決策問題的一門重要理論。馬氏決策是一類可連續(xù)進行觀察的隨機動態(tài)系統(tǒng)的最優(yōu)化決策，它將(確定性)動態(tài)規(guī)劃與馬爾可夫過程相結(jié)合，是隨機離散事件動態(tài)系統(tǒng)惟一的動態(tài)控制方法。

關(guān)于馬氏決策的具體說明可參考百度百科：https://baike.baidu.com/item/%E9%A9%AC%E6%B0%8F%E5%86%B3%E7%AD%96

二、實現(xiàn)代碼

#include<stdio.h>
#include<cstdlib>
#define N 100 
float p[N][N],s[N][N],a[N],b[N];
int o;
void set_TPM()  //輸入轉(zhuǎn)移概率矩陣(Transition Probability Matrix) 
{  int i,j;
 printf("Please input Number of State:");
 scanf("%d",&o);
 for(i=0;i<o;i++)
 for(j=0;j<o;j++)
 {
  printf("Please input state%d,state%d:",i,j);
  scanf("%f",&p[i][j]);
  rewind(stdin);
 }
}
void set_Initial_Prob() //輸入初始概率狀態(tài)（Initial Probability）
{
 int i;
 for(i=0;i<o;i++)
 {
 printf("Please input state%d Initial Prob:",i);
 scanf("%f",&a[i]);
 rewind(stdin);
 }
}
void run_Markov(int count) //Markov主算法
{
 int i,j,k;
 float c[N];
 for(i=0;i<o;i++) c[i]=a[i];
 for(k=0;k<count;k++)
 {
 for(i=0;i<o;i++)
  for(j=0;j<o;j++)
  {
  s[i][j]=p[i][j]*c[i]; 
  }
 for(i=0;i<o;i++)
 { 
  b[i]=0;
  for(j=0;j<o;j++)
  {
  b[i]=b[i]+s[j][i];  
  }
  c[i]=b[i]; 
 }
 }
 for(i=0;i<o;i++) c[i]=0;
}
void print_Result() //輸出周期結(jié)果
{
 int i,j;
 for(i=0;i<o;i++)
 for(j=0;j<o;j++)
 {
  printf(" %f",s[i][j]);
  if(j==2) printf("/n");
 }
 for(i=0;i<o;i++)
 {
  printf(" %f",b[i]);
 }
 printf("/n");
}
main() //主函數(shù)
{
 int a,count,i,j;
 for(count=0;;)
 {
 printf("Create New Project:/n");
 set_TPM();
 set_Initial_Prob();
 for(;;)
 {
  printf("***********************************/n"); //展示選擇菜單
  printf("1.Times periods from initial./n");
  printf("2.Next Period./n");
  printf("3.Create New Porject./n");
  printf("4.Exit./n/n");
  printf("**********************************/n");
  printf("Please input your choose:/n");
  scanf("%d",&a);
  rewind(stdin);
  if(a==3) break;
  switch(a)
  {
  case 1: 
  printf("Input number of time periods from initial:/n");
  scanf("%d",&count);
  rewind(stdin);
  run_Markov(count);
  print_Result();
  break;
  case 2: 
  run_Markov(count++);
  print_Result();
  break;
  case 4: exit(1);
  default: printf("Error choose!!/n");break;
  }
 }
 }
 system("pause");
}

希望本文所述對大家C語言程序設計有所幫助。

您可能感興趣的文章: