资源描述：
《具有自由面格林函数计算程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、%具有自由面格林函数计算程序clearallifk=12;s(1,2)=4;if(ifk==12);%H(a,b)数据输出,一维法+迭代法0.1na=92;%nd=10;%节点个数hd=pi/(nd+1);%ha=0.1;%步长forj=1:na;%a=ha*j;s0,1)=a;s(j,2)=exp(-a)*(-0.5772-log(a));%s(j,3)=s(j,2);%ja=fix(a);%截尾取整函数，朝零方向取整m=fix(a)+15;%w=0;%wa=0;%wi=1;%ae=a/exp(1
2、.0);%fork=1:ja;%wi=wi*ae/k;%wa=wa+exp(k-a)/k*wi;%end;%wi=wi*exp(ja-a);fork=ja+1:m;wi=wi*a/k;wa=wa+wi/k;end;s(j,2)=s(j,2)-wa;s(j,4)=-wa;end;%j=1:y=-0.5*pi-hd/2;forj=1:nd;y=y+hd;gO)=y;end;fork=1:na;a=s(k,1);J=fix(a);forj=1:nd;w=g(j)；X=a*exp(w*i);d=1;ae=a
3、/exp(1.0);wj=O;forj1=1:J;d=d*ae/j1;wj=wj+exp(j1-X)*d/j1*(exp(j1*w*i)-1);end;d=d*exp(J-X);forj1=J+1:J+25;d=d*a/j1;wj=wj+d/j1*(exp(j1*w*i)-1);end;H=s(k,2)*exp(a-X)-i*g(j)*exp(-X)-wj;P(kj)=H;p1((k-1)*nd+j)=H;end;end;n4=10;h0=ha/n4;forj=1:nd;W二g(j);X=-exp(
4、w*i);y=-h0/2;X0)=X；fork=1:n4;y=y+hO;b0,k)=exp(y*X);end;end;n1=na;na=na+100;fork二n1+1:na;a=ha*k;forj=1:nd;X=x(j);H=exp(ha*X)*p(k-1,j);w二0;y=-h0/2;fork1=1:n4;y=y+h0;w=w+b(j,k1)/(a-y);end;H=H+w*hO;P(kj)=H;p1((k-1)*nd+j)=H;end;end;figure(2),plot(real(p1),i
5、mag(p1),'.*)；xlabelfreal*);ylabelfimag');ww=-1;if(ww>0);fid=fopen('oned.txf/w');fprintf(fid/%s',,%a=0.1i,i=1:150;angle=-pi/2+jh,h=pi/21j=1:20')fprintf(fid,*%s','%a,w,realH,imagH');foria=1:na;a二0.們a;forid=1:nd;w=-pi/2+hd*id;fprintf(fid,'%7.2f%7.4f%7
6、.4f%7.4f,,a,wJreal(p(ia,id)),imag(p(ia,id)));end;end;fclose(fid)end;%ww=-1;if(ww>0);end;%if(ifk==9);if(ifk==11);%数值积分法a=1.0;X=-0.26126+0.42632*i;na=20;nd=20;ha=0.1;hd=pi/(nd+1);foria=1:na;a=0.1*ia;forid=1:nd;w=-pi/2+hd*id;X=-exp(w*i);h=0.01;n=fix(a/h
7、);y=-0.5*h;w=0;fork=1:n;y二y+h;w=w+(exp(X*y)-exp(-X*y))/y;end;H=w*h;n1=5550;h1=0.01;y=a-0.5*h1;w=0;fork=1:n1;y=y+h1;w=w+exp(X*y)/y;end;H=H+w*h1;H=exp(a*X)*H;p(aidlrH-P£ia-Lrnd土dllH」end-end-figu「e(3)・P_ot(「ea_(p1)jmag(p1).:)_wwha=(wwvsfidnfopen(・imegra_・
8、tx
