电子文档交易市场
安卓APP | ios版本
电子文档交易市场
安卓APP | ios版本

非线性最优化理论与方法 实验报告

26页
  • 卖家[上传人]:靳**
  • 文档编号:69218189
  • 上传时间:2019-01-13
  • 文档格式:DOCX
  • 文档大小:335.05KB
  • / 26 举报 版权申诉 马上下载
  • 文本预览
  • 下载提示
  • 常见问题
    • 1、XIAN TECHNOLOGICAL UNIVERSITY课程设计报告课程名称 非线性最优化理论与方法 专 业: 数学与应用数学 班 级: 17100306 姓 名: 靳玉梅 学 号: 16100305120 指导教师: 张襄松 成 绩: 2018 年 11 月 28 日西安工业大学课程设计(论文)用纸第 页一、 实验目的(1)掌握计算步骤,并能熟练使用该方法求解线性规划问题。(2)了解算法算法程序实现的过程及方法二、 实验要求用c语言实现单纯形法、0.618法、最速下降法、Newton法、共轭梯度法中的FR方法、变度量法中的BFP算法,求解非线性规划问题,包括详细的接题步骤,运行程序,输出结果等。三、 实验环境Microsoft VisualC+6.0(绿色版)四、 实验内容及步骤、数据记录及处理(一)单纯形法(simplex method)Dantzing在1949年提出的单纯形法(simplex method)是一种行之有效的方法。基本自想是:从一个基本可行解出发,求出使目标函数值下降的另一个基本可行解。因此需要解决三个问题:(1) 求(LP)的初始基本可行解的方法;(2) 判别

      2、一个基本可行解是否为最优解的准则;(3) 从一个基本可行解转换到使目标数值下降的另一个基本可行解的方法。单纯形算法解线性规划的步骤及流程图:算法步骤(1)将线性规划化为标准型(2)用最快的方法确定一个初始基本可行解X(0)。当st均为“”形式时,以松驰变量做初始基本变量最快。(3)求X(0)中非基本变量xj的检验数j。若,则停止运算,X(0)=X*(表示最优解),否则转下一步。(4)由确定xk进基;由确定xl出基,其中alk称为主元素;利用初等变换将alk化为1,并利用alk将同列中其它元素化为0,得新解X(1)。(5)返回(3),直至求得最优解为止。程序图找出初始基本可行解列出初始单纯形表计算校验数已得最优解结束所有i=0?对于某个,j0是否存在Pj=0无最优解结束计算新的单纯形表确定为入基变量确定出基变量2、单纯形算法程序的规格说明各段代码功能描述:(1)、定义程序中使用的变量#include#include#define m 3 /*定义约束条件方程组的个数*/#define n 5 /*定义未知量的个数*/float M=1000000.0;float Amn; /*用于记录方

      3、程组的数目和系数;*/float Cn; /*用于存储目标函数中各个变量的系数*/float bm; /*用于存储常约束条件中的常数*/float CBm; /*用于存储基变量的系数*/float setam; /*存放出基与入基的变化情况*/float deltan; /*存储检验数矩阵*/float xn; /*存储决策变量*/int numm; /*用于存放出基与进基变量的情况*/float ZB=0; /*记录目标函数值*/(2)、定义程序中使用的函数void input();void print();int danchunxing1();int danchunxing2(int a);void danchunxing3(int a,int b);(3)、确定入基变量,对于所有校验数均小于等于0,则当前解为最优解。int danchunxing1()int i,k=0;int flag=0;float max=0;for(i=0;in;i+)if(deltai=0)flag=1;else flag=0;break;if(flag=1)return -1;for(i=0;in;i

      4、+) if(max deltai) max =deltai;k=i;return k;(4)、确定出基变量,如果某个大于0的校验数,对应的列向量中所有元素小于等于0,则线性规划问题无解。int danchunxing2(int a)int i,k,j;int flag=0;float min;k=a;for(i=0;im;i+)if(Aik=0)flag=1;else flag=0;break;if(flag=1)printf(n该线性规划无最优解!n); return -1;for(i=0;i0) setai=bi/Aik; else setai=M;min=M;for(i=0;i=setai) min=setai;j=i;numj=k+1;CBj=Ck;return j;(5)、迭代运算,计算新的单纯形表。void danchunxing3(int p,int q)int i,j,c,l;float temp1,temp2,temp3;c=p;/*行号*/l=q;/*列号*/temp1=Acl;bc=bc/temp1;for(j=0;jn;j+)Acj=Acj/temp1; for

      5、(i=0;im;i+) if(i!=c) if(Ail!=0) temp2=Ail; bi=bi-bc*temp2; for(j=0;jn;j+) Aij=Aij-Acj*temp2; temp3=deltal;for(i=0;in;i+)deltai=deltai-Aci*temp3;(6)、输入函数,输入方程组的系数矩阵、初始基变量的数字代码、方程组右边的值矩阵、目标函数各个变量的系数所构成的系数阵。void print()int i,j=0;printf(n-n);for(i=0;im;i+) printf(%8.2ftX(%d) %8.2f ,CBi,numi,bi); for(j=0;jn;j+) printf(%8.2f ,Aij); printf(n);printf(n-n);printf(ttt);for(i=0;in;i+)printf( %8.2f,deltai);printf(n-n);void input()int i,j; /*循环变量*/int k;printf(请输入方程组的系数矩阵A(%d行%d列):n,m,n);for(i=0;im;i+)for(j=

      6、0;jn;j+)scanf(%f,&Aij);printf(n请输入初始基变量的数字代码num矩阵:n);for(i=0;im;i+)scanf(%d,&numi);printf(n请输入方程组右边的值矩阵b:n);for(i=0;im;i+)scanf(%f,&bi);printf(n请输入目标函数各个变量的系数所构成的系数阵C:n);for(i=0;in;i+)scanf(%f,&Ci);for(i=0;in;i+)deltai=Ci;for(i=0;im;i+) k=numi-1; CBi=Ck;(7)、主函数,调用前面定义的函数。main()int i,j=0;int p,q,temp;input();printf(n-n);printf( tCBtXBtbt);for(i=0;in;i+)printf( X(%d)t,i+1);for(i=0;in;i+)xi=0;printf(n);while(1) q=danchunxing1(); if(q=-1) print(); printf(n所得解已经是最优解!n); printf(n最优解为:n); for(j=0;jm;j+) temp=numj-1; xtemp=bj; for(i=0;in;i+) printf(x%d=%.2f ,i+1,xi); ZB=ZB+xi*Ci; printf(ZB=%.2f,ZB); break; print(); p=danchunxing2(q); printf(np=%d,q=%d,p,q); if(q=-1) break; danchunxing3(p,q);输入:(1)、输入方程组的系数矩阵A(3行5列)(2)、输入初始基变量的数字代码num矩阵(3)、输入方程组右边的值矩阵b(4)、输入目标函数各个变量的系数所构成的系数阵C输出:(1)、输出是否为最优解(2)、输出最优解为多少 例题:P96页例4.2.1求解线性规划问题解:问题有一个明显的可行基B=(P3,P4,P5)对应的初始基本可行解x=(0,0,5,0,21)T,编译得:所以:最优解x=(2.75, 2.25,

      《非线性最优化理论与方法 实验报告》由会员靳**分享,可在线阅读,更多相关《非线性最优化理论与方法 实验报告》请在金锄头文库上搜索。

      点击阅读更多内容
    最新标签
    发车时刻表 长途客运 入党志愿书填写模板精品 庆祝建党101周年多体裁诗歌朗诵素材汇编10篇唯一微庆祝 智能家居系统本科论文 心得感悟 雁楠中学 20230513224122 2022 公安主题党日 部编版四年级第三单元综合性学习课件 机关事务中心2022年全面依法治区工作总结及来年工作安排 入党积极分子自我推荐 世界水日ppt 关于构建更高水平的全民健身公共服务体系的意见 空气单元分析 哈里德课件 2022年乡村振兴驻村工作计划 空气教材分析 五年级下册科学教材分析 退役军人事务局季度工作总结 集装箱房合同 2021年财务报表 2022年继续教育公需课 2022年公需课 2022年日历每月一张 名词性从句在写作中的应用 局域网技术与局域网组建 施工网格 薪资体系 运维实施方案 硫酸安全技术 柔韧训练 既有居住建筑节能改造技术规程 建筑工地疫情防控 大型工程技术风险 磷酸二氢钾 2022年小学三年级语文下册教学总结例文 少儿美术-小花 2022年环保倡议书模板六篇 2022年监理辞职报告精选 2022年畅想未来记叙文精品 企业信息化建设与管理课程实验指导书范本 草房子读后感-第1篇 小数乘整数教学PPT课件人教版五年级数学上册 2022年教师个人工作计划范本-工作计划 国学小名士经典诵读电视大赛观后感诵读经典传承美德 医疗质量管理制度 2 2022年小学体育教师学期工作总结 2022年家长会心得体会集合15篇
    关于金锄头网 - 版权申诉 - 免责声明 - 诚邀英才 - 联系我们
    手机版 | 川公网安备 51140202000112号 | 经营许可证(蜀ICP备13022795号)
    ©2008-2016 by Sichuan Goldhoe Inc. All Rights Reserved.