当前位置:网站首页>GPS高程拟合抗差中误差的求取代码实现
GPS高程拟合抗差中误差的求取代码实现
2022-07-06 09:18:00 【Proletarians】
话不多说,直接上代码了,c代码
// 2.1、使用中位数法,而非算术中位法
qsort(tmp_v_res, ns, sizeof(double), cmpres);
median_v_res =median(tmp_v_res, ns);// 观测值中数 Eq.14
// 2.1.1 求取偏离量,中心化
for (j = 0; j < ns; j++) {
tmp_v_res[j] -= median_v_res;// Eq.15
/*v_res[j] -= median_v_res;*/
if (tmp_v_res[j] < 0) tmp_v_res[j] *= -1;//取绝对值
}
// 2.1.2 使用中位数计算方差因子
qsort(tmp_v_res, ns, sizeof(double), cmpres);
median_v_res = median(tmp_v_res, ns);//观测值中数
sigma1 = median_v_res / 0.6745;// Eq.17
printf("抗差单位权中误差(weighted)sigma1= %f\n", sigma1);
/* 残差比较 */
static int cmpres(const void *p1, const void *p2){
double *q1 = (double *)p1, *q2 = (double *)p2;
double tt = (*q1) - (*q2);
return tt<-0.0 ? -1 : (tt>0.0 ? 1 : 0);
}
/* 求中位数 */
static double median(const double* A,unsigned int n){
if (n>0){
return (n & 1 == 1) ? A[(n - 1) / 2] : (A[n / 2] + A[n / 2 - 1]) / 2;
}
return 0;
}
边栏推荐
猜你喜欢
[Red Treasure Book Notes simplified version] Chapter 12 BOM
ESP8266连接onenet(旧版MQTT方式)
Unity3D基础入门之粒子系统(属性介绍+火焰粒子系统案例制作)
NovAtel 板卡OEM617D配置步骤记录
(3) Introduction to bioinformatics of R language - function, data Frame, simple DNA reading and analysis
MySQL時間、時區、自動填充0的問題
Unity3D,阿里云服务器,平台配置
Redis based distributed ID generator
PR 2021 quick start tutorial, first understanding the Premiere Pro working interface
Fabrication of fairygui simple Backpack
随机推荐
编译原理:源程序的预处理及词法分析程序的设计与实现(含代码)
FairyGUI循环列表
FairyGUI条子家族(滚动条,滑动条,进度条)
MySQL takes up too much memory solution
FairyGUI简单背包的制作
First use of dosbox
[offer78]合并多个有序链表
Pytorch: tensor operation (I) contiguous
Vulnhub target: hacknos_ PLAYER V1.1
MySQL time, time zone, auto fill 0
Guided package method in idea
Unity场景跳转及退出
Mysqldump error1066 error solution
Page performance optimization of video scene
Unity3d camera, the keyboard controls the front and rear left and right up and down movement, and the mouse controls the rotation, zoom in and out
Unity3D基础入门之粒子系统(属性介绍+火焰粒子系统案例制作)
[offer18] delete the node of the linked list
(课设第一套)1-4 消息传递接口 (100 分)(模拟:线程)
[leetcode15] sum of three numbers
Postman 中级使用教程【环境变量、测试脚本、断言、接口文档等】