juzheng.cpp

本源码包括矩阵运算的基本功能

#include "matrix.h"
#ifndef _NO_NAMESPACE
using namespace std;
using namespace math;
#define STD std
#else
#define STD
#endif
#ifndef _NO_TEMPLATE
typedef matrix<double> Matrix;
#else
typedef matrix Matrix;
#endif
#include <stdio.h>
#include <stdlib.h>
void main()
{
	FILE *fp;
    Matrix A(3,5),B(5,5),C(5,3),D(3,3);
	int i,j;
	float f;
		if((fp=fopen("1.txt","r"))==NULL)
	{
	printf("this file can't be open\n");
	exit(0) ;
	}
	for(i=0;i<3;i++)
	{	
		for(j=0;j<5;j++)
		{
			fscanf(fp,"%f",&f);
			A(i,j)=f;
		}
	}
	fclose(fp);
	if((fp=fopen("2.txt","r"))==NULL)
	{
	printf("this file can't be open\n");
	exit(0) ;
	}
	for(i=0;i<5;i++)
	{	
		for(j=0;j<5;j++)
		{
			fscanf(fp,"%f",&f);
			B(i,j)=f;
		}
	}
	fclose(fp);
	if((fp=fopen("3.txt","r"))==NULL)
	{
	printf("this file can't be open\n");
	exit(0) ;
	}
	for(i=0;i<5;i++)
	{	
		for(j=0;j<3;j++)
		{
			fscanf(fp,"%f",&f);
			C(i,j)=f;
		}
	}
	fclose(fp);
	D=A*B*C;
	for(i=0;i<3;i++)
	{
		for(j=0;j<3;j++)
		{
			printf("%f\t",D(i,j));		
		}
			printf("\n");
	}
}