精密星历Lagrange插值法java实现

1、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.pri

2、ntln("请输入文件名:"); 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;a<9;a+) u = 1; for(b = 0;b<9;b+) if(b!= a) u = u*(X0 - Xb)/(Xa-Xb); u = u*Ya; t = t + u; return t; public static vo

3、id readFileByLines(String fileName) File file = new File(fileName); BufferedReader reader = null; try System.out.println("以行为单位读取文件内容，一次读一整行："); reader = new BufferedReader(new FileReader(file); int time1 = new int596; int time11 = new int9; double satelite = new double396; double satelite

4、1 = new double9; double satelite2 = new double9; double satelite3 = new double9; double time2 = new double96; 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

5、(); / 一次读入一行，直到读入null为文件结束 while (tempString = reader.readLine() != null) / 显示行号 if(j<96) String m = tempString.substring(1,4).trim(); if(m.equals(numb) String ss1 = tempString.substring(4,18).trim(); satelite0j = Double.parseDouble(ss1); /System.out.println(satelite0j); /读取所需X坐标 String ss2 = tem

6、pString.substring(19,33).trim(); satelite1j = Double.parseDouble(ss2); /System.out.println(satelite1j); /读取所需Y坐标 String ss3 = tempString.substring(33,47).trim(); satelite2j = Double.parseDouble(ss3); /System.out.println(satelite2j); /读取所需Y坐标 j+; if(line=3) String s = tempString.substring(4,6).trim()

7、; n = Integer.parseInt(s); /System.out.println(n); if(line =23 + i*(n+1) && i<96) s1 = tempString.substring(3,7).trim(); time10i = Integer.parseInt(s1); /System.out.println(time10i); / 读取年 String s2 = tempString.substring(8,10).trim(); time11i = Integer.parseInt(s2); /System.out.println(t

8、ime11i); / 读取月 String s3 = tempString.substring(11,13).trim(); time12i = Integer.parseInt(s3); /System.out.println(time12i); / 读取日 String s4 = tempString.substring(14,16).trim(); time13i = Integer.parseInt(s4); /System.out.println(time13i); / 读取时 String s5 = tempString.substring(17,19).trim(); time1

9、4i = Integer.parseInt(s5); /System.out.println(time14i); / 读取分 String s6 = tempString.substring(20,22).trim(); time2i = Double.parseDouble(s6); /System.out.println(time2i); / 读取秒 i+; line+; reader.close(); Scanner x1 = new Scanner(System.in); Time time = new Time(); System.out.println("请输入要插入的点

10、的时:"); time.hour = x1.nextInt(); if(time.hour<0 | time.hour >24) System.out.println("输入错误,请重新输入："); time.hour = x1.nextInt(); System.out.println("请输入要插入的点的分:"); time.minute = x1.nextInt(); if(time.minute < 0 | time.minute > 60) System.out.println("输入错误,请重新输入

11、："); time.minute = x1.nextInt(); h = time.hour * 3600 + time.minute * 60; for(i=0;i<96;i+) if(time.hour = time13i && time.minute >= time14i && time.minute<time14i+1) for(k=0;k<9;k+) if(i<4) satelite1k = satelite0k; satelite2k = satelite1k; satelite3k = satelite2k;

12、 /System.out.println(satelite1k); time11k = k*15*60; else if(i>91) satelite1k = satelite087+k; satelite2k = satelite187+k; satelite3k = satelite287+k; /System.out.println(satelite1k); time11k = 87*15*60 + k*15*60; else satelite1k = satelite0i-4+k; satelite2k = satelite1i-4+k; satelite3k = satelit

13、e2i-4+k; /System.out.println(satelite1k); time11k = (i-4+k)*15*60; System.out.println(time.hour + "时" + time.minute + "分" + numb + "卫星的X坐标为：" + Lagrange(time11,satelite1,h);System.out.println(time.hour + "时" + time.minute + "分" + numb + "卫星的Y坐标为

14、：" + Lagrange(time11,satelite2,h);System.out.println(time.hour + "时" + time.minute + "分" + numb + "卫星的Z坐标为：" + Lagrange(time11,satelite3,h); catch (IOException e) e.printStackTrace(); finally if (reader != null) try reader.close(); catch (IOException e1) class Timeint year=2010;int month=1;int day=1;in