1

私のコードでは、while ループの前にステップ サイズ "h" を定義しています。ループで使ってみるとなんとなく勝手に変わっていくような気がします。ループ内で定義すると問題ないように見えますが、取得したデータが正しくないように見えるので、問題が関連している可能性があると推測しています。

この場所に印刷しても (コードの printf を参照)、出力は 3 つの値であり、その理由はわかりません。関連していないが間違っていると思われるものを他に見つけた場合は、教えてください。予期しない値が得られていると言ったように (それは私の数式かもしれません)。

int main()
{
    FILE *f1;
    f1 = fopen("Question2 part 3 solution.txt", "w");
   double r0=0.05;
   double dr0=-a*a*r0/sqrt(1+pow(a*a*r0,2)),h=0.01;
   double k[4][3],x[]={dr0,z0,T0,r0},x1[]={0,0,0,0}, s=0;
   int i,j;

   while(s<=1)
    {
        //Runge-Kutta
        for (j=0;j<4;j++)
        {
            for (i=0;i<4;i++)
            {
                if (j==0){k[i][0]=h*System(i,x[0],x[1],x[2],x[3]);}
                if (j==1){k[i][1]=h*System(i,x[0]+k[0][0]/2.0,x[1]+k[1][0]/2.0,x[2]+k[2][0]/2.0,x[3]+k[3][0]/2.0);}
                if (j==2){k[i][2]=h*System(i,x[0]+k[0][1]/2.0,x[1]+k[1][1]/2.0,x[2]+k[2][1]/2.0,x[3]+k[3][1]/2.0);}
                if (j==3){k[i][3]=h*System(i,x[0]+k[0][2],x[1]+k[1][2],x[2]+k[2][2],x[3]+k[3][2]);}
            }
        }

        for (i=0;i<4;i++)
        {
            x[i]=x[i]+(k[i][0]+2.0*k[i][1]+2.0*k[i][2]+k[i][3])/6.0;
        }
printf("%8.3lf",h);
        s+=h;
    }

   fclose(f1);
    return(0);
}


double System(int i,double dr, double z, double T, double r)
    {
    //printf("%e\t%e\t%e\t%e\n",dr,z,T,r);
    if (T==T0 && z==z0 && i==0) {return  (-a*a*dr)*pow(1-dr*dr,3/2)/2.0;}
    if (i==0 && T!=0){return  (-a*a*r*(1-dr*dr)-dr*sqrt(1-dr*dr))/T;}
    if (i==1){return  (-sqrt(1-dr*dr));}
    if (i==2){return  (-a*a*r*dr+sqrt(1-dr*dr));}
    if (i==3){return  (dr);}
    //if (i==3){return  (-m2*l1*l2*B*theta1dt*theta2dt*sin(theta2-theta1)-l2*m2*g*B*sin(theta2));}
    }

前もって感謝します！

4

2 に答える 2

3

の宣言を参照してくださいk。

double k[4][3]

そして、この声明を見てください

k[i][3]=...

ここでは、配列の境界を超えて書き込み、未定義の動作につながります。

于 2013-11-10T10:27:07.707 に答える
3

メモリをオーバーランしています。変数kは として定義されていますが、次の場合double k[4][3]に更新していますk[i][3]j==3

于 2013-11-10T10:27:23.360 に答える