數(shù)值逼近課程設(shè)計(jì)報(bào)告-資料下載頁
2025-03-23 03:20本頁面
【正文】 x對應(yīng)柱高的近似值為：m = 在求解該題的過程中，我分別運(yùn)用了C語言和matlab兩個(gè)軟件進(jìn)行編程實(shí)現(xiàn)，在C語言中，在時(shí)，運(yùn)用了牛頓插值多項(xiàng)式來進(jìn)行逼近原曲線，為了驗(yàn)證程序的準(zhǔn)確性，在這里取x=53來進(jìn)行計(jì)算，[ 525]之間，即證明了程序的準(zhǔn)確性。同理，可以取任意來驗(yàn)證，在這里取x=1880進(jìn)行驗(yàn)證。在matlab中，進(jìn)行了不同區(qū)間的多項(xiàng)式擬合，并且擬合的多項(xiàng)式表達(dá)式見上面的運(yùn)行結(jié)果。比較C語言和matlab程序可知，兩者的結(jié)果較為接近，但是在matlab中可以得到擬合的多項(xiàng)式的表達(dá)式，并且可以可到結(jié)果圖，讓人更直觀的看出結(jié)果。六、參考文獻(xiàn)[1] 譚浩強(qiáng). C語言程序設(shè)計(jì)[M]. 北京：清華大學(xué)出版社，2005.一、目的意義 (1)提高matlab編程能力 (2) matlab中插值函數(shù)以及其他畫圖函數(shù)的運(yùn)用 (3)提高分析程序的能力二、內(nèi)容要求 丘陵高程測定問題 在一丘陵地帶測量高程，x和y方向每隔100米測一個(gè)點(diǎn)，得高程數(shù)據(jù)如下。試擬合一曲面確定合適的模型，并由此找出最高點(diǎn)和該點(diǎn)的高程。x y 100 200 300 400 100 636 697 624 478 200 698 712 630 478 300 680 674 598 412 400 662 626 552 334 三、問題解決的方法與算法方法：用插值法計(jì)算出各個(gè)節(jié)點(diǎn)的值，用matlab中的mesh函數(shù)繪出地貌圖，即可以直觀的觀察。算法：四、計(jì)算程序程序1：function [s,x0,y0]=high1(N) x=[100 200 300 400]。 y=[100 200 300 400]39。 z=[636 697 624 478。698 712 630 478。680 674 598 412。662 626 552 334]。%mesh(x,y,z) %繪原始數(shù)據(jù)圖xx=linspace(100,400,N)。%即為從100到400分成N等分yy=linspace(100,400,N)39。 %加密縱坐標(biāo)數(shù)據(jù)到N個(gè)zh=interp2(x,y,z,xx,yy,39。cubic39。) %插值 三次多項(xiàng)式插值:39。cubic39。mesh(xx,yy,zh) %加密后的地貌圖s=0。 for i=1:N^2 if zh(i)s s=zh(i)。 %循環(huán)找最大值 n=mod(i,N)。%i對N取余 m=(in)/N。 end end x0=100+300/N*m %最大值點(diǎn)對應(yīng)的X值y0=100+300/N*n hold on。stem3(x0,y0,s,39。fill39。)程序2：clc。clear all。x=[100 100 100 100 200 200 200 200 300 300 300 300 400 400 400 400]。 y=[100 200 300 400 100 200 300 400 100 200 300 400 100 200 300 400]。z=[636 697 624 478 698 712 630 478 680 674 598 412 662 626 552 334]。 xi=100:5:400。yi=100:5:400。[X,Y]=meshgrid(xi,yi)。H=griddata(x,y,z,X,Y,39。cubic39。)%二元非等距插值surf(X,Y,H)。view(120,23)。hold on。maxh=vpa(max(max(H)),6)%控制變量計(jì)算結(jié)果的顯示位數(shù)[r,c]=find(H=single(maxh))% 檢索R中最大元素所在的位置（行標(biāo)r和列標(biāo)c）stem3(X(r,c),Y(r,c),maxh,39。fill39。)五、計(jì)算結(jié)果與分析下面為兩個(gè)不同matlab程序的運(yùn)行結(jié)果，程序1是用三次多項(xiàng)式插值編寫函數(shù)實(shí)現(xiàn)的，程序2是用二元非等距插值編寫M文件實(shí)現(xiàn)的。在程序1中，由于存在等分區(qū)間的問題，即給定不同的區(qū)間個(gè)數(shù)，會(huì)產(chǎn)生不同的結(jié)果，，即證明兩個(gè)程序的正確性。程序1運(yùn)行結(jié)果： high(100) x0 = 166y0 = 196ans = 程序1丘陵地貌圖以及最高值位置程序2運(yùn)行結(jié)果：maxh =r = 16c = 20六、參考文獻(xiàn)[1] 譚浩強(qiáng). C語言程序設(shè)計(jì)[M]. 北京：清華大學(xué)出版社，200
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號-1