C++ 单一附合导线平差程序

于 2024-02-01 14:51:35 发布
阅读量871 收藏 18
点赞数 9
文章标签： c++ 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/learning_well/article/details/135972132
版权

一、以下图附合导线为例，图形如下：

二、第一步，读取测量数据，读取界面设计如下：

 读取数据文本文件格式如下：

（1）已知点坐标数据格式：

（2）角度观测值数据格式：

（3）边长观测值数据格式：

 读取文本文件的代码如下：

void TraverseAdjustment::onReadDataKnown()
{
    //打开文件对话框
    QString fileName = QFileDialog::getOpenFileName(this, tr("打开已知点数据"));
    QFile file(fileName);
    if (!file.open(QIODevice::ReadOnly | QIODevice::Text))
        QMessageBox::StandardButton btnValue = QMessageBox::information(this, tr("提示"), tr("打开失败！"));

    //读取文本文件数据
    QTextStream stream(&file);
    pointKnown pk1;
    while (!stream.atEnd())
    {
        QString str = stream.readLine();
        QStringList list=str.split(",");
        pk1.number = list.at(0);
        pk1.x = list.at(1).toDouble();
        pk1.y = list.at(2).toDouble();
        ui.textEdit->append(str);
        vpKnown.push_back(pk1);
    }

    ui.lineEdit->setText(fileName);
}

说明：需要先创建一个字符串变量用于存储文件对话框返回的路径，然后创建一个文件操作对象，用返回路径进行初始化。使用open()函数打开文本文件，随后设置一个读写流对象，用于逐行读取文本，将单行文本分割后，赋值给创建的结构体对象，然后添加到vecctor中。

三、第二步，计算角度闭合差

观测值为左角时，有

{\alpha _{CD}}'=\alpha _{AB}-n*180+\sum\beta _{i}

f_{\beta }={\alpha }'_{CD}-\alpha _{CD}

其中n为导线转折角数，这里如图1的附合导线，n=8.

注意1：观测角度值之和应为各个观测角度（DMS度分秒格式）的十进制度数的总和，因此，需要将原始观测值逐一转换成十进制度，再求和。（原始角度观测数据为度分秒格式的） 

注意2：\alpha _{CD}等方位角的计算，由于atan()等三角函数返回的是弧度值，因此需要将结果转化成十进制度值再输出。利用坐标反算方法，得到已知边的方位角，如下：

\alpha _{AB}=atan(\left | y_{A} -y_{B}\right |/\left | x_{A}-x_{B} \right |)

（1）方位角计算代码如下：

double TraverseAdjustment::CoordinateInverse(double xa, double ya, double xb, double yb)
{
    double arfa_ab = atan(abs(yb - ya) / abs(xb - xa));//返回的是弧度值
    arfa_ab = Radian2Angle(arfa_ab);
    if (ya<yb&&xa>xb)
    {
        arfa_ab=180 - arfa_ab;
    }

    if (ya > yb && xa > xb)
    {
        arfa_ab = 180 + arfa_ab;
    }

    if (ya > yb && xa < xb)
    {
        arfa_ab = 360 - arfa_ab;
    }
    return arfa_ab;//返回十进制度数
}

 （2）闭合差计算代码如下：

