-
大小: 487KB文件類型: .zip金幣: 2下載: 0 次發(fā)布日期: 2023-12-01
- 語言: C/C++
- 標(biāo)簽: C語言??GPS網(wǎng)??間接平差??測(cè)繪??
資源簡(jiǎn)介
GPS網(wǎng)間接平差C語言程序，采用文件讀取方式，處理成果存入txt文件中
代碼片段和文件信息
#include?
#include?
#include?
#define?M?100
#define?DS?10
#define?GS?25
int?fg[M];
double?Jx[M]Jy[M]Jz[M]X[M]Y[M]Z[M]B[M][M]BT[M][M]BTP[M][M]BTPB[M][M]BTPBN[M][M]l[M]x[M]P[M][M]BTPL[M]xx[M]yy[M]zz[M]v[M];
double?danweiquanzhongwuchaxzwc[M]yzwc[M]zzwc[M]dzwc[M];
struct
{
int?qidian;
int?zhongdian;
double?dx;
double?dy;
double?dz;
}GPSjxxl[M];
void?shuzufuchushizhi（）
{
int?ij;
for?（i=0;i
for?（j=0;j
fg[i]=0;
B[i][j]=P[i][j]=0.0;
}
}
}
//數(shù)據(jù)讀取
void?shujuduqu（）
{
int?iabno;
double?t1t2t3;
FILE?*fp;
fp=fopen（“GPS網(wǎng)觀測(cè)數(shù)據(jù)2.txt““r“）;
if?（fp==NULL）
{
printf（“打開文件失敗!\n“）;
exit（0）;
}
//讀取基線向量
for?（i=1;i<=GS;i++）
{
fscanf（fp“%d%d%d%lf%lf%lf“&no&a&b&t1&t2&t3）;
GPSjxxl[i].qidian=a;
GPSjxxl[i].zhongdian=b;
GPSjxxl[i].dx=t1;
GPSjxxl[i].dy=t2;
GPSjxxl[i].dz=t3;
}
//讀取已知點(diǎn)坐標(biāo)
for?（i=1;i<2;i++）
{
fscanf（fp“%d%lf%lf%lf“&a&t1&t2&t3）;
Jx[a]=X[a]=t1;
Jy[a]=Y[a]=t2;
Jz[a]=Z[a]=t3;
fg[a]=1;
}
fclose（fp）;
}
//計(jì)算坐標(biāo)近似值
void?cal_zbjinsizhi（）
{
int?iqdzd;
double?dxdydz;
for?（i=1;i<=GS;i++）
{
qd=GPSjxxl[i].qidian;
zd=GPSjxxl[i].zhongdian;
dx=GPSjxxl[i].dx;
dy=GPSjxxl[i].dy;
dz=GPSjxxl[i].dz;
if?（fg[qd]!=0&&fg[zd]==0）
{
Jx[zd]=Jx[qd]+dx;
Jy[zd]=Jy[qd]+dy;
Jz[zd]=Jz[qd]+dz;
fg[zd]=1;
}
if?（fg[qd]==0&&fg[zd]!=0）
{
Jx[qd]=Jx[zd]-dx;
Jy[qd]=Jy[zd]-dy;
Jz[qd]=Jz[zd]-dz;
fg[qd]=1;
}
}
}
//計(jì)算B矩陣
void?cal_xishu（）
{
int?iqdzdx1x2;
for?（i=1;i<=GS;i++）
{
x1=x2=1;
qd=GPSjxxl[i].qidian;
zd=GPSjxxl[i].zhongdian;
if?（qd==1）
{
x1=0;
}
if?（zd==1）
{
x2=0;
}
B[3*i-2][3*qd-2]=B[3*i-1][3*qd-1]=B[3*i][3*qd]=-x1;
B[3*i-2][3*zd-2]=B[3*i-1][3*zd-1]=B[3*i][3*zd]=x2;
}
}
//計(jì)算常數(shù)矩陣
void?cal_l（）
{
int?iqdzd;
double?dxdydz;
for?（i=1;i<=GS;i++）
{
qd=GPSjxxl[i].qidian;
zd=GPSjxxl[i].zhongdian;
dx=GPSjxxl[i].dx;
dy=GPSjxxl[i].dy;
dz=GPSjxxl[i].dz;
l[3*i-2]=dx-（Jx[zd]-Jx[qd]）;
l[3*i-1]=dy-（Jy[zd]-Jy[qd]）;
l[3*i]=dz-（Jz[zd]-Jz[qd]）;
}
}
//轉(zhuǎn)置矩陣
void?zhuanzhi（）
{
int?ij;
for?（i=1;i<=3*GS;i++）
{
for（j=1;j<=3*DS;j++）
{
BT[j][i]=B[i][j];
}
}
}
//計(jì)算權(quán)陣
void?cal_quanzhen（）
{
int?i;
double?sdxdydzt;
for?（i=1;i<=GS;i++）
{
dx=GPSjxxl[i].dx;
dy=GPSjxxl[i].dy;
dz=GPSjxxl[i].dz;
s=sqrt（（dx*dx）+（dy*dy）+（dz*dz））;
t=sqrt（2+2*s/1000）;
P[3*i-2][3*i-2]=P[3*i-1][3*i-1]=P[3*i][3*i]=t;
}
}
//矩陣乘法
void?juzhenxiangcheng（）??
{
int?ijk;
double?t;
//計(jì)算BTP
for（i=4;i<=3*DS;i++）
{
for（j=1;j<=3*GS;j++）
{
t=0;
for（k=1;k<=3*GS;k++）
{
t+=BT[i][k]*P[k][j];
}
BTP[i][j]=t;
}
}
//計(jì)算BTPB
for（i=4;i<=3*DS;i++）
{
for（j=4;j<=3*DS;j++）
{
t=0;
for（k=1;k<=3*GS;k++）
{
t+=BTP[i][k]*B[k][j];
}
BTPB[i][j]=t;
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2013-04-27?21:00??GPS網(wǎng)\
?????目錄???????????0??2013-04-27?21:00??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\
?????目錄???????????0??2013-04-27?21:00??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\Debug\
?????文件??????229479??2013-01-06?15:43??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\Debug\GPS網(wǎng).exe
?????文件??????234200??2013-01-06?15:43??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\Debug\GPS網(wǎng).ilk
?????文件???????33104??2013-01-06?15:43??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\Debug\GPS網(wǎng).obj
?????文件??????190452??2013-01-06?15:43??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\Debug\GPS網(wǎng).pch
?????文件??????418816??2013-01-06?15:43??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\Debug\GPS網(wǎng).pdb
?????文件???????33792??2013-01-06?15:43??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\Debug\vc60.idb
?????文件???????45056??2013-01-06?15:43??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\Debug\vc60.pdb
?????文件????????9575??2013-01-06?15:44??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\GPS網(wǎng).c
?????文件????????3389??2013-01-06?15:43??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\GPS網(wǎng).dsp
?????文件?????????518??2013-01-06?15:44??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\GPS網(wǎng).dsw
?????文件???????33792??2013-01-06?15:44??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\GPS網(wǎng).ncb
?????文件???????48640??2013-01-06?15:44??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\GPS網(wǎng).opt
?????文件?????????745??2013-01-06?15:43??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\GPS網(wǎng).plg
?????文件??????106744??2013-01-06?15:43??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\GPS網(wǎng)平差數(shù)據(jù)處理成果.txt
?????文件?????????827??2012-12-26?21:36??GPS網(wǎng)\GPS網(wǎng)（校內(nèi)）\GPS網(wǎng)觀測(cè)數(shù)據(jù)2.txt
?????目錄???????????0??2013-04-27?21:00??GPS網(wǎng)\GPS網(wǎng)（校外）\
?????目錄???????????0??2013-04-27?21:00??GPS網(wǎng)\GPS網(wǎng)（校外）\Debug\
?????文件??????229479??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\Debug\GPS網(wǎng).exe
?????文件??????234200??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\Debug\GPS網(wǎng).ilk
?????文件???????33008??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\Debug\GPS網(wǎng).obj
?????文件??????190452??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\Debug\GPS網(wǎng).pch
?????文件??????418816??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\Debug\GPS網(wǎng).pdb
?????文件???????33792??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\Debug\vc60.idb
?????文件???????45056??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\Debug\vc60.pdb
?????文件????????9613??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\GPS網(wǎng).c
?????文件????????3389??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\GPS網(wǎng).dsp
?????文件?????????518??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\GPS網(wǎng).dsw
?????文件???????33792??2013-01-06?15:46??GPS網(wǎng)\GPS網(wǎng)（校外）\GPS網(wǎng).ncb
............此處省略4個(gè)文件信息
評(píng)論
共有 條評(píng)論
川公網(wǎng)安備 51152502000135號(hào)