第二部分与_常用算法 .pptVIP

第二部分常用算法与MATLAB语言编程 1．算法的概念 2. 方程求根 3．矩阵求逆 4．数字积分（欧拉，梯形，龙格库塔） 5．数字积分的精度及稳定性 常用算法 方程求根 矩阵求逆 数值积分 欧拉法 重复上述近似过程，当t=t2时， 则有一般近似公式 如果令 ，称为计算步矩，则 (1) 这就是欧拉法数字积分的递推计算公式。 如果hn是固定不便的，称定步长积分； 如果hn是变化的，称变步长积分。 由欧拉法数字积分的递推计算公式可看出，只要给出方程的初值（t0, y0）以及相应的步距，逐步进行递推就可获得微分方程的近似数字解。 欧拉法的计算是十分简单的，其计算的截断误差正比于h2 ， 由此，要获得高精度解，必须减小步距，但这使得计算次数增加，又由于计算机的字长有限，h减小得过小，将引起由于舍入误差累计而产生的累积误差过大，所以此方法的精度提高有限，实际应用中较少采用。 %%%%欧拉法积分程序 x0=0;y0=0;h=0.0001; xm=5; xx=[]; x=x0;y=y0; while x=xm x=x+h; y=y+fx(x)*h; % x=x+h; x