共轭梯度法matlab程序
function x = cg(A,b)
tol=1e-10;
r = b + A*b;
w = -r;
z = A*w;
s = w'*z;
t = (r'*w)/s;
x = -b + t*w;
for k = 1:numel(b);
r = r - t*z;
if( norm(r) < tol )
return;
end
B = (r'*z)/s;
w = -r + B*w;
z = A*w;
s = w'*z;
t = (r'*w)/s;
x = x + t*w;
end
2019-12-21 19:39:19
304B
1
应用该代码进行优化时，首先应将matlab的路径设置至该文件夹的子目录下，然后输入初值，输入函数，点击回车，即可求解
2019-12-21 19:32:03
995B
1
MATLAB的梯度法，内点法，外点法，罚函数，惩罚函数，线性梯度法，源程序，按照提示输入，可直接运行-MATLAB' s gradient method, interior point method, outside the point of law, penalty function, penalty function, the linear gradient method, source code, follow the prompts to input, can be directly run
2019-12-21 19:30:03
3KB
1
共轭梯度法c++程序，也有黄金分割法和进退搜索法，可以使用，有介绍。/*以下是进退法搜索区间源程序*/
void sb(double *a,double *b,double x[],double p[])
{
double t0,t1,t,h,alpha,f0,f1;
int k=0;
t0=2.5; /*初始值*/
h=1; /*初始步长*/
alpha=2; /*加步系数*/
f0=f(x,p,t0);
t1=t0+h;
f1=f(x,p,t1);
while(1)
{
if(f1
2019-12-21 19:24:32
30KB
1