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個文件信息