matlab讲稿·优化插值拟合.ppt

函数 fseminf 格式 [x,fval,flag,output,lambda] fseminf ‘fun’,x0,ntheta,seminfcon,A,b,Aeq,beq,lb,ub,options 说明： x0为初始估计值； fun为目标函数，其定义方式与前面相同； ntheta为半无限约束的个数； seminfcon用来确定非线性约束向量C和Ceq以及半无限约束的向量K1，K2，…，Kn，通过指定函数柄来使用，如： x fseminf @myfun,x0,ntheta,@myinfcon 先建立非线性约束和半无限约束函数文件，并保存为myinfcon.m： function [C,Ceq,K1,K2,…,Kntheta,S] myinfcon x,S %S为向量w的采样值 % 初始化样本间距 if isnan S 1,1 ， S … % S 有ntheta行2列 end w1 … %计算样本集 w2 … %计算样本集 … wntheta … % 计算样本集 K1 … % 在x和w处的第1个半无限约束值 K2 … %在x和w处的第2个半无限约束值 … Kntheta … %在x和w处的第ntheta个半无限约束值 C … % 在x处计算非线性不等式约束值 Ceq … % 在x处计算非线性等式约束值 如果没有约束，则相应的值取为“[ ]”，如Ceq [] 例6-10 求下面一维情形的最优化问题 sub.to 解：将约束方程化为标准形式： 先建立非线性约束和半无限约束函数文件，并保存为mycon.m： function [C,Ceq,K1,K2,S] mycon X,S % 初始化样本间距： if isnan S 1,1 , S [0.2 0; 0.2 0]; end % 产生样本集： w1 1:S 1,1 :100; w2 1:S 2,1 :100; % 计算半无限约束： K1 sin w1*X 1 .*cos w1*X 2 - 1/1000* w1-50 .^2 -sin w1*X 3 -X 3 -1; K2 sin w2*X 2 .*cos w2*X 1 - 1/1000* w2-50 .^2 -sin w2*X 3 -X 3 -1; % 无非线性约束： C [ ]; Ceq [ ]; % 绘制半无限约束图形 plot w1,K1,-,w2,K2,: ,title Semi-infinite constraints 然后在MATLAB命令窗口或编辑器中建立M文件： fun sum x-0.5 .^2 ; x0 [0.5; 0.2; 0.3]; % Starting guess [x,fval] fseminf ’fun‘,x0,2,@mycon 结果为： x 0.6673 0.3013 0.4023 fval 0.0770 [C,Ceq,K1,K2] mycon x,NaN ; % 利用初始样本间距 max K1 ans -0.0017 max K2 ans -0.0845 例6-11 求下面二维情形的最优化问题 sub.to 初始点为x0 [0.25, 0.25, 0.25] 解：先建立非线性和半无限约束函数文件，并保存为mycon.m： function [C,Ceq,K1,S] mycon X,S % 初始化样本间距： if isnan s 1,1 , s [2 2]; end % 设置样本集 w1x 1:s 1,1 :100; w1y 1:s 1,2 :100; [wx, wy] meshgrid w1x,w1y ; % 计算半无限约束函数值 K1 sin wx*X 1 .*cos wx*X 2 -1/1000* wx-50 .^2 -sin wx*X 3 -X 3 +… sin wy*X 2 .*cos wx*X 1 -1/1000* wy-50 .^2-sin wy*X 3 -X 3 -1.5; % 无非线性约束 C [ ]; Ceq [ ]; %作约束曲面图形 m surf wx,wy,K1,edgecolor,none,facecolor,interp ; camlight headlight title Semi-infinite constraint drawnow 然后在MATLAB命