void TraverseAdjustment::onAngleClosureError()
{
    double arfa_1=0;//角度观测值之和
    for (auto& val : vA)
    {
        double a = val.value;
        a = DMS2Angle(a);
        arfa_1 = arfa_1 + a;
    }

    qDebug() << arfa_1;
    qDebug() << DMS2Angle(1644.2828);

    double arfa_ab;//AB方位角值
    double arfa_cd;//CD方位角值

    arfa_ab = CoordinateInverse(vpKnown.at(0).x, vpKnown.at(0).y, vpKnown.at(1).x, vpKnown.at(1).y);
    arfa_cd = CoordinateInverse(vpKnown.at(2).x, vpKnown.at(2).y, vpKnown.at(3).x, vpKnown.at(3).y);
    qDebug() << arfa_ab;
    qDebug() << arfa_cd;
    qDebug() << DMS2Angle(303.2728);
    qDebug() << DMS2Angle(147.5628);
    Rab = arfa_ab;
    Rcd = arfa_cd;

    n = vA.size();
    double arfa_cd1 = 0;//理论角度值

    arfa_cd1 = arfa_ab - n * 180+arfa_1;

    if (arfa_cd1 > 360)
    {
        int t1 = arfa_cd1 / 360;
        arfa_cd1 = arfa_cd1 - t1 * 360;
    }

    double f_beta = arfa_cd1-arfa_cd;
    f_beta = f_beta * 3600;
    f = f_beta;

    ui.textEdit_2->append("角度闭合差：");
    ui.textEdit_2->append(QString::number(f_beta));
    ui.textEdit_2->append("限差：");
    ui.textEdit_2->append(QString::number(24*sqrt(n)));
    if (f_beta <= 24 * sqrt(n))
    {
        ui.textEdit_2->append("合格");
    }
    else
    {
        ui.textEdit_2->append("不合格");
    }
}

四、第三步，计算方位角

\hat{\beta _{i}}=\beta _{i}+v_{i} 

若角度闭合差附合限差，则将角度闭合差均分至各角度观测值（即求改正数），然后使用改正后的角度观测值计算各测站的方位角信息。观测值为左角时，有如下：

\alpha _{i}=\alpha _{i-1}+\beta _{i}-180

以此推算各测站方位角。

（1）方位角计算代码如下：

void TraverseAdjustment::CalculateAzimuth()
{
    //计算角度观测值改正数
    double t1 = -f / n;
    t1 = t1 / 3600;
    Angle a1;
    for (auto& val : vA)
    {
        a1.value = DMS2Angle(val.value)+t1;
        a1.number = val.number;
        a1.number1 = val.number1;
        a1.number2 = val.number2;
        va.push_back(a1);
    }

    ui.textEdit_2->clear();
    ui.textEdit_2->append("方位角结果：");
    Angle a2; double R1 = Rab;
    for (int i = 0; i < va.size(); i++)
    {
        double R = R1 + va.at(i).value-180;
        if (R > 360)
        {
            int t2 = R / 360;
            R = R - t2 * 360;
        }
        a2.number = va.at(i).number;
        a2.number1 = va.at(i).number1;
        a2.number2 = va.at(i).number2;
        a2.value = R;
        ui.textEdit_2->append(QString::number(R));
        vAzi.push_back(a2);
        R1 = R;
    }
}

五、第四步，计算坐标增量

 注意：这里的cos()函数等三角函数的输入值为弧度，应把方位角转化成护弧度在计算，否则结果错误。

（1）计算坐标增量的代码如下：

void TraverseAdjustment::onCalIncrement()
{
    Side s1; ui.textEdit_2->clear();
    ui.textEdit_2->append("X坐标增量：");
    for (int i = 0; i < vAzi.size(); i++)
    {
        s1.number1 = vAzi.at(i).number;
        s1.number2 = vAzi.at(i).number2;
        for (int j = 0; j < vS.size(); j++)
        {
            if (vAzi.at(i).number == vS.at(j).number1 && vAzi.at(i).number2 == vS.at(j).number2)
            {
                double b1 = vAzi.at(i).value;
                b1 = Angle2Radian(b1);
                s1.value = vS.at(j).value * cos(b1);//cos函数输入值需为弧度
                ui.textEdit_2->append(QString::number(s1.value));
                vXIn.push_back(s1);
                s1.value = vS.at(j).value * sin(b1);//sin函数输入值需为弧度
                vYIn.push_back(s1);
                break;
            }
        }
    }

    ui.textEdit_2->append("Y坐标增量：");
    for (auto& val : vYIn)
    {
        ui.textEdit_2->append(QString::number(val.value));
    }
}

六、第五步，计算坐标改正数

 这里只需按照上述公式计算即可。

（1）计算坐标改正数代码如下：

