我的M文件如下： " f <+~
function [G5]=fun(x); ds UY[X-<6
syms n1 n2 k1 k2 d; OLM}en_L
[subwave,nsub]=textread('subout.txt','%f %f'); #mcU);s
[filmwave,filmtran]=textread('filmtr.dat','%f %f'); xK C{P{:
[p,q]=size(subwave); meNz0v e
for j=1:p; 9J'3b <
wavel(j)=subwave(j,1); tM;S
)S(=
n(j)=n1+n2/((wavel(j))^2); i7(\i2_P
k(j)=k1*exp(k2/wavel(j)); &24z`ZS[w6
N(j)=n(j)-i*k(j); qQ "O;_
phase(j)=2*pi*N(j)*d/wavel(j); jW!)5(B[A
ns(j)=nsub(j,1); `T;Y%"X!
B(j)=cos(phase(j))+ns(j)*i*sin(phase(j))/N(j); 8cyC\Rs
C(j)=i*N(j)*sin(phase(j))+ns(j)*cos(phase(j)); o|0QstS Cl
T(j)=4*ns(j)/((B(j)+C(j))*conj(B(j)+C(j))); K~JXP5`(
ts(j)=1-((1-ns(j))/(1+ns(j)))^2;
@s@ 67\
t1(j)=ts(j)*T(j); 4\%XC
F!
tf(j)=filmtran(j,1)/100; jIs2 R3B
difference(j)=(t1(j)-tf(j))^2; IB+) 2 `
end Ghpk0ia%d
G=symsum(difference(j),1,401); Co#_Cyxg=9
G1=subs(G,n1,sym('x(1)')); *X4$'LSx1
G2=subs(G1,n2,sym('x(2)')); z7P] g
C$\
G3=subs(G2,k1,sym('x(3)')); qEbzF#a-:
G4=subs(G3,k2,sym('x(4)')); >ey- j\_v
G5=subs(G4,d,sym('x(5)')); @|GKNW#
_C9 7G&
CAD:if V
qxe%RYdA'j
在matlab下运行 T{%' "mm;
>> options=optimset('tolfun',1e-8,'tolx',1e-8); /4Lmu+G4
>> [xmin,value,flag,output]=fminsearch(@fun,[1.7,4000,0.005,400,500],options); E[RLBO[*n
Ew kZzVuX
xz$S5tgDQK
结果出错如下： p-t*?p
C
z.7'yJIP#
??? Conversion to double from sym is not possible. `i)&nW)R
Error in ==> fminsearch at 175 .5~ W3 v
<
fv(:,1) = funfcn(x,varargin{:}); [o,S.!W8
?~5J!|r#
哪位高手给看看 g$f ;
