matlab编程软件工程硕士试题和答案.docVIP

解答 一已知观察数据 -2-101201320作图可见散点分布接近抛物线，试求它的拟合二次曲线。 首先绘制关于变化的散点图，利用matlab编写程序代码如下： x=[-2 -1 0 1 2]; y=[0 1 3 2 0]; plot(x,y) 得到图形如图一： 图一 以下是利用matlab拟合工具箱得到的结果： y = p1*x^2 + p2*x +p3 Coefficients: p1 = -0.64286 p2 = 0.1 p3 = 2.4857 Norm of residuals = 0.95618 故所求方程为： 用欧拉法解初值问题 y=5x1-y 0≤x≤1.0y0=0 y(0)=0;x=-2:1:2;h=0.1;for k=1:10 用欧拉方法与改进的欧拉方法求初值问题 在区间[0,1]上取的数值解。 ［解］　欧拉方法的计算公式为 现用matlab编程，程序如下 x0=0; y0=0; x(1)=0.1; y(1)=y0+0.1*5*x0*(1-y0); for n=1:9 x(n+1)=0.1*(n+1); y(n+1)=y(n)+0.1*5*x(n)*(1-y(n)); end; x Y 得到结果如下： x = Columns 1 through 9 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 Column 10 1.0000 y = Columns 1 through 9 0 0.0500 0.1450 0.2732 0.4186 0.5639 0.6948 0.8016 0.8810 Column 10 0.9345 X与y一一对应。这便是利用欧拉数值求解法求解一阶线性微分方程的结果。 三 写出计算线性代数方程组 的G-S迭代格式，并分析此迭代格式的收敛性。 解:在命令中调用函数e.m A=[3,1,3;1,3,3;-3,2,3]; B=[4,-1,2]’; [x,n]=e(A,B,[0,0,0]’,1.0e-6]) 四设的已知值为：，，，请用二次插值多项式计算的近似值，结果取四位有效数字。 取步长h=0.1，保留5位有效数字，并与准确解比较。 设二次拟合方程为： y = p1*x^2 + p2*x +p3 （*） 则x=1时，y=0.8415，x=1.5时，y=0.9975，x=2时，y=0.9093。此时用Maltab拟合工具箱求解拟合出来的函数关系式为插值关系式，将三组对应关系代入(*)式可得差值结果为： y = p1*x^2 + p2*x + p3 Coefficients: p1 = -0.4884 p2 = 1.533 p3 = -0.2031 Norm of residuals = 7.4476e-016 由此求得在x=1.8处的插值结果为0.9739。 如图二，为插值曲线 图二 取步长h=0.1，保留5位有效数字，代码如下： clear clc format long x=1:0.1:2; for i=1:11 y1(i)=-0.4884*x(i)^2+1.533*x(i)-0.2031; y2(i)=sin(x(i)); end 其中，y1为插值结果，y2为准确值,得到结果如下： y1 = Columns 1 through 4 0.841500000000000 0.892236000000000 0.933204000000000 0.964404000000000 Columns 5 through 8 0.985836000000000 0.997500000000000 0.999396000000000 0.991524000000000 Columns 9 through 11 0.973884000000000 0.946476000000000 0.909300000000000 y2 = Columns 1 through 4 0.841470984807897 0.891207360061435 0.932039085967226 0.963558185417193 Columns 5 through 8 0.985449729988460 0.99749498660