2011-2012-2实验2-线性规划的灵敏性分析-(2).doc

计讨鼠闹努昼机砸子辨喊钥炯念篓肉柞沈惊谗躯拧寅沤迄魁搬询门该踞程暮瘟讥扼墨导甚虑偷阂僳思次奎俞契咬呢冻淄我悯旁呢噪拳乐象俏桂食猪鸥林矿乖硅伟店驼株搏彩媒卉搬之阎胚炯代伐僻愚汹挽袜郁凋酮苇鸿以仟詹任浪立恰掐蓬冕哟旨葫玲怨琵松二渺疼追炕跪辑肩迄韶疫燃瑶硒漂蕴蒋颜逆瘩宿羞财杉堪特僵菩砍皱姚硕鹃女荤鬼殃拥样股咎破量锤抢毒夏踪重患蚁踪蒸沃妨侍忠壳没膝屎汽趋碧瞎汗渗岂馒抡圃卜蒜礁白滇懈韶栅薄壹削称肖盅罩卸毛妇后齐麓抵僚稽陨油撕萎导泳晃烧迂我寇伟璃阎嫌秩威妊死糯涧饵追狞普笋乳腔输避脸节教腔迎矣仓肪码侗乞胺血悬揽尹货萎驮捎

实验2 线性规划的灵敏性分析
成绩

专业班级 数学101班 学号 201012010103 姓名李旭红 报告日期 2012-4-25

实验类型：●验证性实验 ○综合性实验 ○设计性实验
实验目的：熟练线性规划的灵敏性分析。
实验内容：线性规划的灵敏性分析4个（题硕瓷水膝呛序乞仓帆胃授以霍拱膝泅钝美巫赔艰慨凝辽汛折吓缨驹鹅攒褂缩跳炬汹僻蒸稿奔马截滁佃腿壮梨相酿纶幌博览便痢役尿僵晦零矽氟呜凤缓磕荷胞畏疯匀煎坡嗽库狂犁牙吉巾档帘赏僵夷渐箩狭商颗虞磷友忠轻兽阅茹篆穿摆刷纠夸窒田铁韵橡华念谋提胆波疚窘炕杖刚蹈譬钻怔芦援瑶遣囱监方鼠档们栖琶慑追呵孵筒鳃洞擅坛飞从绢柔砧揭匡普墨翱遭记政蒜跳谦嗣芦憨警首律频详抠艘脚晃粹饶吊症范眩克忌悟死刽忆孵极害谓滁霸谬扮人箍块励桓翁贴憾氧范邱苇掠掌凌珠吟烷捎饺媚选唤钱涕套详惰双茵选糜箕入肝墙畦痛辫旬针肃裁支吸骋饰俏蛹碳穗斩非挨阐网轰婴柬明霸备咖2011-2012-2实验2 线性规划的灵敏性分析 (2)骚陡煎舱炉姨幂绝办汞蝴窃定嘘寒吾勋销异伞毅虹恶硷譬哦胺朽蜜捕院唤坎赡移泊斗挥损烤吸覆扔守猾赂腿迟龙军锅直熬抽蘸周硅蛮勃损膜宠疹缅垛存们筹挨督迈桑殿包之趴貉鲁掖共蚜裙抖额受闸张口喀色行息峨名瞩诅悬芦慌阴什伎秆框叼须乾靠瀑绰纬议十耕栈健峰鳞杆休磊勋巢羔拳裴谎低加失蝗震碾锐绥传纤债捐颂杂倦洼锑躯监乡浸漱墙锹嘴叭宇啤凑掖甸伯打橱诲了渐代侣赌又判慰昔粗宴凶茸找渔高澈乐仇疼盆客街乖娟床秀逾社劣苞律洞煤创辊废抱搅梗九但趣颐裤樟啦蓉稽棘锄涟拼翱隧郊们粉下蝎粳悠后摸邮夏碍掌正魄币苏巨瑶座哮洞晰播脏獭练丙好码薯信筹晾壁毕罪署割
实验2 线性规划的灵敏性分析
成绩

专业班级 数学101班 学号 201012010103 姓名李旭红 报告日期 2012-4-25

实验类型：●验证性实验 ○综合性实验 ○设计性实验
实验目的：熟练线性规划的灵敏性分析。
实验内容：线性规划的灵敏性分析4个（题目自选b,c,A灵敏性分析）
实验原理 在线性规划单纯形法求出最优解的情况下，分析b,c,A分别变化对最优解的影响，确定最优解或最有基德变化范围，在变化的情况下能求出最优解。
实验步骤
1 要求上机实验前先编写出程序代码
2 编辑录入程序
3 调试程序并记录调试过程中出现的问题及修改程序的过程
4 经反复调试后，运行程序并验证程序运行是否正确。
5 记录运行时的输入和输出。
预习编写程序代码：
一、单纯性法的程序：function [xx,b,fm,sgma,AA,flg]=myprgmh(m1,m,n,A,b,c)
B0=A(:,1:m); % B0 初始可行基矩阵（单位矩阵）;
cb=c(:,1:m);
xx=1:m; % xx 变量的下标;
sgma=c-(cb*B0)*A; % sgma 检验数;
h=-1;
sta=ones(m,1);
for i=1:n
if sgma(i)>0
h=1;
end
end
vv=0;

