| 查看: 437 | 回复: 2 | ||
[求助]
求助，程序怎么改？运行有错误。
|
|
%光振幅函数S1和S2程序% function result=Mie_S12(m,x,u) m=1.5; x=8; u=0.5; nmax=round(2+x+4*x^(1/3)); n=(1:nmax);nu=(n+0.5);z=m.^x;m2=m.*m; sqx=sqrt(0.5*pi./x);sqz=sqrt(0.5*pi./z); bx=besselj(nu,x).*sqx; bz=besselj(nu,z).*sqz; yx=bessely(nu,x).*sqx; hx=bx+i*yx; blx=[sin(x)/x,bx(1:nmax-1)]; blz=[sin(z)/z,bz(1:nmax-1)]; ylx=[-cos(x)/x,yx(1:nmax-1)]; hlx=blx+i*ylx; ax=x.*blx-n.*bx; az=z.*blz-n.*bz; ahx=x.*hlx-n.*hx; an=(m2.*bz.*ax-bx.*az)./(m2.*bz.*ahx-hx.*az); bn=(bz.*ax-bx.*az)./(bz.*ahx-hx.*az); %-1<=u<=1,u=cos(theta) p(1)=1; t(1)=u; p(2)=3*u; t(2)=3*cos(2*(acos(u))); for n1=3:nmax, p1=(2*n1-1)./(n1-1).*p(n1-1).*u; p2=n1./(n1-1).*p(n1-2); p(n1)=p1-p2; t1=n1*u.*p(n1); t2=(n1+1).*p(n1-1); t(n1)=t1-t2; end; ab=Mie_ab(m,x); an=ab(1,  ;bn=ab(2, ;pt=Mie_pt(u,nmax); pin=pt(1, ;tin=pt(2, ;n=(1:nmax); n2=(2*n+1)./(n.*(n+1)); pin=n2.*pin; tin=n2.*tin; S1=(an*pin'+bn*tin'); S2=(an*tin'+bn*pin'); result=[S1;S2]; |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
2楼2015-12-29 23:48:49
3楼2017-02-25 15:12:37
200