void TraverseAdjustment::onCalCorrection()
{
    ui.textEdit_2->clear();
    ui.textEdit_2->append("X增量改正数：");
    //计算坐标增量之和
    double detax = 0;
    double detay = 0;

    for (auto& val : vXIn)
    {
        double b1 = val.value;
        detax = detax + b1;
    }

    for (auto& val : vYIn)
    {
        double b1 = val.value;
        detay = detay + b1;
    }

    //计算坐标增量闭合差
    fx = detax - (vpKnown.at(2).x - vpKnown.at(1).x);
    fy = detay - (vpKnown.at(2).y - vpKnown.at(1).y);

    //计算坐标改正数
    double s1 = 0;//导线边距离之和
    for (auto& val : vS)
    {
        s1 = s1 + val.value;
    }

    Side s2;
    for (int i = 0; i<vXIn.size(); i++)
    {
        double d;
        for (int j=0; j < vS.size(); j++)
        {
            if (vXIn.at(i).number1 == vS.at(j).number1 && vXIn.at(i).number2 == vS.at(j).number2)
            {
                d = vS.at(j).value;
            }
        }
        s2.number1 = vXIn.at(i).number1;
        s2.number2 = vXIn.at(i).number2;
        s2.value = -fx / s1 * d;
        ui.textEdit_2->append(QString::number(s2.value));
        vX.push_back(s2);
    }

    ui.textEdit_2->append("Y增量改正数：");
    for (int i = 0; i < vYIn.size(); i++)
    {
        double d;
        for (int j = 0; j < vS.size(); j++)
        {
            if (vYIn.at(i).number1 == vS.at(j).number1 && vYIn.at(i).number2 == vS.at(j).number2)
            {
                d = vS.at(j).value;
            }
        }
        s2.number1 = vXIn.at(i).number1;
        s2.number2 = vXIn.at(i).number2;
        s2.value = -fy / s1 * d;
        ui.textEdit_2->append(QString::number(s2.value));
        vY.push_back(s2);
    }
}

七、第六步，计算未知点坐标值

按公式计算即可。

（1）计算坐标值代码如下：

void TraverseAdjustment::onCalCoordinate()
{
    ui.textEdit_2->clear();
    ui.textEdit_2->append("改正后X增量：");
    double x0 = vpKnown.at(1).x;
    double y0 = vpKnown.at(1).y;

    for (int i=0;i<vXIn.size();i++)
    {
        vXIn.at(i).value = vXIn.at(i).value + vX.at(i).value;
        ui.textEdit_2->append(QString::number(vXIn.at(i).value));
    }

    ui.textEdit_2->append("改正后Y增量：");
    for (int i = 0; i < vYIn.size(); i++)
    {
        vYIn.at(i).value = vYIn.at(i).value + vY.at(i).value;
        ui.textEdit_2->append(QString::number(vYIn.at(i).value));
    }
    
    ui.textEdit_2->append("X坐标：");
    pointKnown pk1; int t2 = 0;
    for (auto& val : vXIn)
    {
        pk1.number = val.number1;
        pk1.x = x0 + val.value;
        pk1.y = y0 + vYIn.at(t2).value;
        vp.push_back(pk1);
        x0 = pk1.x;
        y0 = pk1.y;
        if (t2 <= vYIn.size()-2)
        {
            ui.textEdit_2->append(QString::number(pk1.x, 'f', 8));
        }
        t2++;
    }

    ui.textEdit_2->append("Y坐标：");
    vp.pop_back();
    for (auto& val : vp)
    {
        ui.textEdit_2->append(QString::number(val.y,'f',8));
    }
}

 八、程序运行结果如下：

 九、上述计算结果均校对过，结果准确。附上完整代码

//.cpp文件
#include "TraverseAdjustment.h"

