matlab线性规划实验报告

合集下载

MATLAB解决线性规划问题

MATLAB解决线性规划问题

运行环境:Windows+MATLAB解决问题:线性规划问题(特定题目)实验简述:MATLAB 可以高效、方便地解决线性规划问题。

线性规划是合理利用、调配资源的一种应用数学的方法。

它的基本思路就是在满足一定的约束条件下,使预定的目标达到最优。

它的研究内容可归纳为两个方面:一是系统的任务已定,如何合理筹划,精细安排,用最少的资源去实现这个任务:二是资源的数量已定,如何利用、分配,使任务完成得最多。

前者是求极小,后者是求极大。

线性规划是在满足企业内、外部的条件下,实现管理目标和极值问题,就是要以尽少的资源输入来实现更多的社会需要的产品的产出。

现在通过专门的数学MATLAB 软件,只要将模型中的目标函数系数、约束条件系数、不等关系输入计算机,就会很快算出结果。

在生活实践中,很多重要的实际问题都是线性的(至少能够用线性函数很好的近似表示),所以我们一般把这些问题化为线性的目标函数和约束条件进行分析,通常将目标函数和约束都是线性表达式的规划问题称为线性规划。

它的一般形式是:),,2,1(0..min 221122222121112121112211n i x b x a x a x a b x a x a x a b x a x a x a t s x c x c x c f i m n mn m m n n n n nn=>=⎪⎪⎩⎪⎪⎨⎧<=+++<=+++<=++++++= 也可以用矩阵形式来表示:0,..min>=<==x b Ax t s x c f T线性规划的可行解是满足约束条件的解;线性规划的最优解是使目标函数达到最优的可行解。

线性规划关于解的情况可以是:1、无可行解,即不存在满足约束条件的解;2、有唯一最优解,即在可行解中有唯一的最有解;3、有无穷最优解,即在可行解中有无穷个解都可使目标函数达到最优;4、有可行解,但由于目标函数值无界而无最优解。

一般求解线性规划的常用方法是单纯形法和改进的单纯形法,这类方法的基本思路是先求得一个可行解,检验是否为最优解;若不是,可用迭代的方法找到另一个更优的可行解,经过有限次迭代后,可以找到可行解中的最优解或者判定无最优解。

线性规划模型 实验报告

线性规划模型 实验报告
可设多个变量 x1~x16
M文件:
c=-[190,90,244,193,261,199,170,110,260,150,280,165,140,80,186,103]';
L=zeros(16,1);
H=[31,52,22,41,10,60,25,33,20,31,8,41,34,59,13,15];
烤制部有10座大烤炉,每座烤炉的容量是每天出140台,每台可容纳10个唐师面包或5个更大的宋赐面包。可以在一台上同时放两种面包,只需注意宋赐面包所占的空间是唐师面包的两倍。
调配部每天可以调配最多8000个唐师面包和5000个宋赐面包。有两个自动调配器分别用于两种面包的调配而不至于发生冲突。
田园公司决定找出这两种面包产品的最佳产量配比,即确定两种面包的日产量,使得在公司面包厂的现有生产条件下利润最高。
L=zeros(5,1);
[x,fmin]=linprog(c,[],[],Aeq,Beq,L,H);
Min=fmin+10400
X=[x(1) x(2) x(3) x(4) x(5)]
求解结果
Min =
7.5500e+003
X =
600.0000 270.0167 220.0167 350.0000 300.0000
所以
从Toronto和Detroit到Chicago运输的货物为 600 和 270
从从Toronto和Detroit到Buffalo运输的货物为 0 和 230
从Chicago运输到NewYork、Phila.和St.louis的货物为分别 220 、350 、300
剩余的将从Buffalo运往。
X =
Columns 1 through 12

用matlab求解线性规划问题

用matlab求解线性规划问题

用m a t l a b求解线性规划问题Company number:【0089WT-8898YT-W8CCB-BUUT-202108】实验四 用M A T L A B 求解线性规划问题一、实验目的: 了解Matlab 的优化工具箱,能利用Matlab 求解线性规划问题。

二、实验内容:线性规划的数学模型有各种不同的形式,其一般形式可以写为:目标函数: n n x f x f x f z +++= 2211m in约束条件: s n sn s s n n b x a x a x a b x a x a x a ≤+++≤+++221111212111这里nn x f x f x f z +++= 2211称为目标函数,j f 称为价值系数,T n f f f f ),,,(21 =称为价值向量,j x 为求解的变量,由系数ij a 组成的矩阵 称为不等式约束矩阵,由系数ij c 组成的矩阵 称为等式约束矩阵,列向量T n b b b b ),,,(21 =和T n d d d d ),,,(21 =为右端向量,条件0≥j x 称为非负约束。

一个向量Tn x x x x ),,,(21 =,满足约束条件,称为可行解或可行点,所有可行点的集合称为可行区域,达到目标函数值最大的可行解称为该线性规划的最优解,相应的目标函数值称为最优目标函数值,简称最优值。

我们这里介绍利用Matlab 来求解线性规划问题的求解。

在Matlab 中有一个专门的函数linprog()来解决这类问题,我们知道,极值有最大和最小两种,但求z 的极大就是求z -的极小,因此在Matlab 中以求极小为标准形式,函数linprog()的具体格式如下:X=linprog(f,A,b)[X,fval,exitflag,ouyput,lamnda]=linprog(f,A,b,Aeq,Beq,LB,UB,X0,options)这里X 是问题的解向量,f 是由目标函数的系数构成的向量,A 是一个矩阵,b 是一个向量,A ,b 和变量x={x1,x2,…,xn}一起,表示了线性规划中不等式约束条件,A ,b 是系数矩阵和右端向量。

MATLAB实验之线性规划问题求解

MATLAB实验之线性规划问题求解

封面作者:PanHongliang仅供个人学习桂林电子科技大学数学与计算科学学院实验报告实验室:实验日期:年月日x附录Ⅱ综合性、设计性实验报告格式桂林电子科技大学数学与计算科学学院综合性、设计性实验报告版权申明本文部分内容,包括文字、图片、以及设计等在网上搜集整理。

版权为潘宏亮个人所有This article includes some parts, including text, pictures, and design. Copyright is Pan Hongliang's personal ownership.用户可将本文的内容或服务用于个人学习、研究或欣赏,以及其他非商业性或非盈利性用途,但同时应遵守著作权法及其他相关法律的规定,不得侵犯本网站及相关权利人的合法权利。

除此以外,将本文任何内容或服务用于其他用途时,须征得本人及相关权利人的书面许可,并支付报酬。

Users may use the contents or services of thisarticle for personal study, research or appreciation, andother non-commercial or non-profit purposes, but at the same time, they shall abide by the provisions of copyright law and other relevant laws, and shall not infringe upon the legitimate rights of this website and its relevant obligees. In addition, when any content or service of this article is used for other purposes, written permission and remuneration shall be obtained from the person concerned and the relevant obligee.转载或引用本文内容必须是以新闻性或资料性公共免费信息为使用目的的合理、善意引用,不得对本文内容原意进行曲解、修改,并自负版权等法律责任。

MATLAB应用实验线性规划

MATLAB应用实验线性规划

A15
500
A14
A13
210
A12
S1~S7 钢管厂 铁路
火车站 公路 管道
450里程(km)
s1
s2 … si

s7
C1j
C1,15
A1
A2 …
Aj
… A15
—单位运费; —在第 厂提供的量; —第 地需要量; 求从运多少钢管到, 可使总运费最少. 决策变量: —从运到的钢管数量
a1
a2

ai
范例
加工奶制品的生产计划 程序
)<
() ?
)<
)<
范例
)
) ) ) .
桶牛奶生产, 桶生产,利润元。
)
) ) ) .

