C++ C++;更改实例变量
我正试图写一个程序来计算一个n体问题。我遇到了一个问题,试图改变我的代码,以便它可以很容易地适应任何数量的机构。函数corr()似乎有问题,在使用函数corr()后,对该函数中类粒子实例的某些变量所做的更改似乎丢失了。我如何解决这个问题C++ C++;更改实例变量,c++,C++,我正试图写一个程序来计算一个n体问题。我遇到了一个问题,试图改变我的代码,以便它可以很容易地适应任何数量的机构。函数corr()似乎有问题,在使用函数corr()后,对该函数中类粒子实例的某些变量所做的更改似乎丢失了。我如何解决这个问题 #include <cstdlib> #include <iostream> #include <cmath> #include <fstream> #define h 1000.0 #define G 6.6
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <fstream>
#define h 1000.0
#define G 6.67384*pow(10.0,-11)
using namespace std;
class particle{
public:
double kx1,kx2,kx3,kx4, kv1, kv2, kv3, kv4;
double ky1, ky2, ky3, ky4, kvy1, kvy2, kvy3, kvy4;
double x,y,vx,vy,m;
double dist(particle aap){
double dx = x - aap.x;
double dy = y - aap.y;
return sqrt(pow(dx,2.0)+pow(dy,2.0));
}
double g(double x1, double y1,particle aap){
return G*aap.m*(aap.x-x1)/pow(dist(aap),3.0);
}
double p(double x1, double y1, particle aap){
return G*aap.m*(aap.y-y1)/pow(dist(aap),3.0);
}
void update(){ //zet het object 1 stap vooruit
x = x + (1/6.0)*(kx1+2*kx2+2*kx3+kx4);
vx = vx + (1/6.0)*(kv1+2*kv2+2*kv3+kv4);
y = y + (1/6.0)*(ky1+2*ky2+2*ky3+ky4);
vy = vy + (1/6.0)*(kvy1+2*kvy2+2*kvy3+kvy4);
}
void create(double x1, double y1, double vx1, double vy1, double m1){
x = x1;
y = y1;
vx = vx1;
vy = vy1;
m =m1;
}
bool operator ==(particle &other){
if(x == other.x && y == other.y && vx == other.vx && vy == other.vy){
return true;
}
}
};
particle zon, maan, aarde;
void set(){
zon.create(1, 1, -2, 1, 2*pow(10.0,30));
aarde.create(1.5*pow(10.0,11), 0, 2, 29780, 6*pow(10.0,24));
maan.create(aarde.x + 1, aarde .y + 3.844399*pow(10.0,8), aarde.vx + -1022.0, aarde.vy + 1, 7.3347*pow(10.0,22));
}
double xforce(double x1, double y1, particle aap){ //kracht in de x-richting
particle bodies[] = {zon, aarde, maan};
double fx;
for (int i = 0; i < 3; i++){
if (bodies[i].x == aap.x && bodies[i].y == aap.y && bodies[i].vx == aap.vx && bodies[i].vy == aap.vy ){;}
else{
fx += aap.g(x1,y1,bodies[i]);
}
}
return fx;
}
double yforce(double x1, double y1, particle aap){ //kracht in de y-richting
particle bodies[] = {zon, aarde, maan};
double fy;
for (int i = 0; i <= 3; i++){
if (bodies[i].x == aap.x && bodies[i].y == aap.y && bodies[i].vx == aap.vx && bodies[i].vy == aap.vy) {;}
else{
fy += aap.p(x1,y1,bodies[i]);
}
}
return fy;
}
void corr(particle body){
body.kx1 = h*body.vx;
body.kv1 = h*xforce(body.x, body.y, body);
body.ky1 = h*body.vy;
body.kvy1 =h*yforce(body.x, body.y, body);
body.kx2 = h*(body.vx + 0.5*body.kv1);
body.kv2 = h*xforce(body.x + 0.5*body.kx1, body.y + 0.5*body.ky1, body);
body.ky2 = h*(body.vy + 0.5*body.kvy1);
body.kvy2 = h*yforce(body.x + 0.5*body.kx1, body.y + 0.5*body.ky1, body);
body.kx3 = h*(body.vx+ 0.5*body.kv2);
body.kv3 = h*xforce(body.x + 0.5*body.kx2, body.y + 0.5*body.ky2, body);
body.ky3 = h*(body.vy+ 0.5*body.kvy2);
body.kvy3 = h*yforce(body.x + 0.5*body.kx2, body.y + 0.5*body.ky2,body);
body.kx4 = h*(body.vx+body.kv3);
body.kv4 = h*xforce(body.x+ body.kx3, body.y + body.ky3, body);
body.ky4 = h*(body.vy + body.kvy3);
body.kvy4 = h*yforce(body.x + body.kx3, body.y + body.ky3, body);
}
void bereken(){
set();
ofstream file;
file.open("3body.txt");
for(int i =0; i <=30000; i++){
corr(maan);
corr(zon);
corr(aarde);
zon.update();
aarde.update();
maan.update();
file << i*h <<" "<< zon.x << " "<< zon.y << " "<< zon.vx<< " "<< zon.vy <<" "<< aarde.x << " " << aarde.y <<" "<< aarde.vx <<" " << aarde.vy <<" "<< maan.x<<" "<<maan.y<<"\n";
}
file.close();
}
int main()
{
set();
bereken();
system("pause");
return 0;
}
#包括
#包括
#包括
#包括
#定义h1000.0
#定义G 6.67384*功率(10.0,-11)
使用名称空间std;
类粒子{
公众:
双kx1、kx2、kx3、kx4、kv1、kv2、kv3、kv4;
双ky1、ky2、ky3、ky4、kvy1、kvy2、kvy3、kvy4;
双x,y,vx,vy,m;
双距离(粒子aap){
双dx=x-aap.x;
双dy=y-aap.y;
返回sqrt(pow(dx,2.0)+pow(dy,2.0));
}
双g(双x1,双y1,粒子aap){
返回G*aap.m*(aap.x-x1)/pow(距离(aap),3.0);
}
双p(双x1,双y1,粒子aap){
返回G*aap.m*(aap.y-y1)/pow(距离(aap),3.0);
}
void update(){//zet het对象1 stap vooruit
x=x+(1/6.0)*(kx1+2*kx2+2*kx3+kx4);
vx=vx+(1/6.0)*(kv1+2*kv2+2*kv3+kv4);
y=y+(1/6.0)*(ky1+2*ky2+2*ky3+ky4);
vy=vy+(1/6.0)*(kvy1+2*kvy2+2*kvy3+kvy4);
}
创建空洞(双x1、双y1、双vx1、双vy1、双m1){
x=x1;
y=y1;
vx=vx1;
vy=vy1;
m=m1;
}
布尔运算符==(粒子和其他){
如果(x==other.x&&y==other.y&&vx==other.vx&&vy==other.vy){
返回true;
}
}
};
阿尔德马恩颗粒区;
空集(){
区域创建(1,1,-2,1,2*pow(10.0,30));
创建(1.5*pow(10.0,11),0,229780,6*pow(10.0,24));
创建(aarde.x+1,aarde.y+3.844399*功率(10.0,8),aarde.vx+-1022.0,aarde.vy+1,7.3347*功率(10.0,22));
}
双xforce(双x1,双y1,粒子aap){//kracht in de x-richting
粒子体[]={zon,aarde,maan};
双外汇;
对于(int i=0;i<3;i++){
如果(bodies[i].x==aap.x&&bodies[i].y==aap.y&&bodies[i].vx==aap.vx&&bodies[i].vy==aap.vy{;}
否则{
fx+=aap.g(x1,y1,body[i]);
}
}
外汇收益;
}
双y力(双x1,双y1,粒子aap){//y-richting中的kracht
粒子体[]={zon,aarde,maan};
双fy;
对于(int i=0;i只需使用引用:
void corr(粒子和物体)
同样的修复也必须应用于其他地方。您当前拥有的是一个代码,它复制您的对象,对其进行一些计算,然后删除此临时对象…通过使用引用(&
),您可以在调用者和被调用者之间有效地“共享”对象
在你的类方法中,你也应该使用引用来传递这个对象,但这一次是出于效率的考虑——传递引用很简单(通常与指针相同),但是复制这个巨大的对象需要时间。Freddie Chopin的建议奏效了,并修复了函数corr()但是,现在我发现函数xforce()和yforce()不能正常工作。它们应该返回某个粒子上的总力,这是它们在第一步中所做的,但之后它们似乎返回了错误的值。有人知道为什么会这样吗?void corr(粒子体)形体参数正在被值传递。我对C++和编程都是相当新的,你能解释得更深入吗?请阅读C++教程,引用是一个基本的特性。非常感谢你的AWNSER,它似乎已经解决了我的问题。不要做他们应该做的事情(即计算某个粒子上的总作用力)。知道怎么回事吗。