试题七[线代方程组误差分析数值积分].docVIP

考试课程 数学实验 2004.6.24下午 班级 姓名 学号 得分 [说明] （1）第一、二、三题的答案直接填在试题纸上； （2）第四题将数学模型、简要解题过程和结果写在试题纸上；卷面空间不够时，请写在背面； （3）除非特别说明，所有计算结果小数点后保留4位数字。 （4）考试时间为120分钟。 一、（10分）已知一组数据 x 0 0.1 0.2 0.3 0.4 0.5 y 0 5.3 1.8 2.6 4.9 6.3 用数值积分求x在（0, 0.5）内y的积分为 1.7750 。根据表中数据，试建立2次多项式回归模型，其回归系数为。在显著性水平为0.05情况下有无异常数据（若有，请指出第几个点为异常数据）？第二个。 数值积分： 梯形公式：x=0:0.1:0.5; y=[0 5.3 1.8 2.6 4.9 6.3]; S1=trapz(x,y) 输出结果： S1 =1.7750 多项式回归： y=[0 5.3 1.8 2.6 4.9 6.3]; x1=0:0.1:0.5; %切削时间 x2=x1.^2; %刀具厚度 n=6; %已知的数据容量 T=[ones(n,1),x1, x2]; %1与自变量组成的输入矩阵 [b,bint,r,rint,s]=regress(y,T); %回归分析程序(α=0.05) b,bint,s, rcoplot(r,rint) 输出结果： b = 1.4821 5.5821 6.6071 二、(10分) 已知常微分方程组初值问题，试用数值方法求__0.897893___(保留小数点后6位数字)。 你用的方法是__龙格—库塔方法 _，其精度为_四阶__。 %待解常微分方程组函数M文件源程序： function dy=ff (x,y) dy=[y(2);-y(1)*sin(x)- exp(x)]; %应用欧拉方法和龙格-库塔方法求解该常微分方程： ts=:0.1:0.4; y0=[1,0]; [x,y]=ode45(@ff, ts,y0); %龙格-库塔方法求数值解 [, y(:,1)] 输出结果： (10分) 已知线性代数方程组Ax=b, 其中 , , , 若方程组右端项有小扰动，试根据误差估计式估计__0.06416__ (保留小数点后5位数字)；若取初值，则用高斯-赛德尔迭代法求解Ax=b时，_(1.07103,0.37900,-0.00072,-0.06084)__，此迭代是否收敛____是_______。 线性代数方程组解的误差分析： 故其误差上限为： A=[10 -7 0 1 ;-3 22 6 2 ;5 -1 51 -1 ;2 1 0 25]; b=[8 5 5 1]; db=[0 0 0 0.1]; d=cond(A)*norm(db)/norm(b) 输出结果： d = 0.061460749095949 A=[10 -7 0 1 ;-3 22 6 2 ;5 -1 51 -1 ;2 1 0 25]; D=diag(diag(A)); %从稀疏矩阵A中提取D L=-tril(A,-1); %从稀疏矩阵A中提取L U=-triu(A,1); %从稀疏矩阵A中提取U b=[8 5 5 1]; %设定方程组右端项向量 x= zeros(4,1); %设定方程组初始向量 m= inv(D-L)*U; n= inv(D-L)*b; %高斯-赛德尔