#include?“QrInv.h“
int?c_qrs（int?ndouble?**?realdouble?**?imgdouble?**R_realdouble**?R_imgdouble?**?Q_realdouble?**?Q_img）
{
	int?ncount?=?0;
	double?**?HMatReal?=?new?double?*[n];
	double?**?HMatImg?=?new?double?*[n];
	double?**?ResMatReal=?new?double?*[n];
	double?**?ResMatImg?=?new?double?*[n];
	for?（?int?i?=?0;?i?	{
		HMatReal[i]?=?new?double[n];
????????HMatImg[i]?=?new?double[n];
		ResMatReal[i]?=?new?double[n];
		ResMatImg[i]?=?new?double[n];
	}

????for?（i?=?0;?i?????{
		for?（int?j?=?0;?j?		{
			if?（i?==?j）
			{
				Q_real[i][j]?=?1;
			}
			else
			{
????????????????Q_real[i][j]?=?0;
			}
			Q_img[i][j]?=?0;
		}
????}
????c_copmatris（n?R_realR_imgrealimg）;
	for?（?i?=?n;?i?>?1;?i--?）????//?總體循環(huán)，進行QR分解
	{
		double?*?a1Real?=?new?double?[i];
		double?*?a1Img?=?new?double[i];
		double?*?e1Real?=?new?double[i];
		double?*?e1Img?=??new?double[i];
		
		double?*?w1Real?=?new?double?[i];
		double?*?w1Img?=?new?double[i];
		ncount++;
		for?（int?j?=?0;?j?		{
			for?（int?k?=?0;?k?			{
				if?（j?==?k）
				{
					HMatReal[j][k]?=?1;
				}
				else
				{
					HMatReal[j][k]?=?0;
				}
				HMatImg[j][k]?=?0;
			}
		}?//?初始化單位矩陣
		for?（?j?=?0;?j?		{
			a1Real[j]?=?R_real[（n?-i）?+?j][n?-?i];??//獲取第一例實數(shù)部分
			a1Img[j]?=?R_img[（n?-i）?+?j][n?-?i];????//獲取第一例虛數(shù)部分
			if?（j?==?0）
			{
				e1Real[j]?=?a1Real[j];???????????//?構造單位向量實部（共軛）
				e1Img[j]?=??a1Img[j];
			}
			else
			{
				e1Real[j]?=?0;???????????//?構造單位向量實部
				e1Img[j]?=?0;?
			}
		????????????????????????????????//?構造單位向量虛部
		}
		double?**?MatReal?=?new?double?*[i];
		double?**?MatImg?=?new?double?*[i];

		for?（j?=?0;?j?		{
			MatReal[j]?=?new?double[i];
			MatImg[j]?=?new?double[i];

		}
		double?norm?=?c_norm（ia1Reala1Img）;??//?求a向量模??
????????
		???????????????????????????????????????//?判斷向量a的第一個元素是否為復數(shù)
		if?（norm?==?0）
		{
			return?-2;??//?矩陣不可逆該列出現(xiàn)全零現(xiàn)象，該矩陣不可逆
		}
		if?（e1Img[0]?!=?0）
		{
			double?eNorm?=?sqrt（e1Real[0]?*?e1Real[0]?+?e1Img[0]?*?e1Img[0]）;
			e1Real[0]?=?（norm?/?eNorm）?*?e1Real[0];
			e1Img[0]?=??（norm?/?eNorm）?*?e1Img[0]?;

		}?else
		{
			e1Real[0]?=?1;
????????????e1Img[0]?=?0;
			c_multvector（i?norm?e1Reale1Img）;???//?實數(shù)與向量乘積
		}
		c_subvector（i?a1Reala1Img?e1Reale1Img）;?//向量相減
????????norm?=?c_norm（ia1Reala1Img）;???//?求相減后的模
		if?（norm?		{
			continue;?????
		}
	????c_dividvector（i?norma1Reala1Img）;???//?除以模
		c_mutvectors（i?a1Reala1Imga1Reala1ImgMatRealMatImg）;??//?向量相乘得到
		for?（j?=?0;?j?		{
			for?（int?k?=?0;?k?			{
				if?（k?==?j）
				{
					MatReal[j][k]?=?1?-?2?*?MatReal[j][k];
				}
				else
				{
					MatReal[j][k]?=??-2?*?MatReal[j][k];	
				}
				MatImg[j][k]?=?-2?*?MatImg[j][k];?
			}
		}
????//?求得H1（MatRealMatImg）?H1?*?A

		for?（j?=?0;?j?