实用最优化方法编程大作业

这是实用最优化方法课程的大作业程序，压缩文件里有6个文档并且是对应于大连理工大学编写的实用最优化方法这本书。

function [xkkykk]=cs（HcAebeAibixk）
%功能：利用约束集法求解一般二次规划问题
%输入：hessen阵H，等式约束系数Ae，不等式约束系数Ai
%输出：最优解xkk，最优函数值ykk
epsilon=1e-9;err=1e-6;
k=0;xkk=xk;n=length（xkk）;maxk=1000;
ne=length（be）;ni=length（bi）;lamk=zeros（ne+ni1）;
index=ones（ni1）;
%初始化
for i=1:ni
	if Ai（i:）*xkk>bi（i）+epsilon
		index（i）=0;
	end
end
while k	Aee=[];
	if ne>0
		Aee=Ae;
	end
	for j=1:ni
		if index（j）>0
			Aee=[Aee;Ai（j:）];
		end
	end
	gk=H*xkk+c;
	[m1n1]=size（Aee）;
	[dklamk]=qsubp（HgkAeezeros（m11））;
	if norm（dk）<=err
		y=0;
		if length（lamk）>ne
			[yjk]=min（lamk（ne+1:length（lamk）））;
		end
		if y>=0
			exitflag=0;
		else
			exitflag=1;
			for i=1:ni
				if index（i）&&（ne+sum（index（1:i）））==jk
					index（i）=0;
					break;
				end
			end
		end
		k=k+1;
	else
		exitflag=1;
		alpha=1.0;tm=1.0;
		for i=1:ni
			if （index（i）==0）&&（Ai（i:）*dk<0）
				tm1=（bi（i）-Ai（i:）*xkk）/（Ai（i:）*dk）;
				if tm1					tm=tm1;ti=i;
				end
			end
		end
		alpha=min（alphatm）;
		xkk=xkk+alpha*dk;
		if tm<1
			index（ti）=1;
		end
	end
	if exitflag==0
		break;
	end
	k=k+1;
end
ykk=fun4（xkk）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-11-08 11:46  优化方法编程大作业\
     目录           0  2012-11-08 11:38  优化方法编程大作业\4\
     文件        1188  2012-11-08 11:38  优化方法编程大作业\4\cs.m
     文件          65  2012-11-08 11:38  优化方法编程大作业\4\fun4.m
     文件         157  2012-11-08 11:38  优化方法编程大作业\4\function42.m
     文件         192  2012-11-08 11:38  优化方法编程大作业\4\qsubp.m
     目录           0  2012-11-08 18:16  优化方法编程大作业\6\
     文件        4608  2012-11-08 18:16  优化方法编程大作业\6\qpsubp.m
     文件        3299  2012-11-08 18:16  优化方法编程大作业\6\sqpm.m
     文件         220  2012-11-08 18:16  优化方法编程大作业\6\v6.m
     文件      217018  2012-10-24 15:21  优化方法编程大作业\homework.pdf
     目录           0  2012-11-08 11:46  优化方法编程大作业\作业1\
     文件         104  2012-11-08 17:45  优化方法编程大作业\作业1\fun1.m
     文件         122  2012-11-08 17:46  优化方法编程大作业\作业1\gfun1.m
     文件         670  2012-11-08 11:46  优化方法编程大作业\作业1\wolfe_powell.m
     目录           0  2012-11-08 11:46  优化方法编程大作业\作业2\
     文件         686  2012-11-08 11:46  优化方法编程大作业\作业2\fr.m
     文件         135  2012-11-08 11:46  优化方法编程大作业\作业2\fun2.m
     文件         123  2012-11-08 11:46  优化方法编程大作业\作业2\gfun2.m
     文件         872  2012-11-08 11:46  优化方法编程大作业\作业2\wolfe_powell1.m
     目录           0  2012-11-08 11:46  优化方法编程大作业\作业3\
     文件         671  2012-11-08 11:46  优化方法编程大作业\作业3\dfp.m
     文件         104  2012-11-08 11:46  优化方法编程大作业\作业3\fun3.m
     文件         105  2012-11-08 11:46  优化方法编程大作业\作业3\gfun3.m
     文件         835  2012-11-08 11:46  优化方法编程大作业\作业3\wolfe_powell2.m
     目录           0  2012-11-08 11:46  优化方法编程大作业\作业4\
     文件         179  2012-11-08 11:46  优化方法编程大作业\作业4\callqpact.asv
     文件         183  2012-11-08 11:46  优化方法编程大作业\作业4\callqpact.m
     文件         177  2012-11-08 11:46  优化方法编程大作业\作业4\callqpact1.asv
     文件         179  2012-11-08 11:46  优化方法编程大作业\作业4\callqpact1.m
     文件        1222  2012-11-08 11:46  优化方法编程大作业\作业4\qpact.m
............此处省略13个文件信息