matlab 符号 约束,非线性优化 问题约束函数带有符号的条件函数怎么写？

本帖最后由 houleilei 于 2020-11-3 13:20 编辑

syms theta z1 z2 a1 a2

z=[0.2,0.4,0.6,0.8];

p=[0.1,0.5,0.6,1.0];for i=1:4

z11(i)=z1-z(i);z22(i)=z2-z(i);

A1(i)=a1./p(i);A2(i)=a2./p(i);Z11(i)=z1(i)/p(i);

if Z11(i)~=0

fun1(i)=sqrt(1+A1(i).^2-2*A1(i).*cos(theta)+Z11(i).^2).*cos(theta);

fun2(i)=sqrt(1+A2(i).^2-2*A2(i).*cos(theta)+Z11(i).^2).*cos(theta);

FP1(i)=int(fun1(i),theta,0,pi);

FP2(i)=int(fun2(i),theta,0,pi);

if Z11(i)==0

if A1(i)>1

fun1(i)=(sqrt(1+A1(i).^2-2*A1(i).*cos(theta)+Z11(i).^2).*cos(theta));

FP1(i)=int(fun1(i),theta,0,pi);

else if A1(i)==1

FP1(i)=-1.042260218;

else if ((0

fun1(i)=((sqrt(1+A1(i).^2-2*A1(i).*cos(theta)+Z11(i).^2)));

FP1(i)=-pi/4*(1+log4)+int(fun1(i),theta0,pi);

end

end

end

end

if Z11(i)==0

if A2(i)>1

fun2(i)=(sqrt(1+A2(i).^2-2*A2(i).*cos(theta)+Z11(i).^2).*cos(theta));

FP2(i)=int(fun2(i),theta,0,pi);

else if A2(i)==1

FP2(i)=-1.042260218;

else if ((0

fun2(i)=((sqrt(1+A2(i).^2-2*A2(i).*cos(theta)+Z11(i).^2)));

FP2(i)=-pi/4*(1+log4)+int(fun2(i),theta0,pi);

end

end

end

end

end

end

FP1,FP2

请问一下大神，怎么表达出Fp1,FP2关于变量 z1,z2 a1,a2在点(z,p)处的函数表达式呢？

因为我最终要做的是一个非线性优化问题，这个Fp1,FP2关于变量 z1,z2 a1,a2在点(z,p)处的函数表达式是非线性优化问题的约束条件。