我的M文件如下： #es9d3 ~\
function [G5]=fun(x); AW,53\ 0
syms n1 n2 k1 k2 d; "70WUx(\t
[subwave,nsub]=textread('subout.txt','%f %f'); mVR P~:+
[filmwave,filmtran]=textread('filmtr.dat','%f %f'); El@(mOu|
[p,q]=size(subwave); u{*SX k
for j=1:p; YJo ["Q
wavel(j)=subwave(j,1); phgm0D 7
n(j)=n1+n2/((wavel(j))^2); VP6ZiQ|
k(j)=k1*exp(k2/wavel(j)); ,%)6jYHR w
N(j)=n(j)-i*k(j); yfm^?G|sW
phase(j)=2*pi*N(j)*d/wavel(j); $5*WLG&AK
ns(j)=nsub(j,1); ~t~- A,1
B(j)=cos(phase(j))+ns(j)*i*sin(phase(j))/N(j); {>F7CT'G6
C(j)=i*N(j)*sin(phase(j))+ns(j)*cos(phase(j)); M{*Lp6h
T(j)=4*ns(j)/((B(j)+C(j))*conj(B(j)+C(j))); uh?>-
]r`
ts(j)=1-((1-ns(j))/(1+ns(j)))^2; U\R}`l
t1(j)=ts(j)*T(j); nG;8:f `
tf(j)=filmtran(j,1)/100; Q*b]_0 Rb
difference(j)=(t1(j)-tf(j))^2; R[;z X(y
end 1u6^z
G=symsum(difference(j),1,401); ;W^o@*i{>
G1=subs(G,n1,sym('x(1)')); Oj^,m .R
G2=subs(G1,n2,sym('x(2)')); D`VM6/iQR
G3=subs(G2,k1,sym('x(3)')); X0n~-m"m
G4=subs(G3,k2,sym('x(4)')); `3h SL R
G5=subs(G4,d,sym('x(5)')); uxzze~_+C
E~_]Lfs)
OdB?_.+ $
dx+hhg \L
在matlab下运行 UNkC L4N
>> options=optimset('tolfun',1e-8,'tolx',1e-8); 7=DjI ~
>> [xmin,value,flag,output]=fminsearch(@fun,[1.7,4000,0.005,400,500],options); ruazOmnn~
EMfdBY5
Yx>"bv
结果出错如下： oD$J0{K6
rhb@FE)Mc
??? Conversion to double from sym is not possible. 4'LB7}WG
Error in ==> fminsearch at 175 )-`;1ca)s
fv(:,1) = funfcn(x,varargin{:}); HdPo O;
N+ y&,N,
哪位高手给看看 m2v'WY5u
