基於C語言實現計算生辰八字五行的範例詳解

Canvas cs ;
string sBarDes[10];   //按鈕
int nBarId[10];
float src[4];  //ClearDraw (cls) clear screen
int ds[3] ;
 int i,j,t,k;    //t = times
string fname;  //filename 加底圖
int sx, sy, x0, y0 ;
int alln,alln1;
int ns,leap;
double nleap,nleapd;  //農曆閏月,0無閏,閏月大小
int n,n0,n1,n2 ;
string s,s0,ss,s5,s6 ;
string s1,s2,s3,s4;
 string s11, s12, s13, s14 ;

int wd ;
string wds,wdss;  //weekday
 double dy0,dy1,dm1,dd1;
 int dy, dm, dd ;
 string dy2, dm2, dd2 ;   //日期輸入dy2
double hh1,mm1,ss1;
  string hh2,mm2,ss2;
int gzy,gzm,gzd,gzh;  //干支 yy mm dd hh
string gzys,gzms,gzds,gzhs;
int gzn1;
int yz,mz; //年柱  月柱
string gzs, bzs ;  //八字
string whss, g12ss,st28ss; //五行 十二建 二十八宿
string tgs, dzs, whs, whs1, whs2 ;  //計算五行
int dzd0,dzd1;  //冬至日

 string jqnames ;      //節氣中文名稱
 string t1, t2, t3 ;       //計算時間串 hh;mm;ss
double jd;            //計算節氣
double juD,tht,yrD;
double sD;    //shouD
double vs,dalt;     //solar

int jqdn ;        //節氣計算
int jqd1,jqd2 ;
 string jqs1,jqs2;
 string jqss1,jqss2;
 string jnum,jname1, jname2 ;
 
main(){
setDisplay(1);
  cs.SetBackground(240,240,240);
  cs.Active();
   sBarDes[0]="輸入日期";
   nBarId[0]=100;
   sBarDes[1]="    ";
   nBarId[1]=101;
   sBarDes[2]="文字顯示";
   nBarId[2]=102;
   sBarDes[3]="圖片顯示";
   nBarId[3]=103;
   sBarDes[4]="退出程式";
   nBarId[4]=104;
   sBarDes[5]="...";
   nBarId[5]=105;
   setToolBarHeight(6);
   setButtonTextSize(14);
   setToolBarBackgroundColor(255,220,220,220);
   setButtonColor(255,240,240,240);
   setButtonTextColor(255,0,0,120);
   setToolBar(100,myToolBarProc,sBarDes,nBarId,6);
  setTitle("生辰八字五行計算   ");
//**********
   jqnames="小寒大寒立春雨水驚蟄春分清明穀雨立夏小滿芒種夏至小暑大暑立秋處暑白露秋分寒露霜降立冬小雪大雪冬至";
  wdss="星期日星期一星期二星期三星期四星期五星期六" ;
whss="海中金爐中火大林木路旁土劍鋒金山頭火洞下水城牆土白臘金楊柳木泉中水屋上土霹雷火松柏木長流水沙中金山下火平地木壁上土金箔金佛燈火天河水大驛土釵釧金桑松木大溪水沙中土天上火石榴木大海水";
g12ss="建除滿平定執破危成收開閉";
st28ss="角亢氐房心尾箕斗牛女虛危室壁奎婁胃昴畢觜參井鬼柳星張翼軫";

      setDisplay (0);
      print " 請選擇：輸入日期 " ; 
      bzwh ( ) ;    //八字五行計算
 while (){} 
 }//main ()

bzwh ( ){        //八字干支表gzs五行表 whs
int n1,n2 ;
    tgs="甲乙丙丁戊己庚辛壬癸" ;
   dzs="子醜寅卯辰巳午未申酉戌亥" ;
   whs1="木木火火土土金金水水" ;
   whs2="水土木木土火火土金金土水" ;
           setDisplay (0);
      for (i=0; i<60; i++){
           n1=i-i/10*10 ;   n2=i-i/12*12 ;
           s1=subString (tgs, n1, 1) ;
           s2=subString (dzs, n2, 1) ;
           s3=subString (whs1, n1, 1) ;
           s4=subString (whs2, n2, 1) ;
               gzs=gzs+s1+s2 ;           //建八字表
               whs=whs+s3+s4 ;    }   //建五行表
}//bzwh ()

