北航数值分析作业三

《北航数值分析作业三》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、《数值分析》第三次上机大作业姓名：学号：一、题目关于x,y,t,u,v,e的下列方程组厂0.5cost+u+v+w-x=2.67t+0.5sinu+v+w-y=1.070.5t+u+cosv+w-x=3.74jt+0.5u+v+sinw-y=0.79以及关于z,t,u的下列二维数表t00.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.10.620.80.78・0.02・0.5■0.66■0.5■0.021.

2、01.50.46・0.26・0.66・0.74•0.5确定了一个二元函数Z二f(x,y)o1.试用数值方法求出f(x,y)在区域D={(x,y)

3、0

4、牛顿迭代法求解非线性方程组，对兀=0.08厂儿=0.5+0.05;-(0

5、得到结果/(£』)，以观察pgy)逼近/(x,y)的效果。其中兀：=0.1厂y；=0・5+0.2厂(心1,2,…,8,)二1,2,..・,5)。三、源程序:#include#include#include#include#include#defineN4〃方程组未知个数#defineM6//z,t,u数表阶数#defineX_Num11^defineY_Num21〃定义数衣人小^defineEPSLle-12//定义阶数，粘度#defineEPSL2le-7usingnamespacestd;typed

6、efvectorEquation(Matinput);〃定义求解非线性方程的函数，同时供Inverse,Zxy函数调用Matlnverse(intrank,Matinput2);〃定义求解逆矩阵的函数doubleAccuracy(vectorX_l/vectorX_2);〃定义求解近似向吊滯度的函数doubleInterpolation(doubleu_l,doublet_l);〃定义分片代数二次插值函数MatCrs(vectorX/vector

7、ble>Y,MatU);〃最小二乘法求解近似衣达式系数MatZxy(vectorXl,vectorYl);//>£义非线性方程组,调用Equation,Accuracy和Interpolation完成求解〃所有的output应该调整，是否调整为输出到文件为好voidoutput(vectorFinall,vectorFinal2zMatFinal3);//定义输Hl函数，输出矩阵voidoutput2(MatXi);doublevector_u[M]={0,0.4,0.8,1.24.6,2};doublevector_t[M]=

8、{0,0.2,0.4,0.6,0.84}；doublemat_z[M][M]={{-0.5,-0.34,0.14,0.94,2.06,3.5},{-0.42,-0.5,-0.26,0.3,1.18,2.38},{-0.18,-0.5,-0.5,-0.18,0.46,1.42},{0.2乙034广0・58广0・5广0丄0.62},{0.78,-0.02,-0.5^0.66,-0.5^0.02},{1.5,0.4