- 1、本文档共14页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验汇报阻尼牛顿法
太原理工大学机械学院机测系课程上机实验报告
课 程 名 称: 机械优化设计
班级?日 期?成绩评定?姓名实验室?老师签名?实验名称用阻尼牛顿法程序解题所用软件C++ DEV实验目的
及
内
容?实验目的:
1.掌握并能够建立最优化基本类型问题的数学模型。
2.掌握最优化方法的基本概念、基本理论和基本方法,奠定最优化的理论基础。
3.能够熟练编制和调试最优化方法的程序,奠定解决实际中的优化问题的基础
实验内容:
理解阻尼牛顿法并编写相关程序求其最优解。
牛顿法程序考核题
实
验
?原
理
步
骤
、实验原理:
实验步骤:
1,画流程图,编写程序;
2,将目标函数代入;
3,编译运行,将结果保存实
验
结
果
及
分
析 *********阻尼牛顿法计算结果**********
++++++一维搜索方法:黄金分割法++++++
初始坐标:
x( 0)=[ 0.0000000, 0.0000000], f( 0)= 16.0000000
迭代轮数 k= 1
x( 1)=[ -1.1249988, 0.7499992], f( 1)= 9.8125000
迭代精度: 0.000009986
迭代轮数 k= 2
x( 2)=[ -1.1250000, 0.7500000], f( 2)= 9.8125000
迭代精度: 0.000000114
*********************
阻尼牛顿法法优化最优点及目标函数值为:
x( * )=[ -1.1250000, 0.7500000], f( * )= 9.8125000
迭代精度: 0.000000114
算法程序实现
/*csssqj.cpp */
#include string.h
#include stdio.h
#include math.h
#include stdlib.h
#include time.h
#define N 2 /*优化设计维数*/
#define EPSIN 0.000001 /*迭代精度*/
#define H_QJ 1.0 /*初始区间搜索步长*/
#define Y_F 1 /*一维搜索方法选择:1——黄金分割法*/
/* 2——二次插值法*/
FILE *fp;
char outname[50]=阻尼牛顿法计算结果.txt; /*计算结果输出文件*/
/*给出初始点坐标*/
void csd_x(double x0[])
{
int i;
for(i=0;iN;i++) /*初始点为坐标原点的情况*/
x0[i]=0.0;
return;
}
/*目标函数*/
double hanshu(double x[])
{
double f;
f=4.0*(x[0]+1.0)*(x[0]+1.0)+2.0*(x[1]-1.0)*(x[1]-1.0)
+x[0]+x[1]+10.0; /*阻尼牛顿法*/
return f;
}
/*函数的梯度*/
void gread(double x[],double gf[])
{
/*阻尼牛顿法*/
gf[0]=8.0*(x[0]+1.0)+1.0;
gf[1]=4.0*(x[1]-1.0)+1.0;
return;
}
/*Heisen矩阵*/
void heisen(double x[],double hei[][N])
{
hei[0][0]=8.0;
hei[0][1]=0.0;
hei[1][0]=0.0;
hei[1][1]=4.0;
return;
}
/*计算f(xk+as)*/
double xkadd(double x[],double d[],double a)
{
int i;
double x1[N];
for(i=0;iN;i++)
x1[i]=
您可能关注的文档
最近下载
- 工程弹塑性力学题库及答案.pdf
- SHS35-39除尘脱硫课程设计.docx
- 人教版美术六年级下册期末测试题+答案.docx
- 北师大版八年级下册数学期末几何压轴题专练(含答案).docx VIP
- 2024首届全国红旗杯班组长大赛题库及答案(第4001-6000题).docx VIP
- 2024首届全国红旗杯班组长大赛题库及答案(2001-4000题).docx VIP
- 2024首届全国红旗杯班组长大赛题库及答案 (1一2000题).docx VIP
- YJK隔震结构设计手册.pdf
- 小学四年级春季学期《道德与法治》(统编版)学习任务单(全册汇总).doc
- 党史新中国史改革开放史社会主义发展史四史学习教育.pptx VIP
文档评论(0)