test (){  //  生辰八字計算，測試
       cs.ClearDraw (0,src);     //清屏
       clearOutput () ;
          cs.SetTextSize (30);
          cs.SetColor (255,0,0,150) ;
    fname="/storage/emulated/0/tutu.jpg";
    cs.DrawBitmap(fname,0,2);  //加底圖
           sx=60 ;   sy=100 ;
    ss="InputDate = "+dy2+"-"+dm2+"-"+dd2;
       cs.DrawText (ss, sx, sy-60);
       print ss ;
      
    ss= "是日為："+wds ;
        print ss ;
        cs.DrawText (ss, sx, sy-30);
    print "   " ;
    print "all days = ",alln,"    ",alln1;
 
   // 計算節氣日，年柱以立春日起算
   // 計算節氣日，月柱以月節氣日起算
//**** calculate solarterm    //show 節氣 
          dy2=intToString(dy);    //show 節氣
          dm2=intToString(dm);
          dy1=stringToDouble(dy2);
          dm1=stringToDouble(dm2);
           dy0=dy1;
          n=dm*2-1;
          caljq();      //節氣計算
          jqd1=(int)(dd1);
          jqs1=subString(jname2,0,2);
          cs.DrawText (jname2,sx,sy+20);
                jqss1=jname2;
          n=dm*2;
          caljq ();      //節氣計算
          jqd2=(int)(dd1);
          jqs2=subString(jname2,0,2);
          cs.DrawText(jname2,sx,sy+50);         
               jqss2=jname2;
          print jqss1 ;
          print jqss2 ;
  
//calculate GZ,gzy,gzm,gzd,alln
      ss= "節氣日 :  " +intToString(jqd1) ;
         print ss ;
         cs.DrawText(ss, sx, sy+90);
         print "  " ;
         print "  " ;
      ss= "八字五行計算如下：" ;
          print ss ;
          print " " ;
          cs.DrawText(ss, sx, sy+150);
  //計算天數，dy+2096，算公元前避免出現負數
    alln1=alln1+dd-1;       //all days
    n=(dy+2096);
    if(dm<2)n=n-1;
    if(dm==2&&dd<jqd1)n=n-1;
    gzy=n-(int)(n/60*60);
        s1=subString(gzs,gzy*2,2);   //年柱
    if(dd<jqd1){mz=1;}else{mz=0;}
    gzn1=(dy+2096)*12+dm-mz;
    gzm=gzn1-(int)(gzn1/60*60);
        s2=subString (gzs,gzm*2,2);  //月柱
    n1=alln1+9;
    gzd=n1-(int)(n1/60*60);
        s3=subString(gzs,gzd*2,2);    //日柱
    gzh=(gzd-(int)(gzd/5*5))*12;
         s4=subString(gzs,gzh*2,2);    //時柱
     dd2=intToString(dd);  //this day $
     if(len(dd2)==1)dd2="0"+dd2;    //format
     s="日期："+dy2+" 年 "+dm2+" 月 "+dd2+ " 日";
         //cs.SetTextSize (27);   
         cs.DrawText(s,sx, sy+200);
         print s ;
     bzs="八字："+s1+" "+s2+" "+s3+" "+s4;   
         cs.DrawText(bzs,sx , sy+230);
         print bzs ;
  string wh ;
        s11=subString (whs,gzy*2,2);
        s12=subString (whs,gzm*2,2);
        s13=subString (whs,gzd*2,2);
        s14=subString (whs,gzh*2,2);
        wh=s11+" "+s12+" "+s13+" "+s14 ;
         wh="五行："+wh ;       
        cs.DrawText(wh, sx, sy+260);    //output
        print wh ;
        print "  " ;

 //八字五行十二建二十八宿
        s=" " ; ss=" " ;    
    s11=subString (whss,(gzy/2)*3,3);
    s12=subString (whss,(gzm/2)*3,3);
    s13=subString (whss,(gzd/2)*3,3);
    s14=subString (whss,(gzh/2)*3,3);
    s="年 : "+s1+s11;       //show 五行
    cs.DrawText (s,sx, sy+320);
         print "年： "+s1+s11 ;
    s="月 : "+s2+s12;
    cs.DrawText (s, sx, sy+350);
         print "月： "+s2+s12 ;

        n=alln1+3;             //print 二十八宿
         n1=n-n/28*28;
         s5=subString (st28ss,n1,1);
     cs.DrawText (ss, sx+190, sy+380);
         n=alln1-(dy-1900)*12-dm+10;    // 12建日
         if (dd>=jqd1)n=n-1;
         n1=n-n/12*12;
         s6=s5+subString(g12ss,n1,1);
      cs.DrawText (s6, sx+200, sy+380);
    s="日 : "+s3+s13 ;
    cs.DrawText (s, sx, sy+380);
         print "日： "+s3+s13+s6 ;
    s="時 : "+s4+s14;
    cs.DrawText (s, sx, sy+410);
          print "時： "+s4+s14 ;

 //** 八字起時 ，十二時辰對照表 *************
         x0=60;
         y0=600;
         cs.SetTextSize(24);
         s="起時 : ";
         cs.DrawText(s,x0,y0-30);
         s="23 - 1              11 - 13";
         cs.DrawText(s,x0,y0-5);
         s=" 1 -  3              13 - 15";
         cs.DrawText(s,x0,y0+20);
         s=" 3 -  5              15 - 17";
         cs.DrawText(s,x0,y0+45);
         s=" 5 -  7              17 - 19";
         cs.DrawText(s,x0,y0+70);
         s=" 7 -  9              19 - 21";
         cs.DrawText(s,x0,y0+95);
         s=" 9 - 11             21 - 23";
         cs.DrawText(s,x0,y0+120);
    for(i=0;i<6;i++){
          s5=subString(gzs,(gzh+i)*2,2);
          s6=subString(gzs,(gzh+6+i)*2,2);
          s=s5+"                  "+s6;
         cs.DrawText(s,x0+80,y0-5+i*25);     }
 string wh1, wh2, h1, h2 ;
    wh1=" ";  wh2=" ";  h1=" ";  h2=" ";
         n=gzh ;   
    for(i=n ; i<n+6; i++){
         whs1=subString (whs, i*2, 2) ;
         whs2=subString (whs, (i+6)*2, 2) ;   
         wh1=wh1+whs1+" " ;
         wh2=wh2+whs2+" " ; 
         s5=subString(gzs, i*2,2);
         s6=subString(gzs, (i+6)*2,2);
         h1=h1+s5+" " ;
         h2=h2+s6+" " ;          }      
          wh1="23 - 11  "+wh1 ;
          wh2="11 - 23  "+wh2 ;
       cs.SetTextSize (24);   
        print "  " ;
        print "起時：" ;
       cs.DrawText (wh1,50,760) ;
       cs.DrawText (wh2,50,790) ;
         print "  " ;
         print "23 - 11  ", h1 ;
         print wh1 ;
         print "  " ;
         print "11 - 23  ", h2 ;
         print wh2 ;
 
      cs.Update () ;
}//test ( )

 
inputDate ()｛ //輸入日期
string sDat[101];
     getDate(ds);    //getDate = today
     pickDate("輸入日期:",ds);       //input date
        dy=ds[0];
        dm=ds[1];
        dd=ds[2];
      dy2=intToString(dy);  //this year
      dm2=intToString(dm);  //this month
      if (dm<10) dm2="0"+dm2 ; 
      dd2=intToString(dd);  //this day 
      if (dd<10) dd2="0"+dd2 ; 
  //print "InputDate = ",dy2+"-",dm2+"-"+dd2;
      getWeekday () ;
}//inputDate()

