function?[]=fminemaxc（ab）
%調用的格式為fminemax1（ab），a=[0.30.69]b=[9.9111.17]
%由于全局變量只有兩個xy所以這里就引用兩個二維矩陣作為參變量
if?（a（1）<0.3）|（a（2）>0.69）|（b（1）<9.91）|（b（2）>11.17）
????disp（‘wrong?input‘）;
else
????runfun（ab）;
end
%以下為子函數
function?runfun（ab）
?y1=9.91;
?y2=10.54;
?y3=11.17;
?Rg=10;
r4=0.0033;
Dis42=0.69;
xtem=size（a（1）:0.01:a（2））;
ytem=size（b（1）:0.01:b（2））;
Utem=zeros（6xtem（2）*ytem（2））;
i=1;
for?xi=a（1）:0.01:a（2）
????for?yj=b（1）:0.01:b（2）
????????
????????D41=sqrt（abs（Dis42-xi）^2+（yj+y1）^2）;
????????D42=sqrt（abs（Dis42-xi）^2+（yj+y2）^2）;
????????D43=sqrt（abs（Dis42-xi）^2+（yj+y3）^2）;
????????
????????d41=sqrt（abs（Dis42-xi）^2+abs（yj-y1）^2）;
????????d42=sqrt（abs（Dis42-xi）^2+abs（yj-y2）^2）;
????????d43=sqrt（abs（Dis42-xi）^2+abs（yj-y3）^2）;
?
????????Z44=60*log（2*yj/r4）;
????????Z41=60*log（D41/d41）;
????????