求对偶单纯形法的源程序,C的或者matlab的均可.
来源:学生作业帮 编辑:百度作业网作业帮 分类:综合作业 时间:2024/06/16 14:21:43
求对偶单纯形法的源程序,C的或者matlab的均可.
![求对偶单纯形法的源程序,C的或者matlab的均可.](/uploads/image/z/18752558-14-8.jpg?t=%E6%B1%82%E5%AF%B9%E5%81%B6%E5%8D%95%E7%BA%AF%E5%BD%A2%E6%B3%95%E7%9A%84%E6%BA%90%E7%A8%8B%E5%BA%8F%2CC%E7%9A%84%E6%88%96%E8%80%85matlab%E7%9A%84%E5%9D%87%E5%8F%AF.)
function [val, X, exitFlag] = simplex(A, b, B, c)
%function [val, X, exitFlag] = simplex(A, b, B, c)
% simplex.m
% 单纯形法求解线性规划
% 函数输入各项为标准型对应值
% val 最优解的值
% X 最优解 1*n
% exitFlag 求解结果
% 1: succeed
% 0: unbounded
% -1: iterate 100 times
% A 系数矩阵 m*n
% b 初始基可行解 m*1
% B 初始基可行解的索引 m*1
% c 价格系数 1*n
[m, n] = size(A);
x = 1:m;
y = 1:n;
d = zeros(m,1);
X = zeros(1,n);
c1 = c;
%计算判别数
for kk = 1:m
c1 = c1 - A(kk,:).*c(B(kk));
end
%最多迭代100次
for count=1:100
[cMin, cIndex] = min(c1);
if (cMin >= 0)
exitFlag = 1;
break;
end
[dMin, dIndex] = minD(b, A(:,cIndex));
if (isinf(dMin))
exitFlag = 0;
break;
end
B(dIndex) = cIndex;
major = A(dIndex, cIndex);
%矩阵A,b,c的迭代更新
A(dIndex,:) = A(dIndex,:)./major;
b(dIndex) = b(dIndex)/major;
temp = b(dIndex);
b(x) = b(x) - A(x,cIndex)*b(dIndex);
b(dIndex) = temp;
temp = A(dIndex,:);
A(x,y) = A(x,y) - A(x,cIndex)*A(dIndex,y);
A(dIndex,:) = temp;
c1 = c1 - c1(cIndex)*A(dIndex,:);
end
if count==100
exitFlag = -1;
end
X(B(x)) = b(x);
val = sum(X.*c);
%function [val, X, exitFlag] = simplex(A, b, B, c)
% simplex.m
% 单纯形法求解线性规划
% 函数输入各项为标准型对应值
% val 最优解的值
% X 最优解 1*n
% exitFlag 求解结果
% 1: succeed
% 0: unbounded
% -1: iterate 100 times
% A 系数矩阵 m*n
% b 初始基可行解 m*1
% B 初始基可行解的索引 m*1
% c 价格系数 1*n
[m, n] = size(A);
x = 1:m;
y = 1:n;
d = zeros(m,1);
X = zeros(1,n);
c1 = c;
%计算判别数
for kk = 1:m
c1 = c1 - A(kk,:).*c(B(kk));
end
%最多迭代100次
for count=1:100
[cMin, cIndex] = min(c1);
if (cMin >= 0)
exitFlag = 1;
break;
end
[dMin, dIndex] = minD(b, A(:,cIndex));
if (isinf(dMin))
exitFlag = 0;
break;
end
B(dIndex) = cIndex;
major = A(dIndex, cIndex);
%矩阵A,b,c的迭代更新
A(dIndex,:) = A(dIndex,:)./major;
b(dIndex) = b(dIndex)/major;
temp = b(dIndex);
b(x) = b(x) - A(x,cIndex)*b(dIndex);
b(dIndex) = temp;
temp = A(dIndex,:);
A(x,y) = A(x,y) - A(x,cIndex)*A(dIndex,y);
A(dIndex,:) = temp;
c1 = c1 - c1(cIndex)*A(dIndex,:);
end
if count==100
exitFlag = -1;
end
X(B(x)) = b(x);
val = sum(X.*c);
求Matlab中regress函数的源程序代码
急求matlab在图像处理方面的源程序
求理想点法matlab源程序,或熵权理想点法matlab源程序代码都可以,有其一就行的
数据结构C语言括号的检验源程序
求MATLAB图像分割源程序 最好简单点带说明的
MATLAB中的源程序注释 解释每一行的意思
求单纯形法的算法?
求sin(pi*t)/(pi*t)的fourier变换,用matlab或者直接算均可,最好有程序或者计算步骤
管理运筹学 单纯形法的灵敏度分析与对偶问题,b1在什么范围内,其对偶价格不变 怎么算啊
求matlab大神给下面一个源程序给个详细的注释,就是程序每一步的操作,速求
求matlab 用Bezier曲线绘制一个兔子的图案~后天就要交了!带源程序代码!
matlab 累积求和,求源程序