其他回答

• 给你一个kernel做参考吧voidJacobidouble*A,double*x,doublelagma,unsignedintn{unsignedinti,j,k=0;unsignedintL=20;doubletemp;//doublemu=1e-12;//精度//doubleFanshu1;//x的一范数double*AA=double*mallocn+1*n+1;fori=1;i<=n;i++{*x+i=1;//迭代初始值设置为1*AA+i*n+i=*A+i*n+i-lagma;//AA=A-lagma*I}while1{fori=1;i<=n;i++{temp=0;//Fanshu1=0;forj=1;j<=n;j++ifj!=itemp-=*AA+i*n+j**x+j;*x+i=temp/*AA+i*n+i;}k++;ifk==L{//printf"Jacobi法求解Ax=b时迭代次数超过最大值";break;}}freeAA;。
  龚小霞2019-11-05 22:19:46
• 根据牛顿拉弗森迭代思想，把F看成x，g-x；x_k+1=x_k-gx_k/g''x_k,这里是矩阵，只要把除号变成求逆就可以了。
  窦道琴2019-11-05 22:05:47