《最优化方法及其应用》,郭科、陈聆、魏友华编7)解线性方程组的共轭梯度法(3人)/*刘森林、武继飞、芦佩*/用C语言将共轭梯度法编写成通用的子程序,然后用你编写的程序求解
# include <stdio.h> # include <math.h> # define N 5 /*目标函数*/ double F(double x[],double p[],double t) { double f=0; int i; for(i=0;i<N;i++) { f+=(x[i]+t*p[i])*(x[i]+t*p[i]); } return f; } /*欧氏范数平方*/ double D(double g[]) { double d=0...
共轭梯度法C语言(西安交大)#include<stdio.h> #include<math.h> #define N 10 /*定义矩阵阶数*/ void main() { int i,j,m,A[N][N],B[N]; double X[N],akv[N],dka[N],rk[N],dk[N],pk,pkk,ak,bk; for(i=0;i<N;i++) /*输入系数矩阵A*/ { for(j=0;j<N;j++) { if(i=...
c语言实现socket 热度: 用C语言实现雅可比迭代法 热度: 共轭梯度法.txt #include #include intGong(double**A,double*b,double*x,int n,intCND=100000); voidprint(double*b,intn); voidprint(double**A,intn); intmain() { constintn=5;inti,j; ...
共轭梯度法c语言 #include<stdio.h> #include<stdlib.h> #include<math.h> #include<string.h> #define N 4 #define epsilon 0.00005 /* 矩阵A和向量b,相乘结果存在Ab里 */ void matrixTimesVec(double A[N][N], double b[N], double Ab[N])...
共轭梯度法C语言源程序 //*** //function:minf=x1*x1+25*x2*x2 //method:gradient //*** #include<stdlib.h> #include<stdio.h> #include<math.h> constfloateps=0.000001; voidmain() { doubled[2],x[2],g[2],temp[2],lamda,alpha,g_org,g_now,t; intn=0; x[0]=2; x[1]=2...
共轭梯度法C语言(西安交大) #include #include #define N 10 /*定义矩阵阶数*/ void main() { int i,j,m,A[N][N],B[N]; double X[N],akv[N],dka[N],rk[N],dk[N],pk,pkk,ak,bk; for(i=0;i<=""> for(j=0;j<n;j++)< p=""> { if(i==j) A[i][j]=-2; else if(...