C#:测绘程序设计的条件平差模型

于 2024-09-07 14:48:42 发布
阅读量148 收藏
点赞数 3
文章标签： c#
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_58362929/article/details/120676786
版权

以此题为例：

 写出条件方程，确定系数矩阵A，矩阵W，矩阵Q，L矩阵。

说明：以下程序代码运行的数据，数据有单位的均按单位mm处理。

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;

namespace _10月8号_条件平差
{
    class Program
    {
    
        //输入函数
         static void  InPut(ref int m,ref int n, ref double [,]A )
        {
            Console.WriteLine("输入矩阵的行数列数");
            string r = Console.ReadLine();
            string [] rs = r.Split(' ');
            m = int.Parse(rs[0]);
            n = int.Parse(rs[1]);
            A = new double[m, n];
            Console.WriteLine("矩阵按行输入");
            Console.WriteLine();
            for (int i=0;i<m;i++)
            {
                string str = Console.ReadLine();
                string[] str1 = str.Split(' ');
                for (int j=0;j<n;j++)
                {
                    A[i, j] = double.Parse(str1[j]);
                }
            }
        }

        //输出函数
        static void OutPut(ref double [,] M)
        {
            int m = M.GetLength(0);
            int n = M.GetLength(1);
            for (int i=0;i<m;i++)
            {
                for (int j=0;j<n;j++)
                {
                    Console.Write("{0,12:f4}", M[i, j]);
                }
                Console.WriteLine();
            }
        }
        //输出函数2
        static void OutPut1(ref double [,] M1)
        {
            int m = M1.GetLength(0);
            int n = M1.GetLength(1);
            for (int i = 0; i < m; i++)
            {
                for (int j = 0; j < n; j++)
                {
                    Console.Write("h{0}={1,12:f4}      ",i+1, M1[i, j]);
                }
                Console.WriteLine();
            }
        }
        //矩阵转置函数
        static void Transpose_Matrix (ref double [,] ab,ref double [,]ba )
        {
            int p = ab.GetLength(0);
            int q = ab.GetLength(1);
            ba = new double [q,p];
            for (int i=0;i<q;i++)
            {
                for (int j = 0; j < p; j++)
                    ba[i, j] = ab[j, i];
            }
        }

        //乘法矩阵