TraverseAdjustment::TraverseAdjustment(QWidget *parent)
    : QWidget(parent)
{
    ui.setupUi(this);
    connect(ui.pushButton, SIGNAL(clicked()), this, SLOT(onReadDataKnown()));
    connect(ui.pushButton_2, SIGNAL(clicked()), this, SLOT(onReadDataAngle()));
    connect(ui.pushButton_3, SIGNAL(clicked()), this, SLOT(onReadDataSide()));
    connect(ui.pushButton_5, SIGNAL(clicked()), this, SLOT(onAngleClosureError()));
    connect(ui.pushButton_4, SIGNAL(clicked()), this, SLOT(CalculateAzimuth()));
    connect(ui.pushButton_6, SIGNAL(clicked()), this, SLOT(onCalIncrement()));
    connect(ui.pushButton_7, SIGNAL(clicked()), this, SLOT(onCalCorrection()));
    connect(ui.pushButton_8, SIGNAL(clicked()), this, SLOT(onCalCoordinate()));
    ui.lineEdit->setReadOnly(true);
    ui.lineEdit_2->setReadOnly(true);
    ui.lineEdit_3->setReadOnly(true);
    ui.textEdit->setReadOnly(true);
    ui.textEdit_2->setReadOnly(true);
    ui.textEdit->append("读入数据展示：");
}

TraverseAdjustment::~TraverseAdjustment()
{}

void TraverseAdjustment::onCalCoordinate()
{
    ui.textEdit_2->clear();
    ui.textEdit_2->append("改正后X增量：");
    double x0 = vpKnown.at(1).x;
    double y0 = vpKnown.at(1).y;

    for (int i=0;i<vXIn.size();i++)
    {
        vXIn.at(i).value = vXIn.at(i).value + vX.at(i).value;
        ui.textEdit_2->append(QString::number(vXIn.at(i).value));
    }

    ui.textEdit_2->append("改正后Y增量：");
    for (int i = 0; i < vYIn.size(); i++)
    {
        vYIn.at(i).value = vYIn.at(i).value + vY.at(i).value;
        ui.textEdit_2->append(QString::number(vYIn.at(i).value));
    }
    
    ui.textEdit_2->append("X坐标：");
    pointKnown pk1; int t2 = 0;
    for (auto& val : vXIn)
    {
        pk1.number = val.number1;
        pk1.x = x0 + val.value;
        pk1.y = y0 + vYIn.at(t2).value;
        vp.push_back(pk1);
        x0 = pk1.x;
        y0 = pk1.y;
        if (t2 <= vYIn.size()-2)
        {
            ui.textEdit_2->append(QString::number(pk1.x, 'f', 8));
        }
        t2++;
    }

    ui.textEdit_2->append("Y坐标：");
    vp.pop_back();
    for (auto& val : vp)
    {
        ui.textEdit_2->append(QString::number(val.y,'f',8));
    }
}

void TraverseAdjustment::onCalCorrection()
{
    ui.textEdit_2->clear();
    ui.textEdit_2->append("X增量改正数：");
    //计算坐标增量之和
    double detax = 0;
    double detay = 0;

    for (auto& val : vXIn)
    {
        double b1 = val.value;
        detax = detax + b1;
    }

    for (auto& val : vYIn)
    {
        double b1 = val.value;
        detay = detay + b1;
    }

    //计算坐标增量闭合差
    fx = detax - (vpKnown.at(2).x - vpKnown.at(1).x);
    fy = detay - (vpKnown.at(2).y - vpKnown.at(1).y);

    //计算坐标改正数
    double s1 = 0;//导线边距离之和
    for (auto& val : vS)
    {
        s1 = s1 + val.value;
    }

    Side s2;
    for (int i = 0; i<vXIn.size(); i++)
    {
        double d;
        for (int j=0; j < vS.size(); j++)
        {
            if (vXIn.at(i).number1 == vS.at(j).number1 && vXIn.at(i).number2 == vS.at(j).number2)
            {
                d = vS.at(j).value;
            }
        }
        s2.number1 = vXIn.at(i).number1;
        s2.number2 = vXIn.at(i).number2;
        s2.value = -fx / s1 * d;
        ui.textEdit_2->append(QString::number(s2.value));
        vX.push_back(s2);
    }

    ui.textEdit_2->append("Y增量改正数：");
    for (int i = 0; i < vYIn.size(); i++)
    {
        double d;
        for (int j = 0; j < vS.size(); j++)
        {
            if (vYIn.at(i).number1 == vS.at(j).number1 && vYIn.at(i).number2 == vS.at(j).number2)
            {
                d = vS.at(j).value;
            }
        }
        s2.number1 = vXIn.at(i).number1;
        s2.number2 = vXIn.at(i).number2;
        s2.value = -fy / s1 * d;
        ui.textEdit_2->append(QString::number(s2.value));
        vY.push_back(s2);
    }
}

