| 查看: 400 | 回复: 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]; |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
2楼2015-12-29 23:48:49
3楼2017-02-25 15:12:37
