[资料]

matlab单纯形法实现（GUI based on Matlab）

2012-3-13 15:26:09 10463 matlab

感觉matlab的GUI学起来比较容易上手（Matlab把布局代码放入.fig文件中，而.m文件中只是记录函数代码，这样减轻了代码量，让人感觉思路更清晰；相对比VC++把布局代码和函数代码都放在一起，让人看着不方便）。下面是我做的单纯形法的GUI。

游客
，如果您要查看本帖隐藏内容请回复

内容清单:

文件清单以及说明：
0.jpg,2.jpg,guidetu,jpg,GUItu.jpg,m.jpg,mgui.jpg,mm,jpg,r.jpg,sta.jpg,std.jpg
是论文中相关插入图片的源文件。
finalsimple.m,LPtostandard.m,originalsimple.m,originalsimpleM.m,twosteps.m
是论文中相关m文件的源代码
SimplicialMethod.m是单纯形法GUI的源代码
SimplicialMethod.fig是GUI的界面文件
SimplicialMethod_stdapp文件夹里装的是独立于Matlab的单纯形法的GUI应用程序，使用方法见该文件夹里的readme.txt
部分代码如下
SimplicialMethod.m 界面代码：

function varargout = SimplicialMethod(varargin)
% SIMPLICIALMETHOD M-file for SimplicialMethod.fig
% SIMPLICIALMETHOD, by itself, creates a new SIMPLICIALMETHOD or raises the existing
% singleton*.
%
% H = SIMPLICIALMETHOD returns the handle to a new SIMPLICIALMETHOD or the handle to
% the existing singleton*.
%
% SIMPLICIALMETHOD('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in SIMPLICIALMETHOD.M with the given input arguments.
%
% SIMPLICIALMETHOD('Property','Value',...) creates a new SIMPLICIALMETHOD or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before SimplicialMethod_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to SimplicialMethod_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help SimplicialMethod
% Last Modified by GUIDE v2.5 29-Apr-2008 21:31:56
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @SimplicialMethod_OpeningFcn, ...
'gui_OutputFcn', @SimplicialMethod_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_maiNFCn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
% --- Executes just before SimplicialMethod is made visible.
function SimplicialMethod_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to SimplicialMethod (see VARARGIN)
% Choose default command line output for SimplicialMethod
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes SimplicialMethod wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% --- Outputs from this function are returned to the command line.
function varargout = SimplicialMethod_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
function input_Callback(hObject, eventdata, handles)
% hObject handle to input (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of input as text
% str2double(get(hObject,'String')) returns contents of input as a double
% --- Executes during object creation, after setting all properties.
function input_CreateFcn(hObject, eventdata, handles)
% hObject handle to input (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
% --- Executes on button press in clear.
function clear_Callback(hObject, eventdata, handles)
% hObject handle to clear (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
set(handles.input,'String','0');
set(handles.shuchu,'String','0');
guidata(hObject,handles);
% --- Executes on button press in M.
function M_Callback(hObject, eventdata, handles)
% hObject handle to M (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
in = get(handles.input,'String');
B = str2num(in);
[x,fval,flag,iteration]=finalsimple(B)
result = sprintf('解的状态:nflag = %dn(0 唯一解，1 无穷多，2 无界解，3 无解)n最优解为:nx = ( %s )n最优值为:nfval = %.2fn单纯形法迭代次数:niteration = %dn',flag,num2str(x),fval,iteration);
set(handles.shuchu,'String',result);
guidata(hObject,handles);
% --- Executes on button press in twosteps.
function twosteps_Callback(hObject, eventdata, handles)
% hObject handle to twosteps (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
in = get(handles.input,'String');
B = str2num(in);
[x,fval1,fval2,flag,iteration]=twosteps(B);
result = sprintf('第一阶段:nfval1 = %.2fn第二阶段:n解的状态:nflag = %dn(0 唯一解，1 无穷多，2 无界解，3 无解)n最优解为:nx = ( %s )n最优值为:nfval = %.2fn单纯形法迭代次数:niteration = %dn',fval1,flag,num2str(x),fval2,iteration);
set(handles.shuchu,'String',result);
guidata(hObject,handles);

复制代码

originalsimple.m [原始单纯形法]

function [x,fval,flag,iteration]=originalsimple(C,A,b,XB)
%原始单纯形法（需直接给出初始的基变量）
%Programmed by Liyang(faruto's Studio~!) BNU MATH
%Email:liyangbnu@mail.bnu.edu.cn QQ:516667408
%last modified 2008.4.27
%求解标准型线性规划:max C*X;s.t. A*X=b (b>=0);X>=0;
%输入：C是n维行向量，A是m*n的系数矩阵，b是m维列向量，XB承装初始基变量的下标
%输出：x最优解（如果有的话），fval最优值，flag解的状态说明，interation求解时的循环次数
%flag 最终解的状态说明：
%flag = 0 LP converged to a solution x
%flag = 1 Inf feasible solutions
%flag = 2 LP is unbounded
%flag = 3 No feasible point was found
[m,n] = size(A);
iteration = 0;
pflag = 1;
while pflag
iteration = iteration +1;
flag = 0;
%%
sigma = zeros(1,n);
for col = 1:n
temp = 0;
for row = 1:m
temp = temp+C(XB(row))*A(row,col);
end
sigma(col) = C(col)-temp;
end
%%
if sigma<=0
x = zeros(1,n);
for row = 1:m
x(XB(row)) = b(row);
end
fval = C*x';
for row = 1:m
for temp = n+1:n+m
if XB(row)==temp && x(temp)~=0
flag = 3;
pflag = 0;
break;
end
end
if flag == 3
break;
end
end
if flag == 3
break;
end
for col = 1:n
tflag = 0;
for row = 1:m
if col == XB(row)
tflag = 1;
break;
end
end
if tflag == 0
if sigma(col) == 0
flag = 1;
pflag = 0;
break;
end
end
end
if flag == 1
break;
end
if flag == 0;
pflag = 0;
fval;
break;
end
else
for col = 1:n
if sigma(col)>0 & A(:,col)<=0
flag = 2;
pflag = 0;
break;
end
end
if flag == 2
break;
end
%%
temp = 0;
for col = 1:n
if sigma(col)>temp
temp = sigma(col);
intobase = col;%入基变量的下标
end
end
theta = zeros(m,1);
for row = 1:m
if A(row,intobase)>0
theta(row) = b(row)/A(row,intobase);
end
end
temp = Inf;
for row = 1:m
if theta(row)>0 & theta(row)
temp = theta(row);
outbase = XB(row);%出基变量的下标
outrow = row;%出基变量在基变量中的位置
end
end
%%
b(outrow,1) = b(outrow,1)/A(outrow,intobase);
A(outrow,:) = A(outrow,:)/A(outrow,intobase);
for row = 1:m
if row ~= outrow
b(row) = b(row) - b(outrow)*A(row,intobase);
A(row,:) = A(row,:) - A(outrow,:)*A(row,intobase);
end
end
%%
for row = 1:m
if XB(row) == outbase;
XB(row) = intobase;
end
end
end
end
%%
if flag == 0
disp('LP converged to a solution x!');
x;
fval;
end
if flag == 1
disp('Inf feasible solutions!');
disp('one of the solutions is:');
end
if flag == 2
disp('LP is unbounded!');
end
if flag == 3
disp('No feasible point was found!');
end

复制代码

评分

参与人数 1	威望 +10	+10	积分 +10	收起理由
vvg	+ 10	+ 10	+ 10	多分享经验和资料

查看全部评分

举报 淘帖0 只看该作者

最近下载过的用户(3)

相关推荐

28条评论 发表评论只看该作者

vvg · 2012-3-13 16:01:49 沙发 Good！很好的资料

赞回复举报

717732 · 2012-7-18 17:21:49 板凳牛人啊，学了几遍单纯形法都不知道怎么用

赞回复举报

程宝威 · 2014-2-23 18:33:32 4^# 呵呵，看看学学

赞回复举报

程宝威 · 2014-2-23 18:34:56 5^# {:12:}{:12:}{:12:}

赞回复举报

hfutfjx · 2014-2-24 16:50:54 6^# 学习学习谢谢分享

赞回复举报

程宝威 · 2014-2-25 15:05:49 7^# {:11:}{:11:}{:11:}{:11:}

赞回复举报

JuneL · 2014-7-11 21:51:00 8^# 牛人啊，thank u for ur information

赞回复举报

432482 · 2014-9-1 21:27:02 9^# 谢谢楼主的分享，非常感谢

赞回复举报

txlapollo · 2014-9-20 00:39:36 10^# 感谢啊！！！！很好的学习资料！！！

赞回复举报

txlapollo · 2014-9-20 00:40:24 11^# 谢谢楼主！很好的学习资料感谢分享！！

赞回复举报

马志涛 · 2014-9-20 07:49:43 12^# 谢谢了

赞回复举报

90大川 · 2015-4-16 13:15:37 13^# 楼主牛逼啊，可能给个联系方式啊

赞回复举报

michaelpeng0429 · 2015-11-11 16:09:41 14^# 好东西！！找了好久了！！！感谢楼主！！

赞回复举报

doublesamuel · 2016-3-9 14:44:29 15^# 感谢楼主分享，多谢多谢

赞回复举报

dgbsdkga · 2016-3-28 21:43:19 16^# 薛素素额苏苏苏苏人员收入也是颜色呀实话实说

赞回复举报

hawkgreen · 2016-5-9 17:10:19 17^# 正在查询相关内容，希望版本和我用的版本兼容

赞回复举报

linchen0519 · 2016-6-7 19:57:17 18^# 看看.学习下多谢楼主

赞回复举报

黑暗公爵 · 2016-6-16 09:49:43 19^# matlab单纯形法实现（GUI based on Matlab）

赞回复举报

GameshZhang · 2016-10-16 20:47:33 20^# 不知道代码效率怎么样

赞回复举报

声明：本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人，不代表电子发烧友网立场。文章及其配图仅供工程师学习之用，如有内容图片侵权或者其他问题，请联系本站作侵删。侵权投诉

12 / 2 页下一页

发资料

精选推荐

Matlab中如何将不同的时间的数据进行同时间对比？

874 浏览 3 评论
使用simulink进行三相短路故障分析时，各参数应该如何设置

1547 浏览 1 评论
想请教一下图中是simulink的什么模块，需要这种三段斜率函数模块但没找到在哪

1619 浏览 1 评论
怎么修改matlab与comsol联合的代码啊

2548 浏览 1 评论
200圆！求助大佬给一份VSG并网和离网模式的simulink仿真

2801 浏览 0 评论

快速回复 返回顶部 返回列表

关注微信公众号

电子发烧友网

电子发烧友论坛

社区合作: 刘勇; 联系电话：15994832713; 邮箱地址：liuyong@huaqiu.com

社区管理: elecfans短短; 微信：elecfans_666; 邮箱：users@huaqiu.com

【优惠升级】华秋PCB首单最高立减100元，SMT免费贴片！！！

返回 Matlab论坛

12 / 2 页下一页

搜索历史

matlab单纯形法实现（GUI based on Matlab）

评分

评论

精选推荐

站长推荐 /7