原料无剩余
“资源” 剩余为零的
种 资
时间无剩余
约束为紧约束(有 效约束)
源 加工能力剩余
结果解释 影子价格 最优解下“资源”增加单位时“效益”的 原)料增加单位, 利润增增量长 时间增加单位, 利润增长 加工能力增长不影响利润
课堂延伸
•大陆航空公司乘务组快速补救计划
运筹学解决方案: 大陆航空公司采用技术,研发了 决策支持系统,生成全局最优 或接近最优的乘务组补救方案.
效果:大陆航空公司自实施此计划以来,在每一次严重的 突发事件中,都能实现乘务组快速补救,取得价值数百万 美元的效益。 公司估计,在年系统帮助它节省约了万美 元。由于该系统,使美国航空业自年月日后开始复苏,从 那时起,五家航空公司已采用了类似的基于运筹学的软 件用于突发事件后航班时刻表、飞行员安排等的补救计 划。
70 30
20 20
S6
110

最优化方法实验

最优化方法实验

《最优化方法》实验报告实验序号:01 实验项目名称:线性规划及MATLAB应用《最优化方法》实验报告实验序号:02 实验项目名称:0.618黄金分割法的应用结果分析:根据以上结果可知,在区间[0,3]上,函数g(x)=x^3-2*x+1的最小值点在x=0.9271处,此时最小值为0。

