精密星历Lagrange插值法java实现浅析.docx

Lagrange算法 程序设计 专业：测绘工程 班级：1班 姓名：陈斯铭  语言：java 实现代码 import java.io.BufferedReader; import java.io.File; import java.io.FileReader; import java.io.IOException; import java.util.Scanner; public class Lagrange { public static void main(String[] args) throws IOException { Scanner sc = new Scanner(System.in); System.out.println(请输入文件名:); String fileName = sc.nextLine(); readFileByLines(fileName); } public static double Lagrange(int[] X,double[] Y,int X0){ int a,b; double t,u; t = 0; for (a = 0;a9;a++){ u = 1; for(b = 0;b9;b++){ if(b!= a){ u = u*(X0 - X[b])/(X[a]-X[b]); } } u = u*Y[a]; t = t + u; } return t; } public static void readFileByLines(String fileName) { File file = new File(fileName); BufferedReader reader = null; try { System.out.println(以行为单位读取文件内容，一次读一整行：); reader = new BufferedReader(new FileReader(file)); int[][] time1 = new int[5][96]; int[] time11 = new int[9]; double[][] satelite = new double[3][96]; double[] satelite1 = new double[9]; double[] satelite2 = new double[9]; double[] satelite3 = new double[9]; double[] time2 = new double[96]; String s1; String tempString = null; int line = 1; int i=0,n = 0,j=0,k=0,h=0; Scanner input = new Scanner(System.in); System.out.println(需要读取卫星的编号：); String numb = input.nextLine(); // 一次读入一行，直到读入null为文件结束 while ((tempString = reader.readLine()) != null) { // 显示行号 if(j96){ String m = tempString.substring(1,4).trim(); if(m.equals(numb)){ String ss1 = tempString.substring(4,18).trim(); satelite[0][j] = Double.parseD