计算物理习题解答(徐文静).docVIP

/** * 计算物理1.1 * @π值计算(通过行列式计算) * @雷闽松*/ public class Evaluatepi { /** * 已知:πκ=π∞+c1/κ+c2/κ^2+c3/κ^3+... * π8 = 3.313708,π16 = 3.182598,π32 =3.151725,π64 = 3.144118 * * 近似:πκ≈π∞+c1/κ+c2/κ^2+c3/κ^3 * * 因此，通过求4个未知量的线性方程组: * c1/8+c2/64+c3/512+π∞ = 3.313708 * c1/16+c2/256+c3/4096+π∞ = 3.182598 * c1/32+c2/1024+c3/32768+π∞ = 3.151725 * c1/64+c2/4096+c3/262144+π∞ = 3.144118 * 的系数行列式即可求得c1,c2,c3和π * * 系数行列式: * { * {1.0/8, 1.0/64, 1.0/512, 1}, * {1.0/16, 1.0/256, 1.0/4096, 1}, * {1.0/32, 1.0/1024, 1.0/32768, 1}, * {1.0/64, 1.0/4096, 1.0/262144,1} * } */ public static void main(String[] args){ double D = envaluate(new data().xxxx); double[][] temp = new double[4][4]; double[] result = new double[4]; for(int i=0;i4;i++){ System.arraycopy(new data().xxxx, 0, temp, 0, 4); temp[0][i] = 3.061467; temp[1][i] = 3.121445; temp[2][i] = 3.136548; temp[3][i] = 3.140331; double Di = envaluate(temp); result[i] = Di/D; } System.out.println(当π8=3.061467,π16=3.121445,π32=3.136548,π64=3.140331时:); System.out.println( c1的值:+result[0]); System.out.println( c2的值:+result[1]); System.out.println( c3的值:+result[2]); System.out.println( π的值:+result[3]); for(int i=0;i4;i++){ System.arraycopy(new data().xxxx, 0, temp, 0, 4); temp[0][i] = 3.313708; temp[1][i] = 3.182598; temp[2][i] = 3.151725; temp[3][i] = 3.144118; double Di = envaluate(temp); result[i] = Di/D; } System.out.println(\n然而当π8=3.313708,π16=3.182598,π32=3.151725,π64=3.144118时:); System.out.println( c1的值:+result[0]); System.out.println( c2的值:+result[1]); System.out.println( c3的值:+result[2]); System.out.println( π的值:+result[3]); System.out.println(\n由此可见第一组数据更合理。); } public static double envaluate(double[][] d){ int len = d.length; tree[] trees = new tree[len]; for(int i=0;ilen;i++){ int sign = i%2==0?1:-1; double a = d[i][0]; double