第二题:P50 例题3.1程序:function [t,f]=golden3(a,b) %黄金分割函数的m文件t2=a+0.382*(b-a);f2=2*(t2)^2-(t2)-1;t1=a+0.618*(b-a); %按照黄金分割点赋值,更准确可直接算f1=2*(t1)^2-(t1)-1;while abs(t1-t2)>0.16; %判定是否满足精度if f1<f2a=t2;t2=t1;f2=f1;t1=a+0.618*(b-a);f1=2*(t1)^2-(t1)-1;elseb=t1;t1=t2;f1=f2;t2=a+0.382*(b-a);f2=2*(t2)^2-(t2)-1;endendt=(t1+t2)/2; %满足条件取区间中间值输出第四题:P64 T3程序:function [t,d]=newtow2(t0)t0=2.5;t=t0-(4*(t0)^3-12*(t0)^2-12*(t0)-16)/(12*(t0)^2-24*(t0)-12);k=1;T(1)=t;while abs(t-t0)>0.000005t0=t;t=t0-(4*(t0)^3-12*(t0)^2-12*(t0)-16)/(12*(t0)^2-24*(t0)-12); k=k+1;T(k)=t;endt1=t0;d=(t1)^4-4*(t1)^3-6*(t1)^2-16*(t1)+4;kTend运行结果:当x(0)=2.5当x(0)=3四.实验小结:1.通过这次实验,加深了对0.618法的理解。

2.在学习0.618法的过程中,又巩固了倒数、求解函数值等相关知识。

MATLAB数学实验报告

MATLAB数学实验报告

MATLAB数学实验报告姓名:李帆班级:机械(硕)21学号:2120104008第一次数学实验报告——线性规划问题一,实验问题1,某饲养场饲养动物出售,设每头动物每天至少需要700g蛋白质,30g矿物质,100mg 维生素。

现有五种饲料可供选择,各种饲料的每千克营养成分含量和单价如下表。

是确定既能满足动物生长的营养需要,游客是费用最省的选用饲料方案。

2,某工厂生产甲、乙、丙三种产品,单位产品所需工时分别为2、3、1个;单位产品所需原料分别为3、1、5公斤;单位产品利润分别为2、3、5元。

工厂每天可利用的工时为12个,可供应的原料为15公斤。

为使总利润为最大,试确定日生产计划和最大利润。

二,问题分析1,1)该题属于采用线性规划的方式求出最优解的数学问题。

该题有以下特点,1.目标函数有线性,是求目标函数的最小值;2.约束条件为线性方程组;3.未知变量都有非负限制。

1,2)求解该类问题的方法有图解法,理论解法和软件解法。

图解法常用于解变量较少的线性规划问题。

理论解法要构建完整的理论体系。

目前用于解线性规划的理论解法有:单纯形法,椭球算法等。

在此,我们采用单纯形法的MATLAB软件解法来求解该问题。

1,3)此题中,要求既要满足动物生长的营养需要,又要使费用最省,则使每种饲料的选用量为变量,以总费用的最小值为所求量,同时每种饲料的使用量要符合营养成分的要求。

1,4)在此,首先确定建立线性规划模型。

