第四章 拟合与插值.ppt

* clc,clear x=[0,0.4,1.2,2,2.8,3.6,4.4,5.2,6,7.2,8,9.2,10.4,11.6,12.4,13.6,14.4,15]; y=[1,0.85,0.29,-0.27,-0.53,-0.4,-0.12,0.17,0.28,0.15,-0.03,-0.15,-0.07,0.059,0.08,0.032,-0.015,-0.02]; f=fittype(a*cos(k*x)*exp(w*x),independent,x,coeff,{a,k,w}); cfun=fit(x,y,f) xi=0:0.1:20; yi=cfun(xi); plot(x,y,r*,xi’,yi’,b-) * 运行结果： General model: cfun(x) = a*cos(k*x)*exp(w*x) Coefficients (with 95% confidence bounds): a = 0.9988 (0.9837, 1.014) k = 1.001 (0.9959, 1.006) w = -0.2067 (-0.2132, -0.2003) * * 练习题 1. 海水温度随着深度的变化而变化，海面温度较高，随着深度的增加温度越来越低，这样也就影响了海水的对流和混合，使得深层海水中的氧气越来越少，根据这一规律可对海水鱼层作一划分，现通过实验观测得一组海水深度h与温度t的数据如下： t 23.5 22.4 20.9 19.1 15.4 11.5 9.5 8.8 8.4 8.1 h 0 1.5 2.5 4.6 8.2 12.5 16.5 20.2 22.5 26.2 (1)找出深度与温度之间的函数t=t(h)。 (2)找出温度变化最快的深度位置。 * 2.在对山区农作物种植面积进行评估时，需对某一区域的地表面积进行估计，现测得一个矩形区域内一些地面点的高程(见附表)，假设所观测的矩形区域内均是陆地且整个地表可近似看做一张光滑曲面。 要求： （1）通过二维插值绘出该观察区域地表的曲面图 （2）估算出观察区域地表面积的近似值 （3）计算从矩形区域一边上任意一点到其对边上任意一点沿直线走向的地表曲线的近似长度（自己在矩形的对边上取两个定点求出）。 * x y 0 100 200 300 400 500 600 700 800 900 1000 1100 1200 0 381.5 397.3 401.5 396.1 379.2 361.5 339.6 320.7 307.8 299.6 293.5 285.5 273.0 100 398.0 418.5 423.2 416 401 380.5 356.3 336.4 325.3 315.5 310.4 306.3 298.5 200 411.2 428.3 436.5 433.5 420.5 398.5 374.5 353.4 338.3 330.2 325.4 323.1 317.0 300 420.2 439.6 448.8 443.1 435.5 415.2 392.3 370.2 354.1 346.5 342.3 338.2 332.5 400 423.8 444.6 456.7 459.2 447.6 429.2 407.3 385.5 369.4 360.3 356.5 352.5 345.2 500 422.1 441.5 459.5 463.2 454.6 439.5 419.5 400.5 384.5 374.5 369.5 364.5 356.4 600 416.2 439.2 454.5 466.5 457.6 446.5 429.5 412.4 398.5 386.3 381.2 375.5 362.1 700 406.3 429.5 442.3 454.6 450.2 448.1 435.2 422.2 409.3 399.5 391.2 380.2 370.5 800 391.5 414.6 431.4 443.9 447.8 440.2 436.5 426.5 416.5 406.3 397.4 382.2 363.5 900 373.5 395.6 413.5 426.6 433.5 434.5 430.5 424.5 416.5 407.6 390.0 388.6 370.0 1000 352.3 374.5 391.6 405.5 414.3 418.5 416.8 413.6 407.5 398.5 389.5 376.5 350.5 附表：区域[0,1200;0,1000]内观