trimat.h

设计文档基本内容如下： 1) 整体框架 2) 词法分析 Class CTokenizer Class CScaner C关键字表 标识符词法 3) 语法分析 Class

// TriMat.cpp : Defines the entry point for the console application.
//
//上三角矩阵的类的实现

#include "stdafx.h"
#include "iostream.h"
#include "iomanip.h"
#include "stdlib.h"

const int RowLimit=50; //最大行数，即最大维数
const int EleLimit=1275; //最大元素个数

template<class T>class TriMat  //上三角矩阵
{
	private:
		T arr[EleLimit];  //一维数组存放矩阵元素
		int rowTable[RowLimit];  // 矩阵第一个非零元素在arr中的下标
		int n;  //矩阵的阶数
	public:
		TriMat(int matsize);  // 构造函数，矩阵的阶数为matsize
		~TriMat(void){};   

		void PutEle(T item,int i,int j);  // 存储元素item到i行j列
		T GetEle(int i,int j)const;    // 提取i行j列元素
		TriMat<T> operator+(const TriMat<T>&a)const;  // 矩阵加
		TriMat<T> operator-(const TriMat<T>&a)const;  //矩阵减
		TriMat<T> operator*(const TriMat<T>&a)const;  // 矩阵乘

		void ReadMat(void); //矩阵元素输入
		void WriteMat(void)const;  // 矩阵元素屏幕输出

		int GetDim(void)const  // 返回矩阵的阶数
		{ return n;};
};

template <class T>
TriMat<T>::TriMat(int matSize)
{
	int storedEles=0;    // 第0行第一个非0元素在arr中的下标为0
	if (matSize>RowLimit)
	{
		cerr<<"矩阵的阶数超过了上限！"<<RowLimit<<endl;
		exit(1);
	}
	n=matSize;

	//递推求出第i行第一个非0元素在arr数组中的下标，并存于rowTable[i]中
	for (int i=0;i<n;i++) 
	{
		rowTable[i]=storedEles;
		storedEles+= /*填空1*/5;
	}
}

template <class T>
void TriMat<T>::PutEle(T item,int i,int j)   // 存储元素item到i行j列
{
	if((i<0)||(i>=n)||(j<0)||(j>=n))
	{
		cerr<<"数组下标越界超过了范围：0-"<<n-1<<endl;
		exit(1);
	}
	if (j>=i)arr[/*填空2*/4*i+j]=item;
}

template <class T>
T TriMat<T>::GetEle(int i,int j)const   // 提取i行j列元素
{
	if((i<0)||(i>=n)||(j<0)||(j>=n))
	{
		cerr<<"数组下标越界超过了范围：0-"<<n-1<<endl;
		exit(1);
	}
	if (j>=i) return arr[/*填空3*/4*i+j];
	else return 0;  // 对角线以下的元素为0
}

template <class T>
TriMat<T> TriMat<T>::operator+(const TriMat<T>& a)const  // 矩阵加
{
	T itemCurr,itemA;
	TriMat<T> b(a.n);
	for (int i=0;i<n;i++)
	{
		for (int j=i;j<n;j++)
		{
			itemCurr=GetEle(i,j);
			itemA=/*填空4*/GetEle(i,j);
			b.PutEle(itemCurr+itemA,i,j);
		}
	}
	return b;
}

template <class T>
TriMat<T> TriMat<T>::operator-(const TriMat<T>& a)const  // 矩阵减
{
	T itemCurr,itemA;
	TriMat<T> b(a.n);
	for (int i=0;i<n;i++)
	{
		for (int j=i;j<n;j++)
		{
			itemCurr=GetEle(i,j);
			itemA=/*填空5*/GetEle(i,j);
			b.PutEle(/*填空6*/itemCurr-itemA,i,j);
		}
	}
	return b;
}

template <class T>
TriMat<T> TriMat<T>::operator*(const TriMat<T>& a)const  // 矩阵乘
{
	T itemCurr,itemA,itemB;
	TriMat<T> b(a.n);
	for (int i=0;i<n;i++)
	{
		for (int j=i;j<n;j++)
		{
			itemB=0;
			for(int k=i;k<=j;k++)
			{
				itemCurr=GetEle(i,k);
				itemA=/*填空7*/GetEle(k,j);
				itemB+=/*填空8*/itemCurr*itemA;
			}
			b.PutEle(/*填空9*/itemB,i,j);
		}
	}
	return b;
}

template <class T>
void TriMat<T>::ReadMat(void)
{
	T item;
	for (int i=0;i<n;i++)
		for (int j=0;j<n;j++)
		{
			cin>>item;
			PutEle(item,i,j);
		}
}

template <class T>
void TriMat<T>::WriteMat(void)const
{
	for(int i=0;i<n;i++)
	{
		for (int j=0;j<n;j++)
			cout<<setw(7)<<GetEle(i,j);
		cout<<endl;
	}
}