数值分析试题（卷）.doc

WORD格式可编辑 专业技术分享 1.设有某实验数据如下： 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 5.1234 5.3053 5.5684 5.9378 6.4270 7.0798 7.9493 9.0253 10.3627 （1）在MATLAB中作图观察离散点的结构，用最小二乘法拟合一个合适的多项式函数； （2）在MATLAB中作出拟合曲线图. 解：（1）在MATLAB中作图观察离散点的结构，用最小二乘法拟合一个合适的多项式函数。具体程序如下： x=[0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9]; y=[5.1234 5.3053 5.5684 5.9378 6.4270 7.0798 7.9493 9.0253 10.3627]; plot(x,y,r+),legend(x,y),xlabel(x),ylabel(y),title((x,y)的散点结构图) 从图中各点可以看到，（x,y）散点结构图的变化趋势与二次多项式很接近，故选取曲线函数为二次多项式函数，拟合多项式的次数为2。 用最小二乘法拟合一个合适的多项式函数，具体程序如下： s=polyfit(x,y,2); P=poly2str(s,t) P = 8.2191 t^2 - 1.8822 t + 5.3139 则曲线函数为发f（x）=8.2192t2-1.8822t+5.3139 在MATLAB中作出拟合曲线图，具体程序如下： x=linspace(0.1,0.9,9); y=8.2191.*x.^2-1.8822.*x+5.3139; plot(x,y),legend(x,y),xlabel(x),ylabel(y),title(拟合曲线y=f(x)) 2. 在MATLAB中用复合Simpson公式编程计算下列积分，并使其结果至少有6位有效数字. 解：（1）所安装的MATLAB中不含有复合辛普森程序，用M文件编辑器编写被积函数的M文件，并以jfSimpson.m命名保存至MATLAB搜索路径下。程序如下： function [I,step] = jfSimpson(f,a,b,type,eps) %type = 1 辛普森公式 %type = 2 复合辛普森公式 if(type==2 nargin==8) eps=1.0e-8; %缺省精度为0.0001 end I=0; switch type case 1, I=((b-a)/6)*(subs(sym(f),findsym(sym(f)),a)+... 4*subs(sym(f),findsym(sym(f)),(a+b)/2)+... subs(sym(f),findsym(sym(f)),b)); step=1; case 2, n=2; h=(b-a)/2; I1=0; I2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h; while abs(I2-I1)eps n=n+1; h=(b-a)/n; I1=I2; I2=0; for i=0:n-1 x=a+h*i; x1=x+h; I2=I2+(h/6)*(subs(sym(f),findsym(sym(f)),x)+... 4*subs(sym(f),findsym(sym(f)),(x+x1)/2)+... subs(sym(f),findsym(sym(f)),x1)); end end I=I2; step=n; end （2）在命令窗口中，调用jfSimpson函数进行求解。命令窗口显示如下： [I,s]=jfSimpson(inline(exp(x.^2)),0,2,2,1.0e-8) I = 16.4526 s = 102 说明：求解积分函数为16.4526。 3. 在MATLAB中试