sobol灵敏度分析matlab_灵敏度分析

实验二、线性规划的灵敏度分析

实验目的

1. 线性规划求解的单纯形法的灵敏度分析的编程实现

2．掌握使用matlab、Lingo、Excel的规划求解功能求解，并利用“敏感性报告”进行分析。

(二)实验内容

课本例1 解的灵敏度分析

(1)：调用单纯形程序：

function [x,z,flg,sgma]=simplexfun(A,A1,b,c,m,n,n1,cb,xx)

% A,b are the matric in A*x=b

% c is the matrix in max z=c*x

% A1 is the matric in simplex table

% m is the numbers of row in A and n is the con number in A

% n1 is the nubers of artificial variables,and artificial variables are default as the last % n1 variables in x.

% cb is the worth coefficient matrix for basic variables

% xx is the index matrix for basic variables

% B1 is the invers matrix for the basic matrix in simplex table.The initial

% matrix is default as the last m con in the matrix A.

x=zeros(n,1);

z=0;

B1=A1(:,n-m+1:n);

sgma1=c-(cb*B1)*A;

[masg,kk]=max(sgma1);

k=kk(1);

flg=0;

ll=0;

while (masg>0)&&(ll<20)

ll=ll+1;

thita=1000+zeros(m,1);

for i=1:m

if A1(i,k)>0

thita(i)=A1(i,k)\b(i);

end

[r8,c8]=find(thita>999);

if sum(c8)

[mith,rr]=min(thita);

r=rr(1);

aa=A1(r,k);

for i=1:m

if i==r

b(r)=b(r)/aa;

for j=1:n

A1(r,j)=A1(r,j)/aa ;

end

for i=1:m

if i~=r

cc=A1(i,k)

b(i)=b(i)-b(r)*cc;

for j=1:n

A1(i,j)=A1(i,j)-A1(r,j)*cc;

end

cb(r)=c(k);

xx(r)=k;

B1=A1(:,n-m+1:n);

sgma1=c-(cb*B1)*A;

[masg,kk]=max(sgma1);

k=kk(1);

thita=100+zeros(m,1);

else

flg=3;

masg=-1;

x='unbound solution';

z='inf';

end

if flg~=3

if n1==0

sgma1=c-(cb*B1)*A

[rc,ccc]=find(sgma1

if sum(rc)==n-m

flg=1;

else

flg=2;

end

x=zeros(n,1);

for i=1:m

x(xx(i))=b(i);

end

z=c*x;

sobol灵敏度分析matlab_灵敏度分析_使用MATLAB编写.doc相关推荐

matlab 高斯迭代代码_用Matlab编写你的第一个有限元_《数值计算与程序设计》系列课程...
本课主要分析了有限元的基本理论, 不同于其他学习资料的是本文从数学的角度重新阐述了有限元分析的基本步骤,如通解与特解的讲解,并以一个偏微分方程实例,从基础插值出发,详细地讲解了有限元数学推导的每一步过 ...
matlabff2函数_罚函数法MATLAB程序.doc
罚函数法MATLAB程序一.进退法.0.618法.Powell法.罚函数法的Matlab程序设计罚函数法(通用)function y ff x,k y -17.86*0.42*x 1 / 0.8+0 ...
分布分析和分组分析_如何通过群组分析对用户进行分组并获得可行的见解
分布分析和分组分析数据分析 (DATA ANALYSIS) Being a regular at a restaurant is great. 乙 eing定期在餐厅是伟大的. When I sta ...
数据可视化分析票房数据报告_票房收入分析和可视化
数据可视化分析票房数据报告 Welcome back to my 100 Days of Data Science Challenge Journey. On day 4 and 5, I work ...
计算机黑屏的原因及解决办法,导致电脑黑屏的两个常见的原因分析与解决办法_电脑故障...
导致电脑黑屏的两个常见的原因分析与解决办法_电脑故障 2017年04月21日阅读 192 电脑黑屏故障的原因有很多种,有时很简单的一个差失就会导致,找到原因后才恍然大悟.下面就是一个电脑黑屏的案例分 ...
r语言 bsda包_使用R语言creditmodel包进行Vintage分析或留存率分析
1 什么是vintage分析? Vintage分析(账龄分析法)被广泛应用于信用卡及信贷行业,这个概念起源于葡萄酒,即不同年份出产的葡萄酒的品质有差异,那么不同时期开户或者放款的资产质量也有差异,其核 ...
家用电器用户行为分析与事件识别_用户行为分析系统
随着互利网时代多元化商业模式和细分市场的快速发展,众多行业面临着高昂的获客成本和高流失率的双重挑战,如何精准定位获客渠道?如何精细化运管,挖掘更大的用户价值?如何洞察客户流失的真实原因?如何减少用户流 ...
abaqus分析用户手册单元卷_作用卷、分析卷、材料卷三件套，让你也能熟练应用Abaqus...
<Abaqus分析用户手册--指定条件.约束与相互作用卷>是"Abaqus用户手册大系"中的一册,包括指定条件.约束与相互作用三个部分.指定条件部分对各种物理过程中涉及 ...
如何分析案件的性质_律师如何综合分析一个案件
当事人找到律师,一般会简单描述下整个案情,然后提供案件的有关材料.律师拿到这些材料后,就需要对案件材料进行综合性分析.综合分析案件并不是一件容易的事情,而是一种综合能力的体现,对律师的整体能力要求比较 ...
鱼骨图分析法实际案例_【管理工具详解】鱼骨图分析法
第一部分鱼骨头分析法一.鱼骨图分析法的由来鱼骨图是由日本管理大师石川馨先生所发明出来的,故又名石川图.鱼骨图是一种发现问题"根本原因"的方法,它也可以称之为"Ish ...

sobol灵敏度分析matlab_灵敏度分析_使用MATLAB编写.doc

sobol灵敏度分析matlab_灵敏度分析_使用MATLAB编写.doc相关推荐

最新文章

热门文章