C语言计算GPS卫星位置
GPS广播星历计算卫星位置和速度-Read
“GPS广播星历计算卫星位置和速度”及“GPS伪距定位”计算试验1．试验内容及上交成果1.1 试验内容应用C语言按预定格式（函数、输入输出变量之名称、类型）编写“GPS广播星历计算卫星位置和速度”函数SatPos_Vel( )、“GPS伪距定位”函数Positioning( )。
将此两个函数组成文件F2.cpp，并包含于文件GPS_Positioning.cpp中。
编译、连接并运行文件GPS_Positioning.cpp，逐一时刻读取广播星历（Ephemeris.dat）、观测时间及伪距、卫星号（Observation.dat）信息，计算WGS84坐标系中观测时刻相应的卫星位置、速度以及载体位置，结果保存于文件Position.dat中。
1.2 上交成果磁盘文件F2.cpp、Position.dat，并存于“学号作者中文姓名”目录中。
2．函数说明2.1 星历文件读取函数void EFileReading(Efile)功能：读取星历文件，给星历数据结构体Ephemeris赋值。
输入变量：EFile 字符串，文件名。
2.2 观测数据读取函数int ObsReading(fp_Obs,Time,Rho,Mark)功能：从文件Observation.dat中读取某一时刻的伪距、卫星号。
读取成功函数值返回“1”，失败返回“－1”（读错，或至文件尾）。
输入变量：fp_Obs 文件指针；输出变量：Time double，时间（秒）；Rho double[12]，伪距（米）；Mark int[12]，卫星号，“－1”表示此通道无卫星、无伪距。
2.3 最小二乘估计函数int LeastSquareEstimation(Y,A,P,m,n, X)功能：最小二乘方法求解观测方程Y=AX+ε，其中观测值方差阵的逆阵为P（也称为权阵），得未知参数X。
成功返回“1”，失败返回“－1”（亏秩）。
输入变量：Y double[m]，观测方程自由项（米）；A double[m×n]，系数阵（无量纲），按第1行第1、2……n元素，第2行……顺序存放；P double[m]，权矩阵对角线元素，0或1（无量纲）；m int，观测值个数；n int，未知参数个数；输出变量：X double[n]，未知参数（米）。
坐标转换源代码--GPS定位程序(C--)
坐标转换源代码--GPS定位程序(C++)GPS数据处理中为了满足不同的需要，处理的数据要进行坐标转换，得到在不同坐标系统下的结果，下面是笛卡尔坐标系，大地坐标系，站心地平坐标系（线型和极坐标形式）之间的转换源代码：头文件：#ifndef _COORDCOVERT_H#define _COORDCOVERT_H#include "stdlib.h"//WGS-84椭球体参数const double a=6378137.0;//长半轴const double flattening=1/298.257223563;//扁率const double delta=0.0000001;typedef struct tagCRDCARTESIAN{double x;double y;double z;}CRDCARTESIAN;typedef CRDCARTESIAN *PCRDCARTESIAN;//笛卡尔坐标系typedef struct tagCRDGEODETIC{double longitude;double latitude;double height;}CRDGEODETIC;typedef CRDGEODETIC *PCRDGEODETIC;//大地坐标系typedef struct tagCRDTOPOCENTRIC{double northing;double easting;double upping;}CRDTOPOCENTRIC;typedef CRDTOPOCENTRIC *PCRDTOPOCENTRIC; //站心地平坐标系（线坐标形式）typedef struct tagCRDTOPOCENTRICPOLAR{ double range;double azimuth;double elevation;}CRDTOPOCENTRICPOLAR;typedef CRDTOPOCENTRICPOLAR *PCRDTOPOCENTRICPOLAR;//站心地平坐标系（极坐标形式）//由笛卡尔坐标转换为大地坐标void CartesianToGeodetic (PCRDGEODETIC pcg, PCRDCARTESIAN pcc, double dSemiMajorAxis, double dFlattening);//pcg：指向所转换出的大地坐标的指针；//pcc：指向待转换的笛卡尔坐标的指针；//dSemiMajorAxis：参考椭球的长半轴；//dFlattening：参考椭球的扁率。
GPS卫星位置的计算(C++程序计算)
GPS卫星位置的计算利用C++编写了一段能计算单一瞬时卫星坐标的程序，在运行程序之前，需做部分准备工作：（1）在F盘下建立一名为“单一卫星广播星历”的txt文件。
（2）从“广播星历.txt”文件中拷贝从卫星PRN号开始的8行数据到“单一卫星广播星历.txt”中（3）在编辑选项中，将全部的“D”替换为“E”。
下面为我所选取的一个广播星历：18 06 8 25 6 0 0.0-2.472363412380E-04-1.023*********E-12 0.000000000000E+001.410000000000E+02-1.721875000000E+01 4.502687555010E-09 1.413760604187E+00-7.990747690201E-07 7.598234573379E-03 1.118145883083E-05 5.153709835052E+034.536000000000E+05-1.303851604462E-08-1.095067942661E-01 1.527369022369E-079.571235745530E-01 1.640000000000E+02-2.656176299285E+00-8.0374********E-09-5.193073455211E-10 1.000000000000E+00 1.389000000000E+03 0.000000000000E+002.000000000000E+00 0.000000000000E+00-1.024*********E-08 1.410000000000E+024.464490000000E+05 4.000000000000E+00程序设计部分：#include<stdio.h>#include<math.h>int main(){int i = 0;double n[50], n0, nn, t, tk, Mk, Ek, Vk, Yk, Gu, Gr, Gi, uk, rk, ik, xk, yk, zk, X, Y, Z, Lk, UT, yy, mm, JD, gpsz;FILE *fp;fp = fopen("F:\\单一卫星广播星历.txt", "r");if (fp == NULL){printf ("文件打开失败!\n");return 0;}while (! feof (fp)){fscanf(fp, "%lf", &n[i]);i++;}n0 = (sqrt(3986005E+8))/pow(n[17], 3);nn = n0 + n[12];/*计算卫星运行的平均角速度*/UT = n[4] + (n[5] / 60) + (n[7] / 3600);/*民用日的时分秒化为实数时*/if (n[1] >= 80)/*广播星历中年只有后两位，化为4位，参考1980年1月6日0点*/ {if (n[1] == 80 && n[2] == 1 && n[3] < 6){n[1] = n[1] + 2000;}n[1] = n[1] + 1900;}else{n[1] = n[1] + 2000;}if (n[2] <= 2){yy = n[1] - 1;mm = n[2] + 12;}if (n[2] > 2){yy = n[1];mm = n[2];}JD = (int)(365.25 * yy) + (int)(30.6001 * (mm + 1)) + n[3] + (UT / 24) + 1720981.5;/*化为儒略日*/gpsz = (int)((JD - 2444244.5) / 7);/*计算GPS周*/t = (JD - 2444244.5 - 7 * gpsz) * 24 * 3600;/*得出GPS秒*/tk = t - n[18];/*tk1为中间值，用以判断tk与正负302400的关系，然后返回到tk上*/ while (tk > 302400 || tk < -302400){if (tk > 302400){tk = tk - 604800;}else{tk = tk + 604800;}}/*计算归化观测时间*/Mk = n[13] + nn * tk;/*观测时刻的卫星平近点角*/Ek = Mk;Ek = Mk + n[15] * sin(Ek);Ek = Mk + n[15] * sin(Ek);/*迭代两次计算观测时刻的偏近点角*/Vk = atan(sqrt(1 - n[15] * n[15]) * sin(Ek)) / (cos(Ek) - n[15]);/*真近点角*/Yk = Vk + n[24];/*升交距角*/Gu = n[14] * cos(2 * Yk) + n[16] * sin(2 * Yk);Gr = n[23] * cos(2 * Yk) + n[11] * sin(2 * Yk);Gi = n[19] * cos(2 * Yk) + n[21] * sin(2 * Yk);/*摄动改正项*/uk = Yk + Gu;rk = n[17] * n[17] * (1 - n[15] * cos(Ek)) + Gr;ik = n[22] + Gi + n[26] * tk;/*经摄动改正后的升交距角、卫星矢径、轨道倾角*/xk = rk * cos(uk);yk = rk * sin(uk);zk = 0;/*卫星在轨道坐标系的坐标*/Lk = n[20] + (n[25] - 7.29211515E-5) * tk - 7.29211515E-5 * n[18];/*观测时刻t的升交点经度*/X = xk * cos(Lk) - yk * cos(ik) * sin(Lk);Y = xk * sin(Lk) + yk * cos(ik) * cos(Lk);Z = yk * sin(ik);/*卫星在WGS-84坐标系的坐标*/printf("该卫星在WGS-84坐标系中的坐标为：\nX = %lf m\nY = %lf m\nZ = %lf m\n", X, Y, Z);fclose(fp);return 0;}计算结果：该卫星在WGS-84坐标系中的坐标为：X = 9223153.692525 mY = 24133486.931401 mZ = 6032585.919385 m。
基于C语言的GPS解析程序
if(x==4||x==7) {}
else dat[x]=ch[y++];
}
break;
case 5: //东西经
dat[0]=ch[0];
break;
case 6: //航速
for(x=0;x<6;x++)
while((RMC_Data[i]!=',')&&(RMC_Data[i]!='*')) //RMC_Data[i]==','或者=='*'时结束循环。
{i++;}
if(j+1!=i)
{
for(j=j+1,k=0;j<i;j++)
{
uchar x,y;
send_byte('d');
switch(ch1)
{
case 0: //时间
for(x=0,y=0;x<8;x++)
{
if(x==2||x==5){time[x]=0x3a;}
else time[x]=ch[y++];
break;
case '*':
i=strlen(RMC_Data); //处理完成，软件结束。
break;
default:break;
}
}
// send_byte('e');
}
void GPRMC_Handle(uchar *ch,uchar ch1) //GPRMC信息，与，之间的子信息处理，存储到相应的结构体中元素中
C语言计算GPS卫星位置
C语言计算GPS卫星位置GPS（全球定位系统）是一种利用卫星定位来确定全球地理位置的技术。
GPS系统由一组卫星和地面接收器组成，可以帮助我们在地球上的任何位置确定自己的准确位置。
在C语言中，我们可以利用一些数学和物理公式来计算GPS卫星位置。
首先，我们需要了解GPS系统的原理。
GPS卫星运行在地球轨道上，同时向地球发送信号。
地面接收器接收到卫星发出的信号，并通过计算信号的距离和时间来确定自身的位置。
计算GPS卫星位置的关键是测量信号的传播时间。
当接收器接收到来自至少4个不同卫星的信号时，我们可以使用以下公式来计算GPS卫星的位置：速度=距离/时间由于信号的传播速度恒定（大约为光速），我们可以知道信号的传播时间等于距离与速度的比。
在C语言中，我们可以使用以下代码计算GPS卫星的位置：```c#include <stdio.h>#include <math.h>typedef structdouble x;double y;double z;} Point3D;Point3D calculateSatellitePosition(double distance, double latitude, double longitude, double altitude)Point3D position;position.x = (altitude + distance) * cos(E) * cos(longitude);position.y = (altitude + distance) * cos(E) * sin(longitude);return position;int maidouble latitude = 37.7749; // 纬度（假设）double longitude = -122.4194; // 经度（假设）double altitude = 0.0; // 海拔（假设）Point3D satellitePosition =calculateSatellitePosition(distance, latitude, longitude, altitude);printf("GPS卫星位置：(x=%.2f, y=%.2f, z=%.2f)\n", satellitePosition.x, satellitePosition.y, satellitePosition.z);return 0;```在这段代码中，我们首先定义了光速的常量。
利用C语言实现手持GPS航点的批量输入
（ ２ ）命名 空间 ，主 要用来 作为 一种 区别其他
Ｘ Ｍ Ｌ文 件 标 签 的 标 志 ； （ ３ ）元 数 据 信 息 ，描 述 该 文 档 所 创 建 的时 间 以 及航点坐标范围：
ｄ ａ ｔ ａ ， ｉ ｍｐ ｌ ｅ ｍｅ ｎ ｔ ｓ ｔ ｈ ｅ ｗａ ｙ ｐ ｏ ｉ ｎ ｔ ｄ ａ ｔ ａ ｂ ａ ｔ ｃ ｈ ｉ ｎ ｐ ｕ ｔ ， ａ ｎ ｄ ｉ ｍｐ ｒ ｏ ｖ ｅ ｓ
ｈ
图 １包含航点 的最小 Ｇ Ｐ Ｘ文件
在这 五个部分 中，Ｘ Ｍ Ｌ 序言 、命 名空间以及元
数据 信 息 可 以 不 用 改 动 ， 直 接 参 照模 板 也 可 以 使 ． Ｅ ｆ ｊ 。
在航 点信息部分 ，这里 只关注航 点部 分，忽略掉航 迹 、航线等 。一 个航 点信 息包括 经纬度 （ Ｗ Ｇ Ｓ 一 ８ ４坐 标） 、航 点名称 、航 点符号、创 建时 间等，在文件末 尾 附上 结 束 符 （ ／ ｇ ｐ ｘ ） 表示文件结束 。
Ｓ一 ８ ４ Ｃ Ｏ Ｏｒ ｄｉ ｎａ ｔ ｅ
等信 息 ；
（ ５ ）结束符，文件的结束标志 。
引 言
目前，手 持 Ｇ Ｐ Ｓ接收机在测绘 、物探 等勘 察工 作 中得 到了广泛的应用 ，已经成为 了野外勘察工 作 的常用工 具。手持 Ｇ Ｐ Ｓ轻巧方便 ，其定位 精度达 到
２ ～５ ｍ 的水 平 ，在 定 位 、导 航 等 方 面 深 受 野 外 勘 察
为野 外勘察 人员的常用工具。本文通过 Ｃ语言编程来处理
Ｗ Ｇ Ｓ ８ ４坐标航 点数据 ，实现了航 点数据 的批 量输 入，提高
了手 持 Ｇ Ｐ Ｓ的使用效率。
式 ）的缩写。Ｇ Ｐ Ｘ交换格 式文件 是 Ｇ Ｐ Ｓ设备 中的地 图 信 息 格 式 ，主 要 用 来 描 写 和 保 存地 理 数 据 信 息 。
c语言gps课程设计
c语言 gps课程设计一、教学目标本课程的目标是让学生掌握C语言在GPS领域的应用。
通过本课程的学习，学生将能够：1.理解GPS系统的基本原理和工作流程。
2.掌握C语言的基本语法和编程技巧。
3.能够使用C语言进行GPS数据的采集、解析和处理。
4.能够利用C语言实现简单的GPS导航功能。
二、教学内容本课程的教学内容主要包括以下几个部分：1.GPS系统的基本原理和工作流程。
2.C语言的基本语法和编程技巧。
3.GPS数据的采集、解析和处理。
4.GPS导航功能的实现。
三、教学方法为了达到本课程的教学目标，我们将采用多种教学方法，包括：1.讲授法：用于讲解GPS系统的基本原理和工作流程，以及C语言的基本语法和编程技巧。
2.案例分析法：通过分析具体的GPS应用案例，让学生了解GPS技术的实际应用。
3.实验法：通过实验让学生亲手操作，掌握GPS数据的采集、解析和处理方法，以及实现简单的导航功能。
四、教学资源为了支持本课程的教学内容和教学方法的实施，我们将准备以下教学资源：1.教材：选用合适的C语言和GPS相关教材，为学生提供理论学习的参考。
2.参考书：提供相关的参考书籍，拓展学生的知识面。
3.多媒体资料：制作课件和视频资料，直观地展示GPS系统和C语言的相关内容。
4.实验设备：准备GPS接收器、计算机等实验设备，让学生进行实际操作。
五、教学评估为了全面、客观地评估学生在C语言GPS课程中的学习成果，我们将采取以下评估方式：1.平时表现：通过学生在课堂上的参与度、提问回答、小组讨论等表现来评估其学习态度和理解程度。
2.作业：布置与课程内容相关的编程作业，评估学生对C语言编程和GPS应用的掌握情况。
3.考试：定期进行理论知识考试和编程实践考试，全面评估学生的知识掌握和实际应用能力。
4.项目报告：要求学生完成一个GPS相关的编程项目，通过项目报告评估学生的综合应用能力。
六、教学安排本课程的教学安排如下：1.教学进度：按照教材和大纲进行，确保覆盖所有重要知识点。
GPS-GPZDA数据解析C语言
这是一种使用没有在标准下预定义的特殊内容的方法。它通常的格式为：$PmmmA,df1,df2,....,[CR][LF]P说明是属性内容，mmm定义为厂商信息代码，A(A-Z）表明信息类型。
NMEA-0183协议定义的语句非常的多，主要有$GPGLL、 $GPGGA、$GPVTG、$GPZDA等。因为本次课题我们需要使用GPS-GPZDA进行数据解析所以我们只介绍GPZDA标准格式的样式。
2.3.2
查询的标准格式为：$ttllQ,sss,[CR][LF]
前面两个字符表示请求者的信息源的识别码，后面两个字符是被查询设备的识别码，最后一个字符说明这是一个查询的信息。紧跟着的字符（SSS）包含了三个字的被查询内容的记忆信息。例如： $CCGPQ,ZDA[CR][LF] “CC”说明计算机正从“GP”这个设备查询ZDA中的内容。GPS将会每隔1秒发送内容直到有新的请求。
char chhour[5]={'\0'};
char chmin[5]={'\0'};
char chsec[8]={'\0'};
char chsum[10]={'\0'};
char afterch[128]={'\0'};
int i,j,n,m; /*定义并初始化一些数组，指针，整型变量等用于后续程序使用*/
通过这次课程设计我们可以了解，如何用GPZDA解析GPS的时间信号，在计算机上输出普通用户所能理解的语句，让用户能够在时间信息这方面更加精确地得到信息，给用户提供更可靠的更高水准的服务。
第一章
GPS全球卫星定位系统由21颗工作卫星和3颗备用卫星组成。它们分布在6个等间距的轨道平面上，每个轨道面上有4颗工作卫星。GPS能覆盖全球，能够连续、实时、隐蔽地定位，一次定位时间仅几秒到十几秒，用户不发射任何电磁信号，只要接受卫星导航信号即可定位，所以可全天候昼夜作业。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
C语言计算G P S 卫星位置1 概述 在用GPS 信号进行导航定位以及制订观测计划时，都必须已知GPS 卫星在空间的瞬间位置。
卫星位置的计算是根据卫星电文所提供的轨道参数按一定的公式计算的。
本节专门讲解观测瞬间GPS 卫星在地固坐标系中坐标的计算方法。
2 卫星位置的计算1. 计算卫星运行的平均角速度n根据开普勒第三定律，卫星运行的平均角速度n0可以用下式计算：式中μ为WGS-84坐标系中的地球引力常数，且μ=3.986005×1014m 3/s 2。
平均角速度n 0加上卫星电文给出的摄动改正数Δn ，便得到卫星运行的平均角速度nn=n 0+Δn (4-12)2. 计算归化时间t k首先对观测时刻t ′作卫星钟差改正t=t ′-Δt然后对观测时刻t 归化到GPS 时系t k =t-t oc (4-13)式中t k 称作相对于参考时刻t oe 的归化时间（读者注意：toc ≠toe ）。
3. 观测时刻卫星平近点角M k 的计算M k =M 0+n tk (4-14)式中M 0是卫星电文给出的参考时刻toe 的平近点角。
4. 计算偏近点角E kE k =M k +esinE k (E k ,M k 以弧度计) （4-15）上述方程可用迭代法进行解算，即先令E k =M k ,代入上式，求出E k 再代入上式计算，因为GPS 卫星轨道的偏心率e 很小，因此收敛快，只需迭代计算两次便可求得偏近点角E k 。
5. 真近点角V k 的计算由于:因此:6．升交距角Φk 的计算ω为卫星电文给出的近地点角距。
7. 摄动改正项δu,δr,δi 的计算δu,δr,δi 分别为升交距角u 的摄动量，卫星矢径r 的摄动量和轨道倾角i 的摄动量。
8. 计算经过摄动改正的升交距角u k 、卫星矢径r k 和轨道倾角i k9. 计算卫星在轨道平面坐标系的坐标卫星在轨道平面直角坐标系（X 轴指向升交点）中的坐标为10. 观测时刻升交点经度Ωk 的计算升交点经度Ωk 等于观测时刻升交点赤经Ω（春分点和升交点之间的角距）与格林泥治视恒星时GAST （春分点和格林尼治起始子午线之间的角距）之差，Ωk =Ω-GAST （4-23）又因为:tk oe Ω+Ω=Ω (4-24)其中Ωoe 为参与时刻t oe 的升交点的赤经；Ω是升交点赤经的变化率，卫星电文每小时更新一次Ω和t oe 。
此外，卫星电文中提供了一周的开始时刻t w 的格林尼治视恒星时GAST w 。
由于地球自转作用，GAST 不断增加，所以:GAST=GAST w +ωe t （4-25）式中ωe ×10-5rad/s 为地球自转的速率；t 为观测时刻。
由式（4-24）和（4-25），得:由（4-13）式，得:其中0oe w GAST Ω=Ω-，o Ω、Ω、oe t 的值可从卫星电文中获取。
11. 计算卫星在地心固定坐标系中的直角坐标把卫星在轨道平面直角坐标系中的坐标进行旋转变换，可得出卫星在地心固定坐标系中的三维坐标：12. 卫星在协议地球坐标系中的坐标计算考虑极移的影响，卫星在协议地球坐标系中的坐标为利用C语言程序实现#include <stdio.h>#include <string.h>#include <stdlib.h>#include <math.h>#define WE 7.292115e-6struct canshu{int prn, nian, yue, ri, shi, fen;//卫星PRN号，年，月，日，时，分double miao;//秒long double adoe, a0, a1, a2, mo, dn, e, ga, pio, io, w, pid, ii, cuc, cus, cue, crs, crc, cis, cic, toe, aodc, wn;/*参数说明： ADOE值，a0 卫星钟偏差, a1 卫星钟漂移, a2 卫星钟频率漂移, M0 平近点角, Δn平运动差, e偏心率, a1/2半长轴的平方根, Ω0 轨道平面升交点经度,i0 倾角, ω近地点角距, *Ω升交点速率, IDot 倾角速率, Cuc Cus 升交角距的摄动改正项, Crc Crs 地心距的摄动改正项, Cic Cis 倾角的摄动改正项,toe 参考历元*/};void wxzbjx(struct canshu *pt){long double a, n0, n, t, tk, toc, mk, ek, vk, fik, uk, rk, ik;long double xk, yk ,zk, lk;long double XK, YK, ZK;int temp;pt->nian = pt->nian + 2000;t = (long double)(((pt ->nian)- 1980) * 365 * 24 * 3600 + (pt ->yue - 1) * 30 * 24 * 3600 + pt ->ri* 24 * 3600 + pt->shi * 3600 + pt ->miao);a = pt ->ga * pt ->ga;n0 = sqrt(WE/(a*a*a));//平均角速度n0n = n0 + pt ->dn;tk = t - pt ->toe;toc = pt ->a0 + pt ->a1 * (t - pt ->toe) + pt ->a2 * (t - pt->toe) * (t - pt ->toe);tk = tk - pt ->toe;mk = pt ->mo + n * tk;ek = mk;for(temp=0;temp<10;temp++){ek=mk+ pt ->e * sin(ek);//利用迭代法求偏近点角ek }vk = 2 * atan(sqrt((1+ pt ->e) / (1 - pt ->e))* (tan(ek)) / 2 );fik = vk + pt ->w;uk = fik + pt ->cuc * cos(2* fik) + pt ->cus * sin(2*fik);rk = pt ->ga * pt ->ga * (1 - pt ->e * ek) + pt ->crc * cos(2* fik) + pt ->crs * sin(2*fik);ik = pt ->io + pt ->cic * cos(2 * fik) + pt ->cis * sin(2* fik) + pt ->ii * tk;xk = rk * cos(uk);yk = rk * sin(uk);zk = 0;lk = pt ->pio + (pt ->pid - WE) * tk - WE * pt->toe;XK = xk * cos(lk) - yk * cos(ik) * sin(lk);YK = xk * sin(lk) + yk * cos(ik) * cos(lk);ZK = yk * sin(ik);printf("\n%d年%d月%d号%d时%.2f秒 %d号卫星的坐标:", pt->nian,pt->yue ,pt->ri , pt->shi ,pt->miao, pt->prn);printf("\nXk = %.9f\nYk= %.9f\nZK = %.9f\n\n", XK, YK, ZK);}int main(void){FILE *fp, *fp1, *fp2;struct canshu a;int i=0, hanhao = 1;long double t emp1, temp2, temp3, temp5, temp4, temp6, temp7;char ch, ch1;if((fp1 = fopen("E:\\星历文件\\guangboxingli2.txt", "r")) == NULL)//请自定义星历文件位置及名称{printf("文件无法打开！");exit(0);}else{if((fp2 = fopen("E:\\星历文件\\guangboxingli2fu.txt", "w")) == NULL){printf("文件无法打开！");exit(0);}else{while((ch1 = fgetc(fp1)) != EOF){if(ch1 == '\n'){i ++;}putchar(ch1);if(i == 15){break;}}while(!feof(fp1)){ch1=fgetc(fp1);if(ch1 == 'D'){ch1 = 'e';}fputc(ch1,fp2);}fclose(fp1);}fclose(fp2);}printf("以上是星历文件的头文件！\n");system("pause");printf("读取文件参数数据\n！");if((fp = fopen("E:\\星历文件\\guangboxingli2fu.txt", "r")) == NULL)//创建计算结果文档{printf("文件无法打开！");exit(0);}while(!feof(fp)){switch(hanhao){case 1:{fscanf(fp,"\n%d%d%d%d%d%d%lf %le %le %le", &a.prn, &a.nian, &a.yue, &a.ri, &a.shi,&a.fen, &a.miao, &a.a0, &a.a1, &a.a2);printf("%d %d %d %d %d %d %lf %le %le %le", a.prn, a.nian, a.yue, a.ri, a.shi,a.fen, a.miao, a.a0, a.a1, a.a2);hanhao++;}case 2:{fscanf(fp,"%le %le %le %le", &a.adoe, &a.crs,&a.dn, &a.mo);printf("\n%le %le %le %le", a.adoe, a.crs, a.dn, a.mo);hanhao++;}case 3:{fscanf(fp,"%le %le %le %le", &a.cue, &a.e, &a.cus, &a.ga);printf("\n%le %le %le %le", a.cue, a.e, a.cus,a.ga);hanhao++;}case 4:{fscanf(fp,"%le %le %le %le", &a.toe, &a.cic,&a.pio, &a.cis);printf("\n%le %le %le %le", a.toe, a.cic, a.pio, a.cis);hanhao++;}case 5:{fscanf(fp,"%le %le %le %le", &a.io, &a.crc, &a.w, &a.pid);printf("\n%le %le %le %le", a.io, a.crc, a.w,a.pid);hanhao++;}case 6:{fscanf(fp,"%le %le %le %le", &a.ii, &temp1, &a.wn, &temp2);printf("\n%le %le %le %le", a.ii, temp1, a.wn, temp2);hanhao++;}case 7:{fscanf(fp,"%le %le %le %le", &temp3, &temp4,&temp5, &a.aodc);printf("\n%le %le %le %le", temp3, temp1, temp5, a.aodc);hanhao++;}case 8:{fscanf(fp,"%le", &temp6);printf("\n%le", temp6);if((ch=fgetc(fp)) != '\n'){fscanf(fp,"%le", &temp7);printf(" %le\n", temp7);}}}wxzbjx(&a);system("pause");hanhao = 1;}fclose(fp);return 0;}。