外点法C++.zip

  • 波波b
    了解作者
  • C/C++
    开发工具
  • 6KB
    文件大小
  • zip
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 1
    下载次数
  • 2019-07-18 04:39
    上传日期
C++实现外点法求解约束问题,其中有使用到单纯性法
外点法C++.zip
  • 外点法C++.doc
    35.5KB
内容介绍
<html xmlns="http://www.w3.org/1999/xhtml"> <head> <meta charset="utf-8"> <meta name="generator" content="pdf2htmlEX"> <meta http-equiv="X-UA-Compatible" content="IE=edge,chrome=1"> <link rel="stylesheet" href="https://static.pudn.com/base/css/base.min.css"> <link rel="stylesheet" href="https://static.pudn.com/base/css/fancy.min.css"> <link rel="stylesheet" href="https://static.pudn.com/prod/directory_preview_static/6275b4b29221806f9d796bef/raw.css"> <script src="https://static.pudn.com/base/js/compatibility.min.js"></script> <script src="https://static.pudn.com/base/js/pdf2htmlEX.min.js"></script> <script> try{ pdf2htmlEX.defaultViewer = new pdf2htmlEX.Viewer({}); }catch(e){} </script> <title></title> </head> <body> <div id="sidebar" style="display: none"> <div id="outline"> </div> </div> <div id="pf1" class="pf w0 h0" data-page-no="1"><div class="pc pc1 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/6275b4b29221806f9d796bef/bg1.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x1 h3 y2 ff1 fs0 fc0 sc0 ls0 ws0">#include&lt;iostream.h&gt;</div><div class="t m0 x1 h3 y3 ff1 fs0 fc0 sc0 ls0 ws0">#include&lt;math.h&gt;</div><div class="t m0 x1 h3 y4 ff1 fs0 fc0 sc0 ls0 ws0">int cons<span class="_ _0"></span>t n=2;</div><div class="t m0 x1 h3 y5 ff1 fs0 fc0 sc0 ls0 ws0">sta<span class="_ _0"></span>&#58903;c double Mk=1;</div><div class="t m0 x1 h4 y6 ff1 fs0 fc0 sc0 ls0 ws0">double fun_original(double x[n]);//<span class="ff2">&#26082;&#26159;&#20027;&#20989;&#25968;&#65292;&#20063;&#26159;<span class="_ _1"> </span></span>F(X,Mk)<span class="ff2">&#20043;&#19968;</span></div><div class="t m0 x1 h4 y7 ff1 fs0 fc0 sc0 ls0 ws0">double g_fun(double x[n]);//<span class="ff2">&#32422;&#26463;&#20989;&#25968;</span></div><div class="t m0 x1 h4 y8 ff1 fs0 fc0 sc0 ls0 ws0">double h_fun(double x[n]);;//<span class="ff2">&#31561;&#24335;&#20989;&#25968;</span></div><div class="t m0 x1 h3 y9 ff1 fs0 fc0 sc0 ls0 ws0">double a_e_f(double x[n]);</div><div class="t m0 x1 h4 ya ff1 fs0 fc0 sc0 ls0 ws0">double fun(double x[n]);//<span class="ff2">&#35753;&#21333;&#32431;&#24418;&#27861;&#27714;&#35299;&#27492;&#20989;&#25968;&#26368;&#23567;&#20540;</span></div><div class="t m0 x1 h4 yb ff1 fs0 fc0 sc0 ls0 ws0">void<span class="_ _2"></span> <span class="_ _2"></span>compare(double<span class="_ _2"></span> <span class="_ _2"></span>x[n+1][n],double<span class="_ _3"></span> <span class="_ _2"></span>LGH[3][n],double<span class="_ _2"></span> <span class="_ _2"></span>LGH_tap[3]);//<span class="_ _4"> </span><span class="ff2">&#22312;<span class="_ _4"></span>&#35843;<span class="_ _4"></span>&#29992;<span class="_ _4"></span>&#20043;<span class="_ _4"></span>&#21069;<span class="_ _5"> </span></span>LGH<span class="_ _5"> </span><span class="ff2">&#20808;<span class="_ _4"></span>&#32473;</span></div><div class="t m0 x1 h4 yc ff1 fs0 fc0 sc0 ls0 ws0">x[n+1]<span class="ff2">&#36171;&#20540;</span></div><div class="t m0 x1 h3 yd ff1 fs0 fc0 sc0 ls0 ws0">void zx(double x[n+1][n],double LGH[<span class="_ _0"></span>3][n],double LGH_ta<span class="_ _0"></span>p[3],double xn_1[n],double xn_2[n]);</div><div class="t m0 x1 h4 ye ff1 fs0 fc0 sc0 ls0 ws0">double H(double x[n+1][n],double LGH[3][n],double c);//<span class="ff2">&#21028;&#21035;&#20934;&#21017;</span></div><div class="t m0 x1 h3 yf ff1 fs0 fc0 sc0 ls0 ws0">void danchun(double original_point[n]);</div><div class="t m0 x1 h3 y10 ff1 fs0 fc0 sc0 ls0 ws0">void main()</div><div class="t m0 x1 h3 y11 ff1 fs0 fc0 sc0 ls0 ws0">{</div><div class="t m0 x2 h3 y12 ff1 fs0 fc0 sc0 ls0 ws0">int k=1;</div><div class="t m0 x2 h3 y13 ff1 fs0 fc0 sc0 ls0 ws0">int i;</div><div class="t m0 x2 h3 y14 ff1 fs0 fc0 sc0 ls0 ws0">double r=10;</div><div class="t m0 x2 h3 y15 ff1 fs0 fc0 sc0 ls0 ws0">double c=0.000001;</div><div class="t m0 x2 h4 y16 ff1 fs0 fc0 sc0 ls0 ws0">double x[n]={0,0.5};//<span class="ff2">&#32473;&#23450;&#30340;&#21021;&#22987;&#28857;&#22312;&#26576;&#20010;&#22686;&#24191;&#30446;&#26631;&#20989;&#25968;&#30340;&#27714;&#35299;&#22495;&#20043;&#20013;</span></div><div class="t m0 x2 h3 y17 ff1 fs0 fc0 sc0 ls0 ws0">do</div><div class="t m0 x2 h3 y18 ff1 fs0 fc0 sc0 ls0 ws0">{</div><div class="t m0 x3 h3 y19 ff1 fs0 fc0 sc0 ls0 ws0">danchun(x);</div><div class="t m0 x3 h4 y1a ff1 fs0 fc0 sc0 ls0 ws0">//ln(x2)&gt;=0<span class="_"> </span><span class="ff2">&#26102;<span class="_ _6"> </span></span>F(X,Mk)=-X1+X2;<span class="ff2">&#26080;&#26497;&#23567;&#28857;&#65292;&#19981;&#29992;&#31649;&#36825;&#20010;&#22686;&#24191;&#30446;&#26631;&#20989;&#25968;</span></div><div class="t m0 x3 h3 y1b ff1 fs0 fc0 sc0 ls0 ws0">if(Mk*a_e_f(x)&lt;c)</div><div class="t m0 x4 h3 y1c ff1 fs0 fc0 sc0 ls0 ws0">break;</div><div class="t m0 x3 h3 y1d ff1 fs0 fc0 sc0 ls0 ws0">Mk=r*Mk;</div><div class="t m0 x2 h3 y1e ff1 fs0 fc0 sc0 ls0 ws0">}while(1);</div><div class="t m0 x2 h4 y1f ff1 fs0 fc0 sc0 ls0 ws0">cout&lt;&lt;"<span class="ff2">&#29992;&#22806;&#28857;&#27861;&#27714;&#35299;&#65306;</span>"&lt;&lt;endl;</div><div class="t m0 x2 h3 y20 ff1 fs0 fc0 sc0 ls0 ws0">cout&lt;&lt;"minf(X)=-X1+X2"&lt;&lt;endl;</div><div class="t m0 x2 h3 y21 ff1 fs0 fc0 sc0 ls0 ws0">cout&lt;&lt;" g1(X)=lnX2&gt;=0"&lt;&lt;endl&lt;&lt;"s.t<span class="_ _0"></span>."&lt;&lt;endl;;</div><div class="t m0 x2 h3 y22 ff1 fs0 fc0 sc0 ls0 ws0">cout&lt;&lt;" g2(X)=X1+X2-1=0"&lt;&lt;endl;</div><div class="t m0 x2 h4 y23 ff1 fs0 fc0 sc0 ls0 ws0">cout&lt;&lt;"<span class="ff2">&#32467;&#26524;&#26159;&#65306;</span>"&lt;&lt;endl;</div><div class="t m0 x2 h3 y24 ff1 fs0 fc0 sc0 ls0 ws0">cout&lt;&lt;"Xk="&lt;&lt;endl;</div><div class="t m0 x2 h3 y25 ff1 fs0 fc0 sc0 ls0 ws0">for(i=0;i&lt;n;i++)</div><div class="t m0 x3 h3 y26 ff1 fs0 fc0 sc0 ls0 ws0">cout&lt;&lt;x[i]&lt;&lt;" ";</div><div class="t m0 x2 h3 y27 ff1 fs0 fc0 sc0 ls0 ws0">cout&lt;&lt;endl;</div><div class="t m0 x2 h3 y28 ff1 fs0 fc0 sc0 ls0 ws0">cout&lt;&lt;"f(Xk)="&lt;&lt;fun_original(x)&lt;&lt;endl;</div><div class="t m0 x2 h3 y29 ff1 fs0 fc0 sc0 ls0 ws0">//c<span class="_ _0"></span>out&lt;&lt;"Mk="&lt;&lt;Mk&lt;&lt;endl;</div><div class="t m0 x1 h3 y2a ff1 fs0 fc0 sc0 ls0 ws0">}</div><div class="t m0 x1 h4 y2b ff1 fs0 fc0 sc0 ls0 ws0">double fun_original(double x[n])//<span class="ff2">&#26082;&#26159;&#20027;&#20989;&#25968;&#65292;&#20063;&#26159;<span class="_ _1"> </span></span>F(X,Mk)<span class="ff2">&#20043;&#19968;</span></div><div class="t m0 x1 h3 y2c ff1 fs0 fc0 sc0 ls0 ws0">{</div><div class="t m0 x2 h3 y2d ff1 fs0 fc0 sc0 ls0 ws0">double y;</div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div> </body> </html>
评论
    相关推荐
    • 1.25外点法终稿.rar
      vc编程实现bp神经网络,比较实用!!还 有一些别的程序也在里面,希望能有用!!
    • structe.rar
      数据结构的资料,自己整理的,可用性一般,下载不好,请自己删掉!
    • GaussDB_100_1.0.1-DATABASE-REDHAT-64bit.tar.gz
      guassdb100在redhat上安装包,单机部署的包,安装步骤请看我的文中介绍,经过大量实验搭建总结出来的文档
    • SIM800C_MQTT.rar
      使用SIM800C模块,使用MQTT协议,连接中国移动onenet平台,能实现数据的订阅、发布、存储等
    • 卷积神经网络
      这是卷积神经网络的一个实际用例,已经调试好了,能够在matlab上成功运行,适合从事卷积神经网络(CNN)研究的人员学习使用。
    • android从bootloader到launcher启动流程整理
      讲述android 开机流程 从boot rom---bootloader---init--zygote---systemserver---ams 并附上自己整理的每个流程的流程图 ,清晰熟悉android 启动流程
    • 分数阶混沌系统:分数阶混沌系统的数值解。-matlab开发
      该工具箱包含可用于模拟一些著名的分数阶混沌系统的函数,例如: - 陈的系统, - Arneodo的系统, - Genesio-Tesi 的系统, - 洛伦兹系统, - 牛顿-莱普尼克系统, - 罗斯勒的系统, - Lotka-Volterra系统, - 达芬的系统, - 范德波尔的振荡器, -伏打的系统- 陆氏系统, - 刘的系统, - Chua的系统, - 金融系统, - 3 细胞 CNN。 这些函数以数值方式计算描述混沌系统的分数阶非线性微分方程的解。 每个函数返回总模拟时间的状态轨迹(吸引器)。 更多详情请看书: Ivo Petras,分数阶非线性系统:建模、分析和仿真,Springer,系列:非线性物理科学,2011,ISBN 978-3-642-18100-9。 http://www.springer.com/engineering/control/book/978-3-
    • matlabpam代码-dsp-library:用于光通信的MATLABDSP函数库
      matlab pam代码 DSP库 该项目包含几种不同的功能,可将DSP算法应用于光通信。 这些功能适用于相干和非相干(PAM,DMT)光通信。 大多数功能是相互独立的。 因此,该代码中的功能可以轻松使用,并与其他DSP功能结合使用。 用法 功能列表以及简短说明在文件中。 输入和输出参数的描述在每个函数的标题中。 讯号 通常,输入(和输出)信号在第一维度上具有时间(例如,列向量),而第二维度用于一次管理多个信号(例如,不同的极化,不同的参数等)。 之所以选择这种约定,是因为MATLAB通过将列保留在内存的连续部分中来存储矩阵,因此,这种约定比其他方法(行向量)要快。 参数 大多数功能使用参数结构作为输入参数。 函数中使用的参数的描述通常在标头中,而coherent-dsp函数的默认参数在文件中。 参考 执照 此代码在下发布。
    • 有关多目标跟踪的PHD滤波的一些资料
      一些关于多目标跟踪的新的文献,主要是有关概率假设密度(PHD)的
    • 基于python开发的全国新工商采集工具 v1.2版本
      这个软件是通过scrapy爬虫框架结合代理IP池再加上request模拟请求技术以及验证码识别技术,可以做到日更新采集全国新工商信息。采集的数据自动存储在mysql数据库表里,可下载全量1.8亿多企业工商基本信息和36维度的详细信息.支持sql和excel导出数据包格式。