《数值分析第四章外推法计算数值微分MATLAB计算实验报告(共3页).docx》由会员分享，可在线阅读，更多相关《数值分析第四章外推法计算数值微分MATLAB计算实验报告(共3页).docx（3页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上数值分析MATLAB计算实验报告姓名 班级 学号 一、实验名称用MATLAB编程实现数值微分的外推法计算。二、实验目的1掌握数值微分和定义和外推法的计算过程；2了解数值微分外推法的计算方法并且编写出与其算法对应的MATLAB程序代码；3体会利用MATLAB软件进行数值计算。三、实验内容 用外推法计算f(x)=x2e-x在x=0.5的导数。四、算法描述1.命名函数。2.如果输入未知数少于四个，默认精度10-33.描述T表矩阵坐标4.依次赋值计算 T表第一列5.根据数值微分计算公式求出T表矩阵的值6.若达到精度则运算结束，若未达到循环计算7.输出T表，得出的值就是导数值五
2、、实验结果六、实验结果分析此实验通过MATLAB实现外推法数值微分计算，得到相应的数据，方便对数据进行分析。从结果可以看出，当步长h=0.025时用中点微分公式只有3位有效数字，外推一次达到5位有效数字，外推两次达到9位有效数字。七、附录（程序）function g=waituifa(fname,x,h,e)if nargine i=i+1; G(i+1,1)=(feval(fname,x+h/2i)-feval(fname,x-h/2i)/(2*h/2i);for j=1:i G(i+1,j+1)=(4j)*G(i+1,j)-G(i,j)/(4j-1);endendGg=G(i+1,i+1);专心-专注-专业