开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 分形几何 > 查看源码
fractus4_03_b.cpp
资源名称：fenxin.rar　[点击查看]
上传用户：hzcygd
上传日期：2022-04-30
资源大小：109k
文件大小：5k
源码类别：

分形几何

开发平台：

Visual C++

fractus4_03_b.cpp：源码内容
  1. //---------------------------------------------------------------------------
  3. #include <vcl.h>
  4. #pragma hdrstop
  6. #include "fractus4_03_b.h"
  7. //---------------------------------------------------------------------------
  8. #pragma package(smart_init)
  9. #pragma resource "*.dfm"
  10. TForm1 *Form1;
  11. #include "glib.h"
  12. const int KL=100;
  13. const double BOX=0.01;
  14. const int STEP=4;
  15. int sx,sy,dx=192,dy=192;
  16. double x_min=-2.0,x_max=2.0,y_min=-2.0,y_max=2.0;
  17. double xc,yc,x0,y0;
  18. double a1,b1,a2,b2;
  19. double delta=0.0001;TRect rect;
  20. int colors[16]={0x00000000,0x00FF0000,0x0000FF00,0x00FF8C00,
  21.                 0x000000FF,0x00FF00FF,0x0000FFFF,0x00A0A0A0,
  22.                 0x009BFFFF,0x009BFFFF,0x00A0FF00,0x00FFFF00,
  23.                 0x006E6EFF,0x00FF9BFF,0x009BFFFF,0x00FFFFFF};
  24. //---------------------------------------------------------------------------
  25. __fastcall TForm1::TForm1(TComponent* Owner)
  26.         : TForm(Owner)
  27. {
  28. }
  29. //---------------------------------------------------------------------------
  30. //Zc为复数的朱利亚集合
  31. void __fastcall coord(void)
  32.  {a1=(x_max-x_min)/(2*dx);
  33.   b1=(x_max+x_min)/2;
  34.   a2=(y_max-y_min)/(2*dy);
  35.   b2=(y_max+y_min)/2;
  36.  }
  37. //---------------------------------------------------------------------------
  39. void __fastcall plot(int col)
  40.  {putpixel(sx,-sy,colors[col]);
  41.   putpixel(-sx,sy,colors[col]);
  42.  }
  43. //---------------------------------------------------------------------------
  45. void julia(void)
  46.  {int k;double s;
  47.   static double x[KL+1],y[KL+1];
  48.   x[0]=x0;y[0]=y0;
  49.   for(k=1;k<=KL;k++)
  50.    {x[k]=x[k-1]*x[k-1]-y[k-1]*y[k-1]+xc;
  51.     y[k]=2*x[k-1]*y[k-1]+yc;
  52.     if((x[k]*x[k]+y[k]*y[k])>=4.0)
  53.       return;
  54.    }
  55.   for(k=1;k<=13;k++)
  56.    {if(fabs(x[KL]-x[KL-k])<BOX
  57.      &&fabs(y[KL]-y[KL-k])<BOX)
  58.      {plot(k);return;}
  59.    }
  60.   plot(15);
  61.  }
  62. //---------------------------------------------------------------------------
  64. void __fastcall TForm1::FormCreate(TObject *Sender)
  65. {
  66.   Image1->Height=0;Image1->Top=0;
  67.   Image1->Width=Screen->Width;Image1->Height=Screen->Height;
  68.   rect=Rect(0,0,Image1->Width,Image1->Height);//清屏
  69.   Image1->Canvas->CopyMode=cmBlackness;
  70.   Image1->Canvas->CopyRect(rect,Image1->Canvas,rect);
  71. }
  72. //---------------------------------------------------------------------------
  74. void __fastcall TForm1::N1Click(TObject *Sender)
  75. {
  76.   Image1->Canvas->CopyRect(rect,Image1->Canvas,rect);
  77.   set0(320,200);
  78.   coord();
  79.   Image1->Canvas->Pen->Color=colors[7];
  80.   Image1->Canvas->Pen->Style=psSolid;
  81.   Image1->Canvas->Brush->Style=bsSolid;
  82.   Image1->Canvas->Brush->Color=clBlack;
  83.   Image1->Canvas->Rectangle(-dx,-dy,dx,dy);
  84.   Image1->Canvas->Pen->Style=psDot;
  85.   for(sx=-dx+dx/4;sx<dx;sx+=dx/4)
  86.    {Form1->Image1->Canvas->MoveTo(sx,-dy);
  87.     Form1->Image1->Canvas->LineTo(sx,dy);
  88.    }
  89.   for(sy=-dy+dy/4;sy<dy;sy+=dy/4)
  90.    {Form1->Image1->Canvas->MoveTo(-dx,sy);
  91.     Form1->Image1->Canvas->LineTo(dx,sy);
  92.    }
  93.   for(sx=0;sx<=dx;sx+=STEP)
  94.    {x0=a1*sx+b1;
  95.     for(sy=-dy;sy<=dy;sy+=STEP)
  96.      {y0=a2*sy+b2;
  97.       julia();
  98.      }
  99.    }
  100. }
  101. //---------------------------------------------------------------------------
  103. void __fastcall TForm1::N2Click(TObject *Sender)
  104. {
  105.   if(SavePictureDialog1->Execute())
  106.     Image1->Picture->SaveToFile(SavePictureDialog1->FileName);
  107. }
  108. //---------------------------------------------------------------------------
  110. void __fastcall TForm1::N3Click(TObject *Sender)
  111. {
  112.   Close();        
  113. }
  114. //---------------------------------------------------------------------------
  116. void __fastcall TForm1::N4_1Click(TObject *Sender)
  117. {
  118.   xc=0;yc=0.6;
  119. }
  120. //---------------------------------------------------------------------------
  122. void __fastcall TForm1::N4_2Click(TObject *Sender)
  123. {
  124.   xc=-1.1;yc=0.2;
  125. }
  126. //---------------------------------------------------------------------------
  128. void __fastcall TForm1::N4_3Click(TObject *Sender)
  129. {
  130.   xc=-0.1;yc=0.77;
  131. }
  132. //---------------------------------------------------------------------------
  134. void __fastcall TForm1::N4_4Click(TObject *Sender)
  135. {
  136.   xc=0.3;yc=0.5;
  137. }
  138. //---------------------------------------------------------------------------
  140. void __fastcall TForm1::N4_5Click(TObject *Sender)
  141. {
  142.   xc=-0.52;yc=0.55;
  143. }
  144. //---------------------------------------------------------------------------
  146. void __fastcall TForm1::N4_6Click(TObject *Sender)
  147. {
  148.   xc=-1.15;yc=0.25;
  149. }
  150. //---------------------------------------------------------------------------
  152. void __fastcall TForm1::N4_7Click(TObject *Sender)
  153. {
  154.   xc=0.12;yc=0.61;        
  155. }
  156. //---------------------------------------------------------------------------
  158. void __fastcall TForm1::N4_8Click(TObject *Sender)
  159. {
  160.   xc=-0.37;yc=0.61;        
  161. }
  162. //---------------------------------------------------------------------------