基于MATLAB的数值分析中的相关算法实现.ppt

基于MATLAB的《数值分析》有关算法实现 ;目的;;;;例：已知数据如下：;运行结果： 1 -1.0143 1.8863 0.91272 -0.73137 0.79518 0.80039 -0.51667 0.63198 0.74522 -0.40292 0.3151 各区间三次样条插值函数图形如下：;（2）求第二边界问题，在matlab命令窗口输入： [yy2 b2 c2 d2]=myspline3(x1,y1,xx1,2,0,0); disp([b2 c2(1:end-1,1) d2]); plot(x1,y1,bo,xx1,yy2,--k) 运行结果： 0.96966 0 -6.2652 0.92267 -0.93977 1.8813 0.79923 -0.43181 -0.46 0.74245 -0.51461 2.1442 各区间的三次样条插值函数图像为：;;例：用二分法求方程 在[0.5,0.6]间的解，要求精度为 ，输出各二分点的值和计算次数。;;例：用牛顿法求解 在 附近的根，根的准确值为1要求精度为 ;三、常微分方程数值解;例：用梯形法求解下列初值问题的数值解 精度??? ，n=10，列出各节点数值解，并与精确解比较 解：将方程函数写入函数mmyfunc.m中，然后再命令窗口中输入： y=ladder(0,1,1,10,1e-4)，可得到各节点数值解。 近似解和精确解的对比表格如下：;;例：用改进的欧拉法求解 精度为 ，n=10，列出各节点数值解，并与精确解比较。 解：将方程函数写入函数mmyfuno.m中，然后再命令窗口中输入： myenler(0,1,0,10,1e-4)，可得到各节点数值解。 近似解和精确解的对比表格如下：;四、雅克比迭代法;例：求解方程组 精度为 ，最大迭代次数为10，初始向量为 . 解：在matlab命令窗口中输入： A=[8 -3 2;4 11 -1;6 3 12]; b=[20;33;36]; X0=[0 0 0]; jkb(A,b,X0, 0.001,10) 运行结果为： 系数矩阵A不是严格对角占优的，此雅可比迭代不一定收敛 k = 8 雅可比迭代收敛,此方程组近似解X如下： ans = 3.0003 1.9999 0.9997;总结; 敬请各位评委老师批评指正 谢 谢!