开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图象 > 查看源码
StdAfx.cpp
资源名称：DLT.rar　[点击查看]
上传用户：szacenet
上传日期：2022-06-20
资源大小：317k
文件大小：1k
源码类别：

图形图象

开发平台：

Visual C++

StdAfx.cpp：源码内容
  1. // stdafx.cpp : source file that includes just the standard includes
  2. // DLT.pch will be the pre-compiled header
  3. // stdafx.obj will contain the pre-compiled type information
  5. #include "stdafx.h"
  6. //矩阵转置
  7. void transpose(double *m1,double *m2,int m,int n) 
  8. { int i,j;
  9. for(i=0;i<m;i++)
  10. for(j=0;j<n;j++)
  11. m2[j*m+i]=m1[i*n+j];
  12. return;
  13. }
  15. /*正定矩阵求逆*/
  16. void inv(double *a,int n)
  17. { 
  18. int i,j,k;
  19.     for(k=0;k<n;k++)
  20. {
  21.         for(i=0;i<n;i++)
  22. {
  23.             if(i!=k)
  24. *(a+i*n+k)=-*(a+i*n+k)/(*(a+k*n+k));
  25. }
  26.         *(a+k*n+k)=1/(*(a+k*n+k));
  27.         for(i=0;i<n;i++)
  28. {
  29.             if(i!=k)
  30. {
  31.                 for(j=0;j<n;j++)
  32. {
  33.                     if(j!=k)
  34. *(a+i*n+j)+=*(a+k*n+j)* *(a+i*n+k);
  35. }
  36. }
  37. }
  38.         for(j=0;j<n;j++)
  39. {
  40.             if(j!=k)
  41. *(a+k*n+j)*=*(a+k*n+k);
  42. }
  43. }
  44. }
  45. //矩阵相乘
  46. void mult(double *m1,double *m2,double *result,int i_1,int j_12,int j_2) 
  47. { 
  48. int i,j,k; 
  49. for(i=0;i<i_1;i++)
  50.         for(j=0;j<j_2;j++)
  51. {
  52.             result[i*j_2+j]=0.0;
  53.             for(k=0;k<j_12;k++)
  54. result[i*j_2+j]+=m1[i*j_12+k]*m2[j+k*j_2];
  55. }
  56. return;
  57. }