设饲料i选用量为xi公斤,i=1,2,3,4,5.则有模型:Minz=0.2x1+0.7x2+0.4x3+0.3x4+0.8x5s.t. {3x1+2x2+6x4+18x5>=700;x1+0.5x2+0.2x3+2x4+0.5x5>=300.5x1+x2+0.2x3+2x4+0.8x5>=100Xj>=0,j=1,2,3,4,5解之得:x1=x2=x3=0X4=39.74359X5=25.14603Zmin=32.435902, 1)该问题与第一题分析步骤相似,故只在此写出其线性规划模型Z=2x+3y+5z2x+3y+z<=123x+y+5z<=15三,程序设计流程图第一题:c=[0.2,0.7,0.4,0.3,0.8]A=[3,2,1,6,18;1,0.5,0.2,2,0.5;0.5,1,0.2,2,0.8;1,0,0,0,0;0,1 ,0,0,0;0,0,1,0,0;0,0,0,1,0;0,0,0,0,1]b=[700,30,100,0,0,0,0,0][x,fval]=linprog(c,-A,-b)c =0.2000 0.7000 0.4000 0.3000 0.8000A =3.0000 2.0000 1.0000 6.0000 18.00001.0000 0.5000 0.20002.0000 0.50000.5000 1.0000 0.2000 2.0000 0.80001.0000 0 0 0 00 1.0000 0 0 00 0 1.0000 0 00 0 0 1.0000 00 0 0 0 1.0000b =700 30 100 0 0 0 0 0Optimization terminated.x =0.0000-0.00000.000039.743625.6410fval =32.4359第二题c=[-2 -3 -5]A=[2 3 1;3 1 5]b=[12;15]lb=[0 0 0][x,Z,exitflag,output]=linprog(c,A,b,[],[],lb,[])将上述程序输入matlab。

第十一章 MATLAB求解线性规划问题-2014120809203072

第十一章 MATLAB求解线性规划问题-2014120809203072

3.完成实验报告。

设生产A产品 x1 件,生产B产品 x2 件, z 为所获利 润,我们将问题归结为如下的线性规划问题:
min {(7000x1 10000 x2 )}
8 x1 6 x2 380 s.t. 4 x1 8 x2 300 4 x 6 x 220 2 1
接着写出MATLAB程序如下: clear f=-[7000,10000]; A=[8,6;4,8;4,6]; b=[380,300,220]; [X,fval]=linprog(f,A,b)

这里 z f1 x1 f 2 x2 f n xn 称为目标函数, f j 称为 价值系数,f ( f1 , f 2 ,, f n )T 称为价值向量, x j 为求 解的变量,由系数aij 组成的矩阵
a11 A a s1 a1n asn

解决上述问题的MATLAB程序为: clear f=-[5,4,6]; A=[1,-2,1;3,2,4;3,2,0]; b=[20,42,30]; LB=[0;0;0]; [X,fval,exitflag,output,lambda]=linprog(f,A,b,[],[ ],LB)

求解线性规划问题已有一些成熟的方法,我们这里 介绍利用MATLAB来求解线性规划问题的求解。
在MATLAB中有一个专门的函数linprog()来解决这 类问题,我们知道,极值有最大和最小两种,但求 z的极大就是求-z的极小,因此在MATLAB中以求 极小为标准形式,函数linprog()的具体格式如下: X=linprog(f,A,b) [X,fval,exitflag,output,lamnda]=linprog(f,A,b,Ae q,Beq,LB,UB,X0,options)

数学建模 matlab求解线性规划实验报告

数学建模 matlab求解线性规划实验报告

实验三 线性规划程序: linprogc=[-0.4 -0.28 -0.32 -0.72 -0.64 -0.6];A=[0.01 0.01 0.01 0.03 0.03 0.03;0.02 0 0 0.05 0 0;0 0.02 0 0 0.05 0;0 0 0.03 0 0 0.08]; b=[850;700;100;900]; Aeq=[]; beq=[];vlb=[0;0;0;0;0;0]; vub=[];[x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub)Exam5:function f=fun3(x);f=-x(1)-2*x(2)+(1/2)*x(1)^2+(1/2)*x(2)^2实验目的2、掌握用数学软件包求解线性规划问题。

1、了解线性规划的基本内容。

