山东大学离散数学期末试题答案.doc

山东大学离散数学期末试题答案.doc

ID:56718499

大小:345.50 KB

页数:9页

时间:2020-07-06

山东大学离散数学期末试题答案.doc_第1页
山东大学离散数学期末试题答案.doc_第2页
山东大学离散数学期末试题答案.doc_第3页
山东大学离散数学期末试题答案.doc_第4页
山东大学离散数学期末试题答案.doc_第5页
资源描述:

《山东大学离散数学期末试题答案.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、数学建模作业姓名:王士彬学院:计算机科学与技术班级:2014级计科2班学号:2014001300701.在区域x[-2,2],y[-2,3]内绘制函数z=exp^(-x2-y2)曲面图及等值线图。解:曲面图如下:>>x=-2:0.5:2;>>y=-2:0.5:3;>>[X,Y]=meshgrid(x,y);>>Z=exp(-X.^2-``Y.^2);>>mesh(X,Y,Z)>>等值线图如下:>>x=-2:0.5:2;>>y=-2:0.5:3;>>[X,Y]=meshgrid(x,y);>>Z=exp(-X.^2-Y

2、.^2);>>mesh(X,Y,Z)>>surf(X,Y,Z)>>surf(X,Y,Z)>>contour(X,Y,Z)>>2.已知一组观测数据,如表1所示.(1)试用差值方法绘制出x[-2,4.9]区间内的光滑曲线,并比较各种差值算法的优劣.(2)试用最小二乘多项式拟合的方法拟合表中的数据,选择一个能较好拟合数据点的多项式的阶次,给出相应多项式的系数和偏差平方和.(3)若表中数据满足正态分布函数.试用最小二乘非线性拟合的方法求出分布参数值,并利用锁求参数值绘制拟合曲线,观察拟合效果.解:(1)分别用最领近插值,分段

3、线性插值(缺省值),分段三次样条插值,保形分段三次插值方法绘制在x[-2,4.9]的光滑曲线,图形如下:样条插值效果最好,其次线性插值,最近点插值效果最差,在这里效果好像不太明显。最近点插值优点就是速度快,线性插值速度稍微慢一点,但效果好不少。所以线性插值是个不错的折中方法。样条插值,它的目的是试图让插值的曲线显得更平滑,为了这个目的,它们不得不利用到周围若干范围内的点,不过计算显然要比前两种大许多。MATLAB文件如下:>>x0=-2:0.3:4.9;>>y0=[0.102890.117410.131580.144

4、830.156560.166220.173320.177500.17853...0.176350.171090.163020.152550.14020.126550.112190.097680.08353...0.070150.058760.046870.037290.029140.02236];>>cx=-2:0.3:4.9;>>y1=interp1(cx,y0,cx,'nearest');>>y2=interp1(cx,y0,cx,'linear');>>y3=interp1(cx,y0,cx,'spline')

5、;>>y4=interp1(cx,y0,cx,'cubic');>>subplot(2,2,1),plot(cx,y0,'o',cx,y1,'-r'),title('NearestInterpolant');>>subplot(2,2,2),plot(cx,y0,'o',cx,y1,'-k'),title('LinearInterpolant');>>subplot(2,2,3),plot(cx,y0,'o',cx,y1,'-b'),title('SplineInterpolant');>>subplot(2,2,4

6、),plot(cx,y0,'o',cx,y1,'-k'),title('CubicInterpolant');>>subplot(2,2,1),plot(cx,y0,'o',cx,y1,'-r'),title('NearestInterpolant');(2),从图形可以看出曲线函数遵从幂函数的形式,设幂函数形式为:可化为即把非线性函数转化为线性函数,原线性函数形式为由此我们可以得出p(x)等价于lny;x等价于lnx;,我们可以先求出。求一个线性多项式使之在最小二乘准则下拟合这些观测值,问题即化为求使E()=利用多

7、元函数极值原理可知,若目标函数E()的极小值存在,一定有,用MATLAB工具我们可以求得最后的结果。>>log(x0);>>log(y0);>>x0=log(x0);>>y0=log(y0);>>n=length(x0);>>a=sum(x0);>>b=sum(y0);>>c=sum(x0.*y0);>>d=sum(x0.^2);>>a0=(d*b-c*a)*(n*d-a^2);>>a1=(n*c-a*b)/(n*d-a^2);>>a0,a1a0=-2.5891e+05-1.7515e+06ia1=0.1045-0.

8、3558i即系数a0为-2.5891e+05-1.7515e+06i,a1为0.1045-0.3558i其相应多项式的系数和偏差平方和.我们可以求出E=-7.2019e+13+2.1767e+13i其MATLAB文件如下:>>Y=a1*x0+a0;>>e=Y-y0;>>E=sum(e.^2)E=-7.2019e+13+2.1767e+13i即其

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。