(整理)二维波动方程第一类吸收边界条件c++实现代码.

1、精品文档 #include stdafx.h #include #include #include #include using namespace std; const double pi=4*atan(1.0); double freq=45; double sb=7.45; double t1=2*pi/(sb*4); double source(double t) /double t2=0.0; 子波 if(t=t1) return (sin(sb*4*t-pi/2)+1)/10; else double tep=0.0; return tep; /return (1-2*pi*pi*

2、freq*freq*t*t)*exp(-pi*pi*freq*freq*t*t)+1);/Ricker void update_Vn(double upt,double lowt,double upx1,double lowx1) int i,j,m; const int Csize=300; double deg=0; double stepx1=abs(upx1-lowx1)/(Csize-1); /double te=sqrt(static_cast(3.0/8.0); double stept=sqrt(static_cast(1.0/2.0)*stepx1/2.0;/ int tn=

3、static_cast(upt/stept); double r=stept/stepx1; double *u_current,*u_old,*u_past; u_current=new double *Csize; u_old=new double*Csize; u_past=new double*Csize; for(i=0;iCsize;i+) u_currenti=new double Csize; u_oldi=new doubleCsize; u_pasti=new doubleCsize; for(i=0;iCsize;i+) for(j=0;jCsize;j+) u_curr

4、entij=0; u_oldij=0; 精品文档 u_pastij=0; double ckCsizeCsize; int flag=0; for(j=0;jCsize;j+) for(i=0;iCsize;i+) if(ji) ckij=4; else ckij=1; string str; coutstr; ofstream fout(str); if(!fout) coutn 不能打开文件 strendl;exit(1); m=0; double f0=2.0/(stept*30); double t0=4.0/f0; while(m1500iCsize;i+) for(j=0;jps曰

5、芳*jCsize;j+) fout u_currentji; /fout ckij; foutendl; foutendl; for(i=0;iCsize;i+) for(j=0;jCsize;j+) u_pastij=u_oldij; u_oldij=u_currentij; u_currentij=0; / m+; fout.close(); for(i=0;iCsize;i+) delete u_currenti,u_oldi,u_pasti; delete u_current,u_old,u_past; void main ()/主函数 double upt=17.8*t1*10; double lowt=0,upx1=3.0,lowx1=-3.0; / const int Csize=1