《山东大学数值计算实验报告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的画图依旧感 觉很神奇。