《山东大学数值计算实验报告7》由会员分享，可在线阅读，更多相关《山东大学数值计算实验报告7（6页珍藏版）》请在金锄头文库上搜索。
1、山东大学软件学院数值计算课程实验报告学号：姓名：班级：5班实验题目：1、 在二次曲面z=x2 + 2y2与圆柱面任1严+护=1交线上取10个不同的点1.1画散点图（三维图，plot3（）1.2用平面拟合这些点2、例题3.2 （图3.1）的21个点（选做）1.1画出散点图1.2用一次多项式拟合。给出拟合函数，确定矛盾方程组的系数矩阵、右端 项。正规方程的方法得到结果，并画图显示该拟合函数1.3用二次多项式拟合，画图1.4用双曲线拟合（选做）实验学时：2实验日期：2017年11月21日实验目的：最小二乘法数据拟合硬件环境：PC机软件环境：MATLAB实验步骤与内容：1在二次曲面和圆柱面上取以下十个
2、点：x02/52/51118/1318/138/58/52y0-4/54/5-11-12/1312/13-4/54/50z036/2536/2533612/169612/16996/2596/254使用plot函数画出散点图：代码：x=0,2/5,2/5,1,1,18/13,18/13,8/5,8/5,2;y=0,-4/5,4/5,-1,1,-12/13,12/13,-4/5,4/5,0;z=0,36/25,36/25,3,3,612/169,612/169,96/25,96/25,4;plot3(x,y,z,*)运彳丁结果：y=0,-4/5,4/5,-1,1,-12/13,12/13,-4/
3、5,4/5,0; z=0,36/25,36/25,3,3,612/169,612/169,96/25,96/25,4; x=x;y=y；z = z;A=ones(10,3);A(:,1)=x;A(：,2)=y; %得到右端项ans = Az %直接解方程得出结果运彳丁结果：ans =2. 0472-0. 00000. 5756从结果我们可以看出拟合的结果是：2.0472x+0.5756二z2.1利用plot函数画出散点图代码如下：t=0.0,0.5,1.0,1.5,2.0,2.5,3.0,3.5,4.0,4.5,5.0,5.5,6.0,6.5,7. 0,7.5,8.0,8.5,9.0,9.5,
4、10.0;y=2.9,2.7,4.8,5.3,7.1,7.6,7.7,7.6,9.4,9.0,9.6,10.0,10.2,9.7, 8.3,8.4,9.0,8.3,6.6,6.7,4.1;plot(t,y,*)运彳丁结果：Figure 1 X文件(E)编辑旧 查看(V)插入(D 工具CD 桌面(D)窗口酗 帮助(H)鼻10987654日除 气鼻卩相建 目 旦米米米米 米米米米米0123456789102.2用一次多项式拟合(恕我直言，这个和上个实验差不多)拟合函数为： f(t,x)=x1+x2t代码如下：t=0.0,0.5,1.0,1.5,2.0,2.5,3.0,3.5,4.0,4.5,5.0
5、,5.5,6.0,6.5,7. 0,7.5,8.0,8.5,9.0,9.5,10.0;y=2.9,2.7,4.8,5.3,7.1,7.6,7.7,7.6,9.4,9.0,9.6,10.0,10.2,9.7, 8.3,8.4,9.0,8.3,6.6,6.7,4.1;plot(t,y,*)t=t;y=y；A = ones(21,2); %个六行两列全是1的矩阵A(:,2) = t; %将A的第二列替换成tans = Ay%直接解方程 plot(t,y, *) %得到散点图hold on; %不覆盖原图，继续在这个图上画 syms m y1ezplot(y1=ans(1)+ans(2)*m,-1,10) hold on;运彳丁结果：解方程结果：5. 95110. 2860图像:线hold on; 运彳丁结果： 解方程结果:023斗568910图像:# Figure 1 X文件旧编辑(E)査看(V)插入(D 工具CD 桌面窗口團帮助止DD曰fe故纵的驭詔-復 S 口10-111111 米*111来来米8氣迸6/斗*a20/fy2 =(6013236363349989 m)/225179981368524819266/8855结论分析与体会：现在觉得拟合也挺简单的，感觉这个二次多项式拟合的超级好看！ Matlab的画图依旧感 觉很神奇。
