我的M文件如下： E {6~oZ#L
function [G5]=fun(x); 9nu3+.&P
syms n1 n2 k1 k2 d; oJ#,XMKga
[subwave,nsub]=textread('subout.txt','%f %f'); |t$Ma'P
[filmwave,filmtran]=textread('filmtr.dat','%f %f'); Zmbfq8K
[p,q]=size(subwave); 4|A>b})H
for j=1:p; </uO e.l>Q
wavel(j)=subwave(j,1); t1E[uu ,V8
n(j)=n1+n2/((wavel(j))^2); .^J2.>.
k(j)=k1*exp(k2/wavel(j)); T843 ":
N(j)=n(j)-i*k(j); #plY\0 E@
phase(j)=2*pi*N(j)*d/wavel(j); 0"~i ^
ns(j)=nsub(j,1); L(GjZAP
B(j)=cos(phase(j))+ns(j)*i*sin(phase(j))/N(j); 9@Cv5L?p\
C(j)=i*N(j)*sin(phase(j))+ns(j)*cos(phase(j)); 8* Jw0mSw
T(j)=4*ns(j)/((B(j)+C(j))*conj(B(j)+C(j))); 1+?^0%AC
ts(j)=1-((1-ns(j))/(1+ns(j)))^2; fnx-s{c?
t1(j)=ts(j)*T(j); SkU9ON
tf(j)=filmtran(j,1)/100; U]a*uF~h
difference(j)=(t1(j)-tf(j))^2; UGK*G y
end R $zH]
G=symsum(difference(j),1,401); <w3!!+oK"
G1=subs(G,n1,sym('x(1)')); g^zs,4pPU<
G2=subs(G1,n2,sym('x(2)')); 0p89: I *0
G3=subs(G2,k1,sym('x(3)')); yCCrK@{oo
G4=subs(G3,k2,sym('x(4)')); j 65<8svl
G5=subs(G4,d,sym('x(5)')); Q_l'o3
J90v!p-
`({T]@]V
K\|FQ^#UYm
在matlab下运行 H#ClIh?'b
>> options=optimset('tolfun',1e-8,'tolx',1e-8); <$JaWL
>> [xmin,value,flag,output]=fminsearch(@fun,[1.7,4000,0.005,400,500],options); A_KW(;50
0Won9P
Y3Q9=u*5
结果出错如下： a MTu-hA
\6\<~UX^
??? Conversion to double from sym is not possible. B7imV@<
Error in ==> fminsearch at 175 X1~1&:V,<
fv(:,1) = funfcn(x,varargin{:}); o=`C<}
y+7PwBo%e
哪位高手给看看 Nz)l<S9>