例1 max 6543216.064.072.032.028.04.0x x x x x x z +++++= 85003.003.003.001.001.001.0..654321≤+++++x x x x x x t s 70005.002.041≤+x x 10005.002.052≤+x x 90008.003.063≤+x x 6,2,10 =≥j x jx0=[1;1];A=[2 3 ;1 4]; b=[6;5]; Aeq=[];beq=[]; VLB=[0;0]; VUB=[];[x,fval]=fmincon('fun3',x0,A,b,Aeq,beq,VLB,VUB)书 求下列非线性规划2221232212322123212223123min 8020..2023,,0x x x x x x x x x s t x x x x x x x +++⎧-+≥⎪++≤⎪⎪--+=⎨⎪+=⎪⎪≥⎩在Matlab 2013软件中输入如下程序: (i )编写M 文件fun1.m 定义目标函数function f=fun1(x); f=sum(x.^2)+8;(ii )编写M 文件fun2.m 定义非线性约束条件 function [g,h]=fun2(x); g=[-x(1)^2+x(2)-x(3)^2x(1)+x(2)^2+x(3)^3-20]; %非线性不等式约束 h=[-x(1)-x(2)^2+2x(2)+2*x(3)^2-3]; %非线性等式约束 (iii )编写主程序文件example2.m 如下:options=optimset('largescale','off');[x,y]=fmincon('fun1',rand(3,1),[],[],[],[],zeros(3,1),[], ... 'fun2', options)就可以求得当1230.5522 1.2033,,0.9478x x x ===时,最小值y =10.6511。

matlab实验报告

matlab实验报告

matlab实验报告实验报告:Matlab实验分析1. 实验目的本实验旨在通过Matlab软件完成一系列数值计算和数据分析的任务,包括绘制曲线、解方程、矩阵运算等,以加深对Matlab软件的理解和掌握。

2. 实验内容2.1 绘制函数曲线首先,我们通过在Matlab中输入函数的表达式来绘制函数曲线。

例如,我们可以输入y = sin(x)来绘制正弦函数的曲线。

另外,我们还可以设置曲线的颜色、线型和坐标轴范围等。

2.2 解方程接下来,我们使用Matlab来解方程。

对于一元方程,我们可以使用solve函数来求出方程的解。

例如,我们输入syms x; solve(x^2 - 2*x - 8)来解方程x^2 - 2x - 8 = 0。

而对于多元方程组,我们可以使用solve函数的向量输入形式来求解。

例如,我们输入syms x y; solve(x^2 + y^2 - 1, x - y - 1)来求解方程组x^2 + y^2 - 1 = 0和x - y - 1 = 0的解。

2.3 矩阵运算Matlab也可以进行矩阵运算。

我们可以使用矩阵相乘、相加和取逆等运算。

例如,我们可以输入A = [1 2; 3 4]和B = [5 6;7 8]来定义两个矩阵,然后使用A * B来计算它们的乘积。

3. 实验结果与分析在本实验中,我们成功完成了绘制函数曲线、解方程和矩阵运算等任务。

通过Matlab软件,我们可以快速、准确地进行数值计算和数据分析。

使用Matlab的高级函数和工具箱,我们可以更方便地处理复杂的数值计算和数据分析问题。

4. 实验总结通过本次实验,我们进一步加深了对Matlab软件的理解和掌握。

Matlab提供了丰富的函数库和工具箱,适用于各种不同的数值计算和数据分析任务。

在日常科研和工程实践中,Matlab是一个非常强大和方便的工具,可以帮助我们更高效地完成任务。

数学建模实验报告之线性规划

数学建模实验报告之线性规划

数学模型实验报告——线性规划专业:数学与应用数学L081姓名: XXX 学号: 08L1002106姓名: XXX 学号: 08L1002109姓名: XXX 学号: 08L1002112数学模型实验报告(线性规划)一、 实验目的:1、了解线性规划的基本内容。

2、掌握用数学软件包求解线性规划问题。

二、实验内容:1、用MATLAB 优化工具箱解线性规划 ;2、两个例题;3、实验作业。

