我的M文件如下： 5D>cbzP@  
function [G5]=fun(x); Rk}=SB-  
syms n1 n2 k1 k2 d; Y{L|ja%9?  
[subwave,nsub]=textread('subout.txt','%f %f'); xI8v'[3  
[filmwave,filmtran]=textread('filmtr.dat','%f %f'); 2G$-:4B  
[p,q]=size(subwave); sNJ?Z"5k1h  
for j=1:p; Jo9c|\4  
    wavel(j)=subwave(j,1); d{DBG}/Yg  
    n(j)=n1+n2/((wavel(j))^2); \{1Vjo  
    k(j)=k1*exp(k2/wavel(j)); Ps~)l#gue  
    N(j)=n(j)-i*k(j); w^/"j_p@  
    phase(j)=2*pi*N(j)*d/wavel(j); M+lI,j+  
    ns(j)=nsub(j,1); +Q!Kj7EU/  
    B(j)=cos(phase(j))+ns(j)*i*sin(phase(j))/N(j); [H>/N7v19*  
    C(j)=i*N(j)*sin(phase(j))+ns(j)*cos(phase(j)); kLbo |p"cT  
    T(j)=4*ns(j)/((B(j)+C(j))*conj(B(j)+C(j))); <"`P;,S  
    ts(j)=1-((1-ns(j))/(1+ns(j)))^2; )==Qo/N