void TraverseAdjustment::onCalIncrement()
{
    Side s1; ui.textEdit_2->clear();
    ui.textEdit_2->append("X坐标增量：");
    for (int i = 0; i < vAzi.size(); i++)
    {
        s1.number1 = vAzi.at(i).number;
        s1.number2 = vAzi.at(i).number2;
        for (int j = 0; j < vS.size(); j++)
        {
            if (vAzi.at(i).number == vS.at(j).number1 && vAzi.at(i).number2 == vS.at(j).number2)
            {
                double b1 = vAzi.at(i).value;
                b1 = Angle2Radian(b1);
                s1.value = vS.at(j).value * cos(b1);//cos函数输入值需为弧度
                ui.textEdit_2->append(QString::number(s1.value));
                vXIn.push_back(s1);
                s1.value = vS.at(j).value * sin(b1);//sin函数输入值需为弧度
                vYIn.push_back(s1);
                break;
            }
        }
    }

    ui.textEdit_2->append("Y坐标增量：");
    for (auto& val : vYIn)
    {
        ui.textEdit_2->append(QString::number(val.value));
    }
}

void TraverseAdjustment::onReadDataKnown()
{
    //打开文件对话框
    QString fileName = QFileDialog::getOpenFileName(this, tr("打开已知点数据"));
    QFile file(fileName);
    if (!file.open(QIODevice::ReadOnly | QIODevice::Text))
        QMessageBox::StandardButton btnValue = QMessageBox::information(this, tr("提示"), tr("打开失败！"));

    //读取文本文件数据
    QTextStream stream(&file);
    pointKnown pk1;
    while (!stream.atEnd())
    {
        QString str = stream.readLine();
        QStringList list=str.split(",");
        pk1.number = list.at(0);
        pk1.x = list.at(1).toDouble();
        pk1.y = list.at(2).toDouble();
        ui.textEdit->append(str);
        vpKnown.push_back(pk1);
    }

    ui.lineEdit->setText(fileName);
}

void TraverseAdjustment::onReadDataAngle()
{
    //打开文件对话框
    QString fileName = QFileDialog::getOpenFileName(this, tr("打开角度观测数据"));
    QFile file(fileName);
    if (!file.open(QIODevice::ReadOnly | QIODevice::Text))
        QMessageBox::StandardButton btnValue = QMessageBox::information(this, tr("提示"), tr("打开失败！"));

    //读取数据
    QTextStream stream(&file);
    Angle va1;
    while (!stream.atEnd())
    {
        QString str = stream.readLine();
        QStringList list = str.split(",");
        va1.number = list.at(0);
        va1.number1 = list.at(1);
        va1.number2 = list.at(2);
        va1.value = list.at(3).toDouble();
        ui.textEdit->append(str);
        vA.push_back(va1);
    }
    ui.lineEdit_2->setText(fileName);
}

void TraverseAdjustment::onReadDataSide()
{
    //打开文件对话框
    QString fileName = QFileDialog::getOpenFileName(this, tr("打开边长观测数据"));
    QFile file(fileName);
    if (!file.open(QIODevice::ReadOnly | QIODevice::Text))
        QMessageBox::StandardButton btnValue = QMessageBox::information(this, tr("提示"), tr("打开失败！"));

    //读取数据
    QTextStream stream(&file);
    Side vs1;
    while (!stream.atEnd())
    {
        QString str = stream.readLine();
        QStringList list = str.split(",");
        vs1.number1 = list.at(0);
        vs1.number2 = list.at(1);
        vs1.value = list.at(2).toDouble();
        ui.textEdit->append(str);
        vS.push_back(vs1);
    }
    ui.lineEdit_3->setText(fileName);
}