三、内容分析:(一)用MATLAB 优化工具箱解线性规划1、模型: min z=cXb AX t s ≤..命令:x=linprog (c ,A ,b )2、模型: min z=cXb AX t s ≤..beq X Aeq =⋅命令:x=linprog (c ,A ,b ,Aeq, beq ) 注意:若没有不等式:b AX ≤ 存在,则令A=[ ],b=[ ].3、模型:min z=cX b AX t s ≤..beq X Aeq =⋅VLB ≤X ≤VUB命令:[1] x=linprog (c ,A ,b ,Aeq,beq, VLB ,VUB )[2] x=linprog (c ,A ,b ,Aeq,beq, VLB ,VUB, X 0)注意:[1] 若没有等式约束: beq X Aeq =⋅, 则令Aeq=[ ], beq=[ ]. [2]其中X 0表示初始点4、命令:[x,fval]=linprog(…) 返回最优解x及x处的目标函数值fval.例1 max 6543216.064.072.032.028.04.0x x x x x x z +++++=85003.003.003.001.001.001.0..654321≤+++++x x x x x x t s70005.002.041≤+x x 10005.002.052≤+x x 90008.003.063≤+x x 6,2,10=≥j x j解 :编写M 文件程序如下:c=[-0.4 -0.28 -0.32 -0.72 -0.64 -0.6]; A=[0.01 0.01 0.01 0.03 0.03 0.03;0.02 0 0 0.05 0 0; 0 0.02 0 0 0.05 0; 0 0 0.03 0 0 0.08]; b=[850;700;100;900]; Aeq=[]; beq=[];vlb=[0;0;0;0;0;0]; vub=[];[x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub)例2321436m in x x x z ++= 120..321=++x x x t s301≥x 5002≤≤x 203≥x解:编写M 文件程序如下: c=[6 3 4]; A=[0 1 0]; b=[50];Aeq=[1 1 1]; beq=[120]; vlb=[30,0,20];vub=[];[x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub)运行结果如下:Optimization terminated. (最优解为) x =1.0e+004 * 3.5000 0.5000 3.0000 0.0000 0.0000 0.0000 fval =-2.5000e+004(二)例题例1:任务分配问题:某车间有甲、乙两台机床,可用于加工三种工件。

matlab实验一实验报告

matlab实验一实验报告

matlab实验一实验报告实验一:Matlab实验报告引言:Matlab是一种强大的数学软件工具,广泛应用于科学计算、数据分析和工程设计等领域。

本实验旨在通过使用Matlab解决实际问题,探索其功能和应用。

一、实验目的本次实验的主要目的是熟悉Matlab的基本操作和常用函数,了解其在科学计算中的应用。

二、实验内容1. 数值计算在Matlab中,我们可以进行各种数值计算,包括基本的加减乘除运算,以及更复杂的矩阵运算和方程求解。

通过编写相应的代码,我们可以实现这些功能。

例如,我们可以使用Matlab计算两个矩阵的乘积,并输出结果。

代码如下:```matlabA = [1 2; 3 4];B = [5 6; 7 8];C = A * B;disp(C);```2. 数据可视化Matlab还提供了强大的数据可视化功能,可以将数据以图表的形式展示出来,更直观地观察数据的规律和趋势。

例如,我们可以使用Matlab绘制一个简单的折线图,来展示某个物体在不同时间下的位置变化。

代码如下:```matlabt = 0:0.1:10;x = sin(t);plot(t, x);xlabel('Time');ylabel('Position');title('Position vs. Time');```3. 图像处理Matlab还可以进行图像处理,包括图像的读取、处理和保存等操作。

我们可以通过Matlab对图像进行增强、滤波、分割等处理,以及进行图像的压缩和重建。

例如,我们可以使用Matlab读取一张图片,并对其进行灰度化处理。

代码如下:```matlabimg = imread('image.jpg');gray_img = rgb2gray(img);imshow(gray_img);```三、实验结果与分析在本次实验中,我们成功完成了数值计算、数据可视化和图像处理等任务。

数学实验——线性规划

数学实验——线性规划

实验5 线性规划分1 黄浩 43一、实验目的1.掌握用MATLAB工具箱求解线性规划的方法2.练习建立实际问题的线性规划模型二、实验内容1.《数学实验》第二版(问题6)问题叙述:某银行经理计划用一笔资金进行有价证券的投资,可供购进的证券以及其信用等级、到期年限、收益如下表所示。

