梯度下降法C语言源程序梯度下降法C语言源程序
//****************************** //function: min f=x1*x1+25*x2*x2 //method:gradient
//****************************** #include
#include
#include
float eps = 0.000001;
void main()
{
double x[2],p[2],g[2];
double t,gradient,f;
int n = 0; ...
梯度下降法C语言源程序
//****************************** //function: min f=x1*x1+25*x2*x2 //method:gradient
//****************************** #include
#include
#include
float eps = 0.000001;
void main()
{
double x[2],p[2],g[2];
double t,gradient,f;
int n = 0;
x[0] = 2;
x[1] = 2;
g[0] = 2*x[0];
g[1] = 50*x[1];
while (sqrt(g[0]*g[0]+g[1]*g[1])>eps)
{
n++;
f = x[0]*x[0]+25*x[1]*x[1];
gradient = sqrt(g[0]*g[0]+g[1]*g[1]);
printf("Iteration number: %d\n",n);
printf("x[0] = %f; x[1]=%lf\n",x[0],x[1]);
printf("gradient: %lf\n",gradient);
printf("funciton value: %lf\n",f);
p[0] = -2*x[0];
p[1] = -50*x[1];
printf("g[0] = %f;g[1] =%f\n",g[0],g[1]);
t = -(2*x[0]*p[0]+50*x[1]*p[1])/(2*p[0]*p[0]+50*p[1]*p[1]);
printf("t = %lf\n",t);
printf("\n");
x[0] = x[0]+t*p[0];
x[1] = x[1]+t*p[1];
g[0] = 2*x[0];
g[1] = 50*x[1];
}
}
本文档为【梯度下降法C语言源程序】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。