        static void Multiplication_Matrix(ref double[,] a, ref double[,] b, ref double[,] ab2)
        {
            ab2 = new double[a.GetLength(0), b.GetLength(1)];
            for (int i = 0; i < a.GetLength(0); i++)
            {
                for (int j = 0; j < b.GetLength(1); j++)
                {
                    for (int k = 0; k < a.GetLength(1); k++)
                    {
                        ab2[i, j] += a[i, k] * b[k, j];
                    }
                }
            }
        }
 
       
        //逆矩阵函数
        static void Inverse(ref double[,] Array, ref double[,] NI)
        {
            int m = Array.GetLength(0);
            int n = Array.GetLength(1);
            double[,] array = new double[2 * m + 1, 2 * n + 1];
            for (int k = 0; k < 2 * m + 1; k++)  //初始化数组
            {
                for (int t = 0; t < 2 * n + 1; t++)
                {
                    array[k, t] = 0.00000000;
                }
            }
            for (int i = 0; i < m; i++)
            {
                for (int j = 0; j < n; j++)
                {
                    array[i, j] = Array[i, j];
                }
            }
            for (int k = 0; k < m; k++)
            {
                for (int t = n; t <= 2 * n; t++)
                {
                    if ((t - k) == m)
                    {
                        array[k, t] = 1.0;
                    }
                    else
                    {
                        array[k, t] = 0;
                    }
                }
            }
            //得到逆矩阵
            for (int k = 0; k < m; k++)
            {
                if (array[k, k] != 1)
                {
                    double bs = array[k, k];
                    array[k, k] = 1;
                    for (int p = k + 1; p < 2 * n; p++)
                    {
                        array[k, p] /= bs;
                    }
                }
                for (int q = 0; q < m; q++)
                {
                    if (q != k)
                    {
                        double bs = array[q, k];
                        for (int p = 0; p < 2 * n; p++)
                        {
                            array[q, p] -= bs * array[k, p];
                        }
                    }
                    else
                    {
                        continue;
                    }
                }
            }
            NI = new double[m, n];
            for (int x = 0; x < m; x++)
            {
                for (int y = n; y < 2 * n; y++)
                {
                    NI[x, y - n] = array[x, y];
                }
            }
        }
        static void Matrix_plus(ref double [,] a1,ref double [,] a2,ref double [,] plus)
        {
            int m1 = a1.GetLength(0);
            int n1 = a1.GetLength(1);
            int m2 = a2.GetLength(0);
            int n2 = a2.GetLength(1);
            plus = new double[m1, n1];
            if (m1 == m2 & n1 == n2)
            {
                for (int i = 0; i < m1; i++)
                {
                    for (int j = 0; j < n1; j++)
                        plus[i, j] = a1[i, j] + a2[i, j];
                }
            }
            else
                Console.WriteLine("两矩阵不能相加");
        }
        static void Main(string[] args)
        {
            Console.WriteLine("------------输入系数矩阵A[i,j]------------");
            int m = 0, n = 0;
            double [,] A={};
            double[,] M = { };
            InPut(ref m,  ref n,  ref A);
            double[,] AT = { };
            double[,] ab = { };
            double[,] ba = { };
            double[,] A1 = { };
            A1 = A;//系数矩阵
            ab = A;
            Transpose_Matrix(ref ab, ref ba);
            AT = ba;//系数矩阵的转置
            Console.WriteLine("------------矩阵A的转置为------------");
            Console.WriteLine();
            M = AT;
            OutPut(ref M);
            Console.WriteLine("------------输入协因数矩阵Q[i,j]------------");
            Console.WriteLine();
            InPut(ref m,ref n,ref A);
            double[,] Q = { };
            Q = A;//协因数矩阵
            Console.WriteLine("------------Q[i,j]------------");
            Console.WriteLine();
            M=Q;
            OutPut(ref M);
            Console .WriteLine("------------输入W矩阵------------");
            Console .WriteLine();
            double [,] W={};
            InPut(ref m,ref n ,ref A);
            Console .WriteLine("------------W[i,j]------------");
            Console .WriteLine();
            W=A;
            M=W;
            OutPut(ref M);
            Console.WriteLine("------------按单位为毫米输入L[i,j]矩阵------------");
            Console .WriteLine();
            InPut (ref m, ref n, ref A);
            double[,] L = { };
            L = A;//l矩阵
            Console .WriteLine("------------L[i,j]------------");
            Console .WriteLine();
            M=L;
            OutPut(ref M);

            //计算A1Q
            double[,] a = { };
            double[,] b = { };
            double[,] ab2 = { };
            a = A1;
            b = Q;
            Multiplication_Matrix(ref a, ref b, ref ab2);
            double[,] A1Q = { };
            A1Q = ab2;
            double[,] AQAT = { };
            a = A1Q;
            b = AT;
            Multiplication_Matrix(ref a, ref b, ref ab2);
            AQAT = ab2;  //Naa矩阵
            Console.WriteLine("------------Naa[i,j]------------");
            Console.WriteLine();
            M = AQAT;
            OutPut(ref M);

            //计算-K矩阵.............................由于乘以-1过于繁琐，所以把直接计算-K简化过程
            double[,] Arry = { };
            double[,] NI = { };
            Arry = AQAT;
            Inverse(ref Arry,ref NI);
            a=NI;
            b=W;
            Multiplication_Matrix(ref a, ref b, ref ab2);
            double[,] K = { };
            K = ab2;//-K矩阵
            Console.WriteLine("------------/-K[i,j]/------------");
            Console.WriteLine();
            M=K;
            OutPut(ref M);

            //计算-V矩阵
            double[,] QAT = { };
            a = Q;
            b = AT;
            Multiplication_Matrix(ref a, ref b, ref ab2);
            QAT = ab2;
            double[,] v = { };
            a = QAT;
            b = K;
            Multiplication_Matrix(ref a, ref b, ref ab2);
            v = ab2;//-V[i,j]矩阵
            Console.WriteLine("------------/-v[i,j]/------------");
            Console.WriteLine();
            M = v;
            OutPut(ref M);

            //计算最终改正后的高差
            double [,] a1 = { };
            double [,] a2 = { };
            double [,] plus = { };
            double[,] V = { };
            a1 = L;
            a2 = v;
            Matrix_plus(ref a1, ref a2, ref plus);
            V = plus;
            Console.WriteLine("------------高差的改正后矩阵------------");
            Console.WriteLine();
            double [,] M1={};
            M1=V;
            OutPut1(ref M1);
        }
    }
}

运行结果如下：

