EasyIFS分形学演示，命令行版

happy886rr

EasyIFS分形学演示工具是ZH@EASYX.CN用easyx库写的作品，体积比较大，并且不属于控制台程序，我用GDI库改写了该工具，使其可以在批处理中批量调用，并修复原版个别漏洞。体积也由112KB缩小到13KB。
下载附件：http://bcn.bathome.net/s/tool/index.html?key=EasyIFS
图片均为外链图

1. /*
   
2.         CONSOLE FRACTAL IFS, AUTHOR: ZH@EASYX.CN, MODIFYED BY HAPPY.
   
3.         EASYIFS.EXE
   
4. */
   
5. 
   
6. #include    <string>
   
7. #include    <vector>
   
8. #include   <stdio.h>
   
9. #include <windows.h>
   
10. #include   <conio.h>
    
11. #include    <time.h>
    
12. #include      <io.h>
    
13. 
    
14. using namespace std;
    
15. 
    
16. // 对VC6.0、VS2010、GCC作兼容性编译
    
17. #if defined (__GNUC__) || (_MSC_VER == 1200)
    
18. extern "C" HWND WINAPI GetConsoleWindow();
    
19. #endif
    
20. 
    
21. // 宏名变量定义
    
22. #define F_EXIST 4
    
23. #define E_WITH  640
    
24. #define E_HIGH  480
    
25. 
    
26. // 定义帮助说明
    
27. #define HELP_INFORMATION "\
    
28. Usage: easyifs [file]\n\
    
29. \n\
    
30.     The probability of requiring each formula and equ 1\n\
    
31.     IFS configuration file:\n\
    
32. \n\
    
33.     [Main]\n\
    
34.     name = IFS name\n\
    
35.     iterLimit = number of iterations\n\
    
36.     color = drawing color\n\
    
37.     minX = the lower limit of x\n\
    
38.     maxX = the upper limit of x\n\
    
39.     minY = the lower limit of y\n\
    
40.     maxY = the upper limit of y\n\
    
41.     condition = the number of formulas\n\
    
42.     [Conditionx]\n\
    
43.     p = probability\n"
    
44. 
    
45. 
    
46. // IFS 全局参数
    
47. struct MAIN
    
48. {
    
49.         string name;             // ifs 名称
    
50.         int iterLimit;           // 迭代次数
    
51.         int color;               // 绘图颜色
    
52.         double scale;            // 相对绘图窗口的缩放比例
    
53.         int offsetX, offsetY;    // 相对绘图窗口的偏移量
    
54. };
    
55. 
    
56. // 每个IFS公式的参数
    
57. struct IFS
    
58. {
    
59.         int p;                   // 概率
    
60.         double a, b, c, d, e, f; // IFS 公式参数
    
61. };
    
62. 
    
63. // 全局变量
    
64. MAIN g_main;
    
65. vector<IFS> g_ifs;
    
66. 
    
67. // 从配置文件中获取 double 类型数据
    
68. double GetPrivateProfileDouble(LPCTSTR lpAppName, LPCTSTR lpKeyName, double fDefault, LPCTSTR lpFileName)
    
69. {
    
70.         TCHAR d[50];
    
71.         GetPrivateProfileString(lpAppName, lpKeyName, NULL, d, 50, lpFileName);
    
72.         return (d[0]!=0)?atof(d):fDefault;
    
73. }
    
74. 
    
75. // 根据配置文件初始化 IFS 系统
    
76. int initargs(int argc, char** argv)
    