void TraverseAdjustment::CalculateAzimuth()
{
    //计算角度观测值改正数
    double t1 = -f / n;
    t1 = t1 / 3600;
    Angle a1;
    for (auto& val : vA)
    {
        a1.value = DMS2Angle(val.value)+t1;
        a1.number = val.number;
        a1.number1 = val.number1;
        a1.number2 = val.number2;
        va.push_back(a1);
    }

    ui.textEdit_2->clear();
    ui.textEdit_2->append("方位角结果：");
    Angle a2; double R1 = Rab;
    for (int i = 0; i < va.size(); i++)
    {
        double R = R1 + va.at(i).value-180;
        if (R > 360)
        {
            int t2 = R / 360;
            R = R - t2 * 360;
        }
        a2.number = va.at(i).number;
        a2.number1 = va.at(i).number1;
        a2.number2 = va.at(i).number2;
        a2.value = R;
        ui.textEdit_2->append(QString::number(R));
        vAzi.push_back(a2);
        R1 = R;
    }
}

double TraverseAdjustment::Angle2Radian(double a)
{
    a = a * 0.017453293;
    return a;
}

double TraverseAdjustment::Radian2Angle(double a)
{
    a = a * 57.29578;
    return a;
}

double TraverseAdjustment::DMS2Angle(double a)
{
    int deg; double min, sed;
    deg = int(a);
    min = int((a - deg) * 100);
    sed = ((a - deg) * 100 - min) * 100;
    return deg + min / 60 + sed / 3600;
}

void TraverseAdjustment::onAngleClosureError()
{
    double arfa_1=0;//角度观测值之和
    for (auto& val : vA)
    {
        double a = val.value;
        a = DMS2Angle(a);
        arfa_1 = arfa_1 + a;
    }

    qDebug() << arfa_1;
    qDebug() << DMS2Angle(1644.2828);

    double arfa_ab;//BA方位角值
    double arfa_cd;//CD方位角值

    arfa_ab = CoordinateInverse(vpKnown.at(0).x, vpKnown.at(0).y, vpKnown.at(1).x, vpKnown.at(1).y);
    arfa_cd = CoordinateInverse(vpKnown.at(2).x, vpKnown.at(2).y, vpKnown.at(3).x, vpKnown.at(3).y);
    qDebug() << arfa_ab;
    qDebug() << arfa_cd;
    qDebug() << DMS2Angle(303.2728);
    qDebug() << DMS2Angle(147.5628);
    Rab = arfa_ab;
    Rcd = arfa_cd;

    n = vA.size();
    double arfa_cd1 = 0;//理论角度值

    arfa_cd1 = arfa_ab - n * 180+arfa_1;

    if (arfa_cd1 > 360)
    {
        int t1 = arfa_cd1 / 360;
        arfa_cd1 = arfa_cd1 - t1 * 360;
    }

    double f_beta = arfa_cd1-arfa_cd;
    f_beta = f_beta * 3600;
    f = f_beta;

    ui.textEdit_2->append("角度闭合差：");
    ui.textEdit_2->append(QString::number(f_beta));
    ui.textEdit_2->append("限差：");
    ui.textEdit_2->append(QString::number(24*sqrt(n)));
    if (f_beta <= 24 * sqrt(n))
    {
        ui.textEdit_2->append("合格");
    }
    else
    {
        ui.textEdit_2->append("不合格");
    }
}

double TraverseAdjustment::CoordinateInverse(double xa, double ya, double xb, double yb)
{
    double arfa_ab = atan(abs(yb - ya) / abs(xb - xa));//返回的是弧度值
    arfa_ab = Radian2Angle(arfa_ab);
    if (ya<yb&&xa>xb)
    {
        arfa_ab=180 - arfa_ab;
    }

    if (ya > yb && xa > xb)
    {
        arfa_ab = 180 + arfa_ab;
    }

    if (ya > yb && xa < xb)
    {
        arfa_ab = 360 - arfa_ab;
    }
    return arfa_ab;//返回十进制度数
}
//.h文件
#pragma once

#include <QtWidgets/QWidget>
#include "ui_TraverseAdjustment.h"
#include<QFileDialog>
#include<QFile>
#include<QMessageBox>
#include<QTextStream>
#include<QStringList>
#pragma execution_character_set("UTF-8")
#include"math.h"
#include<qDebug>


