河南城建學(xué)院matlab上機(jī)實(shí)驗(yàn)答案(編輯修改稿)
2024-07-22 00:53 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 = +003 y4=interp1(x,y,)y4 = 1696 y5=spline(x,y,)y5 = 1682習(xí)題（1）用函數(shù)roots求方程x2x1=0的根 roots([1 1 1])ans = （2）y=sinx，0≤x≤2π，在n個(gè)節(jié)點(diǎn)（n不要太大，如取5~11）上用分段線性和三次樣條插值方法，計(jì)算m個(gè)插值點(diǎn)（m可取50~100）的函數(shù)值。通過數(shù)值和圖形輸出，將兩種插值結(jié)果與精度進(jìn)行比較。適當(dāng)增加n，再作比較。 x=linspace(0,2*pi,8),y=sin(x)x = 0 y = 0 xi=linspace(0,2*pi,100)。y0=sin(xi)。y1=interp1(x,y,xi)。y2=interp1(x,y,xi,39。spline39。)。 plot(xi,y0,39。*39。,xi,y1,39。.39。,xi,y2) e1=y1y0。e2=y2y0。 plot(xi,e1) plot(xi,e2)（3）大氣壓強(qiáng)p隨高度x變化的理論公式為p=+5007756，為驗(yàn)證這一公式，測(cè)得某地大氣壓強(qiáng)隨高度變化的一組數(shù)據(jù)如表所示。試用插值法和擬合法進(jìn)行計(jì)算并繪圖，看那種方法較為合理，且總誤差最小。高度/m0300600100015002000壓強(qiáng)/Pa插值法： x=[0 300 600 1000 1500 2000]。 p=[ ]。 xi=linspace(0,2000)。p0=*exp((xi+500)/7756)。 p1=interp1(x,p,xi,39。spline39。)。 plot(xi,p0,39。*39。,xi,p1) e1=p1p0。 e=sum(e1.^2)e = 擬合法： x=[0 300 600 1000 1500 2000]。 p=[ ]。 P=log10(p)P = p1=polyfit(x,P,1)p1 = b=p1(1)/,a=10.^p1(2)b = a = xi=linspace(0,2000)。p0=*exp((xi+500)/7756)。 p2=polyval(p1,xi)。P2=10.^p2。 e2=P2p0。e=sum(e2.^2)e = 四 Matlab數(shù)值運(yùn)算數(shù)值微積分練習(xí)：瑞士地圖如圖所示，為了算出其國(guó)土面積，首先對(duì)地圖作如下測(cè)量：以由西向東方向?yàn)閄軸，由南到北方向?yàn)閅軸，選擇方便的原點(diǎn)，并將從最西邊界點(diǎn)到最東邊界點(diǎn)在X軸上的區(qū)間適當(dāng)劃分為若干段，在每個(gè)分點(diǎn)的Y方向測(cè)出南邊界點(diǎn)和北邊界點(diǎn)的Y坐標(biāo)Y1和Y2，根據(jù)地圖比例尺知道18mm相當(dāng)于40km，試由測(cè)量數(shù)據(jù)計(jì)算瑞士國(guó)土近似面積，與其精確值41228km2比較。X7133448566191Y14445475050383030343634414546Y24459707293100110110110117118116118118X96101104118142146150157158Y143373328326555545250666668Y2121124121121121116122838182868568 x=[7,13,34,48,56,61,91,96,101,104,118,142,146,150,157,158]。 y1=[44,45,47,50,50,38,30,30,34,36,34,41,45,46,43,37,33,28,32,65,55,54,52,50,66,66,68]。 y2=[44,59,70,72,93,100,110,110,110,117,118,116,118,118,121,124,121,121,121,116,122,83,81,82,86,85,68]。 X=x./18*40。Y1=y1./18*40。Y2=y2./18*40。 t1=trapz(X,Y1),t2=trapz(X,Y2), t=t2t1t1 = +004t2 = +004t = +004 expt=t41228expt = +003習(xí)題（4）利用梯形法和辛普森法求定積分的值，并對(duì)結(jié)果進(jìn)行比較。如果積分區(qū)間改為5~5結(jié)果有何不同？梯形積分中改變自變量x的維數(shù)，結(jié)果有何不同？12π33ex22dx x=linspace(3,3)。y=exp(x.^2/2)。 t=(1/2*pi)*trapz(x,y)t = q=(1/2*pi)*quad(39。exp(x.^2/2)39。,3,3)q =
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1