//求取DEM模型的程序

#include?“stdafx.h“
#include?“iostream.h“
#include?“math.h“
#include?“stdio.h“
#include?“stdlib.h“

int?main（）
{
	int?duknown（double?a[20][3]）;
????int?duunknown（double?a[2][2]）;
	int?inverse（double?a[]?int?n）;
	double?xyh[20][3]={0}a[20][6]={0}n[6][6]={0}nn[36]={0}u[6]={0}r[6]={0}v[20]={0};
	/*設(shè)置各類計(jì)算過(guò)程中所要用到的矩陣，依次為已知20組測(cè)量值，系數(shù)矩陣A，
	法方程中的N矩陣，N矩陣的一維數(shù)組表示，法方程中的U矩陣，最終的參數(shù)解R矩陣改正數(shù)矩陣V*/
	double?xy[2][2]={0}h[2]={0};??//未知點(diǎn)坐標(biāo)，所要求的未知點(diǎn)高程
	double?miu=0sm=0;
	int?ijk;
	duknown（xyh）;
	duunknown（xy）;
	for（i=0;i<20;i++）???//由已知數(shù)據(jù)得出系數(shù)矩陣A
	{
		a[i][0]=1;
		a[i][1]=xyh[i][0];
		a[i][2]=xyh[i][1];
		a[i][3]=xyh[i][0]*xyh[i][0];
		a[i][4]=xyh[i][1]*xyh[i][0];
		a[i][5]=xyh[i][1]*xyh[i][1];
	}
	for（i=0;i<6;i++）???//由矩陣A求出N矩陣，此時(shí)各觀測(cè)量為獨(dú)立等精度，故權(quán)矩陣為單位陣
		for（j=0;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件?????????420??2013-06-10?12:57??最小二乘平差C++程序\Known.txt
?????文件?????????262??2013-06-23?20:02??最小二乘平差C++程序\Result.txt
?????文件??????????29??2013-06-10?18:07??最小二乘平差C++程序\UnKnown.txt
?????文件????????4252??2013-06-30?18:28??最小二乘平差C++程序\thefirst.cpp
?????目錄???????????0??2013-06-30?20:34??最小二乘平差C++程序\