| 查看: 330 | 回复: 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]; |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有31人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
回收溶剂求助
已经有6人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
2楼2015-12-29 23:48:49
3楼2017-02-25 15:12:37