我的M文件如下： *5_8\7d  
function [G5]=fun(x); XvkI+c  
syms n1 n2 k1 k2 d; c]!D`FA*K  
[subwave,nsub]=textread('subout.txt','%f %f'); X)iWb(@k"7  
[filmwave,filmtran]=textread('filmtr.dat','%f %f'); .6 ?>t!&W  
[p,q]=size(subwave); a.&#dxgW[  
for j=1:p; VfA5r`^  
    wavel(j)=subwave(j,1); NuS|X   
    n(j)=n1+n2/((wavel(j))^2); r|\'9"@  
    k(j)=k1*exp(k2/wavel(j)); CNRU"I+jU  
    N(j)=n(j)-i*k(j); &;h~JS=  
    phase(j)=2*pi*N(j)*d/wavel(j); EfBVu  
    ns(j)=nsub(j,1); BUA6(  
    B(j)=cos(phase(j))+ns(j)*i*sin(phase(j))/N(j); @FdCbPl$  
    C(j)=i*N(j)*sin(phase(j))+ns(j)*cos(phase(j)); !ga(L3vf  
    T(j)=4*ns(j)/((B(j)+C(j))*conj(B(j)+C(j))); ['OCw {<  
    ts(j)=1-((1-ns(j))/(1+ns(j)))^2; @U /3iDB\  
    t1(j)=ts(j)*T(j); Bpm5dT;  
    tf(j)=filmtran(j,1)/100; F2PLy q  
    difference(j)=(t1(j)-tf(j))^2; 7=p-A_X  
end ADv"_bB:h  
G=symsum(difference(j),1,401); OBw`!G*w  
G1=subs(G,n1,sym('x(1)')); [ 4Y `O  
G2=subs(G1,n2,sym('x(2)')); 97]a-)SA  
G3=subs(G2,k1,sym('x(3)')); C/CfjRzd  
G4=subs(G3,k2,sym('x(4)')); >xZhK63C/  
G5=subs(G4,d,sym('x(5)')); m]=oaj@9  
(XG[_  
ueE?"Hk  
Z?~d']XD  
在matlab下運(yùn)行 Zu0;/_rN  
>> options=optimset('tolfun',1e-8,'tolx',1e-8); mj(&`HRs4  
>> [xmin,value,flag,output]=fminsearch(@fun,[1.7,4000,0.005,400,500],options); |:e|~sism  
^0s\/qyqm  
(+gL#/u  
結(jié)果出錯(cuò)如下： s*)41\V0  
Oa}V>a  
??? Conversion to double from sym is not possible. zOiY0`=  
Error in ==> fminsearch at 175 yK&*,J |  
fv(:,1) = funfcn(x,varargin{:}); H7yg9zFT N  
y3KcM#[  
哪位高手給看看 n(Q\',C  

哪位高人告訴我，錯(cuò)誤在哪里呀。M文件復(fù)制到命令行，執(zhí)行，一切正常，就是在求解多元函數(shù)0值時(shí)，不知道哪里錯(cuò)了。 oTF^<
I-C  
`\/toddUh[  
43J\8WBn@  
萬(wàn)分感謝呀