计算方法实验.

1计算 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 实验指导姓名______________学号______________院系______________专业______________哈尔滨工业大学2计算方法实验指导根据实际问题建立的数学模型，一般不能求出所谓的解析解，必须针对数学模型的特点确定适当的计算方法，编制出计算机能够执行的计算程序，输入计算机，进行调试，完成运算，如果计算结果存在问题或不知是否正确，还需要重新确定新的计算方法，再编制出计算程序，输入计算机，重新调试，完成运算，直至获得正确的计算结果，这就是数值计算的全部过程。学生在学习“计算方法”和“高级语言”等课程时普遍存在的问题是：只会套用教科书中的 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 程序进行数值计算，很少有人能够独立地将学过的数值算法编制成计算机程序，至于灵活应用已经掌握的算法求解综合性较大的课题，则更是困难的事情。编写《计算方法实验指导》的目的是：突出数值计算程序结构化的思想。提高学生的编程能力，加深对“计算方法”课程内容的理解和掌握，为”计算方法“课程的教学服务，进一步奠定从事数值计算工作的基础。具体地1.根据“计算方法”课程内容的特点，给出五个典型算法的分析流程，学生可以利用所掌握的“高级语言”顺利地编制出计算机程序，上机实习，完成实验环节的教学要求。32.所有的计算实习题目都经过任课教师逐一检验，准确无误。3.充分利用循环的思想、迭代的思想，给出算法结构描述和程序语言的对应关系，有利于学生编制相应的程序。4.结合实习题目，提出实验要求，要求学生按 规范 编程规范下载gsp规范下载钢格栅规范下载警徽规范下载建设厅规范下载 格式写出相应的实验报告，实验报告成绩记入期末总成绩。需要提醒学生：不能简单地套用现成的标准程序完成实验题目，应当把重点放在对算法的理解、程序的优化设计、上机调试和计算结果分析上，否则就失去实验课的目的啦。5.五个具体的实验题目是：实验题目1拉格朗日(Lagrange)插值实验题目2龙贝格(Romberg)积分法实验题目3四阶龙格—库塔(Runge—Kutta)方法实验题目4牛顿(Newton)迭代法实验题目5高斯(Gauss)列主元消去法要求必须完成其中三个（如果全部完成更好）。4实验题目1拉格朗日(Lagrange)插值方法概要：给定平面上个不同的数据点，，，1n(,())kkxfx0,1,,knijxx；则满足条件ij，()()nkkPxfx0,1,,kn的次拉格朗日插值多项式n0()()()nnkkkPxfxlx是存在唯一的。若，且 函数 excel方差函数excelsd函数已知函数     2 f x m x mx m      2 1 4 2拉格朗日函数pdf函数公式下载 充分光滑，[,],0,1,,kxabkn()fx则当时，有误差估计式[,]xab，(1)01()()()()()()(1)!nnnffxPxxxxxxxn[,]ab拉格朗日插值算法实验实验目的：利用拉格朗日插值多项式求的近似值()nPx()fx输入：个数据点，；插值点1n(,())kkxfx0,1,,knx输出：在插值点的近似值()fxx()nPx程序流程：1置；0.0y0k2当时，做2.1—2.4kn52.1置；1.0l2.2对，置0,1,,1,1,,jkkn()/()jkjllxxxx2.3置()kyylfx2.4置1kk3输出,xy4停机问题1拉格朗日插值多项式的次数越大越好吗？n考虑下面两个拉格朗日插值问题：（1）设，，考虑等距节点的拉格朗日插值多21()1fxx[5,5]x项式，即将区间进行等分，记，，()nPx[5,5]n10.0hn5.0kxkh，构造，利用拉格朗日插值多项式作为的0,1,,kn()nPx()nPx()fx近似值。分别取，，，同时计算在，5n10n20n()nPx0.75x，，，处的函数值。1.75x2.75x3.75x4.75x（2）设，，考虑等距节点的拉格朗日插值多项()xfxe[1,1]x式，即将区间进行等分，记，，()nPx[1,1]n2.0hn1.0kxkh，构造，利用拉格朗日插值多项式作为的0,1,,kn()nPx()nPx()fx近似值。分别取，，，同时计算在，5n10n20n()nPx0.95x6，，处的函数值。0.05x0.05x0.95x问题2插值区间越小越好吗？考虑下面两个拉格朗日插值问题：（1）设，，考虑等距节点的拉格朗日插值多21()1fxx[1,1]x项式，即将区间进行等分，记，，()nPx[1,1]n2.0hn1.0kxkh，构造，利用拉格朗日插值多项式作为的0,1,,kn()nPx()nPx()fx近似值。分别取，，，同时计算在，5n10n20n()nPx0.95x，，处的函数值。0.05x0.05x0.95x（2）设，，考虑等距节点的拉格朗日插值多项()xfxe[5,5]x式，即将区间进行等分，记，，()nPx[5,5]n2.0hn1.0kxkh，构造，利用拉格朗日插值多项式作为的0,1,,kn()nPx()nPx()fx近似值。分别取，，，同时计算在，5n10n20n()nPx4.75x，，处的函数值。0.25x0.25x4.75x问题3在区间考虑拉格朗日插值问题，为了使得插值误[1,1]差较小，应如何选取插值节点？考虑下面两个拉格朗日插值问题：7（1）设，，考虑非等距节点的拉格朗日插值21()1fxx[1,1]x多项式，记，，构造，利用拉格()nPx(21)cos2(1)kkxn0,1,,kn()nPx朗日插值多项式作为的近似值。分别取，，()nPx()fx5n10n，同时计算在，，，处的函20n()nPx0.95x0.05x0.05x0.95x数值。（2）设，，考虑非等距节点的拉格朗日插值多()xfxe[1,1]x项式，记，，构造，利用拉格朗()nPx(21)cos2(1)kkxn0,1,,kn()nPx日插值多项式作为的近似值。分别取，，，()nPx()fx5n10n20n同时计算在，，，处的函数值。()nPx0.95x0.05x0.05x0.95x问题4考虑拉格朗日插值问题，内插比外推更可靠吗？考虑下面两个拉格朗日插值问题：（1）设，关于以，，为节点的拉格朗()fxx01x14x29x日插值多项式，利用拉格朗日插值多项式作为的近2()Px2()Px()fx似值。同时计算在，，，处的函数值。2()Px5x50x115x185x（2）设，关于以，，为节点的拉格()fxx036x149x264x朗日插值多项式，利用拉格朗日插值多项式作为的2()Px2()Px()fx8近似值。同时计算在，，，处的函数值。2()Px5x50x115x185x（3）设，关于以，，为节点的拉()fxx0100x1121x2144x格朗日插值多项式，利用拉格朗日插值多项式作为2()Px2()Px()fx的近似值。同时计算在，，，处的函数2()Px5x50x115x185x值。（4）设，关于以，，为节点的拉()fxx0169x1196x2225x格朗日插值多项式，利用拉格朗日插值多项式作为2()Px2()Px()fx的近似值。同时计算在，，，处的函数2()Px5x50x115x185x值。思考题：1.对实验1存在的问题，应如何解决？2.对实验2存在的问题的回答，试加以说明3.对实验3存在的问题的回答，试加以说明4.如何理解插值问题中的内插和外推？写出实验报告9实验题目2龙贝格(Romberg)积分法方法概要：利用复化梯形求积公式、复化辛普生求积公式、复化柯特斯求积公式的误差估计式计算积分。记，()bafxdxbahn，，其计算公式：kxakh0,1,,kn111[()()]2nnkkkThfxfx21111()222nnnkkTThfxh21(4)3nnnSTT21(16)15nnnCSS21(64)63nnnRCC10一般地，利用龙贝格算法计算积分，要输出所谓的数表T1214218421TTSTSCTSCR龙贝格(Romberg)积分法实验实验目的：利用龙贝格(Romberg)积分法计算积分()bafxdx输入：,,,abN输出：龙贝格数表T程序流程：1置，bahn1m11[()()]2Thfafb2输出1T3对，做3.1—3,52,3,,iN3.1置12iii置211111(())222iikTThfakh输出2T3.2置2211(4)3STT11输出2S3.3对，置1m2211(16)15CSS输出，转3.62C3.4对，置2m2211(64)63RCC输出，转3.62R3.5对，置3m21tolRR如果，则停机，否则转3.6tol3.6置，，，，，12RR12CC12SS12TT2hh1mm4停机问题1：利用龙贝格(Romberg)积分法计算积分（1），120exxdx610（2），31esinxxdx610（3），12041dxx610（4），1011dxx610问题2：被积函数无界，如何处理？（1），10sinxdxx610提示：0sin(0)lim1xxfx12（2），10cos1xdxx610提示：引进变换1tx（3），10cosxdxx610提示：利用等式，第一个积分值111000cos1cos1xxdxdxdxxxx等于2，第二个积分，利用；0cos1(0)lim0xxfx也可以考虑利用分部积分1100cos2cos()xdxxdxx（4），121sin1xxdxx610提示：利用第一类Gauss-Chebyshev求积公式问题3：积分区间无限，如何处理？（1），2e,xdx610提示：利用作近似21010exdx（2），11(1)dxxx610提示：利用变换1tx（3），23ecosxxdx610提示：Gauss-Hermite求积公式13（4），20esinxxdx610提示：Gauss-Lagurre求积公式思考题：1.输入的参数有什么意义？N2.在实验1中二分次数和精度的关系如何？3.在实验2中给出的提示具有普遍性吗？存在其它的方法吗？试加以说明。4.在实验3中给出的提示具有普遍性吗？存在其它的方法吗？试加以说明。写出实验报告14实验题目3四阶龙格—库塔(Runge—Kutta)方法方法概要：给定常微分方程初值问题d(,),d()yfxyaxbxbayahN记，，利用四阶龙格—库塔方法nxanh0,1,,nN112234311234(,)(,)22(,)22(,)1(22)6nnnnnnnnnnKhfxyKhKhfxyKhKhfxyKhfxhyKyyKKKK0,1,,1nN15可逐次求出微分方程初值问题的数值解，。ny1,2,,nN四阶龙格—库塔(Runge—Kutta)方法实验实验目的：利用四阶龙格—库塔(Runge—Kutta)方法求解微分方程初值问题输入：,,,abN输出：初值问题的数值解，。,nnxy0,1,2,,nN程序流程：1置00,,baxayhN2对，做2.1—2.41,2,,nN2.1置100(,)Khfxy1200(,)22KhKhfxy2300(,)22KhKhfxy4003(,)KhfxhyK2.2置10xxh1012341(22)6yyKKKK2.3输出11,xy162.4置0101,xxyy3停机问题1（1）d,01,5,10,20d(0)1yxyxNxbayhN准确解1yx（2）2d,01,5,10,20d(0)1yyxNxbayhN准确解11yx问题2（1）2d2e,13,5,10,20d(1)0xyyxxNxxbayhN17准确解2(ee)xyx（2）2d1(),13,5,10,20d(1)2yyyxNxxbayhN准确解212xyx问题3（1）2d20()2,01,5,10,20d1(0)3yyxxxNxbayhN准确解2201e3xyx（2）d2020sincos,01,5,10,20d(0)1yyxxxNxbayhN准确解1820esinxyx（3）d20(esin)e(sincos),01d(0)05,10,20xxyyxxxxxbayNhN准确解esinxyx思考题：1.对实验1，数值解和解析解相同吗？为什么？试加以说明。2.对实验2，越大越精确吗？试加以说明。N3.对实验3，较小会出现什么现象？试加以说明N写出实验报告19实验题目4牛顿(Newton)迭代法方法概要：求非线性方程的根，牛顿迭代法计算公式()0fx*x01()()nnnnxfxxxfx0,1,n一般地，牛顿迭代法具有局部收敛性，为保证迭代收敛，要求，对充分小的，。如果，，0*(,)Ox2()[,]fxCab*()0fx，那么，对充分小的，当时，由牛顿迭代法*()0fx0*(,)Ox计算出的收敛于，且收敛速度是2阶的；如果，{}nx*x()[,]mfxCab，，那么，对充分小**(1)*()()()0mfxfxfx()*()0(1)mfxm的，当时，由牛顿迭代法计算出的收敛于，且0*(,)Ox{}nx*x收敛速度是1阶的；20牛顿(Newton)迭代法实验实验目的：利用牛顿迭代法求的根()0fx输入：初值，精度，最大迭代次数12,N输出：方程根的近似值或计算失败标志()0fx*x程序流程：1置1n2当时，做2.1—2.4nN2.1置,0()Ffx0()DFfx如果，输出；停机1F0x如果，输出失败标志；停机2DF2.2置10/xxFDF2.3置10Tolxx如果，输出；停机1Tol1x2.4置，1nn01xx3输出失败标志4停机问题1：（1），，，，cos0xx6110421010N2100.7853981634x（2），，，，esin0xx6110421010N00.6x问题2：（1），，，，e0xx6110421010N00.5x（2），，，，222ee0xxxx6110421020N00.5x问题3：（1）由下面的递推公式可以生成勒让德（Legendre）多项式0121()1()231()()(),022nnnPxPxxnnPxxPxPxnnn①试确定和234(),(),()PxPxPx5()Px②确定，求得所有零点，精度6()Px6()Px6100.9324695142，0.6612093865，0.2386191861（2）由下面的递推公式可以生成切比雪夫勒让德（Chebyshev）多项式0121()1()()2()(),0nnnTxTxxTxxTxTxn①试确定和234(),(),()TxTxTx5()Tx22②确定，求得所有零点，精度6()Tx6()Tx610，21cos()2(1)jjxn0,1,,jn（3）由下面的递推公式可以生成拉盖尔（Laguerre）多项式01221()1()()(23)()(1)(),0nnnLxLxxLxnxLxnLxn①试确定和234(),(),()LxLxLx5()Lx②求得所有零点，精度5()Lx6100.2635603197，1.4134030591，3.5964257710，7.0858100059，12.6408008443（4）由下面的递推公式可以生成埃尔米特（Hermite）多项式0121()1()()2()2(1)(),0nnnHxHxxHxxHxnHxn①试确定和234(),(),()HxHxHx5()Hx②确定，求得所有零点，精度6()Hx6()Hx6102.3506049737，1.3358490740，0.4360774119思考题：231.对实验1确定初值的原则是什么？实际计算中应如何解决？2.对实验2如何解释在计算中出现的现象？试加以说明3.对实验3存在的问题的回答，试加以说明写出实验报告实验题目5相对高斯(Gauss)列主元消去法方法概要：高斯（Gauss）列主元消去法：对给定的阶线性方程组，nAxb首先进行列主元消元过程，然后进行回代过程，最后得到解或确定该线性方程组是奇异的。如果系数矩阵的元素按绝对值在数量级方面相差很大，那么，在进行列主元消元过程前，先把系数矩阵的元素进行行平衡：系数矩阵的每行元素和相应的右端向量元素同除以该行元素绝对值最大的元素。这就是所谓的平衡技术。然后再进行列主元消元过程。如果真正进行运算去确定相对主元，则称为显式相对Gauss24列主元消去法；如果不进行运算，也能确定相对主元，则称为隐式相对Gauss列主元消去法。显式相对Gauss列主元消去法：对给定的阶线性方程组n，首先进行列主元消元过程，在消元过程中利用显式平衡技Axb术，然后进行回代过程，最后得到解或确定该线性方程组是奇异的。隐式相对Gauss列主元消去法：对给定的阶线性方程组n，首先进行列主元消元过程，在消元过程中利用隐式平衡技Axb术，然后进行回代过程，最后得到解或确定该线性方程组是奇异的。实验目的：利用Gauss列主元消去法、显式相对Gauss列主元消去法、隐式相对Gauss列主元消去法求解线性方程组。Axb输入：；，，nijaib,1,2,,ijn输出：线性方程组的近似解，Axbix1,2,,in程序流程：一、Gauss列主元消去法251对，做1.1—1.3，消元过程1,2,,1kn1.1寻找最小的正整数，和。如果pkpnmaxpkjkkjnaa，输出奇异标志，停机；0pka1.2如果，那么交换两行；pk,pk1.3对，记，计算1,,ikn/ikikkkmaa1,,1,,1,,ijijkjikiikikaaamiknjknbbbmikn2.如果输出奇异标志，停机；0nna3.置，回代过程nnnnabx/4.对，置1,2,,1nkkknkjjkjkkaxabx/)(1二、显式相对Gauss列主元消去法1.对，做1.1—1.4，消元过程1,2,,1kn1.1对，计算，如果，输出奇异,1,,ikknmaxiijkjnsa0is标志，停机；计算，；/ijijiaas,1,,jkkn1.2寻找最小的正整数，和，如果pkpnmaxpkjkkjnaa，输出奇异标志，停机；0pka261.3如果，那么交换两行；pk,pk1.4对，记，计算1,,ikn/ikikkkmaa1,,1,,1,,ijijkjikiikikaaamiknjknbbbmikn2.如果输出奇异标志，停机；0nna3.置，回代过程nnnnabx/4.对，置1,2,,1nkkknkjjkjkkaxabx/)(1三、隐式相对Gauss列主元消去法1.对，做1.1—1.3消元过程1,2,,1kn1.1对，计算，如果，输出奇,1,,ikknmaxijkkjnsa0is异标志，停机；寻找最小的正整数，和pkpn；/max/pkpikikinasas1.2如果，那么交换两行；pk,pk1.3对，记，计算1,,ikn/ikikkkmaa271,,1,,1,,ijijkjikiikikaaamiknjknbbbmikn2.如果输出奇异标志，停机；0nna3.置，回代过程nnnnabx/4.对，置1,2,,1nkkknkjjkjkkaxabx/)(1问题1实验题目：（1）2499.19989.01262.12951.13927.02786.04002.01784.00643.03781.01920.03645.01129.04015.03872.02246.02943.03678.01234.04096.04321xxxx（2）43.22368.11008.12287.22617.177260.4600260.4601.132590.6700590.67810.98860.9000860.9001.1364321xxxx（3）420/31960/5760/7712/257/16/15/14/16/15/14/13/15/14/13/12/14/13/12/114321xxxx（4）3133233210957910685657787104321xxxx实验题目的准确结果：28（1）；1234(,,,)(1,1,1,1)TTxxxxx（2）；1234(,,,)(1,1,1,1)TTxxxxx（3）；1234(,,,)(1,1,1,1)TTxxxxx（4）。1234(,,,)(1,1,1,1)TTxxxxx问题2（1）123419730520680413646.871.347.452.011.788.676.410.880225.11.455.906.1336.56.60xxxx（2）12340.53980.71610.55540.29820.20580.52570.69240.35650.62550.05030.64650.81870.18720.12910.10700.58140.94000.77790.40420.1859xxxx（3）1231012131102131157xxx（4）12342422171025410915xxx思考题：1.计算实验1、实验2的各个题目说明：对什么类型的线性方程组三种方法是一致的？2.用三种方法计算实验1、实验2的各个题目，哪种方法最29好？试加以说明。3.综合上述两种结论，总结三种方法的关系，试加以说明。写出实验报告30说明1本课程给出五类实验题目，供学生选用，要求必须完成其中三个实验。2实验课的目的是为了让学生深入理解和掌握“计算方法”课程的基本内容，同时有助于培养学生的上机调试程序进行数值计算的动手能力，进一步提高利用数值方法求解数学问题、分析计算结果、选择算法的综合能力。3为了顺利完成实验教学的规定内容，建议学生按下面方法准备实验、进行实验、写出实验报告：（1）应明确实验的目的，清楚实验的内容，包括算法和误差分析；（2）写出内容摘要，包括算法的理论基础和对算法的初步认识；（3）上机前，编写好计算程序；（4）上机调试程序要做到快速、准确；（5） 记录 混凝土 养护记录下载土方回填监理旁站记录免费下载集备记录下载集备记录下载集备记录下载 计算结果要做到真实、准确；（6）课后，认真写好实验报告，包括对算法的新认识和体会，要特别注意对计算结果的分析和讨论，当然包括对计算结果的误差分析。31实验报告一题目（摘要）Lagrange插值给定平面上个不同的数据点，，，1n(,())kkxfx0,1,,knijxx；则满足条件ij，()()nkkPxfx0,1,,kn的次拉格朗日插值多项式n320()()()nnkkkPxfxlx是存在唯一的。若，且函数充分光滑，[,],0,1,,kxabkn()fx则当时，有误差估计式[,]xab，(1)01()()()()()()(1)!nnnffxPxxxxxxxn[,]ab前言：（目的和意义）目的：利用拉格朗日插值多项式求的近似值()nPx()fx意义：33数学原理34程序设计流程本实验采用CodeBlocks的C文件编写。Lagrange插值源程序：#include#include#includefloatlagrange(float*x,float*y,floatxx,intn){inti,j;float*a,yy=0.0;a=(float*)malloc(n*sizeof(float));for(i=0;i<=n-1;i++){a[i]=y[i];for(j=0;j<=n-1;j++)if(j!=i)35a[i]*=(xx-x[j])/(x[i]-x[j]);yy+=a[i];}free(a);returnyy;}main(){inti,n;floatx[20],y[20],xx,yy;printf("Inputn:");scanf("%d",&n);if(n>=20){printf("Error!Thevalueofnmustin(0,20).");getch();return1;}36if(n<=0){printf("Error!Thevalueofnmustin(0,20).");getch();return1;}for(i=0;i<=n-1;i++){printf("x[%d]:",i);scanf("%f",&x[i]);}printf("\n");for(i=0;i<=n-1;i++){printf("y[%d]:",i);scanf("%f",&y[i]);}printf("\n");37printf("Inputxx:");scanf("%f",&xx);yy=lagrange(x,y,xx,n);printf("x=%f,y=%f\n",xx,yy);getch();}实验结果、结论与讨论问题1（1）38N=5时，程序运行如下：N=10时，程序运行如下：39实验报告二题目（摘要）40前言：（目的和意义）数学原理4142程序设计流程43实验结果、结论与讨论44实验报告三题目（摘要）前言：（目的和意义）45数学原理46程序设计流程47实验结果、结论与讨论48实验报告四题目（摘要）0.75前言：（目的和意义）49数学原理50程序设计流程51实验结果、结论与讨论52