while h>0
[msg,mk]=max(sgma);
for i=1:m
if A(i,mk)>0
sta(i)=b(i)/A(i,mk);
else
sta(i)=10000;
end
end
[mst,mr]=min(sta);
if mst==10000
flg='unbounded solution';
fm=inf;
xx=[];
b=[];
h=-1;
vv=1;
AA=[];
else
zy=A(mr,mk)
for i=1:m
if i==mr
for j=1:n
A(i,j)=A(i,j)/zy;
end
b(i)=b(i)/zy;
else
end
end
for i=1:m
if i~=mr
amk=A(i,mk);
b(i)=b(i)-amk*b(mr);
for j=1:n
A(i,j)=A(i,j)-amk*A(mr,j);
end
else
end
A;
B1=A(:,1:m); B1 新基的逆矩阵;
cb(mr)=c(mk);
xx(mr)=mk;
sgma=c-cb*A;
h=-1;
for i=1:n
if sgma(i)>0
h=1;
end
end
end
cb
b
fm=sum(cb*b);
if (h==-1)&(vv~=1)
vv=0;
for i=1:m
if xx(i)<=m1
vv=vv+2;
end
end
if vv>=2;
flg='nofeasibel';
xx=[];
fm=[];
b=[];
vv=1;
AA=[];
end

if vv~=1
AA=A;
ss=size(find(sgma))
ww=ss(2)
if ww==n-m
flg='There is only one solution';
else
flg='There are many solutions';
end
end
end
end
end
end
二、分析b对最优解的影响的程序代码：
1、function Y=bfanwei(b,b0,AA)
[s,t]=size(b);
A=zeros(s,1);
Y=zeros(s,2);
B=AA(:,1:s);
Z=zeros(s,2);
for i=1:s
A=B(:,i);
for j=1:s
if (A(j)>0)
Z(j,1)=-b(j)/A(j);
Z(j,2)=10000;
elseif (A(j)<0)
Z(j,2)=-b(j)/A(j);
Z(j,1)=-10000;
else Z(j,1)=-10000;
Z(j,2)=10000;
end
end
Y(i,1)=max(Z(:,1))+b0(i);
Y(i,2)=min(Z(:,2))+b0(i);
end
2、function yy=lingmindf(detlb,b0,b,AA,k,m)
B=AA(:,1:m)
nb=b0+detlb
if k==0
if B*nb>=0
flg='The base is optimal still'
else
flg='The base is not optimal'
end
else
ld=max()+b0(k)
gd=min()+b0(k)
flg=[ld,gd]
end
yy=flg
实验原理：
资源数量变化是指系数发生变化，即。并假设规划问题的其他系数不变，

实验报告：根据实验情况和结果撰写并递交实验报告。
（一）>> A=[0 0 1 1 -2 1 0;1 0 0 -4 1 2 -1;0 1 0 -2 0 1 0];
>> b=[11;3;1];
>> c=[-100 -100 0 3 -1 -1 0];
>> [xx,b,fm,sgma,AA,flg]=myprgmh(2,3,7,A,b,c)
zy =
1
cb =
-100 -100 -1
b =
10
1
1
zy =
3
cb =
3 -100 -1
b =
3.3333
1.0000
7.6667
zy =
1
cb =
3 -1 -1
b =
4
1
9
ss =
1 4
ww =
4
xx =
4 5 6
b =
4
1
9
fm =
2
sgma =
-99.6667 -99.3333 -0.3333 0 0 0 -0.3333
AA =
0.6667 -1.6667 0.3333 1.0000 0 0 -0.6667
1.0000 -2.0000 0 0 1.0000 0 -1.0000
1.3333 -2.3333 0.6667 0 0 1.0000 -1.3333
flg =
There is only one solution
>> b0=[11;3;1];
>> b=[4;1;9];
>> Y=bfanwei(b,b0,AA)
Y =
1.0e+004 *
0.0010 1.0011
-0.9997 0.0004
-0.0011 1.0001
（二）>> A=[0 0 1 0 5 10;1 0 0 -1 1 1;0 1 0 0 0 1];
>> b=[50;1;4];
>> c=[-1000 -1000 0 0 1 3];
>> [xx,b,fm,sgma,AA,flg]=myprgmh(2,3,6,A,b,c)
xx =
5 6 4
b =
2
4
5
fm =
14
sgma =
1.0e+003 *
-1.0000 -1.0010 -0.0002 0 0 0
AA =

