| 查看: 376 | 回复: 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]; |
回复此楼» 猜你喜欢
青C资助名额大幅增加！
已经有12人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有6人回复
-
有带发论文的吗
已经有3人回复
-
导师各种操作恶心咋办
已经有13人回复
-
本人最近太闲了，谁有问题可以提，每天会统一回复
已经有12人回复
-
护理论文 晋升
已经有5人回复
-
求助大佬sci投稿哪个好中
已经有4人回复
-
评审感受-评审感受-评审感受
已经有14人回复
-
这年头没有找到涵评专家，还有中面上的可能吗
已经有6人回复
-
2026博士申请求助
已经有10人回复
2楼2015-12-29 23:48:49
3楼2017-02-25 15:12:37