fjdon.cpp

自己写的有关计算矩阵的小程序

#include "iostream.h"
using namespace std;

void main()
{
	float **a,x,result;
	int n;
	float k,b;
	void smallest(float **a,int n,float *k,float *b);
	cout<<"输入点数:";
	cin>>n;
	a=new float *[n];
	for(int i=0;i<n;i++)
	{
		a[i]=new float[2];
		cout<<"输入各个点:"; 
		for(i=0;i<5;i++)
		{
			for(int j=0;j<2;j++)
			{
				cin>>a[i][j];
				smallest(a,n,&k,&b);
				cout<<"输入自变量的值:";
				cin>>x;
				result=k*x+b;
				cout<<"拟合曲线方程:Y="; 
				cout<<"运算结果是:";
			}
		}
	}
}

void smallest(float **a,int n,float *k,float *b)
{
	float can[2][3]={n,0,0,0,0,0};
	for(int i=0;i<n;i++)
	{
		can[0][1]+=a[i][0];
		can[0][2]+=a[i][1];
		can[1][1]+=a[i][0]*a[i][0];
		can[1][2]+=a[i][0]*a[i][1];
	}
	can[1][0]=can[0][1];
	for(i=2;i>=0;i--)
	{
		can[1][i]=can[1][i]-can[0][i]*can[1][0]/can[0][0];
		*b=can[1][2]/can[1][1];
		*k=(can[0][2]-can[0][1]**b)/can[0][0];
	}
}