| fanyushan |
2010-07-19 14:12 |
fminsearch应用出错
我的M文件如下： >,w\lf9
function [G5]=fun(x); <3SFP3^:
syms n1 n2 k1 k2 d; VYl_ U?D
[subwave,nsub]=textread('subout.txt','%f %f'); >:Rt>po8|w
[filmwave,filmtran]=textread('filmtr.dat','%f %f'); I3}HNGvU
[p,q]=size(subwave); SeRK7Q&_
for j=1:p; <r_P?
lZW
wavel(j)=subwave(j,1); I 2t-D1X
n(j)=n1+n2/((wavel(j))^2); ! # ,-
k(j)=k1*exp(k2/wavel(j)); W,9k0t
N(j)=n(j)-i*k(j); .RQ Xxw
phase(j)=2*pi*N(j)*d/wavel(j); kbIY%\QSO
ns(j)=nsub(j,1); *([0"
B(j)=cos(phase(j))+ns(j)*i*sin(phase(j))/N(j); \j2 :
6]Hm
C(j)=i*N(j)*sin(phase(j))+ns(j)*cos(phase(j)); 6i9Q ,4~
T(j)=4*ns(j)/((B(j)+C(j))*conj(B(j)+C(j))); p" hm.=,
ts(j)=1-((1-ns(j))/(1+ns(j)))^2; `@f hge
t1(j)=ts(j)*T(j); bl:a&<F
tf(j)=filmtran(j,1)/100; |g7nh[
difference(j)=(t1(j)-tf(j))^2; /I`cS%U
end r)9i1rI+
G=symsum(difference(j),1,401); K ,lK\^y
G1=subs(G,n1,sym('x(1)'));
m+72C]9
G2=subs(G1,n2,sym('x(2)')); #K<=xP
G3=subs(G2,k1,sym('x(3)')); h8iaJqqvJ
G4=subs(G3,k2,sym('x(4)')); C;58z 5*,
G5=subs(G4,d,sym('x(5)')); 0,vj,ic*WX
I&'S2=s
F- !}dzO
zn2"swhq\V
在matlab下运行 75hFyh;u
>> options=optimset('tolfun',1e-8,'tolx',1e-8); nU? Xc(Xy
>> [xmin,value,flag,output]=fminsearch(@fun,[1.7,4000,0.005,400,500],options); WV,?G e
4t +/
8r`VbgI&
结果出错如下： |k7ts&2
?v ht~5'
??? Conversion to double from sym is not possible. ~L \( /[
Error in ==> fminsearch at 175 JhMrm%
fv(:,1) = funfcn(x,varargin{:}); ySr091Q
0btmao-
哪位高手给看看 oSu|Y n
|
|