struct pointKnown
{
    QString number;
    double x;
    double y;
};

struct Angle 
{
    QString number;
    QString number1;
    QString number2;
    double value;
};

struct Side
{
    QString number1;
    QString number2;
    double value;
};


class TraverseAdjustment : public QWidget
{
    Q_OBJECT

public:
    TraverseAdjustment(QWidget *parent = nullptr);
    ~TraverseAdjustment();

public:
    std::vector<pointKnown> vpKnown;//已知点坐标
    std::vector<Angle> vA;//角度观测值
    std::vector<Side> vS;//边长观测值
    double f;//角度闭合差
    int n;//导线边数
    std::vector<Angle> va;//添加改正数后的角度值
    std::vector<Angle> vAzi;//方位角值
    double Rab, Rcd;//已知起始、终止边对应方位角
    std::vector<Side> vXIn;//坐标增量X
    std::vector<Side> vYIn;//坐标增量Y
    double fx, fy;//坐标增量闭合差
    std::vector<Side> vX;//X增量改正数
    std::vector<Side> vY;//Y增量改正数
    std::vector<pointKnown> vp;//各测点最终坐标值

public slots:
    void onReadDataKnown();
    void onReadDataAngle();
    void onReadDataSide();
    void onAngleClosureError();
    void CalculateAzimuth();
    void onCalIncrement();
    void onCalCorrection();
    void onCalCoordinate();

public:
    double Angle2Radian(double a);
    double DMS2Angle(double a);
    double Radian2Angle(double a);
    double CoordinateInverse(double xa, double ya, double xb, double yb);

private:
    Ui::TraverseAdjustmentClass ui;
};
create_right
关注
  • 9
    点赞
  • 踩
  • 18
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
水准网平差C++语言实现
12-01
为解决水准网平差，利用间接平差模型，采用最小二乘进行精度评定，可直接利用VS2019以及其他版本运行运行。程序能较好的计算出平差结果
c#编写闭合导线简易平差程序
风风疯风的博客
08-01 5112
  本程序是使用c#语言编写的winform程序，使用VS2012开发，可实现闭合导线的观测数据导入、平差、结果显示、数据超限检查、平差结果导出到Excel表格、导线路线示意图绘制与导出。 
