电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

北航硕士数值分析大作业第三题

北航硕士数值分析大作业第三题_第1页
1/25
北航硕士数值分析大作业第三题_第2页
2/25
北航硕士数值分析大作业第三题_第3页
3/25
1 数值分析 A 计算实习题目 第三题 一、算法设计方案 题目分析如下: 若要进行达到一定精度的曲面拟合,则要得到在各点 ( ,)ijx y的函数值( ,)ijf x y,其中010,020ij; 如果得到各点的函数值,需要根据给出的, ,z t u 二维数表进行二次插值; 如果进行二次插值,则需要根据( ,)ijx y通过题中给出的非线性方程组解出插值点出对应的 ,t u 值,因此要解决整个问题,需要从解非线性方程组开始。 1.1 New ton 法解非线性方程组 New ton 法求方程组的解*x(对应未知数为( , , , )t u v w ),可采用如下的算法: 1) 在*x 附近选取初值(0)x,给定精度水平0 ; 2) 计算( )()kF x和( )'()kF x,其中有 ( )( )( )( )( )1( )( )( )( )( )2( )( )( )( )( )3( )( )( )( )( )4()0.5cos()2.67()0.5sin()1.07()0.5cos()3.74()0.5sin()0.79kkkkkikkkkkjkkkkkikkkkkjF xtuvwxF xtuvwyF xtuvwxF xtuvwy 2 ( )( )( )( )( )0.5sin()11110.5cos()11'()0.51sin()110.51cos()kkkkktuF xvw  3) 通过Doolittle 法求解如下的线性方程组 ( )( )( )'()()kkkF xxF x  4) 若( )( )/kkxx,则取*( )kxx,并停止计算,否则进入下一步; 5) 计算(1)( )( )kkkxxx  ,并回到第(2)步。 如此反复迭代,即可得到满足精度的非线性方程组的解,即得到各个 ( ,)ijt u的值。(其中010,020ij) 1 .2 分片二次插值 根据得到的各个( ,)ijt u值和给出的, ,z t u 二维数表,即可进行分片二次插值,得到插值点出的( , )f x y 值。具体算法如下: 对( , )x y ,若满足 ,2222,2222iiiihhxxxinyyyjm  则应选择(,)krxy(1, ,1;1, ,1)kii irjj j 为插值节点,相应的插值多项式为 112211( , )( ) ( ) (,)ijkrkrk irjpx yl x l y f xy   其中 1111( )(1, ,1)( )(1, ,1)itkt iktt kjtktjktt rxxlxkii ixxxxlyrjj jxx    3 如果12hxx或12nhxx,则在22p的表达式中取1i 或 1in ;如果12yy或12myy,则取1j  或 1jm 。 在此题中,t 对...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

北航硕士数值分析大作业第三题

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部