//?Approach1.cpp?:?此文件包含?“main“?函數。程序執行將在此處開始并結束。
//


#include???

using?namespace?std;


int?main（）
{
	//變量定義
	int?n?t?a;
	double?u_?=?1.96;
	double?c0;
	double?umax;

	cout?<	cout?<	cout?<	cout?<	cout?<	cout?<	//輸入t、n
	cout?<	cin?>>?n;
	cout?<	cin?>>?t;
	cout?<	cout?<	cout?<

	Matrix?B?=?Matrix（n?t）;
	Matrix?BT?=?Matrix（t?n）;
	Matrix?P_input?=?Matrix（n?1）;
	Matrix?P?=?Matrix（n?n）;
	Matrix?L?=?Matrix（n?1）;
	Matrix?d?=?Matrix（n?1）;
	Matrix?X0?=?Matrix（t?1）;
	Matrix?l?=?Matrix（n?1）;
	Matrix?BX0?=?Matrix（n?1）;
	Matrix?L_BX0?=?Matrix（n?1）;
	Matrix?Nbb?=?Matrix（t?t）;
	Matrix?invNbb?=?Matrix（t?t）;
	Matrix?Qvv?=?Matrix（n?n）;
	Matrix?Q?=?Matrix（n?n）;
	Matrix?eye?=?Matrix（n?n）;
	Matrix?BinvNbb?=?Matrix（n?t）;
	Matrix?BinvNbbBT?=?Matrix（n?n）;
	Matrix?x?=?Matrix（t?1）;
	Matrix?invNbbBT?=?Matrix（t?n）;
	Matrix?invNbbBTP?=?Matrix（t?n）;
	Matrix?V?=?Matrix（n?1）;
	Matrix?vtpv?=?Matrix（1?1）;
	Matrix?vtp?=?Matrix（1?n）;
	Matrix?u?=?Matrix（n?1）;
	Matrix?diag_Qvv?=?Matrix（n?1）;
	Matrix?sqrt_diag_Qvv?=?Matrix（n?1）;
	Matrix?uc0?=?Matrix（n?1）;
	Matrix?ru?=?Matrix（n?1）;
	Matrix?find_error?=?Matrix（1?2）;
	Matrix?L_?=?Matrix（n?1）;






	//輸入B
	input（B?“請輸入B矩陣“?“B.txt“）;
	BT?=?B.T（）;
	//輸入P
	input（P_input?“請輸入P矩陣“?“P.txt“）;
	P?=?diag2m（P_input）;
	cout?<	//輸入L
	input（L?“請輸入L向量“?“L.txt“）;
	//輸入d
	input（d?“請輸入d向量“?“d.txt“）;
	//輸入X0
	input（X0?“請輸入X0向量“?“X0.txt“）;
	//計算l
	BX0?=?B?*?X0;
	//cout?<	//cout?<	L_BX0?=?L-BX0;
	//cout?<	//cout?<	l?=?（L_BX0?-?d）*1000;
	cout?<	cout?<	//計算Nbb、Qxx
	Nbb?=?（BT*P）*B;
	cout?<	cout?<	cout?<	invNbb?=?Nbb.inverse（）;
	cout?<	//計算Qvv
	eye?=?eye.eye（）;
	Q?=?inv_diag（P）;
	//cout?<	//cout?<	BinvNbb?=?B?*?invNbb;
	//cout?<	//cout?<	BinvNbbBT?=?BinvNbb?*?BT;
	//cout?<	//cout?<	Qvv?=?Q-?BinvNbbBT;
	//cout?<	//cout?<	//求x
	invNbbBT?=?invNbb?*?BT;
	invNbbBTP?=?invNbbBT?*?P;
	x?=?invNbbBTP?*?l;
	cout?<	cout?<	//求V
	V?=?B?*?x?-?l;
	cout?<	cout?<	vtp?=?V.T（）*P;
	vtpv?=?vtp*V;
	//求單位權中誤差
	c0?=?sqrt（vtpv.det（）?/?（n?-?t））;
	cout?<	//求檢測量u
	diag_Qvv=ext_diag（Qvv）;
	//cout?<	//cout?<	sqrt_diag_Qvv?=?Msqrt（diag_Qvv）;
	//cout?<