资源描述:
《[推荐]数值计算方法实验指导书》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、数值计算方法实验指导书陕西科技大学信息与计算科学教研室目录实验1代数插值实验(一)1实验2代数插值实验(二)4实验3数据拟合实验7实验4数值积分实验9实验5解非线性方程实验11实验6解线性方程组实验13实验7解非线性方程组实验15实验8解常微分方程组一导弹跟踪问题仿真实验17实验0Matlab基础实验专业班级学号姓名成绩.实验类型:•验证性实验O综合性实验O设计性实验实验目的:熟悉MATLAB系统的启动、退出、演示系统、帮助系统、MATLAB的运行环境,掌握Matlab矩阵及其运算基础知识实验内容:矩阵的建立、矩阵的拆分、特殊矩阵
2、、矩阵求逆,求行列式的值,矩阵算符及其表达式表示和计算、常用数学函数、眩图函数应用实验报告:根据实验情况和结果撰写并递交实验报告。实验1代数插值实验(一)实验类型:•验证性实验O综合性实验O设计性实验实验目的:进一•步熟练掌握Lagrange插值算法,提高编程能力和解决插值问题的实践技能。实验内容:依照Lagrange算法编写基于N+!个点的{(心,儿•)}:=()的N次Lagrange插值多项式系数计算的程序实验原理已知N+!个点的{(%儿)}驚NN计算Lagrange插值基函数厶=(口(兀-◎))/(!!(无-厂)(氐=0,・
3、・・,N)j=0y=0N计算N次Lagrange插值多项式函数P(x)=》)t厶仲(兀)"0实验说明需要建立两个数组X,Y分别存放(州,…,xj,(儿,…,儿)需要考虑多项式乘积运算的算法或者函数。实验步骤1要求上机实验前先编写出程序代码2编辑录入程序3调试程序并记录调试过程中岀现的问题及修改程序的过程4经反复调试后,运行程序并验证程序运行是否止确。5记录运行时的输入和输岀。实验总结实验报告:根据实验情况和结果撰写并递交实验报告。参考程序(MATLAB程序)function[C,L]=lagran(X,Y)%Input・Xisav
4、ectorthatcontainsalistofabscissas%・Yisavectorthatcontainsalistofordinates%Output-Cisamatrixthatcontainsthecoefficentsof%theLagrangeinterpolatorypolynomial%・LisamatrixthatcontainstheLagrange%coefficientpolynomials%NUMERICALMETHODS:MatlabPrograms%(c)2004byJohnH.Mathewsa
5、ndKurtisD.Fink%ComplementarySoftwaretoaccompanythetextbook:%NUMERICALMETHODS:UsingMatlab,FourthEdition%ISBN:043-065248-2%Prentice-HallPub.Inc.%OneLakeStreet%UpperSaddleRiver,NJ07458w=length(X);n=vv-l;L=zeros(vv,vv);%FormtheLagrangecoefficientpolynomialsfork=l:n+lV=l;f
6、orj=l:n+lifk〜=jV=conv(V,poly(X(j)))/(X(k)-X(j));endendL(k,:)=V;end%DeterminethecoefficientsoftheLagrangeinterpolator%polynomialC=Y*L;实验2代数插值实验(二)实验类型:•验证性实验O综合性实验O设计性实验实验目的:进一步熟练掌握spline(样条)插值算法,提高编程能力和解决插值问题的实践技能。实验内容:编写程序用计算机求解三次压紧样条曲线,经过点(0,0.0),(1,0.5),(2,2.0)和(3,
7、1.5),而且一阶导数边界条件S5(0)=0.2和S,(3)=1实验原理若已知N+!个点的{(无,儿)肚。及其一阶导数的边界条件S,(a)=d。和S,(b)二厶,则存在唯一的三次样条曲线。求解下列线性方程组绚-3(〃o-S©o)uk(k=23、…,N-2)u-3(S'(xQ-d3(㊁九+2力])“+/i1m2hs+2叽、4-hk)mk+hkmk+l3hN_2mN-2+(2饥—2Sk(X)=%0+%(X一心)+%2(兀一兀J+%3(兀一耳)'rd「5+%)实验说明需要建立两个数组X,Y分别存放(坷,…,£),(儿,…,儿),建立两个
8、变量dxo>dxn分别存放左右端点的一阶导数dxo=S*(x0),dxn=S^xn)实验步骤1要求上机实验前先编写出程序代码2编辑录入程序3调试程序并记录调试过程屮出现的问题及修改程序的过程4经反复调试后,运行程序并验证程序运行是否止确。5记录运行