当前位置:网站首页>Example nonlinear integer programming
Example nonlinear integer programming
2022-07-02 17:26:00 【__ Rain】
example 1
The title is as follows
Using Monte Carlo method ( Random sampling method )
Although calculating the optimal solution is not realistic , But the application of probability theory can prove , In the case of a certain amount of calculation , A satisfactory solution can be obtained
The implementation of this algorithm is very simple , Easy to understand , Maintain the maximum value of random number of violent cycle
1) First write M M M file m e n g t e . m mengte.m mengte.m Define the objective function f f f and Constraint vector function g g g, The procedure is as follows
be aware g g g The essence of 5 × 1 5\times1 5×1 The column vector
function [f,g]=mengte(x);
% Input parameters x, Return value [f,g], The name of the function mengte
f=x(1)^2+x(2)^2+3*x(3)^2+4*x(4)^2+2*x(5)^2-8*x(1)-2*x(2)-3*x(3)...
-x(4)-2*x(5);
g=[sum(x)-400
x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)-800
2*x(1)+x(2)+6*x(3)-200
x(3)+x(4)+5*x(5)-200];
end
2) To write M M M file w o r k . m work.m work.m Answer the question
rand('state', sum(clock));
Max=-10000000;% The maximum value of the function is initialized to negative infinity
tic% Show run time
for i=1:10^6
x=99*rand(5,1);% Generate 5*1 Random real matrix of
x1=floor(x);x2=ceil(x);% integer
[f,g]=mengte(x1);
if sum(g<=0)==4% Input parameter is x1 when , If all inequalities are satisfied
if Max<=f% Maintain the maximum value of the function
Max=f;
ans=x1;
end
end
[f,g]=mengte(x2);% repeat
if sum(g<=0)==4
if Max<=f
Max=f;
ans=x2;
end
end
end
ans,Max
toc
lingo solution
model:
sets:
row/1..4/:b;
col/1..5/:c1,c2,x;
link(row,col):a;
endsets
data:
c1=1,1,3,4,2;
c2=-8,-2,-3,-1,-2;
a=1 1 1 1 1
1 2 2 1 6
2 1 6 0 0
0 0 1 1 5;
b=400,800,200,200;
enddata
max=@sum(col:c1*x^2+c2*x);
@for(row(i):@sum(col(j):a(i,j)*x(j))<b(i));
@for(col:@gin(x));
@for(col:@bnd(0,x,99));
end
Example 2
Explain :
1) To write M M M file f u n 1. m fun1.m fun1.m Define the objective function
function f=fun1(x);
f=sum(x.^2)+8;
2) To write M M M file f u n 2. m fun2.m fun2.m Define nonlinear constraints
function [g,h]=fun2(x);
g=[-x(1)^2+x(2)-x(3)^2% Nonlinear inequality constraints
x(1)+x(2)^2+x(3)^3-20];
h=[-x(1)-x(2)^2+2% Nonlinear equality constraints
x(2)+2*x(3)^2-3];
3) Write main program
options=optimset('largescale', 'off');
[x,y]=fmincon('fun1',rand(3,1),[],[],[],[],...
zeros(3,1),[],'fun2',options)
Can be found x 1 = 0.5522 , x 2 = 1.2033 , x 3 = 0.9478 x_1=0.5522,x_2=1.2033,x_3=0.9478 x1=0.5522,x2=1.2033,x3=0.9478 when , minimum value y = 10.6511 y=10.6511 y=10.6511.
边栏推荐
- 海思Hi3798MV100机顶盒芯片介绍[通俗易懂]
- One year is worth ten years
- 人生的开始
- ceph 原理
- Does digicert SSL certificate support Chinese domain name application?
- Microservice architecture practice: using Jenkins to realize automatic construction
- Win10系统使用pip安装juypter notebook过程记录(安装在系统盘以外的盘)
- JS delete substring in string
- Nexus簡介及小白使用IDEA打包上傳到Nexus3私服詳細教程
- 2022 interview questions
猜你喜欢
线性规划例题 投资的收益与风险
体验居家办公完成项目有感 | 社区征文
Experience home office, feel the completion of the project | community essay solicitation
ETH数据集下载及相关问题
福元医药上交所上市:市值105亿 胡柏藩身价超40亿
How to transfer business data with BorgWarner through EDI?
畅玩集团冲刺港股:年营收2.89亿 刘辉有53.46%投票权
ThreadLocal
剑指 Offer 27. 二叉树的镜像
Blog theme "text" summer fresh Special Edition
随机推荐
移动应用性能工具探索之路
选择 SAP Spartacus 作为 SAP Commerce Cloud Storefront 实现框架的五个理由
Amazon cloud technology community builder application window opens
对接保时捷及3PL EDI案例
Blog theme "text" summer fresh Special Edition
Flutter: 动作反馈
Idea2021.1 installation tutorial
OpenHarmony如何启动FA(本地和远程)
什么是敏捷开发流程
The difference of message mechanism between MFC and QT
Understand one article: four types of data index system
牛客 JS3 分隔符
Schoolbag novel multithreaded crawler [easy to understand]
What if the default browser cannot be set?
超卓航科上市:募资9亿市值超60亿 成襄阳首家科创板企业
牛客 JS3 分隔符
Ocio V2 reverse LUT
Eye of depth (III) -- determinant of matrix
Soul, a social meta universe platform, rushed to Hong Kong stocks: Tencent is a shareholder with an annual revenue of 1.28 billion
Eth data set download and related problems