按照规定,市政证券的收益可以免税,其他证券的收益需按50%的税率纳税。

此外还有如下限制:(1).政府及代办机构的证券总共至少要购进400万元;(2).所购证券的平均信用等级不超过1.4(信用等级数字越小,信用程度越高);(3).所购证券的平均到期年限不超过5年I.若该经理有1000万元资金,该如何投资?II.如果能够以2.75%的利率借到不超过100万元资金,该经理应如何操作?III.在1000万元资金情况下,若证券A的税前收益增加为4.5%,投资应否改变?若证券C的税前收益减少为4.8%,投资应否改变?模型转换及实验过程:I.设经理对于上述五种证券A 、B 、C 、D 、E 的投资额分别为:x 1、x 2、x 3、x 4、x 5(万元),全部到期后的总收益为z 万元。

由题目中的已知条件,可以列出约束条件为:{ x 2+x 3+x 4≥4002x 1+2x 2+x 3+x 4+5x 5≤1.4(x 1+x 2+x 3+x 4+x 5)9x 1+15x 2+4x 3+3x 4+2x 5≤5(x 1+x 2+x 3+x 4+x 5)x 1+x 2+x 3+x 4+x 5≤1000}而决策变量x =(x 1,x 2,x 3,x 4,x 5)T 的上下界约束为:x i ∈[0,1000]目标函数z =0.043x 1+0.027x 2+0.025x 3+0.022x 4+0.045x 5 将上述条件转变为matlab 的要求形式:使用matlab 解上述的线性规划问题(程序见四.1),并整理成表格:得出结论:当经理对A 、B 、C 、D 、E 五种证券分别投资218.18、0、736.36、0、45.45万元时,在全部收回时可得到29.836万元的税后收益,而且这种投资方式所得收益是最大的。

matlab数学实验六

matlab数学实验六

s.t.
0.02 x1 0.05 x4 700 0.02 x2 0.05 x5 100
0.03 x3 0.08 x6 900
xj 0 j 1,2, 6
解 用命令3,编写M文件xxgh1.m如下: c=[-0.4 -0.28 -0.32 -0.72 -0.64 -0.6]; A=[0.01 0.01 0.01 0.03 0.03 0.03;0.02 0 0 0.05 0 0;0 0.02 0 0 0.05 0;0 0 0.03 0 0 0.08]; b=[850;700;100;900]; Aeq=[]; beq=[]; vlb=[0;0;0;0;0;0]; vub=[];
也用于求模型3,其中x0表示初始值。 5、命令:[x,fval]=linprog(…)
返回最优解x及x处的目标函数值fval.
9
例 1 max
z 0.4 x1 0.28 x2 0.32 x3 0.72 x4 0.64 x5 0.6 x6 0.01x1 0.01x2 0.01x3 0.03 x4 0.03 x5 0.03 x6 850
10
[x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub) 结果: x= 1.0e+004 * 3.5000 0.5000 3.0000 0.0000 0.0000 0.0000
fval =
-2.5000e+004
11
例2 min f 10 x11 5x12 6 x13 4 x 21 8 x 22 12 x 23 s.t x11 x12 x13 60
x1 x2 x x n
c c1, c2 , , cn

matlab线性规划实验报告

matlab线性规划实验报告
见附录见附录见附录xmax??xmax??xmax??1300000?1300000?1300000?200000?200000?200000?maxf??maxf??maxf??49000e04?49000e04?49000e04?其中其中其中xmaxxmaxxmax表示当代求函数取最大值时表示当代求函数取最大值时表示当代求函数取最大值时x1x2x1x2x1x2的取值的取值的取值maxfma180];
lb=[0,0];
[xmax maxf]=linprog(f,A,b,[],[],lb);
xmax
maxf=-maxf
作业1
上机题实验报告
(一)实验内容:
调用MATLAB中的线性规划函数讲解讲义中的线性规划例题。
(二)实验代码:
见附录
(三)实验运行结果:
xmax =
130.0000
20.0000
maxf =
4.9000e+04
其中,xmax表示当代求函数取最大值时,x1,x2的取值,maxf表示代求函数最大值。
(四)实验总结:
调用MATLAB中linprog函数,求解线性规划问题,注意问题:
1、linprog函数用来求解代求函数的最小值;
2、linprog函数各参数意义。
附录:
clc
clear
%%线性规划函数linprog,求最小值,A,b为不等式约束,lb是取值上限
f=[-300,-500];
A=[1,1;1,2;0,3];

