1、12014 矩阵实验报告矩阵与数值分析实验报告 任课教师： 张宏伟 教学班号： 02 院 系： 电信（计算机应用技术） 学生姓名： 学生学号： 1. 方程在 x=3.0 附近有根，试写出其三种不同的等价形;%迭代 3 while(epxE%标记 2 end; plot(y); title(迭代 1 误差下降曲线); 迭代公式 1 收敛结果:x=2 迭代 1 误差下降曲线32.92.82.72.62.52.42.32.22.12051015202530 迭代公式 1 误差变化曲线 2x 10138 迭代 2 误差下降曲线 3.532.521.510.50051015 迭代公式 2 误差变化曲线
2、迭代公式 3 收敛结果：x=6.答： （1）复化梯形公式代码： clear; syms x sum1 sum2 a1 b1; fun=(2/3)*(x)*exp(x ); sum1=0;%积分结果变量 1 ，并将计算结果与精确解进行比较： （2）sum2=0;%积分结果变量 2 n=1;%迭代计数变量 epx=1;%阶段误差 E=0.5e-8;%精度 a=1;%积分下限 b=2;%积分上限 while(epxE) h=(b-a)/n; for i=0:(n-1)%for 循环求解一次 N 点积分结果 a1=subs(fun,x,(a+i*h); b1=subs(fun,x,(a+(i+1)*h
3、); t=(a1+b1)*h)/2; sum2=sum2+t; end; epx=abs(sum1-sum2);%计算阶段误差 sum1=sum2;%使用 sum1 暂存上次的计算结果 sum2=0; n=n*2; end; disp(积分结果为：), vpa(sum1) 积分结果为： ans = 54.598157943101583100542484316975 真值约为：354.59815003314423 分析：由结果可知计算结果含有 8 位有效数字，已满足精度要求。使用数据点数为 N=8192。 复化辛普森公式代码： clear; 百度搜索“就爱阅读”,专业资料、生活学习,尽在就爱阅读网 ,您的在线图书馆!
