米散射程序 - 360文档中心

合集下载

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

matlab程序

clc;

clear

n=4+1i*(-2);

re=4;

im=-2;

l=100;

fid1=fopen('C:\Users\USER\Desktop\Î¢²¨\y1.txt','w');

fid2=fopen('C:\Users\USER\Desktop\Î¢²¨\y2.txt','w');

for x=0.1:0.15:15

m1=cos(x)-1i*sin(x);

m2=sin(x)+1i*cos(x);

w(1)=1/x*m2-m1;

A0=(sin(re*x)*cos(re*x)+1i*sinh(abs(im)*x)*cosh(abs(im)*x))/(power(s in(re*x),2)+power(sinh(abs(im)*x),2));

A(1)=-1/(n*x)+1/(1/(n*x)-A0);

a(1)=((A(1)/n+1/x)*real(w(1))-real(m2))/((A(1)/n+1/x)*w(1)-m2); b(1)=((n*A(1)+1/x)*real(w(1))-real(m2))/((n*A(1)+1/x)*w(1)-m2); y1(1)=3*((abs(a(1)))^2+(abs(b(1)))^2);

y2(1)=3*real(a(1)+b(1));

for j=2:l

if (j==2)

w(2)=(3/x)*w(1)-m2;

else

w(j)=((2*j-1)/x)*w(j-1)-w(j-2);

end

A(j)=-j/(n*x)+1/(j/(n*x)-A(j-1));

a(j)=((A(j)/n+j/x)*real(w(j))-real(w(j-1)))/((A(j)/n+j/x)*w(j)-w(j-1 ));

b(j)=((n*A(j)+j/x)*real(w(j))-real(w(j-1)))/((n*A(j)+j/x)*w(j)-w(j-1 ));

y1(j)=y1(j-1)+(2*j+1)*(abs(a(j))^2+abs(b(j))^2);

y2(j)=y2(j-1)+(2*j+1)*real(a(j)+b(j));

end

if (y1(j)-y1(j-1)<10^(-12) && y2(j)-y2(j-1)<10^(-12))

fprintf(fid1,'%f\n',2/(x^2)*y1(j));

fprintf(fid2,'%f\n',2/(x^2)*y2(j));

end

fclose(fid1);

fclose(fid2);