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

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

分形几何

开发平台：

Visual C++

fractus4_03.cpp：源码内容
  1. //---------------------------------------------------------------------------
  3. #include <vcl.h>
  4. #pragma hdrstop
  6. #include "fractus4_03.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=1;
  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.   b1=(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.   putpixel(-sx,-sy,colors[col]);
  43.   putpixel(-sx,sy,colors[col]);
  44.  }
  45. //---------------------------------------------------------------------------
  47. void julia(void)
  48.  {int k;double s;
  49.   static double x[KL+1],y[KL+1];
  50.   x[0]=x0;y[0]=y0;
  51.   for(k=1;k<=KL;k++)
  52.    {x[k]=x[k-1]*x[k-1]-y[k-1]*y[k-1]+xc;
  53.     y[k]=2*x[k-1]*y[k-1]+yc;
  54.     if((x[k]*x[k]+y[k]*y[k])>=4.0)
  55.       return;
  56.    }
  57.   for(k=1;k<=13;k++)
  58.    {if(fabs(x[KL]-x[KL-k])<BOX
  59.      &&fabs(y[KL]-y[KL-k])<BOX)
  60.      {plot(k);return;}
  61.    }
  62.   plot(15);
  63.  }
  64. //---------------------------------------------------------------------------
  66. void __fastcall TForm1::FormCreate(TObject *Sender)
  67. {
  68.   Image1->Height=0;Image1->Top=0;
  69.   Image1->Width=Screen->Width;Image1->Height=Screen->Height;
  70.   rect=Rect(0,0,Image1->Width,Image1->Height);//清屏
  71.   Image1->Canvas->CopyMode=cmBlackness;
  72.   Image1->Canvas->CopyRect(rect,Image1->Canvas,rect);
  73. }
  74. //---------------------------------------------------------------------------
  76. void __fastcall TForm1::N1Click(TObject *Sender)
  77. {
  78.   Image1->Canvas->CopyRect(rect,Image1->Canvas,rect);
  79.   set0(320,200);
  80.   coord();
  81.   Image1->Canvas->Pen->Color=colors[7];
  82.   Image1->Canvas->Pen->Style=psSolid;
  83.   Image1->Canvas->Brush->Style=bsSolid;
  84.   Image1->Canvas->Brush->Color=clBlack;
  85.   Image1->Canvas->Rectangle(-dx,-dy,dx,dy);
  86.   Image1->Canvas->Pen->Style=psDot;
  87.   for(sx=-dx+dx/4;sx<dx;sx+=dx/4)
  88.    {Form1->Image1->Canvas->MoveTo(sx,-dy);
  89.     Form1->Image1->Canvas->LineTo(sx,dy);
  90.    }
  91.   for(sy=-dy+dy/4;sy<dy;sy+=dy/4)
  92.    {Form1->Image1->Canvas->MoveTo(-dx,sy);
  93.     Form1->Image1->Canvas->LineTo(dx,sy);
  94.    }
  95.   for(sx=0;sx<=dx;sx+=STEP)
  96.    {x0=a1*sx+b1;
  97.     for(sy=0;sy<=dy;sy+=STEP)
  98.      {y0=a2*sy+b2;
  99.       julia();
  100.      }
  101.    }
  102. }
  103. //---------------------------------------------------------------------------
  105. void __fastcall TForm1::N2Click(TObject *Sender)
  106. {
  107.   if(SavePictureDialog1->Execute())
  108.     Image1->Picture->SaveToFile(SavePictureDialog1->FileName);
  109. }
  110. //---------------------------------------------------------------------------
  112. void __fastcall TForm1::N3Click(TObject *Sender)
  113. {
  114.   Close();        
  115. }
  116. //---------------------------------------------------------------------------
  118. void __fastcall TForm1::N4_1Click(TObject *Sender)
  119. {
  120.   xc=-0.6;yc=0;
  121. }
  122. //---------------------------------------------------------------------------
  124. void __fastcall TForm1::N4_2Click(TObject *Sender)
  125. {
  126.   xc=-1.1;yc=0;
  127. }
  128. //---------------------------------------------------------------------------
  130. void __fastcall TForm1::N4_3Click(TObject *Sender)
  131. {
  132.   xc=-1.35;yc=0;
  133. }
  134. //---------------------------------------------------------------------------
  136. void __fastcall TForm1::N4_4Click(TObject *Sender)
  137. {
  138.   xc=-1.76;yc=0;
  139. }
  140. //---------------------------------------------------------------------------
  142. void __fastcall TForm1::N4_5Click(TObject *Sender)
  143. {
  144.   xc=-1.65;yc=0;x0=0.2;y0=0;
  145. }
  146. //---------------------------------------------------------------------------
  148. void __fastcall TForm1::N4_6Click(TObject *Sender)
  149. {
  150.   xc=-1.76;yc=0;x0=0.1;y0=0.04;
  151. }
  152. //---------------------------------------------------------------------------