仅供参考，如有错误多多谅解，敬请斧正。 

十万毫升.
关注
测绘程序设计——坐标正反算(C#)
m0_58301791的博客
01-10 2405
若已知A点坐标（200，200）B点坐标（250，250），Xab=50、Yab=50，则tan a=1；若已知A点坐标（200，200）AB两点的距离100，AB两点的坐标方位角为45°；坐标方位角正算是由A点的坐标AB两点间的距离和方位角，来推算B点的坐标。坐标方位角反算是由两点的坐标，来推算AB两点的坐标方位角。若Xab>0,Yab0，Yab>0所以Aab=45°。若Xab
平差课程设计报告（水准网间接平差程序和代码）
12-30
测绘工程专业的平差课程报告，是关于间接平差水准网的，，，里面有程序和报告。。。。
C#水准网间接平差.zip
04-26
1、根据观测值程序自动列B矩阵； 2、自动计算参数估计值 3、程序实现自动定权生成P矩阵； 4、计算平差值及精度评价
间接平差计算代码
07-07
里面包括了间接平差设计步骤以及相关代码，矩阵计算代码
C#:测绘程序设计的间接平差模型
最新发布
m0_58362929的博客
09-07 202
以此题为例 先写出误差方程，确定系数矩阵B[i,j]，l[i,j],P[i,j],带入程序运算 using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.Threading.Tasks; namespace 间接平差 { class Program { //输入函数 static void InPut(re.
利用C#编写一个水准测量近似平差程序
Zj1638的博客
08-01 6648
利用 C#编写一个水准测量近似平差程序
水准网条件平差程序设计与实现.zip
04-10
《水准网条件平差程序设计与实现》是一个关于利用C#编程语言开发测量数据处理软件的项目。在本文中，我们将深入探讨这个项目的背景、核心概念以及实现过程。 水准网条件平差是测绘学中一个重要的计算方法，主要用于...
平差计算C代码
12-16
3环14节点16段环状管网平差计算C程序，课程联系，可供参考
间接平差编程
11-15
平差解算，完整代码。利用vc编程完成间接平差解算！
一些常用的测量程序（c#）
03-31
交会，坐标正反算，坐标转换，三角高程平差，纵断面计算
工程测绘程序
02-14
工程测量5800正反算及隧道超欠挖最终版程序
基于C#的导线平差程序设计
05-05
为解决NET平台下导线数据处理问题,文中以C#为开发语言,使用近似平差模型,开发导线平差程序。使用文献[9]附合导线观测数据,将程序处理结果和文献处理结果对比表明:1程序计算结果中方位角闭合差、方位角闭合差限差、角度改正数、导线全长相对闭合差均与原表一致;2、3、4号点计算坐标或与原表一致或相差0.001 m。2程序中针对平差解算引入并拓展了四舍五入算法,将弧度保留八位有效数字,较好地解决了平差过程中由于数据迭代运算导致的数据丢失问题。3以.NET平台为基础开发的平差程序可以直接与.NET平台下GIS软件的二次开发交互;程序实现了数据录入、平差解算、详细计算过程及计算结果输出等功能。
ControlNetAdjustment_导线网平差_C++_测绘程序_
10-01
"ControlNetAdjustment_导线网平差_C++_测绘程序_"这个项目是利用C++编程语言实现的一种针对导线网进行平差计算的软件工具。导线网平差通常应用于大地测量、工程测量和地籍测量等，它涉及到多个控制点间的观测数据...
c c#附和导线平差程序设计实习报告
10-20
《C#实现附和导线平差程序设计实习报告》 在现代测绘技术的快速发展下，测量数据处理的方法日益多样化，其中附和导线平差是常见的一种。本报告主要探讨如何利用C#编程语言来设计一个附和导线近似平差程序，满足高效...
【原创】C#写的水准网平差程序
weixin_30364147的博客
08-13 3924
Matrix类详见http://hi.baidu.com/aragon2/blog/item/1edbba124ab4beecc2ce796e.html?timeStamp=1311046817866， 一个基本的水准平差控制台程序。 1: using System; 2: using System.Collections.Generic; 3:...
C#详解：2023年附和导线平差算法与应用
3. 建立观测方程：基于观测数据和特定的平差条件（如定向角、距离和坐标），构建数学模型，形成平差方程。 4. 求解平差方程：运用数值方法，如高斯消元法或最小二乘法，解决这些方程，找到平差的最佳解。 5. 计算...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值