北航数值分析大作业第八题.doc

《北航数值分析大作业第八题.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、word格式北京航空航天大学数值分析大作业八学院名称自动化专业方向控制工程学号学生姓名许阳教师孙玉泉日期2014年11月26日....word格式一．题目关于x,y,t,u,v,w的方程组(A.3)(A.3)以及关于z,t,u的二维数表（见表A-1）确定了一个二元函数z=f(x,y)。表A-1二维数表tzu00.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.

2、78-0.02-0.5-0.66-0.5-0.021.01.50.46-0.26-0.66-0.74-0.51.试用数值方法求出f(x,y)在区域上的近似表达式要求p(x,y)以最小的k值达到以下的精度其中。2.计算(i=1,2,…,8;j=1,2,…,5)的值，以观察....word格式p(x,y)逼近f(x,y)的效果，其中。二．算法设计（一）总体思路1.题目要求对f(x,y)进行拟合，可选用乘积型最小二乘拟合。与的数表由方程组与表A-1得到。2.与1使用相同方法求得，由计算得出的p(x,y)直接带入求得。（二）算法实现1.与的数表的获得对区域上的f(x,y

3、)值可由方程组及二维数表得到。将区域D上的分别回代于方程组(A.3)，成为关于t,u,v,w的4元非线性方程组，解出每个对应的t,u。再通过表A-1进行插值近似，得到相应的z值。对应的z即为D区域上对应的。从而得到与的数表。(1)4元非线性方程组求解代入(A.3)后，原方程组变为关于t,u,v,w的4元非线性方程组。观察到方程组中方程形式较为简单，易于对变量t,u,v,w求偏导数，故而选用Newton法对方程组求解。计算方程组矩阵为：....word格式计算方程组偏导数矩阵为：迭代公式为：，k=0,1,2,…,n其中为线性方程组的解。取为迭代终止条件。由表A-1

4、观察到t,u基本在（0,2）上，于是选取为迭代初值。通过以上方法求得与对应的。(2)分片二元双二次代数插值为保证代数插值的收敛性，应采用分片低次插值。故此使用分片双二次代数插值。给定如满足如下关系式：....word格式，，则选择为插值节点，相应插值多项式为其中如果，则上式取m=1或m=4；如果或，则取n=1或n=4。得到表达式后，直接带入，得到的值即为与对应的。2.乘积型最小二乘曲面拟合使用乘积型最小二乘拟合，根据k值不用，有基函数矩阵如下：，数表矩阵如下：记C=[]，则系数的表达式矩阵为：通过求解如下线性方程，即可得到系数矩阵C。....word格式3.计算

5、(i=1,2,…,8;j=1,2,…,5)的值的计算与相同。将代入原方程组，求解响应进行分片双二次插值求得。的计算则可以直接将代入所求p(x,y)。....word格式三．Matlab源程序及结果牛顿法解非线性方程组子程序：function[t,u]=newt(x,y)t=1;u=1;v=1;w=1;ep=1e-12;k=1;N=100;while(k

6、5*u+v+sin(w)-y-0.79;dF=[-0.5*sin(t)111;10.5*cos(u)11;0.51-sin(v)1;10.51cos(w)];deltaX=ones(4,1);deltaX=dF^-1*(-1)*F;ifmax(abs(deltaX))/abs(x)

7、+1;end分片双二次代数插值子程序：functionf=p22(t,u)z=[-0.5-0.340.140.942.063.5;-0.42-0.5-0.260.31.182.38;-0.18-0.5-0.5-0.180.461.42;0.22-0.34-0.58-0.5-0.10.62;0.78-0.02-0.5-0.66-0.5-0.02;1.50.46-0.26-0.66-0.74-0.5];if(t<=0.3)....word格式i=1;endif(t>0.3&&t<=0.5)i=2;endif(t>0.5&&t<=0.7)i=3;endif(t>0.7

8、)i=4;endifu<