?? pen2.m
字號(hào):
% Penalty II function% ------------------- %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function [fvec,J]=pen2(n,m,x,option)% Dimensions -> n=variable, m=2*n% Problem no. 24 % Standard starting point -> x=(1/2,......,1/2)% Minima -> f=9.37629...10^(-6) if n=4% f=2.93660...10^(-4) if n=10 % % 11/21/94 by Madhu Lamba %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function [fvec,J] = pen2(n,m,x,option) J=zeros(m,n);if (option==1 | option==3) fvec(1)=x(1)-0.2;end;if (option==2 | option==3) J(1,1)=1;end;if(n>=2) for i=2:n y(i)=exp(i/10)+exp((i-1)/10); if (option==1 | option==3) fvec(i)=sqrt(1.e-5)*(exp(x(i)/10)+exp(x(i-1)/10)-y(i)); end; if (option==2 | option==3) J(i,i) = sqrt(1.e-5)*exp(x(i)/10)*(1/10); J(i,i-1) = sqrt(1.e-5)*exp(x(i-1)/10)*(1/10); end; end; for i=n+1:(2*n-1) if (option==1 | option==3) fvec(i)=sqrt(1.e-5)*(exp(x(i-n+1)/10)-exp(-1/10)); end; if (option==2 | option==3) J(i,i-n+1) = sqrt(1.e-5)*exp(x(i-n+1)/10)*(1/10); end; end;end;if (option==1 | option==3) sum=0; for j=1:n sum=sum+(n-j+1)*x(j)^2; end; fvec(2*n)=sum-1; end;if (option==2 | option==3) for j=1:n J(m,j) = (n-j+1)*2*x(j); end;end;if (option==1 | option==3) fvec=fvec';end;%
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -