代码搜索结果

x=linspace(-6,6,100); y=1./(x.^2+1); x1=linspace(-6,6,21); y1=1./(x1.^2+1); plot(x,y,x1,y1,x1,y1,'o','LineWidth',1.5), gtext('n=20'),

load xyz xi=0:50:5600; yi=0:50:4800; zi=interp2(x,y,z,xi,yi','nearest'); surfc(xi,yi,zi) rotate3d xlabel('X'),ylabel('Y'),zlabel('Z') hold on; plot3(2000,2800,1100,'c.',0,800,650,'r.',4000,2

function y=lagr1(x0,y0,x) n=length(x0); m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j));

clear x0=[-1;1]; A=[];b=[]; Aeq=[1 1];beq=[0]; vlb=[];vub=[]; [x,fval,exitflag,output]=fmincon('fun4',x0,A,b,Aeq,beq,vlb,vub,'mycon')

clear %x0=[3 5 0 7 0 1 0 0 4 0 6 10 5 1 2 7]'; %x0=[ 3.0000 5.0000 0.0707 7.0000 0 0.9293 0 0 3.9293 0 6.0000 10.0707 6.3875 4.3943

function f=liaoch(x) a=[1.25 8.75 0.5 5.75 3 7.25]; b=[1.25 0.75 4.75 5 6.5 7.75]; d=[3 5 4 7 6 11]; e=[20 20]; f1=0; for i=1:6 s(i)=sqrt((x(13)-a(i))^2+(x(14)-b(i))^2); f1=s(i)*x(i

H=[1 -1; -1 2]; c=[-2 ;-6]; A=[1 1; -1 2]; b=[2;2]; Aeq=[]; beq=[]; VLB=[0;0];VUB=[]; [x,z]=quadprog(H,c,A,b,Aeq,beq,VLB,VUB)

x0=[3;2.5]; VLB=[0 0];VUB=[5 10]; [x,fval,exitflag,output]=fmincon('fun',x0,[],[],[],[],VLB,VUB,'mycon2')

c=[-2 -1]; A=[6 5;6 -5]; b=[40.25 9.75]; vlb=[2 2]; vub=[4 3]; x=linprog(c,A,b,[],[],vlb,vub) f=c*x g1=25-x(1)^2-x(2)^2 g2=7-x(1)^2+x(2)^2