getWeekday () {

//傳人dy dm dd , 輸出星期wds , 天數alln1
//計算月首日星期幾 alln=... +1+1 ;
//計算當日或輸入 是日 alln=... +1+dd ;
int ty, tm ;
      ty=dy;
      tm=dm;
      if (dm<3){
          ty=dy-1;
          tm=dm+12;   }
  //first day=1 this month
  //iWeek=(d+2*m+3*(m+1)/5+y+y/4-y/100+y/400) 
  //mod 7   基姆拉爾森公式
      alln=2*tm+3*(tm+1)/5+ty+(int)(ty/4);
      alln=alln-(int)(ty/100)+(int)(ty/400)+1+dd;
      wd=alln-alln/7*7;      //mod (alln,7)
      wds=subString(wdss,wd*3,3);
 string yds;     //計算天數
 string ydds;
 double tms;
 yds="000031059090120151181212243273304334";

ydds=subString(yds,(dm-1)*3,3);
  tms=stringToDouble(ydds);
  alln1=(dy-1900)*365+((dy-1900-1)/4)+(int)(tms)+1;
  if ((dy-dy/4*4)==0&&dm>2)alln1=alln1+1;
      // return weekday = wds , alln1
     test () ;
}//getWeekday ()

myToolBarProc(int nBtn,int nContext){
      if(nBtn==100){//input date
            setDisplay (0);
            inputDate ( ) ;
                   }  
      if(nBtn==101){//test  input
             inputDate ( ) ;
                    }
      if(nBtn==102){//關閉圖板文字列印輸出
            setDisplay (0);
               }
      if(nBtn==103){//顯示圖板
           setDisplay (1);
               }
      if(nBtn==104){//退出程式
         clearOutput();
         cs.ClearDraw (0,src);
         setDisplay (0);
          exit (0);
               }
     if(nBtn==105){//@testdraw
          setDisplay (1);
                }
}//myToolBar ()

caljq(){//calculate solar ** return jname2

  //傳入dy0, n節氣序數，輸出節氣時間
  juD=dy0*(365.2423112-0.000000000000064
     *(dy0-100)*(dy0-100)- 0.00000003047
     *(dy0-100))+15.218427*n+1721050.71301;
tht=0.0003*dy0-0.372781384-0.2617913325*n;
yrD=(1.945*sin(tht)-0.01206*sin(2*tht))
        *(1.048994-0.00002583 *dy0);
sD=-0.0018*sin(2.313908653*dy0
        -0.439822951-3.0443*n);
  //'vs = juD          '** 平氣
vs = (juD + yrD + sD);  //   '** 定氣
dalt=-15+(juD-2382148)*(juD-2382148)/41048480;
//print dalt;
dalt = dalt/86400;
//print dalt;
 jd = vs - 0.5 - dalt-0.025;     
int gg;
gg=n-(int)(n/2*2);

jdtoGD ();
jqs1=subString (jqnames,(n-1)*2,2);
//print "jqs1 = ",jqs1;
s=jqs1+" " +jname2;
//print jname2;           
}//caljq() ****************************

//calculate jd to GD ** return jname2  *****
jdtoGD(){
double a,b,c,d,e;
double F;
double allss;
F=jd-(int)(jd);
//print "   ";
//print " JD =  ",jd;
//print " 時分秒 日小數 =  ",F;
 a=(int)(jd+0.5);
 b=a+1537;
 c=(int)((b-122.1)/365.25);
 d=(int)(365.25*c);
 e=(int)((b-d)/30.6001);
dd1=b-d-(int)(30.6001*e);
dm1=e-1-(int)((e/14)*12);
dy1=c-4715-(int)((7+dm1)/10);
//print a," ",b,"  ",c,"  ",d,"  ",e;
dy2=intToString (dy);
dm2=intToString (dm);
dd2=doubleToString (dd1);
if (len(dm2)==1)dm2="0"+dm2;
if (dd1<10)dd2="0"+dd2;
dd2=subString (dd2,0,2);
//print dy2+" 年 "+dm2+" 月 "+dd2+" 日 ";
if(n==24)jqdn=(int)(dd1);
if (n==24&&dm1<3)jqdn=jqdn+1;

   //日allss 的小數轉為時分秒
  allss=(int)((jd-a)*86400+43200.5);
//print "allss = ", allss;
  hh1=(int)(allss/3600);
  mm1=(int)((allss-hh1*3600)/60);
  ss1=(int)(allss-hh1*3600-mm1*60);
   if(ss1>=60){
        ss1=ss1-60;
        mm1=mm1+1;}
   if(mm1>=60){
        mm1=mm1-60;
        hh1=hh1+1;}
//print "JD  轉為 GD，計算結果：」;
    hh2=doubleToString(hh1);
    mm2=doubleToString(mm1);
    ss2=doubleToString(ss1);
    if(hh1<10){
        hh2="0"+doubleToString(hh1);}
    if(mm1<10){
        mm2="0"+doubleToString(mm1);}
    if(ss1<10){
        ss2="0"+doubleToString(ss1);}
    hh2=subString (hh2,0,2);
    mm2=subString (mm2,0,2);
    ss2=subString (ss2,0,2);
    jname1=subString(jqnames,(n-1)*2,2);
    //print "節氣 :  ",jname1;
     jname2=jname1+"： "+dm2+"-"+dd2+"  "+hh2+":"+mm2+":"+ss2;
     //return jname2;
}//jdtoGD() **************


//**** END ****