北航数值分析大作业3(2)

end do

k=1

write(30,\收敛过程为：')\

do while(1>0)

allocate(b(11,k+1),g(21,k+1),bt(k+1,11),gt(k+1,21),btb(k+1,k+1),gtg(k+1,k+1),c(k+1,k+1)) allocate(btuu(k+1,21),uu(11,21),btbni(k+1,k+1),gtgni(k+1,k+1),bniu(k+1,21),ggni(21,k+1)) do i=1,k+1 do j=1,11

b(j,i)=x(j)**(i-1) end do end do do i=1,k+1 do j=1,21

g(j,i)=y(j)**(i-1)

end do end do

bt=.t.b !调用IMSL库函数，进行曲面拟合的矩阵计算 gt=.t.g btb=bt.x.b gtg=gt.x.g btbni=.i.btb gtgni=.i.gtg do i=1,11 do j=1,21 uu(i,j)=f(i,j) end do end do

btuu=bt.x.uu

bniu=btbni.x.btuu ggni=g.x.gtgni c=bniu.x.ggni

jingdu=0.0 do i=1,11

do j=1,21

jingdu=jingdu+(f(i,j)-pxy(x(i),y(j),c))**2 !调用多项式表达式PXY函数，计算拟合精度 end do

end do

write(*,\write(30,\if(jingdu<=0.1e-6) exit

deallocate(b,g,bt,gt,btb,gtg,c)

deallocate(btuu,uu,btbni,gtgni,bniu,ggni) k=k+1

end do

write(30,\的系数C(R,S)为：')\do i=1,k+1

write(*,\write(30,\end do

do i=1,8 xn(i)=0.1*i end do do j=1,5

yn(j)=0.5+0.2*j end do do i=1,8

do j=1,5

call qiugeng(xn(i),yn(j),t(i,j),u(i,j)) !对X*,Y*进行新的非线性方程组计算，得到新的一组T,U end do end do do i=1,8

do j=1,5

f(i,j)=p22(u(i,j),t(i,j),u0,t0,h,tao,z) !插值计算新的F(X*,Y*) end do end do

write(30,\的值为：')\do i=1,8

do j=1,5

write(*,\

e20.12,

e20.12)\

xn(i),yn(j),pxy(xn(i),yn(j),c),f(i,j)

!

输

出

X*,Y*,F(X*,Y*),P(X*,Y*)，观察逼近程度

write(30,\ end do end do

stop

end program !主程序结束

三．输出结果：

x(i),y(j),f(x,y)的值为：

0.00 0.50 0.446504069242E+00 0.00 0.55 0.324683310597E+00 0.00 0.60 0.210159730631E+00 0.00 0.65 0.103043643055E+00 0.00 0.70 0.340192524001E-02

0.00 0.75 -0.887357886680E-01 0.00 0.00 0.00 0.00

0.80 -0.173371611924E+00 0.85 -0.250534594048E+00 0.90 -0.320276491493E+00 0.95 -0.382668052743E+00

0.00 1.00 -0.437795745638E+00 0.00 1.05 -0.485758921243E+00 0.00 1.10 -0.526667236102E+00 0.00 1.15 -0.560638463024E+00 0.00 1.20 -0.587796524654E+00 0.00 1.25 -0.608269768328E+00 0.00 1.30 -0.622189446194E+00 0.00 1.35 -0.629688384281E+00 0.00 1.40 -0.630899769389E+00 0.00 1.45 -0.625956164380E+00 0.00 1.50 -0.614988550047E+00 0.08 0.50 0.638015234587E+00 0.08 0.55 0.506611795876E+00 0.08 0.60 0.382176408075E+00 0.08 0.65 0.264863525561E+00 0.08 0.70 0.154780230633E+00 0.08 0.08 0.08 0.08

0.75 0.519927097057E-01 0.80 -0.434680179481E-01 0.85 -0.131601038093E+00 0.90 -0.212431072587E+00

0.08 0.95 -0.286004537561E+00 0.08 1.00 -0.352386059629E+00 0.08 1.05 -0.411655437858E+00 0.08 1.10 -0.463904893928E+00 0.08 1.15 -0.509236708669E+00 0.08 1.20 -0.547761104059E+00 0.08 1.25 -0.579594377202E+00 0.08 1.30 -0.604857251213E+00 0.08 1.35 -0.623673426166E+00 0.08 1.40 -0.636168257103E+00 0.08 1.45 -0.642467668455E+00 0.08 1.50 -0.642697116743E+00 0.16 0.50 0.840081396883E+00 0.16 0.55 0.699764166271E+00 0.16 0.60 0.566061444836E+00 0.16 0.65 0.439171614444E+00 0.16 0.70 0.319242167288E+00 0.16 0.75 0.206376218657E+00 0.16 0.80 0.100638546919E+00

0.16 0.85 0.206075984654E-02 0.16 0.16 0.16 0.16

0.90 -0.893540080136E-01 0.95 -0.173626954709E+00 1.00 -0.250799943998E+00 1.05 -0.320932252750E+00

0.16 1.10 -0.384097719188E+00 0.16 1.15 -0.440382160815E+00 0.16 1.20 -0.489881139368E+00 0.16 1.25 -0.532697954784E+00 0.16 1.30 -0.568941871351E+00 0.16 1.35 -0.598726545062E+00 0.16 1.40 -0.622168637157E+00 0.16 1.45 -0.639386546614E+00 0.16 1.50 -0.650499364497E+00 0.24 0.50 0.105151509248E+01 0.24 0.55 0.902927427685E+00 0.24 0.60 0.760580262813E+00 0.24 0.65 0.624715195878E+00 0.24 0.70 0.495519776775E+00 0.24 0.75 0.373134067685E+00 0.24 0.80 0.257656776698E+00 0.24 0.24 0.24 0.24

0.85 0.149150588766E+00 0.90 0.476470043700E-01 0.95 -0.468493081747E-01 1.00 -0.134356747670E+00

0.24 1.05 -0.214913334079E+00 0.24 1.10 -0.288573687009E+00 0.24 1.15 -0.355406352189E+00 0.24 1.20 -0.415491385199E+00 0.24 1.25 -0.468918240380E+00 0.24 1.30 -0.515783875727E+00 0.24 1.35 -0.556191070573E+00 0.24 1.40 -0.590246929357E+00 0.24 1.45 -0.618061557834E+00 0.24 1.50 -0.639746851961E+00 0.32 0.50 0.127124675843E+01 0.32 0.55 0.111500201761E+01 0.32 0.60 0.964607722308E+00 0.32 0.65 0.820347363242E+00 0.32 0.70 0.682447673529E+00 0.32 0.75 0.551085227171E+00 0.32 0.80 0.426392408546E+00 0.32 0.85 0.308463021403E+00 0.32 0.90 0.197357157387E+00

0.32 0.95 0.931056491500E-01 0.32 0.32 0.32 0.32

1.00 -0.428596469252E-02 1.05 -0.948339139383E-01 1.10 -0.178572979137E+00 1.15 -0.255553768875E+00

0.32 1.20 -0.325840141101E+00 0.32 1.25 -0.389506980643E+00 0.32 1.30 -0.446638198547E+00 0.32 1.35 -0.497324947822E+00 0.32 1.40 -0.541664031390E+00 0.32 1.45 -0.579756487909E+00 0.32 1.50 -0.611706299458E+00 0.40 0.50 0.149832107231E+01 0.40 0.55 0.133499864071E+01 0.40 0.60 0.117712512422E+01 0.40 0.65 0.102502405047E+01 0.40 0.70 0.878960016740E+00 0.40 0.75 0.739145103470E+00 0.40 0.80 0.605744887693E+00 0.40 0.85 0.478883880984E+00 0.40 0.90 0.358650648829E+00 0.40 0.40 0.40 0.40

0.95 0.245102261159E+00 1.00 0.138268376673E+00 1.05 0.381548905914E-01 1.10 -0.552527979027E-01

0.40 1.15 -0.141986870438E+00 0.40 1.20 -0.222094430862E+00 0.40 1.25 -0.295635227198E+00 0.40 1.30 -0.362679507525E+00 0.40 1.35 -0.423306161793E+00 0.40 1.40 -0.477601035622E+00 0.40 1.45 -0.525655428548E+00 0.40 1.50 -0.567564753485E+00 0.48 0.50 0.173189277935E+01 0.48 0.55 0.156203460147E+01 0.48 0.60 0.139721693052E+01 0.48 0.65 0.123780101007E+01 0.48 0.70 0.108408753012E+01 0.48 0.75 0.936322767079E+00 0.48 0.80 0.794704459103E+00 0.48 0.85 0.659387211130E+00 0.48 0.90 0.530487603091E+00 0.48 0.95 0.408088704433E+00 0.48 1.00 0.292244222111E+00

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库北航数值分析大作业3(2)在线全文阅读。