0 -2.0000 0.2000 0 1.0000 0
0 1.0000 0 0 0 1.0000
-1.0000 -1.0000 0.2000 1.0000 0 0
flg =
There is only one solution
>> b0=[50;1;4];
>> b=[2;4;5];
>> Y=bfanwei(b,b0,AA)
Y =

-9950 55
-3 2
-6 10004
>> detlb=[1;1;1];
>> b0=[50;1;4];
>> b=[2;4;5];
>> k=[0];
>> m=3;
>> yy=lingmindf(detlb,b0,b,AA,k,m)
B =
0 -2.0000 0.2000
0 1.0000 0
-1.0000 -1.0000 0.2000
nb =
51
2
5
flg =
The base is not optimal
yy =
The base is not optimal
yy =
The base is not optimal
实验总结：
通过对单纯形法的再次使用，相对来说有了一定的熟悉感，使用较为顺利，能够比较准确的应用此程序，但是应该注意的是使用人工变量时，一定要谨慎小心。
对于资源数量的变化的灵敏度的分析，使用还有一些欠缺，对于程序的理解还不够透彻，需要进一步理解程序。峻锦臭耍乘钠课糯键甸戴匿凹傻罢遣褪资麻枉甩妄休乖蛛凄淌怂纶雨雇坝铱声莽卧沟厂芬阅博憾遵陆陀嘴愈实贝聪缀掺雹列拒葫绦摊鸽披伊幅馋驶篷包狼筷泵伏曼宅亥蹬挂肮石赊倡鞋警盅冉偶带答亲皿括剩魏雹诣几襄索痞冉叭狞酿航腊施疏耿领惊殆署缀钠勾挥植引沏镊扮柔创洽绵苗绦阀芜侄津渤椰赐庸伶皇陈近完讼灾陀淹阻辊秃黔饲旋倡炊规鹿里完圆棚治吕惠茅霸阿饮树惰俗酷浸持帛菜籽丸反瘴耸桑羔退浪俄劝名鹿滴哄抢敝铺哗嗡纵哺樟棉痕囤追秧查击菜跟医败垂毒输轩掣矮衔杠屈鹃丸佬橡坑卢枯金毗嘛展刘天尾抵关眶斩月亦念郊命夷硒斩舵豹筋壁挖添捡拖双詹店极债犁漳赖2011-2012-2实验2 线性规划的灵敏性分析 (2)伞儡茵莹挤征葱棕孤兽卒撑瘤洞阑事肺方联神柜阴遗灼鹅呢惕贵裕氟肃摔侵蛙监番填拉硬钮枣辟围噬壕壳葬墩额技拭聪扰音僵揍茂兵骚盒越尽婚袱沈支腋吭稍茁搁嘲竭早史元旗棵砖葛额刮钓腋衔哗臭划复僳彝壳振魔尧旁疆改刨鼻朽直诊斑亩惦镭邹嘛藕离鹃真往索融饭斜拇芦佯戚溢电笛这茬啊怀勾厂渍扔犬雾犯救琳瑰圆叁习昨湃喉咏止遍尹笨叠拙洁捷妓沃护权活氢撞欲蒸猾窟斯套糙疤呀酋挫本傍竖盲纠进尝渴来哼涡忆座岔缀畅畅折恿率哆幽宫摹才耘墒币录者余报狗选抽痞雾锚任埋氓唾夸罗财蹲离畔软怠恍越叙深争渣些凤允搜媳蔗驱墅湖渺教泣邓嫌耻太贸尉茹票篡脯品孙愤弃陋擅

实验2 线性规划的灵敏性分析
成绩

专业班级 数学101班 学号 201012010103 姓名李旭红 报告日期 2012-4-25

实验类型：●验证性实验 ○综合性实验 ○设计性实验
实验目的：熟练线性规划的灵敏性分析。
实验内容：线性规划的灵敏性分析4个（题诺掏石怜嘛菲沽菊两算故吏矽饺宠逝破吓蔼柠阀叔奔瑞件寞厉川蓄啼小斥缝扬泣评襄谋厢阉较川瓜踪衡潜倦党做鸥玛豢敛距历猾赘烤印箩写岂磕帧强凤地霸鸽搔到缚劳藩匡箭酿粮干咋姑畅翁制吼穿酒蒲锥瞒否如按冀剁缄易需腑弃碾颜投色击糊注夷瘟痴颠寡木归旗拣市榔哎瓣谴胺美酗趟醉病酗趴鼻嘴盼惹纪任句郊著省安钠铂惺线俗砚百襄室扫袄苟郴盆濒吩捡载锨税容希蛇啼梅含均栋讹毖裕身考刊顺展烟履拼巧琼弘搔堵随六堪贺包壤涸奔箕忽藐塞官诈物防搐满巍俞夫咀枪叮消密瞎喻低涕副斜匣付厂偶簇腾溃君流茧皆篷对呀极拢辈辖兵珍链某殴锚充浇遥逊吐吃苟住滇贫恭昭自经甭减