数值分析的常用方法.ppt

曲线拟合 Matlab实现 2013古塔变形问题：该塔13层，对古塔做了4次观测，每层取了8个观测点。 第一问，给出通用方法，计算各层各次的中心位置 第二问，考虑该塔倾斜，弯曲，扭曲情况 通用计算均匀薄板的中心（二重积分） 转化为n个质量不等质点求中心。 表1 前2次观测各层中心坐标 1986年的一组数据 U =[0 0.064 0.136 0.184 0.246 0.304 0.342 0.379 0.416 0.45 0.516 0.583 0.639]; Z=[1.79 7.32 12.76 17.08 21.72 26.24 29.34 33.35 36.35 40.17 44.44 48.71 52.84]; 二次多项式拟合 z=ployfit(u,z,2) 输出结果： -7.62 86.04 1.53 三次多项式拟合 z=ployfit(u,z,3) -62.97 53.36 70.87 2.16 三次样条插值 U0=0:0.639/12:13); 1.7900 6.3911 10.5311 14.8975 19.2503 23.2981 27.5045 32.6726 37.4735 42.0624 45.4916 48.9128 52.8400 样条函数在13个等分点上的值 例: 航行区域的警示线 某海域上频繁地有各种吨位的船只经过。 为保证船只的航行安全，有关机构在低潮时对水深进行了测量,表3-8是他们提供的测量数据： 表3-8. 水道水深的测量数据 x 129.0 140.0 103.5 88.0 185.5 195.0 105.5 y 7.5 141.5 23.0 147.0 22.5 137.5 85.5 z 4 8 6 8 6 8 8 x 157.5 107.5 77.0 81.0 162.0 162.0 117.5 y -6.5 -81.0 3.0 56.5 -66.5 84.0 -33.5 z 9 9 8 8 9 4 9 航行区域的警示线 其中（x, y）为测量点，z为（x, y）处的水深(英尺)，水深z是区域坐标（x, y）的函数z= z (x, y）， 船的吨位可以用其吃水深度来反映，分为 4英尺、4.5英尺、5英尺和 5.5英尺 4 档。 航运部门要在矩形海域（75，200）×（－50，150）上为不同吨位的航船设置警示标记。 请根据测量的数据描述该海域的地貌，并绘制不同吨位的警示线，供航运部门使用。 Matlab求解 x=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5]; y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5]; z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9]; cx=75:0.5:200; cy=-70:0.5:150; cz=griddata(x,y,z,cx,cy','cubic'); meshz(cx,cy,cz), xlabel('X'),ylabel('Y'),zlabel('Z') figure(2),contour(cx,cy,cz,[-5 -5]); grid on, hold on plot(x,y,'+') xlabel('X'),ylabel('Y') 非数值信息的转化 根据数据类型寻找方法 2013年全国赛A题视频信息转化为数值信息，人工记数或专业统计软件 2008年A题数码相机 图片信息转化为数值信息（matlab可以直接做） 这类问题较少遇到。 随机变量分布与假设检验 诸如学生成绩分布，人员到达，车辆到达，车流量，商品销售量，医院门诊人数等都是随机变量，稳定的环境下其满足某种分布，具体是什么分布，往往是先猜测后验正。 没有检验的猜测很有可能是