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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 分形几何 > 查看源码
FractalView.cpp
资源名称：Fractal.rar　[点击查看]
上传用户：dxkjfz
上传日期：2021-06-06
资源大小：1930k
文件大小：15k
源码类别：

分形几何

开发平台：

Visual C++

FractalView.cpp：源码内容
  1. // FractalView.cpp : implementation of the CFractalView class
  2. //
  4. #include "stdafx.h"
  5. #include "Fractal.h"
  7. #include "FractalDoc.h"
  8. #include "FractalView.h"
  9. #include <math.h>
  11. #ifdef _DEBUG
  12. #define new DEBUG_NEW
  13. #undef THIS_FILE
  14. static char THIS_FILE[] = __FILE__;
  15. #endif
  17. #define TREE_TIMER 10
  18. #define SIERPINSKI_TIMER 11
  19. #define FLOWER_TIMER 12
  20. #define CROWN_TIMER 13
  22. const double PI = 3.11415926;
  23. const double PI_tree=3.11415926/180;
  25. /////////////////////////////////////////////////////////////////////////////
  26. // CFractalView
  28. IMPLEMENT_DYNCREATE(CFractalView, CView)
  30. BEGIN_MESSAGE_MAP(CFractalView, CView)
  31. //{{AFX_MSG_MAP(CFractalView)
  32. ON_COMMAND(ID_DRAWTREE, OnDrawtree)
  33. ON_WM_TIMER()
  34. ON_COMMAND(ID_DRAWSTOP, OnDrawstop)
  35. ON_COMMAND(ID_DRAWSIERPINSKI, OnDrawsierpinski)
  36. ON_COMMAND(ID_DRAWFLOWER, OnDrawflower)
  37. ON_COMMAND(ID_DRAWCROWN, OnDrawcrown)
  38. ON_COMMAND(ID_DRAWNEW, OnDrawnew)
  39. //}}AFX_MSG_MAP
  40. // Standard printing commands
  41. ON_COMMAND(ID_FILE_PRINT, CView::OnFilePrint)
  42. ON_COMMAND(ID_FILE_PRINT_DIRECT, CView::OnFilePrint)
  43. ON_COMMAND(ID_FILE_PRINT_PREVIEW, CView::OnFilePrintPreview)
  44. END_MESSAGE_MAP()
  46. /////////////////////////////////////////////////////////////////////////////
  47. // CFractalView construction/destruction
  49. CFractalView::CFractalView()
  50. {
  51. // TODO: add construction code here
  52. m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
  53. m_bgColor = RGB(0, 0, 0);
  54. m_fgColor = RGB(0, 255, 0);
  55. m_dir = TRUE;
  57. tree.m_bRunning = FALSE;
  58. tree.m_Interval = 10;
  59. sier.m_bRunning = FALSE;
  60. sier.m_Interval = 10;
  61. flower.m_bRunning = FALSE;
  62. flower.m_Interval = 10;
  64. crown.m_bRunning = FALSE;
  65. crown.m_Interval = 1;
  66. crown.m_C = 9;
  67. crown.x=5;
  68. crown.y=0;  //与x,y无关
  69. crown.m_pColor= RGB(0,255,0);
  70. crown.stop=false;
  71. crown.ondraw=false;
  72. crown.cxmin=-2;
  73. crown.cxmax=(float)0.8;
  74. crown.cymin=(float)-0.3;
  75. crown.cymax=(float)0.3;
  76.     crown.cxspeed=(float)0.2;
  77. crown.cyspeed=0.25;
  78.     crown.cx=crown.cxmin;
  79.     crown.cy=0;
  80. crown.cxdirect=1;
  81. crown.cydirect=-1;
  82. crown.cnt = 0;
  83. crown.count = 0;
  84. }
  86. CFractalView::~CFractalView()
  87. {
  88. }
  90. BOOL CFractalView::PreCreateWindow(CREATESTRUCT& cs)
  91. {
  92. // TODO: Modify the Window class or styles here by modifying
  93. //  the CREATESTRUCT cs
  95. return CView::PreCreateWindow(cs);
  96. }
  98. /////////////////////////////////////////////////////////////////////////////
  99. // CFractalView drawing
  101. void CFractalView::OnDraw(CDC* pDC)
  102. {
  103. CFractalDoc* pDoc = GetDocument();
  104. ASSERT_VALID(pDoc);
  105. // TODO: add draw code for native data here
  106. static BOOL bDrawflag = true;
  107. //SetIcon(m_hIcon, TRUE); // Set big icon
  108. //SetIcon(m_hIcon, FALSE); // Set small icon
  109. if (bDrawflag){
  110. AfxGetMainWnd()->SetWindowText("分形算法演示");
  112. /************************************************************************/
  113. //获取窗口大小
  114. GetClientRect(&m_Rect);
  115. //得到物理画布
  116. m_DC.Attach(::GetDC(*this));
  117. //创建内存画布
  118. m_MemDC.CreateCompatibleDC(&m_DC);
  119. m_Bitmap.CreateCompatibleBitmap(&m_DC, m_Rect.Width(), m_Rect.Height());
  120. m_MemDC.SelectObject(&m_Bitmap);
  121. //设置画笔颜色
  122. m_Pen.CreatePen(PS_SOLID, 1, m_fgColor);
  123. m_MemDC.SelectObject(&m_Pen);
  124. //填充背景颜色
  125. m_MemDC.FillSolidRect(&m_Rect, m_bgColor);
  126. m_DC.BitBlt(0, 0, m_Rect.Width(), m_Rect.Height(), &m_MemDC, 0, 0, SRCCOPY);
  127. bDrawflag = FALSE;
  128. }
  129. /************************************************************************/
  130. /*
  131. if (IsIconic())
  132. {
  133. CPaintDC dc(this); // device context for painting
  135. //SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);
  137. // Center icon in client rectangle
  138. int cxIcon = GetSystemMetrics(SM_CXICON);
  139. int cyIcon = GetSystemMetrics(SM_CYICON);
  140. CRect rect;
  141. GetClientRect(&rect);
  142. int x = (rect.Width() - cxIcon + 1) / 2;
  143. int y = (rect.Height() - cyIcon + 1) / 2;
  145. // Draw the icon
  146. dc.DrawIcon(x, y, m_hIcon);
  147. }
  148. else
  149. {
  150. m_DC.BitBlt(0, 0, m_Rect.Width(), m_Rect.Height(), &m_MemDC, 0, 0, SRCCOPY);
  152. //CFractalView::OnDraw(pDC);
  153. }*/
  154. }
  156. /////////////////////////////////////////////////////////////////////////////
  157. // CFractalView printing
  159. BOOL CFractalView::OnPreparePrinting(CPrintInfo* pInfo)
  160. {
  161. // default preparation
  162. return DoPreparePrinting(pInfo);
  163. }
  165. void CFractalView::OnBeginPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
  166. {
  167. // TODO: add extra initialization before printing
  168. }
  170. void CFractalView::OnEndPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
  171. {
  172. // TODO: add cleanup after printing
  173. }
  175. /////////////////////////////////////////////////////////////////////////////
  176. // CFractalView diagnostics
  178. #ifdef _DEBUG
  179. void CFractalView::AssertValid() const
  180. {
  181. CView::AssertValid();
  182. }
  184. void CFractalView::Dump(CDumpContext& dc) const
  185. {
  186. CView::Dump(dc);
  187. }
  189. CFractalDoc* CFractalView::GetDocument() // non-debug version is inline
  190. {
  191. ASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CFractalDoc)));
  192. return (CFractalDoc*)m_pDocument;
  193. }
  194. #endif //_DEBUG
  196. /////////////////////////////////////////////////////////////////////////////
  197. // CFractalView message handlers
  198. void CFractalView::OnDrawtree() 
  199. {
  200. // TODO: Add your command handler code here
  201. AfxGetMainWnd()->SetWindowText("摇曳的分形树");
  202. if (tree.m_bRunning) {
  203. KillTimer(TREE_TIMER);
  204. tree.m_bRunning = FALSE;
  205. } else {
  206. OnDrawstop();
  207. m_MemDC.FillSolidRect(&m_Rect, m_bgColor);
  208. SetTimer(TREE_TIMER, tree.m_Interval, NULL);
  209. tree.m_bRunning = TRUE;
  210. }
  212. //CDialog::OnLButtonDown(nFlags, point);
  213. }
  214. void CFractalView::OnDrawsierpinski() 
  215. {
  216. // TODO: Add your command handler code here
  217. AfxGetMainWnd()->SetWindowText("摇摆的Sierpinski三角形");
  218. if (sier.m_bRunning) {
  219. KillTimer(SIERPINSKI_TIMER);
  220. sier.m_bRunning = FALSE;
  221. } else {
  222. OnDrawstop();
  223. m_MemDC.FillSolidRect(&m_Rect, m_bgColor);
  224. SetTimer(SIERPINSKI_TIMER, sier.m_Interval, NULL);
  225. sier.m_bRunning = TRUE;
  226. }
  227. }
  228. void CFractalView::OnDrawflower() 
  229. {
  230. // TODO: Add your command handler code here
  231. AfxGetMainWnd()->SetWindowText("旋转万花筒");
  232. if (flower.m_bRunning) {
  233. KillTimer(FLOWER_TIMER);
  234. flower.m_bRunning = FALSE;
  235. } else {
  236. OnDrawstop();
  237. m_MemDC.FillSolidRect(&m_Rect, m_bgColor);
  238. SetTimer(FLOWER_TIMER, flower.m_Interval, NULL);
  239. flower.m_bRunning = TRUE;
  240. }
  241. }
  242. void CFractalView::OnDrawcrown() 
  243. {
  244. // TODO: Add your command handler code here
  245. AfxGetMainWnd()->SetWindowText("王冠（Julia集）");
  246. if (crown.m_bRunning) {
  247. KillTimer(CROWN_TIMER);
  248. crown.m_bRunning = FALSE;
  249. } else {
  250. OnDrawstop();
  251. m_MemDC.FillSolidRect(&m_Rect, m_bgColor);
  252. SetTimer(CROWN_TIMER, crown.m_Interval, NULL);
  253. crown.m_bRunning = TRUE;
  254. }
  255. }
  256. void CFractalView::OnDrawstop() 
  257. {
  258. // TODO: Add your command handler code here
  259. KillTimer(TREE_TIMER);
  260. tree.m_bRunning = FALSE;
  261. KillTimer(SIERPINSKI_TIMER);
  262. sier.m_bRunning = FALSE;
  263. KillTimer(FLOWER_TIMER);
  264. flower.m_bRunning = FALSE;
  265. KillTimer(CROWN_TIMER);
  266. crown.m_bRunning = FALSE;
  267. }
  269. void CFractalView::OnTimer(UINT nIDEvent) 
  270. {
  271. // TODO: Add your message handler code here and/or call default
  273. static float tree_k = 9;
  274. static float sier_k = 0;
  275. static float flower_k = 0;
  277. //在内存画布中清屏
  278. // if (CROWN_TIMER!=nIDEvent)
  279. {
  280. m_MemDC.FillSolidRect(&m_Rect, m_bgColor);
  281. }
  283. switch(nIDEvent){
  284. /* 分形树定时器 */
  285. case TREE_TIMER:
  286. int x, y;
  288. x = m_Rect.Width()/2;
  289. y = m_Rect.Height()-50;
  291. //在内存画布中绘制树
  292. drawLeaf(x, y, m_Rect.Height()/5, 270, 50, tree_k);
  294. if (m_dir) {
  295. tree_k += 0.5;
  296. if (tree_k > 18)
  297. m_dir = FALSE;
  298. } else {
  299. tree_k -= 0.5;
  300. if (tree_k < -18)
  301. m_dir = TRUE;
  302. }
  303. break;
  304. /* Sierpinski三角形定时器 */
  305. case SIERPINSKI_TIMER:
  306. if (m_dir) {
  307. sier_ifs(&m_MemDC, sier_k+=0.01);
  308. if (sier_k > 0.5)
  309. m_dir = FALSE;
  310. } else {
  311. sier_ifs(&m_MemDC, sier_k-=0.01);
  312. if (sier_k<0)
  313. m_dir =TRUE;
  314. }
  315. break;
  316. /* 万花筒定时器 */
  317. case FLOWER_TIMER:
  318. if (m_dir) {
  319. flower_ifs(&m_MemDC, flower_k+=5);
  320. if (flower_k > 500)
  321. m_dir = FALSE;
  322. } else {
  323. flower_ifs(&m_MemDC, flower_k-=5);
  324. if (flower_k<-500)
  325. m_dir =TRUE;
  326. }
  327. break;
  328. /* 王冠定时器 */
  329. case CROWN_TIMER:
  330. drawCrown();
  331. break;
  332. default:
  333. break;
  334. }
  336. //把树拷贝到物理画布中
  337. m_DC.BitBlt(0, 0, m_Rect.Width(), m_Rect.Height(), &m_MemDC, 0, 0, SRCCOPY);
  339. CView::OnTimer(nIDEvent);
  340. }
  342. void CFractalView::drawLeaf(double x, double y, double L, double A,float B,float C){
  343. double s1, s2, s3;
  344. s1 = 2;
  345. s2 = 3;
  346. s3 = 1.3;
  347. int x1,y1,x1L,y1L,x1R,y1R,x2,y2,x2R,y2R,x2L,y2L;
  348. if(L > s1)
  349. {
  350. x2 = (int)(x + L * cos(A * PI_tree));
  351. y2 = (int)(y + L * sin(A * PI_tree));
  352. x2R = (int)(x2 + L / s2 * cos((A + B) * PI_tree));
  353. y2R = (int)(y2 + L / s2 * sin((A + B) * PI_tree));
  354. x2L = (int)(x2 + L / s2 * cos((A - B) * PI_tree));
  355. y2L = (int)(y2 + L / s2 * sin((A - B) * PI_tree));
  357. x1 = (int)(x + L / s2 * cos(A * PI_tree));
  358. y1 = (int)(y + L / s2 * sin(A * PI_tree));
  359. x1L = (int)(x1 + L / s2 * cos((A - B) * PI_tree));
  360. y1L = (int)(y1 + L / s2 * sin((A - B) * PI_tree));
  361. x1R = (int)(x1 + L / s2 * cos((A + B) * PI_tree));
  362. y1R = (int)(y1 + L / s2 * sin((A + B) * PI_tree));
  365. m_MemDC.MoveTo((int)x,(int)y);
  366. m_MemDC.LineTo(x2,y2);
  367. m_MemDC.LineTo(x2R,y2R);
  368. m_MemDC.MoveTo(x2,y2);
  369. m_MemDC.LineTo(x2L,y2L);
  370. m_MemDC.MoveTo(x1,y1);
  371. m_MemDC.LineTo(x1L,y1L);
  372. m_MemDC.MoveTo(x1,y1);
  373. m_MemDC.LineTo(x1R,y1R);
  375. drawLeaf(x2, y2, L / s3, A + C, B, C);
  376. drawLeaf(x2R, y2R, L / s2, A + B, B, C);
  377. drawLeaf(x2L, y2L, L / s2, A - B, B, C);
  378. drawLeaf(x1L, y1L, L / s2, A - B, B, C);
  379. drawLeaf(x1R, y1R, L / s2, A + B, B, C);
  380. }  
  381. }
  383. void CFractalView::sier_ifs(CDC *pDC, float k)
  384. {
  385. float x=0; // 仿射变换中的自变量
  386. float y=0;
  387. float newx, newy; //仿射变换产生的新点
  388. float a, b, c, d, e, f; //仿射变幻中的系数
  389. long n=10000;  // 迭代次数
  390. float R; //随机变量
  391. float m[7][7]={0};    // 存放IFS码
  393. //IFS码赋值
  394. m[0][0] = 0.5; m[0][1] = 0; m[0][2] = 0;  m[0][3] = k;    m[0][4] = 0;    m[0][5] = 0;   m[0][6] = 0.333;
  395. m[1][0] = 0.5; m[1][1] = 0; m[1][2] = 0;  m[1][3] = 0.5;  m[1][4] =0.5;   m[1][5] = 0;   m[1][6] = 0.333;
  396. m[2][0] = k;   m[2][1] = 0; m[2][2] = 0;  m[2][3] = 0.5;  m[2][4] = 0.25; m[2][5] = 0.5; m[2][6] = 0.334;
  398. //循环迭代，在不同的概率空间下，赋不同的IFS码值
  399. while (n > 0) {
  400.         R = (float)rand()/RAND_MAX;
  401. if (R <= m[0][6]) {
  402. a = m[0][0]; b = m[0][1]; c = m[0][2]; d = m[0][3]; e = m[0][4]; f = m[0][5];
  403. } else if (R <= m[0][6] + m[1][6]) {
  404. a = m[1][0]; b = m[1][1]; c = m[1][2]; d = m[1][3]; e = m[1][4]; f = m[1][5];
  405. } else  {
  406. a = m[2][0]; b = m[2][1]; c = m[2][2]; d = m[2][3]; e = m[2][4]; f = m[2][5];
  407. } 
  408. newx = (a * x) + (b * y) + e;
  409. newy = (c * x) + (d * y) + f;
  410. x = newx; y = newy;
  411. pDC->SetPixelV((1000 + 7000 * x)/15, (6500 - 6000 * y)/15, RGB(x * 500 * R, R * 100, y * 2000 * R));
  412. n--;
  413. }
  414. }
  416. void CFractalView::flower_ifs(CDC *pDC, float k){
  417. float x=0; // '仿射变换中的自变量
  418. float y=0;
  419. float newx, newy; //'仿射变换产生的新点
  420. float a, b, c, d, e, f; //'仿射变幻中的系数
  421. long n=10000;  // '迭代次数
  422. float R; //'随机变量
  423. float m[7][7]={0};    // '存放IFS码
  425. //'IFS码赋值
  426. m[0][0] = 0.2;  m[0][1] = 0.2;  m[0][2] = 0;  m[0][3] = 0;  m[0][4] = 0.7;  m[0][5] = 0;     m[0][6] = 0.2;
  427. m[1][0] = 0.2;  m[1][1] = 0.2;  m[1][2] = 0;  m[1][3] = 0;  m[1][4] =-0.7;  m[1][5] = 0;     m[1][6] = 0.2;
  428. m[2][0] = 0.2;  m[2][1] = 0.2;  m[2][2] = 0;  m[2][3] = 0;  m[2][4] = 0;    m[2][5] = 0.7;   m[2][6] = 0.2;
  429. m[3][0] = 0.2;  m[3][1] = 0.2;  m[3][2] = 0;  m[3][3] = 0;  m[3][4] = 0;    m[3][5] = -0.7;  m[3][6] = 0.2;
  430. m[4][0] = 0.85; m[4][1] = 0.85; m[4][2] = k;  m[4][3] = k;  m[4][4] = 0;    m[4][5] = 0;     m[4][6] = 0.2;
  432. //'循环迭代，在不同的概率空间下，赋不同的IFS码值
  433. while (n > 0) {
  434.         R = (float)rand()/RAND_MAX;
  435. if (R <= m[0][6]) 
  436. {
  437. a = m[0][0]; b = m[0][1]; c = m[0][2]; d = m[0][3]; e = m[0][4]; f = m[0][5];
  438. } else if (R <= m[0][6] + m[1][6]) 
  439. {
  440. a = m[1][0]; b = m[1][1]; c = m[1][2]; d = m[1][3]; e = m[1][4]; f = m[1][5];
  441. } else if (R <= m[0][6] + m[1][6] + m[2][6]) 
  442. {
  443. a = m[2][0]; b = m[2][1]; c = m[2][2]; d = m[2][3]; e = m[2][4]; f = m[2][5];
  444. } else if (R <= m[0][6] + m[1][6] + m[2][6] + m[3][6]) 
  445. {
  446. a = m[3][0]; b = m[3][1]; c = m[3][2]; d = m[3][3]; e = m[3][4]; f = m[3][5];
  447. } else 
  448. {
  449. a = m[4][0]; b = m[4][1]; c = m[4][2]; d = m[4][3]; e = m[4][4]; f = m[4][5];
  450. }
  451. newx = (a * x*cos(c/180)) - (b * y*sin(d/180)) + e;
  452. newy = (a * x*sin(c/180)) + (b * y*cos(d/180)) + f;
  453. x = newx; y = newy;
  454. pDC->SetPixelV((300 + 200 * x), (220 - 200 * y), RGB(x * 500 * R, R * 100, y * 2000 * R));
  455. n--;
  456.   }
  457. }
  459. void CFractalView::drawCrown()
  460. {
  461. static float k;
  463.     crown.count++;
  464. for(int i=0;i<12000;i++)
  465. {   
  467. k=(float)rand();
  468. crown.rnd=(float)(k/RAND_MAX);
  470. crown.wx=crown.x-crown.cx;
  471. crown.wy=crown.y-crown.cy;
  472. if(crown.wx==0){
  473. crown.theta=PI/2;
  474. }
  475. else if(crown.wx>0){
  476. crown.theta=atan(crown.wy/crown.wx);
  477. }
  478. else if(crown.wx<0){
  479. crown.theta=PI-atan(crown.wy/crown.wx); 
  480. }
  482. crown.theta=crown.theta/2;
  483. crown.r=sqrt(crown.wx*crown.wx+crown.wy*crown.wy);
  485. if(crown.rnd<0.5){
  486. crown.r=sqrt(crown.r);
  487. }else{
  488. crown.r=-sqrt(crown.r);
  489. }
  490. crown.x=crown.r*cos(crown.theta);
  491. crown.y=crown.r*sin(crown.theta);
  493. m_MemDC.SetPixel((int)(2000*crown.x/15+3500/15+0.5),(int)(2000*crown.y/15+3500/15+0.5),0x00ff00);
  494. //Invalidate(false);
  496. if(crown.count==2) 
  497. {
  498. crown.count=0;
  499. }
  500. }
  501. crown.cx=crown.cx+0.002*crown.cxdirect;
  503. if(crown.cx>crown.cxmax && crown.cxdirect==1)
  504. {
  505. crown.cxdirect=-1;
  506. crown.cnt++;
  507. }
  508. else if(crown.cx<crown.cxmin && crown.cxdirect==-1)
  509. {
  510. crown.cxdirect=1;
  511. crown.cnt++;
  512. }
  513. if(crown.cnt%3==0 &&crown.cnt!=0)
  514. {
  515. crown.cy=crown.cy+0.003*crown.cydirect;
  516. }
  517. if(crown.cy>crown.cymax && crown.cydirect==1)
  518. {
  519. crown.cydirect=-1;
  520. }
  521. else if(crown.cy<crown.cymin && crown.cydirect==-1)
  522. {
  523. crown.cydirect=1;
  524. }
  525. }
  527. void CFractalView::OnDrawnew() 
  528. {
  529. // TODO: Add your command handler code here
  530. m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
  531. m_bgColor = RGB(0, 0, 0);
  532. m_fgColor = RGB(0, 255, 0);
  534. tree.m_Interval = 10;
  535. sier.m_Interval = 10;
  536. flower.m_Interval = 10;
  538. crown.m_Interval = 1;
  539. crown.m_C = 9;
  540. crown.x=5;
  541. crown.y=0;  //与x,y无关
  542. crown.m_pColor= RGB(0,255,0);
  543. crown.stop=false;
  544. crown.ondraw=false;
  545. crown.cxmin=-2;
  546. crown.cxmax=(float)0.8;
  547. crown.cymin=(float)-0.3;
  548. crown.cymax=(float)0.3;
  549.     crown.cxspeed=(float)0.2;
  550. crown.cyspeed=0.25;
  551.     crown.cx=crown.cxmin;
  552.     crown.cy=0;
  553. crown.cxdirect=1;
  554. crown.cydirect=-1;
  555. crown.cnt = 0;
  556. crown.count = 0;
  557. }