最小二乘法程序

原创 2007年09月17日 18:34:00
#include <math.h>
#include <stdio.h>
//////////////////////////////////////////////////////////////////////////////////////////
//矩阵结构体
struct Matrix
{
 int m,n;//m为行数,n为列数
 double **pm;//指向矩阵二维数组的指针
};
//初始化矩阵mt,并置矩阵的行为m,列为n
void InitMatrix(struct Matrix *mt,int m,int n)
{
 int i;
 //指定矩阵的行和列
 mt->m=m;
 mt->n=n;
 //为矩阵分配内存
 mt->pm=new double *[m];
 for(i=0;i<m;i++)
 {
  mt->pm[i]=new double[n];
 }
}
//用0初始化矩阵mt,并置矩阵的行为m,列为n
void InitMatrix0(struct Matrix *mt,int m,int n)
{
 int i,j;
    InitMatrix(mt,m,n);
 for(i=0;i<m;i++)
  for(j=0;j<n;j++) mt->pm[i][j]=0.0;
}
//销毁矩阵mt
void DestroyMatrix(struct Matrix *mt)
{
 int i;
 //释放矩阵内存
 for(i=0;i<mt->m;i++)
 {
  delete []mt->pm[i];
 }
 delete []mt->pm;
}
//矩阵相乘mt3=mt1*mt2
//成功返回1,失败返回0
int MatrixMul(Matrix *mt1,Matrix *mt2,Matrix *mt3)
{
 int i,j,k;
 double s;
 //检查行列号是否匹配
 if(mt1->n!=mt2->m||mt1->m!=mt3->m||mt2->n!=mt3->n) return 0;
 for(i=0;i<mt1->m;i++)
  for(j=0;j<mt2->n;j++)
  {
            s=0.0;
   for(k=0;k<mt1->n;k++) s=s+mt1->pm[i][k]*mt2->pm[k][j];
   mt3->pm[i][j]=s;
  }
  return 1;
}
//矩阵转置mt2=T(mt1)
//成功返回1,失败返回0
int MatrixTran(Matrix *mt1,Matrix *mt2)
{
 int i,j;
 //检查行列号是否匹配
 if(mt1->m!=mt2->n||mt1->n!=mt2->m) return 0;
 for(i=0;i<mt1->m;i++)
  for(j=0;j<mt1->n;j++) mt2->pm[j][i]=mt1->pm[i][j];
  return 1;
}
//控制台显示矩阵mt
void ConsoleShowMatrix(Matrix *mt)
{
 int i,j;
 for(i=0;i<mt->m;i++)
 {
  printf("/n");
  for(j=0;j<mt->n;j++) printf("%2.4f ",mt->pm[i][j]);
 }
 printf("/n");
}
//////////////////////////////////////////////////////////////////////////////////////////
//Guass列主元素消元法求解方程Ax=b,a=(A,b)
int Guassl(double **a,double *x,int n)
{
 int i,j,k,numl,*h,t;
    double *l,max;
 l=new double[n];
 h=new int[n];
 for(i=0;i<n;i++) h[i]=i;//行标
 for(i=1;i<n;i++)
 {
  max=fabs(a[h[i-1]][i-1]);
  numl=i-1;
  //列元的最大值
  for(j=i;j<n;j++)
  {
   if(fabs(a[h[j]][i-1])>max)
   {
    numl=h[j];
    max=fabs(a[h[j]][i-1]);
   }
  }
  if(max<0.00000000001) return 0;
  //交换行号
  if(numl>i-1)
  {
   t=h[i];
   h[i]=h[numl];
   h[numl]=t;
  }
  for(j=i;j<n;j++) l[j]=a[h[j]][i-1]/a[h[i-1]][i-1];
  for(j=i;j<n;j++)
            for(k=i;k<n+1;k++) a[h[j]][k]=a[h[j]][k]-l[j]*a[h[i-1]][k];
 }
 
 for(i=n-1;i>=0;i--)
 {
        x[i]=a[h[i]][n];
        for(j=i+1;j<n;j++) x[i]=x[i]-a[h[i]][j]*x[j];
        x[i]=x[i]/a[h[i]][i];
 }
 
 //清除临时数组内存
 delete []l;
 delete []h;
 return 1;
}
///////////////////////////////////////////////////////////////////////////////////////////
//最小二乘法求解矩阵Ax=b
int MinMul(Matrix A,Matrix b,double *x)
{
 int i,j;
 if(b.n!=1) return 0;
 if(A.m!=b.m) return 0;
 Matrix TranA;//定义A的转置
 InitMatrix0(&TranA,A.n,A.m);
    MatrixTran(&A,&TranA);
 Matrix TranA_A;//定义A的转置与A的乘积矩阵
 InitMatrix0(&TranA_A,A.n,A.n);
    MatrixMul(&TranA,&A,&TranA_A);//A的转置与A的乘积
 Matrix TranA_b;//定义A的转置与b的乘积矩阵
 InitMatrix0(&TranA_b,A.n,1);
 MatrixMul(&TranA,&b,&TranA_b);//A的转置与b的乘积
 DestroyMatrix(&TranA);//释放A的转置的内存
 Matrix TranA_A_b;//定义增广矩阵
 InitMatrix0(&TranA_A_b,TranA_A.m,TranA_A.m+1);
 //增广矩阵赋值
 for(i=0;i<TranA_A_b.m;i++)
 {
  for(j=0;j<TranA_A_b.m;j++) TranA_A_b.pm[i][j]=TranA_A.pm[i][j];
        TranA_A_b.pm[i][TranA_A_b.m]=TranA_b.pm[i][0];
 }
 DestroyMatrix(&TranA_A);
 DestroyMatrix(&TranA_b);
 //Guass列主消元法求解
 Guassl(TranA_A_b.pm,x,TranA_A_b.m);
 DestroyMatrix(&TranA_A_b);
 return 1;
}
int MinMul(double **A,double *b,int m,int n,double *x)
{
 int r,i;
 Matrix Al,bl;
 Al.pm=new double *[m];
 Al.m=m;
 Al.n=n;
 InitMatrix(&bl,m,1);
 for(i=0;i<m;i++)
 {
  Al.pm[i]=A[i];
  bl.pm[i][0]=b[i];
 }
 r=MinMul(Al,bl,x);
 delete []Al.pm;
 DestroyMatrix(&bl);
 return r;

最小二乘法编程实现方法

1  基于MATLAB编程的最小二乘法 例如,已知某一热敏电容的温度T和电容C的实测数据如表所示: 热敏电容的温度T和电容C的实测数据 T/℃ 50 60 70 80 90 100 10...
  • chanxiaogang
  • chanxiaogang
  • 2015年04月24日 11:50
  • 5893

最小二乘法C语言的实现

1.实验目的: 进一步熟悉曲线拟合的最小二乘法。 掌握编程语言字符处理程序的设计和调试技术。   2.实验要求: 输入:已知点的数目以及各点坐标 。 输出:根据最小二乘法原理以及各点坐标求...
  • huangshanchun
  • huangshanchun
  • 2015年12月26日 14:53
  • 3552

用C语言实现最小二乘法算法

用C语言实现最小二乘法算法本文博客链接:http://blog.csdn.net/jdh99,作者:jdh,转载请注明. 环境:主机:WIN8开发环境:MINGW说明:参考维基百科最小二乘法资料:...
  • jdh99
  • jdh99
  • 2015年03月09日 23:53
  • 10007

最小二乘法的C语言实现

1. 前言 最近断断续续看了一些数学书,有高等数学,也有初等数学,有科普的,有微积分的。 有时候,觉得数学才是世界上最美的东西,但有时候又觉得数学很高冷,不接地气。 不过,前段时间工作中用到...
  • shine_journey
  • shine_journey
  • 2017年06月07日 13:07
  • 1142

MATLAB实现最小二乘法

最小二乘法         最小二乘法(又称最小平方法)是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。         利用最小二乘法可以简便地求得未知的数据,并使得这些求得的数...
  • zengxiantao1994
  • zengxiantao1994
  • 2017年04月17日 15:10
  • 13058

最小二乘法及C语言实现

我们以最简单的一元线性模型来解释最小二乘法。什么是一元线性模型呢? 监督学习中,如果预测的变量是离散的,我们称其为分类(如决策树,支持向量机等),如果预测的变量是连续的,我们称其为回归。回归分析中,如...
  • MySniper11
  • MySniper11
  • 2013年12月08日 15:16
  • 10589

最小二乘法及其应用

基于 MATLAB 实验平台的 正规方程 最小二乘问题求解 了解最小二乘法的基本原理; 了解超定方程组的最小二乘解法; 使用最小二乘法解决实际问题;...
  • zstu_wangrui
  • zstu_wangrui
  • 2015年12月10日 10:47
  • 974

最小二乘法曲线拟合 C语言实现

简单思路如下: 1,采用目标函数对多项式系数求偏导,得到最优值条件,组成一个方程组; 2,方程组的解法采用行列式变换(两次变换:普通行列式——三角行列式——对角行列式——求解),行列式的求解算法上...
  • andylao62
  • andylao62
  • 2014年04月26日 11:55
  • 15927

数值作业:最小二乘法进行线性拟合之C语言代码

作图法求出直线的斜率a和截据b,可以确定这条直线所对应的经验公式,但用作图法拟合直线时,由于作图连线有较大的随意性,尤其在测量数据比较分散时,对同一组测量数据,不同的人去处理,所得结果有差异,因此是一...
  • Chen_dSir
  • Chen_dSir
  • 2017年04月22日 13:02
  • 3015

最小二乘法求AR模型

AR(Autoregressive)模型(自回归模型):用同一变量之前的表现情况来预测该变量现在或未来的表现情况,这种预测方法只与变量自己有关,而与其他变量无关,所以称作是自回归。 数学定义模型:假定...
  • u014557232
  • u014557232
  • 2016年03月26日 18:03
  • 3363
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:最小二乘法程序
举报原因:
原因补充:

(最多只允许输入30个字)