利用C#编写一个附和闭合导线平差程序
Zj1638的博客
07-06 6264
利用C#编写一个附和闭合导线平差程序
导线网平差程序系统设计—CSU测绘实习
SEREND1P17Y的博客
06-14 3030
导线网平差程序系统设计，控制网平差程序对野外控制网观测数据进行平差数据处理，目的是根据最小二乘原理，消除网中的各种几何矛盾，求出全网各待定点元素如未知点的平面坐标。❤来自CSU测绘实习
导线网平差程序系统设计
weixin_57056661的博客
08-08 1794
原始数据 界面设计 点击画图按钮后弹出的对话框 类的设计 Data类：存放点号、坐标、未知点误差椭圆元素 Gangle类：存储起测站、照准点、角度观测值 Length类：存储起测站、照准点、边长观测值 support类：存放各种函数和矩阵 Eigen类：用于矩阵运算 源代码 <Data.h> #pragma once #include <afxext.h> #include <math.h> #include <local....
附和导线平差程序（by C#）
QBigBangQ的博客
12-04 4824
本文主要内容包括本人编写附和导线平差程序的过程及附和导线平差程序源代码，在此记录，以供有兴趣的朋友参考学习。
1.附合导线近似平差_导线平差_附合导线_测量平差导线_conditionbqx_
10-03
在测量领域，导线平差是一项至关重要的工作，它涉及到空间几何定位的精度评估与修正。"附合导线近似平差"是其中的一种方法，适用于闭合或附合导线的处理，以提高测量数据的精确度。在本案例中，我们将深入探讨这一...
附合导线平差程序设计报告.doc
11-30
【附合导线平差程序设计】是测绘工程领域中的一种数据处理技术，它涉及到测量数据的计算和分析。在本文档中，作者罗正材使用C++编程语言设计了一个程序，该程序能够处理导线网的测量数据，进行平差计算，并将结果...
附合导线平差计算及CASS.dat文件导出.zip
04-18
在测绘领域，附合导线平差计算是一项基础且重要的工作。它涉及到测量学中的误差理论和数据处理技术。平差计算是通过一系列数学方法，对测量数据进行修正，以获得更加精确的测量结果。本资源"附合导线平差计算及CASS....
导线网平差程序设计 源代码
03-10
导线网平差程序设计 源代码导线网平差程序设计 源代码导线网平差程序设计 源代码导线网平差程序设计 源代码导线网平差程序设计 源代码导线网平差程序设计 源代码导线网平差程序设计 源代码
导线网平差程序
04-13
导线网平差程序，参数约束平差，经典自由网平差，秩亏自由网拟稳平差
C语言水准网间接平差
04-27
C语言水准网间接平差，国家二等水准网间接平差，处理报告存入txt文件
测绘导线网平差程序
09-19
漂亮完美的导线网平差程序。实现导线网的坐标概算，平差，精度评定，导线网不变形绘制，误差椭圆绘制以及误差椭圆缩放。输出的成功该有的都有。
经典平差C++版（测绘专业课设）
04-25
该平差课设包括条件平差，间接平差，附有参数的条件平差和简介平差，这些将是测绘本科生课设必经之路，之后会有专门的一遍博客讲述它的用法
导线网测量平差C++
09-13
本项目是基于C++编程语言，利用MFC（Microsoft Foundation Classes）对话框结构实现的导线网平差程序。下面将详细讲解这一领域的知识及其在C++和MFC中的应用。 1. **导线网测量基础**： - 导线网是由一系列测点和...
闭合水准网平差
yitian_00的博客
09-09 735
c++水准网间接平差法平差环境：vs2017+qt+eigen-3.3.8矩阵库
利用C#编写一个水准测量近似平差程序
Zj1638的博客
08-01 6521
利用 C#编写一个水准测量近似平差程序
【测绘程序设计】——附合导线近似平差
C_xxy's blog
01-05 1699
测绘程序设计——附合导线近似平差源代码（C++/MFC版，可定制）分享。
c++导线网间接平差程序设计
最新发布
01-03
导线网间接平差程序设计是指对电力系统中导线网的电流进行间接平差的一种计算方法。其设计包括了以下几个步骤。 首先，需要设定导线网的拓扑结构和导线参数。通过电力系统的拓扑图，确定导线网中各个节点之间的连接关系，即导线的连接方式。同时，还需要获取导线的相关参数，包括电阻、电抗、电纳等信息。这些信息可以通过实际测量或者电力系统模拟软件来获取。 接下来，进行导线网中的电流计算。根据欧姆定律和基尔霍夫电流定律，可以得到导线网中各个节点的电流方程。然后，利用迭代法或者高斯-赛德尔法等计算方法，逐步求解导线网中各节点的电流值。在计算过程中，需要考虑导线的损耗和电抗等参数，以得到较为准确的电流计算结果。 最后，进行导线网的间接平差计算。根据电流计算得到的结果，结合导线的参数信息，可以计算出导线网中各个节点的电压、功率等参数。通过对电流和电压的平差计算，可以进一步获得导线网中各个节点的电感、电阻、电容等参数的估计值。 在导线网间接平差程序设计中，需要考虑的因素还包括计算的精度和收敛性等问题。为了提高计算的精度，可以采用高阶数值计算方法或者增加迭代次数等手段。为了保证计算的收敛性，可以设置适当的收敛准则，或者使用松弛因子等技术。同时，还需要对计算结果进行合理性检验，通过与实际测量数据的对比，验证计算结果的准确性和可靠性。 总之，导线网间接平差程序设计是一项复杂而重要的工作，它为电力系统的稳定运行提供了重要的支持和保障。通过合理的设计和高效的计算方法，可以提高导线网的运行效率和电能质量，同时也为电力系统的规划和优化提供了有力的工具和决策依据。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值