Matlab求解线性规划问题

Matlab求解线性规划问题

封面作者:Pan Hongliang仅供个人学习注:上机作业文件夹以自己的班级姓名学号命名,文件夹包括如下上机报告和Matlab程序。

上机报告模板如下:版权申明本文部分内容,包括文字、图片、以及设计等在网上搜集整理。

版权为潘宏亮个人所有This article includes some parts, including text, pictures, and design. Copyright is Pan Hongliang's personal ownership.用户可将本文的内容或服务用于个人学习、研究或欣赏,以及其他非商业性或非盈利性用途,但同时应遵守著作权法及其他相关法律的规定,不得侵犯本网站及相关权利人的合法权利。

除此以外,将本文任何内容或服务用于其他用途时,须征得本人及相关权利人的书面许可,并支付报酬。

Users may use the contents or services of this article for personal study, research or appreciation, and othernon-commercial or non-profit purposes, but at the same time, they shall abide by the provisions of copyright law and other relevant laws, and shall not infringe upon the legitimate rights of this website and its relevant obligees. In addition, when any content or service of this article is used for other purposes, written permission and remuneration shall be obtained from the person concerned and the relevant obligee.转载或引用本文内容必须是以新闻性或资料性公共免费信息为使用目的的合理、善意引用,不得对本文内容原意进行曲解、修改,并自负版权等法律责任。

数学实验报告——利用MALTAB进行线性规划

数学实验报告——利用MALTAB进行线性规划

数学实验报告——利用M A L T A B进行线性规划(总12页)--本页仅作为文档封面,使用时请直接删除即可----内页可以根据需求调整合适字体及大小--数学实验报告——利用MALTAB进行线性规划实验六线性规划一、债券投资㈠问题描述给定可供购进的证券以及其信用等级、到期年限、收益。

市政证券的收益可以免税,其他证券的收益需按50%的税率纳税。

此外还有以下限制:(1) 政府及代办机构的证券总共至少要购进400万元;(2) 所购证券的平均信用等级不超过1.4;(3) 所购证券的平均到期年限不超过5年。

1、若经理有1000万元资金,应如何进行投资?2、如果能以2.75%的利率借到不超过100万元资金,该经理应如何操作?3、在1000万元资金情况下,若证券A的税前收益增加为4.5%,投资应够改变?若证券C的税前收益减少为4.8%,投资应否改变?㈡简要分析本题是一个比较简单的线性规划+扰动分析问题,对所谓问题进行建模,可以得到线性规划如下:设分别购入A、B、C、D、E五种证券a、b、c、d、e万元。

于是对第1问有b+c+d≥4002a+2b+c+d+5e≤1.4(a+b+c+d+e)9a+15b+4c+3d+2e≤5(a+b+c+d+e)a+b+c+d+e≤100a,b,c,d,e≥0max f=4.3%a+5.4%×50%b+5%×50%c+4.4%×50%d+4.5%e对第2问,增设a1,b1,c1,d1,e1分别表示用借来的资金购买证券的金额,于是规划变为b+c+d+b1+c1+d1≥4002(a+a1)+2(b+b1)+c+c1+d+d1+5(e+e1)≤1.4(a+b+c+d+e+a1+b1+c1+d1+e1) 9(a+a1)+15(b+b1)+4(c+c1)+3(d+d1)+2(e+e1)≤5(a+b+c+d+e+a1+b1+c1+d1+e1)a+b+c+d+e≤1000a1+b1+c1+d1+e1≤100a,b,c,d,e,a1,b1,c1,d1,e1≥0max f=4.3%(a+a1)+5.4%×50%(b+b1)+5%×50%(c+c1)+4.4%×50%(d+d1)+4.5%(e+e1)−2.75%(a1+b1+c1+d1+e1)对第三问,仅需将第一问中的规划做一点修改即可。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
相关文档
最新文档