数值分析三-.docxVIP

数值分析作业三题目关于的下列方程组以及关于的下列二维数表00.40.81.21.620-0.5-0.340.140.942.063.50.2-0.42-0.5-0.260.31.182.380.4-0.18-0.5-0.5-0.180.461.420.60.22-0.34-0.58-0.5-0.50.620.80.78-0.02-0.5-0.66-0.5-0.021.01.50.46-0.26-0.66-0.74-0.5确定了一个二元函数。试用数值方法求出在区域上的一个近似表达式,要求最小的值达到以下精度：，其中计算的值，以观察逼近的效果，其中算法设计方案由已知条件：可以确定所有的与。将每一对代入六元非线性方程组中，用牛顿法迭代求解出与对应的一对值。将1得到的值应用分片二次插值，可以得到z的插值结果，这就是对应的函数值。利用曲面拟合的办法，求解针对不同的的中的系数 ，具体流程如下：对于给定的计算针对给定的基函数求对应的矩阵，以及矩阵，其中m+1分别代表给定的拟合点的个数，n+1代表给定的拟合点的个数。矩阵U可以表示为。系数矩阵可以表示为，其中B和G都是范德蒙行列式型的矩阵，记，。并有：，建立线性方程组：，求解该线性方程组得到的所有列向量，从而求解得到矩阵。，建立线性方程组：，为的列向量，求解该线性方程组得到D的所有列向量，从而求解得到矩阵。根据上述和可求解得到。计算在各点（）处的值，计算该情况下的精度:。若，则此时的满足要求，将，以及对应的输出。利用相同的方法计算，（），然后计算二者的差值。源代码!!!-******************************************************************!!!-数值分析计算实习作业3!!!-使用Fortran 95语言编写，采用f90自由格式编码!!!-****************************************************************** program main implicit noneinteger nnum,mnum!nnum,mnum用来表示要用于模拟的点x和y的个数parameter(nnum=10,mnum=20) double precision x(0:nnum),y(0:mnum)!用于拟合的x，y值double precision t(0:nnum,0:mnum),u(0:nnum,0:mnum) !对应于拟合x，y的t，u值double precision z(0:nnum,0:mnum) !对应于拟合x,y的z的值(分片二次插值得到)double precision p(0:nnum,0:mnum)!对应于拟合x,y的p的值integer i,j,k,n,mdouble precision sigma double precision,allocatable:: crs(:,:)double precision,allocatable:: x1(:),y1(:),z1(:,:)!用于x*，y*相关值的计算double precision,allocatable:: p1(:,:),t1(:,:),u1(:,:) open(10,file=data.txt,status=unknown) !对文件的操作do i=0,nnum,1 !xi,yi初始化 x(i)=0.08d0*iend dodo j=0,mnum,1y(j)=0.5d0+0.05d0*jend do !!!求解非线性方程组do i=0,nnum,1do j=0,mnum,1call solve_xy(x(i),y(j),t(i,j),u(i,j))end doend do !!!分片二次代数插值do i=0,nnum,1do j=0,mnum,1call chazhi(t(i,j),u(i,j),z(i,j))end doend do write(10,*) xi,yi以及f(xi，yj):do i=0,nnum,1do j=0,mnum,1 !将计算结果写入文档write(10,(x,i2,=,f7.4,y,i2,=,f7.4, f(x,i2,,y,i2,)=,e23.15)) i,x(i),j,y(j),i,j,z(i,j) end doend do !!!进行曲面拟合k=0do while(.True.) allocate(crs(0:k,0:k))call qumiannihe(k,k,nnum,mnum,x,y,z,crs,p,sigma) write(10,(对应k=,i3, 的精度值：,e23.15)) k,sigmaif(sigma.LE.1e-7) thenwrite(10,(当前k=,i3,