我的M文件如下： Hw[u Sv8
function [G5]=fun(x); nE)?P*$3Z
syms n1 n2 k1 k2 d; x7 jE
Ns )
[subwave,nsub]=textread('subout.txt','%f %f'); H
`Fe |6I&
[filmwave,filmtran]=textread('filmtr.dat','%f %f'); [_ESR/&N
[p,q]=size(subwave); 4)ez0[i$X
for j=1:p; %{s<h6{R
wavel(j)=subwave(j,1); ) 2C`;\/:
n(j)=n1+n2/((wavel(j))^2); %g DMz7$~
k(j)=k1*exp(k2/wavel(j)); ~=i9]%g ?
N(j)=n(j)-i*k(j); 5
rkIK
phase(j)=2*pi*N(j)*d/wavel(j); ?)k;.<6
ns(j)=nsub(j,1); LDHuf<`
B(j)=cos(phase(j))+ns(j)*i*sin(phase(j))/N(j); wN8-M e
C(j)=i*N(j)*sin(phase(j))+ns(j)*cos(phase(j)); 2ku \R 7
T(j)=4*ns(j)/((B(j)+C(j))*conj(B(j)+C(j))); -L(F:
ts(j)=1-((1-ns(j))/(1+ns(j)))^2; ZFh2v]|!
t1(j)=ts(j)*T(j); 8@ck" LUzD
tf(j)=filmtran(j,1)/100; !T02@e /
difference(j)=(t1(j)-tf(j))^2; Au08k}h<G
end yhg^1l|t,
G=symsum(difference(j),1,401); ,
}O>,AU
G1=subs(G,n1,sym('x(1)')); Y ## ftQ
G2=subs(G1,n2,sym('x(2)')); zl\mBSBx"
G3=subs(G2,k1,sym('x(3)')); v)>R)bzqe
G4=subs(G3,k2,sym('x(4)')); B$"CoLC7+
G5=subs(G4,d,sym('x(5)')); j-@3jFu
|13UJ
vR
~itrM3^"w
u {maE ,
在matlab下运行 ] Ec\!,54u
>> options=optimset('tolfun',1e-8,'tolx',1e-8); 6VpT*,2d~
>> [xmin,value,flag,output]=fminsearch(@fun,[1.7,4000,0.005,400,500],options); [f,; +Ze
8R}CvzI
, cD(s(6+
结果出错如下： E$d Pu
l -us j%\
??? Conversion to double from sym is not possible. =B_vQJF2
Error in ==> fminsearch at 175 `)$'1, ]u
fv(:,1) = funfcn(x,varargin{:}); p = _K P9
2bwf(
哪位高手给看看 zts%oIgV
