资源描述：
《数值分析论文2》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、数值分析在水文地质中的应用摘要：本文通过运用数值分析屮线性方程组的直接解法，解决水文地质屮具体的问题，本文将地下水的流动的情况通过数学模型将英演示出来，再运用MATLAB求出地下水的各个参数。关键词:地下水；追赶法MATLABo数值分析是研究各种数学问题求解的数值计算方法，许多实际问题都需耍运用数值分析的各种算法來求解，同时联系计算机各种软件来实现解答。在水文地质中，地卜•水的流动很难描述，通过地卜•水的数值模拟将河流描述，运用数值分析的方法运用MATLAB实现。2实际问题描述考察通过x=0和x二L处的长月.直的河流为界的承压含水层，如下图，该含水层均质各向同
2、性，顶底板水平，上覆弱透水层，垂向补给强度为W(x),两河流边界的水位分别为「和W2,冃不随时间变化。首先，沿河流的方向取单宽作为计算区，并对计算区进行剖分，即江河间距L剖分成N等分，则空间步长为△x=L/No其次，在网格分割线上任取一点作为节点,节点编号由左向右依次为0」,……i,……No任一节点i的坐标为iAx,水位为H”已知节点0的水位为甲i，节点N的水位为甲2。L=800m,屮i=10m,x=8,则共有9个节点,是待求值。W2=5m,W=0.004m/d,T=100m2/d.若取△x=100m即N二L/A编号依次为0,1,……8,其中节点1,2,……7
3、的水头图2-5水文地质模型图。一水文地质平面图#一水文地质剖面图)k一-渗透系数;加一含水层厚度3数学模型的建立建立数学模型:护hT—+W(x)=O(OHi、Hj+i表zji,贝(I02Hg-2H,+H屮(Ax)2+0(2)这里将舍去余项O(Ar2),并以乩表示
4、节点i的水头Hi的近似值，贝用d2H_Hi-x-lHi+HM成立。带入数学模型得:WtAx2T(1=1,2，……N-1)定解条件为H°=0hn=022H-H2yy.Ax2T-H}+2比一比IV2Ax2T可列出方程：■H科_3+2HN_2—H冷上式出N-1个线性方程构成,~Hn-2+2Hnt=+02且未知量冇N・1个，故方程可写为:H、W2^c22-1-12-1-12-1将数值帯入得:_2-1_<510.41-12-10.4<>=<-12-1%0.4-12_—5.4kJ该式子的系数矩阵是正定的，故该题运用追赶法。追赶法的公式：一般形式为：C]a2b2C2aN-2
5、“N-2aN-5一2久-1A=—(i=1,2,……0Ci=bi一aQj(i=2,3,//=/i=ax(i=2,3,Hh2N-2)..…N-i)..…N—1)一2"n-1.利用下式可求出未知量Hi即:HN—=Yn-—H・r・伏H』=N_2,N_3,……2,1）4MATLAB编程及求解结果Clc%清空无用代码dispC用追赶法法求解〃［A］［N］二［门〃'）c二input（'输入系数矩阵＜/+1）:%输入对角线以上的矩阵b=input（'输入系数矩阵A（注）：’）；％输入对角线上的矩阵a二input('输入系数矩阵A(i+l,z)/);%输入对角线以下的矩
6、阵d二input(输入列矩阵b(i):');n=input('输入n);1(l)=b(l);%把b矩阵的第一个值赋给矩阵1y(l)=d(l)/l(l);%初始化矩阵yfori二2:nr(i-l)=c(i-l)/l(i-l);1(i)二b(i)-a(iT)*r(iT);y(i)=(d(i)-a(i-l)*y(i-l))/l(i);end%for循环结束disp('r1y');disp([r,1,y]');%显示r,l,y的值h(n)=y(n);forj=n-l：-1:1h(j)二y(j)-r(j)*h(j+l);end%for循环结束dispC解：');disp
7、(h);%输岀最后的结果节点节点坐标三对角方程组的系数递推系数节点水头IXiaibiCidiaiPiYiHi110002-110.42-1/25.210.7752200-12-10.43/2-2/311.2/311.153300-12-10.44/3-3/43.111.1254400-12-10.45/4-4/52.810.75500-12-10.46/5-5/68/39.8756600-12-10.47/6-6/755.2/21&657700-1205.48/707.0257.0255结论通过运用数值分析的解法，同时结合MATLAB软件，将地下水流的状态模拟
8、出来了，对地下水的研究有很重要的意义，
