共轭梯度法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语言实现
共轭梯度法程序源代码#include<stdio.h>#include<math.h>#defineN10#defineepspow(10,-6)doublef(doublex[],doublep[],doublet){doubles;s=pow(x[0]+t*p[0],2)+25*pow(x[..
文档分类: 幼儿/小学教育--教育管理 文档标签: 共轭梯度法C语言40西安交大41 系统标签: 交大共轭梯度maxtrprintfakv 共轭梯度法C语言(西安交大)#includestdio.h#includemath.h#define10/*定义矩阵阶数*/voidmain()doubleX[N],akv[N],dka[N],rk[N],dk[N],pk,pkk,ak,bk;for(i=0;iN;i++)/*输入系...
共轭梯度法的C实现西安交通大学实验报告 课程名称:数值分析上机实验 实验名称:共轭梯度法 学院:___数学学院___ 班级 姓名: 学号: 实验日期 2015 年 05 月 26 日 自评成绩:97 一、实验目的 (1)熟练掌握改进平方根法和共轭梯度法的迭代过程 (2)尝试使用自己熟悉的计算机语言解决数学中的问题 (3)通过上机实验...
文档标签: 共轭梯度法C语言实现 系统标签: 共轭 梯度 voidprint intn newdouble double 共轭梯度法.txt#include#includeint Gong(double**A,double*b,double*x,int n,int CND=100000);void print(double*b,int n);void print(double**A,int n);int main(){ const int n=5;int i,j; double**A=...
# 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程序代码
共轭梯度法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...
由线性无关方向组的定义,找不到一组非零实数 \alpha_i(i=1,\ldots,m) 使得\alpha_1d_1+\cdots+\alpha_md_m=0 ,则 d_i,\ldots,d_m 为线性无关方向组。 共轭梯度法描述 以下需要证明,以一个 Q- 共轭方向组为更新方向依次更新,可以最终找到全剧最小值。 A.按某一方向进行更新的最优步长 在...