| 查看: 981 | 回复: 6 | ||
WANGGUANG123
|
[求助]
microsoft visual 2013已停止正常工作 已有2人参与
|
|
#include<iostream> #include<fstream> #include<math.h> #include<time.h> #include<stdio.h> #include<stdlib.h> #define N 10 using namespace std; double Min(double array[], int n) { double min = array[1]; if (n > 2) { for (int i = 1; i < n; ++i) { if (array < min) min = array; } } return min; } void main() { //输出文件流 输出到data.txt ofstream outfilex,outfiley,outfilez; outfilex.open("x.txt"  ;outfiley.open("y.txt" ;outfilez.open("z.txt" ;srand((unsigned)time(NULL)); double x[N], y[N], z[N];//定义1000个粒子的三个坐标数组 double r;//半径大小 double step, R;//步长和容限 cout << "输入半径r:" << " "; cin >> r; cout << endl << "输入移动步长step:" << " "; cin >> step; cout << endl << "输入容量大小R:" << " "; cin >> R; cout << endl; for (int i = 0; i < N; i++)//将每个粒子的坐标先全定义成000 { x = 0; y = 0; z = 0; } int a = 2; while (1) { //int hh = 0; x[a] = (rand() % 20000 - 10000) / 10000.0*R;//定义下一个粒子的随机坐标 y[a] = (rand() % 20000 - 10000) / 10000.0*R; z[a] = (rand() % 20000 - 10000) / 10000.0*R; while (1) { double distance[N]; for (int c = 1; c < a; c++) { distance[c] = (x[a] - x[c])*(x[a] - x[c]) + (y[a] - y[c])*(y[a] - y[c]) +(z[a] - z[c])*(z[a] - z[c]);//求两粒子间的距离 } if (Min(distance,a)<4.0* r*r)//如果和之前所有粒子的距离小于球心距，则此粒子生成成功，产生下一个粒子判定 { cout << a << " "; cout << Min(distance, a)<<" "<<x[a]<<" "<<y[a]<<" "<<z[a]<<endl; a++; break; } else//否则在该粒子位置处随机运动给定步长，进行下一次判定 { double theta = (rand() % 1000) / 1000.0*6.28; double fy = (rand() % 1000) / 1000.0*6.28; x[a] = x[a] + step*sin(theta)*cos(fy); y[a] = y[a] + step*sin(theta)*sin(fy); z[a] = z[a] + step*cos(theta); if (abs(x[a]) >= R || abs(y[a]) >= R || abs(z[a]) >= R) break; } } if (a == N+1)//到1000个粒子，退出 break; } system("pause" ;cout << "x y z坐标分别为" << endl; for (int j = 1; j < N+1 ; j++) { cout << j<< " " << x[j] << " " << y[j] << " " << z[j] << endl; //system("pause" ;//输出到文件 outfilex << x[j] << endl; outfiley << y[j] << endl; outfilez<<z[j] << endl; } //关闭输出文件流 outfilex.close(); outfiley.close(); outfilez.close(); system("pause" ;} 结果能运行出来，但是在关闭运行框的时候会报错 |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Microsoft Word
已经有6人回复
- Microsoft Office大全
已经有129人回复
- 2010microsoft office
已经有63人回复
- microsoft developer studio
已经有38人回复
zzy滴号
金虫 (小有名气)
- 应助: 18 (小学生)
- 金币: 5450.9
- 散金: 367
- 红花: 1
- 帖子: 282
- 在线: 897.8小时
- 虫号: 1540568
- 注册: 2011-12-17
- 性别: GG
- 专业: 理论和计算化学
2楼2015-11-15 10:48:04
lvweizhong
铜虫 (著名写手)
- 应助: 57 (初中生)
- 金币: 3181.3
- 红花: 10
- 沙发: 1
- 帖子: 2064
- 在线: 167.6小时
- 虫号: 3683385
- 注册: 2015-02-09
- 专业: 数论
3楼2015-11-15 18:32:41
lvweizhong
铜虫 (著名写手)
- 应助: 57 (初中生)
- 金币: 3181.3
- 红花: 10
- 沙发: 1
- 帖子: 2064
- 在线: 167.6小时
- 虫号: 3683385
- 注册: 2015-02-09
- 专业: 数论
4楼2015-11-15 18:33:39
zzy滴号
金虫 (小有名气)
- 应助: 18 (小学生)
- 金币: 5450.9
- 散金: 367
- 红花: 1
- 帖子: 282
- 在线: 897.8小时
- 虫号: 1540568
- 注册: 2011-12-17
- 性别: GG
- 专业: 理论和计算化学
5楼2015-11-15 21:57:42
WANGGUANG123 
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 695.4
- 散金: 195
- 红花: 1
- 帖子: 525
- 在线: 39小时
- 虫号: 3138440
- 注册: 2014-04-15
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2015-11-17 07:19:09
zzy滴号
金虫 (小有名气)
- 应助: 18 (小学生)
- 金币: 5450.9
- 散金: 367
- 红花: 1
- 帖子: 282
- 在线: 897.8小时
- 虫号: 1540568
- 注册: 2011-12-17
- 性别: GG
- 专业: 理论和计算化学
7楼2015-11-17 20:03:01
