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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > matlab例程 > 查看源码
c7-2.cpp
资源名称：数值分析课程设计.zip　[点击查看]
上传用户：zhdd911129
上传日期：2007-05-11
资源大小：722k
文件大小：1k
源码类别：

matlab例程

开发平台：

Matlab

c7-2.cpp：源码内容
  1. //C7
  2. //Solve a linear system,using Jacobi Iteration
  4. const int N=4;
  5. #include<iostream.h>
  7. int ReachTOL(double x1[],double x2[],double tol,int n)
  8. {
  9. int i;
  10. double error=0;
  11. for (i=0;i<n;i++)
  12. {
  13. error+=(x1[i]-x2[i])*(x1[i]-x2[i]);
  14. }
  15. if(error<tol) 
  16. return 1;
  17. else
  18. return 0;
  19. }
  21. void JI(double a[][N],double b[],double x0[],double tol,int max,int n)
  22. {
  23. double x1[N];
  24. int k,i,j;
  25. for (i=0;i<n;i++)
  26. {
  27. x1[i]=0;
  28. }
  29. for (k=0;k<max;k++)
  30. {
  31. for (i=0;i<n;i++)
  32. {
  33. for (j=0;j<n;j++)
  34. {
  35.     if(j!=i)
  36. {
  37.     x1[i]+=a[i][j]*x0[j];
  38. } 
  39. }
  40.     x1[i]=((-1)*x1[i]+b[i])/a[i][i];
  41. }
  42. if(ReachTOL(x0,x1,tol,n))
  43. {
  44. for (i=0;i<n;i++)
  45. { 
  46.         cout<<"x"<<i<<" is: "<<x1[i]<<endl;
  47. } 
  48. cout<<endl;
  49. return;
  50. }
  51. for (i=0;i<n;i++)
  52. {
  53.     x0[i]=x1[i];
  54. x1[i]=0;
  55. }
  56. }
  57. cout<<"Exceed the maximium iteration time"<<endl;
  58. for (i=0;i<n;i++)
  59. { 
  60.          cout<<"x"<<i<<" is: "<<x1[i];
  61. }
  62. cout<<endl;
  63. }
  65. void main(void)
  66. {
  67. double a[N][N]={{78,-2,-12,-14},   //为克服C语言中输入浮点数的漏洞，这里采用先输大数后输小数的方式
  68. {-2,86,-4,6},
  69. {-12,-4,72,-8},
  70. {-14,6,-8,74}};
  71. double b[N]={76,8,112,68};
  72. double x0[N]={0,0,0,0};
  73. double tol=1e-2;
  74. int max=100,i,j;
  75. for(i=0;i<N;i++)
  76. {
  77. for(j=0;j<N;j++)
  78. {
  79. a[i][j]/=100;
  80. }
  81. }
  82. for (i=0;i<N;i++)
  83. {
  84. b[i]/=100;
  85. }
  86.     JI(a,b,x0,tol,max,N);
  87. }