fanyushan |
2010-07-19 14:12 |
fminsearch应用出错
我的M文件如下: T :S{3 function [G5]=fun(x); `;qv} syms n1 n2 k1 k2 d; ) .-(-6=R [subwave,nsub]=textread('subout.txt','%f %f'); O mh&)|Iql [filmwave,filmtran]=textread('filmtr.dat','%f %f'); D
,ZNh1xt [p,q]=size(subwave); 3zA=q[C for j=1:p; gqdB!l4 wavel(j)=subwave(j,1); u4#~
i0@ n(j)=n1+n2/((wavel(j))^2); `?o1cf A
k(j)=k1*exp(k2/wavel(j)); _ K+V?-= N(j)=n(j)-i*k(j); &]yJCzo] phase(j)=2*pi*N(j)*d/wavel(j); ~R?dDL ns(j)=nsub(j,1); E"1;i B(j)=cos(phase(j))+ns(j)*i*sin(phase(j))/N(j); T3t~=b>&L C(j)=i*N(j)*sin(phase(j))+ns(j)*cos(phase(j)); 1yaIV+_y/ T(j)=4*ns(j)/((B(j)+C(j))*conj(B(j)+C(j))); 9 QCpXy ts(j)=1-((1-ns(j))/(1+ns(j)))^2; -KwL9J4u t1(j)=ts(j)*T(j); d+m6-4[_k tf(j)=filmtran(j,1)/100; J ]Gc difference(j)=(t1(j)-tf(j))^2; So\| Ye end aC8,Y$>?E` G=symsum(difference(j),1,401); k`7.p,;}U G1=subs(G,n1,sym('x(1)')); :YJ7J4 G2=subs(G1,n2,sym('x(2)')); af?\kBm G3=subs(G2,k1,sym('x(3)')); `/wq3+ ? G4=subs(G3,k2,sym('x(4)')); >uchF8)e| G5=subs(G4,d,sym('x(5)')); 4'{hI;&a& 2.Eu+*UC itC *Z6^ ;x3 ]4^ 在matlab下运行 \;B$hT7z* >> options=optimset('tolfun',1e-8,'tolx',1e-8); 3dxnh,]&@ >> [xmin,value,flag,output]=fminsearch(@fun,[1.7,4000,0.005,400,500],options); 1^60I#Vr@ Z}K.^\S9 o9T@uWh+ 结果出错如下: H7tQ# y800(z ??? Conversion to double from sym is not possible. }* B qi7E> Error in ==> fminsearch at 175 1IF'>* fv(:,1) = funfcn(x,varargin{:}); dp
UdFuU" @!&Jgg53G 哪位高手给看看 }[KDE{,V
|
|