77. {
    
78.         // 如果未指定参数，显示帮助信息
    
79.         if (argc <= 1)
    
80.         {
    
81.                 fprintf(stdout, HELP_INFORMATION);
    
82.                 exit(1);
    
83.         }
    
84. 
    
85.         // 定义变量
    
86.         double ix, ax, iy, ay, sx, sy;
    
87.         // 配置文件的文件名
    
88.         char tmpFpath[MAX_PATH*2];
    
89.         if(strrchr(argv[1], ':')==NULL)
    
90.         {
    
91.                 strcpy(tmpFpath, ".\");
    
92.                 strcat(tmpFpath, argv[1]);
    
93.         }
    
94.         else
    
95.         {
    
96.                 strcpy(tmpFpath, argv[1]);
    
97.         }
    
98. 
    
99.         char* filename=tmpFpath;
    
100. 
     
101.         // 如果参数指定的文件不存在或禁止访问，退出系统
     
102.         if (_access(filename, F_EXIST) != 0)
     
103.         {
     
104.                 fprintf(stdout, "The file "%s" does not exist or is disabled.\n", argv[1]);
     
105.                 exit(2);
     
106.         }
     
107. 
     
108.         // 获取 IFS 名称
     
109.         char tmpName[MAX_PATH*2];
     
110.         GetPrivateProfileStringA("main", "name", "noname", tmpName, 50, filename);
     
111.         g_main.name = tmpName;
     
112. 
     
113.         // 获取迭代次数
     
114.         g_main.iterLimit = GetPrivateProfileIntA("main", "iterLimit", 0, filename);
     
115. 
     
116.         // 获取绘图颜色
     
117.         g_main.color = GetPrivateProfileIntA("main", "color", 0xff00, filename);
     
118. 
     
119.         // 获取 x 方向上的缩放比例
     
120.         ix = GetPrivateProfileDouble("main", "minX", 0, filename);
     
121.         ax = GetPrivateProfileDouble("main", "maxX", 0, filename);
     
122.         sx = E_WITH / (ax - ix);
     
123. 
     
124.         // 获取 y 方向上的缩放比例
     
125.         iy = GetPrivateProfileDouble("main", "minY", 0, filename);
     
126.         ay = GetPrivateProfileDouble("main", "maxY", 0, filename);
     
127.         sy = E_HIGH / (ay - iy);
     
128. 
     
129.         // 根据绘图窗口尺寸，确定恰当的缩放比例及偏移量
     
130.         if (sx > sy)
     
131.         {
     
132.                 g_main.scale = sy;
     
133.                 g_main.offsetX = (int)(-ix * g_main.scale + (E_WITH - (ax - ix) * g_main.scale) / 2);
     
134.                 g_main.offsetY = (int)(-iy * g_main.scale);
     
135.         }
     
136.         else
     
137.         {
     
138.                 g_main.scale = sx;
     
139.                 g_main.offsetX = (int)(-ix * g_main.scale);
     
140.                 g_main.offsetY = (int)(-iy * g_main.scale + (E_HIGH - (ay - iy) * g_main.scale) / 2);
     
141.         }
     
142. 
     
143.         // 获取公式数量
     
144.         int n;
     
145.         n = GetPrivateProfileIntA("main", "condition", 0, filename);
     
146. 
     
147.         // 获取每一个公式的常数项及概率
     
148.         IFS ifs;
     
149.         int sump = 0;
     
150.         char app[] = "condition?";
     
151.         for(int i=0; i<n; i++)
     
152.         {
     
153.                 app[9] = ('1' + i);
     
154.                 ifs.a = GetPrivateProfileDouble(app, "a", 0, filename);
     
155.                 ifs.b = GetPrivateProfileDouble(app, "b", 0, filename);
     
156.                 ifs.c = GetPrivateProfileDouble(app, "c", 0, filename);
     
157.                 ifs.d = GetPrivateProfileDouble(app, "d", 0, filename);
     
158.                 ifs.e = GetPrivateProfileDouble(app, "e", 0, filename);
     
159.                 ifs.f = GetPrivateProfileDouble(app, "f", 0, filename);
     
160.                 ifs.p = int(GetPrivateProfileDouble(app, "p", 0, filename) * 1000000 + 0.5);
     
161.                 ifs.p += sump;
     
162.                 sump = ifs.p;
     
163. 
     
164.                 g_ifs.push_back(ifs);
     
165.         }
     
166. 
     
167.         // 如果各公式的概率和不等于 1，返回错误信息
     
168.         if (sump != 1000000)
     
169.         {
     
170.                 fprintf(stdout, "IFS configuration file in the probability of the formula and not equal to 1, please check the configuration file.\n");
     
171.                 exit(3);
     
172.         }
     
173. 
     
174.         return 0;
     
175. }
     
176. 
     
177. // 主函数
     
178. int main(int argc, char* argv[])
     
179. {
     
180.         // 初始化IFS配置
     
181.         initargs(argc, argv);
     
182.         srand(unsigned(time(NULL)));
     
183. 
     
184.         // 初始化图形窗口
     
185.         HWND hCMD=GetConsoleWindow();
     
186.         //刷新窗口
     
187.         InvalidateRect(hCMD, NULL, TRUE);
     
188.        
     
189.         HDC  hDC=GetDC(hCMD);
     
190. 
     
191.         // 创建变量
     
192.         double x = 0, y = 0, tx;
     
193.         int p, k;
     
194. 
     
195.         //设置背景透明
     
196.         //SetBkMode(hDC,TRANSPARENT);
     
197.         //设置字体颜色
     
198.         //SetTextColor(hDC, RGB(255,255,255));
     
199.         // 显示名称
     
200.         //TextOutA(hDC, 0,0, g_main.name.c_str(), strlen(g_main.name.c_str()));
     
201. 
     
202.         // 迭代求解
     
203.         for(int i=0; i<g_main.iterLimit; i++)
     
204.         {
     
205.                 // 生成概率
     
206.                 p = int(double(rand()) / RAND_MAX * 1000000 + 0.5);
     
207. 
     
208.                 // 根据概率获取选用的公式 k
     
209.                 for(k=0; k< (int)g_ifs.size(); k++)
     
210.                         if (p <= g_ifs[k].p)        break;
     
211. 
     
212.                 // 根据公式 k 迭代
     
213.                 tx = g_ifs[k].a * x + g_ifs[k].b * y + g_ifs[k].e;
     
214.                 y  = g_ifs[k].c * x + g_ifs[k].d * y + g_ifs[k].f;
     
215.                 x  = tx;
     
216. 
     
217.                 // 画点（转换到屏幕坐标系）
     
218.                 SetPixel(hDC, int(x * g_main.scale) + g_main.offsetX, E_HIGH - (int(y * g_main.scale) + g_main.offsetY),  g_main.color);
     
219.         }
     
220.         return 0;
     
221. }
     

复制代码