高精度數(shù)值積分公式的構(gòu)造及其應(yīng)用數(shù)學(xué)專業(yè)畢業(yè)設(shè)計(jì)畢業(yè)論文-預(yù)覽頁(yè)
【正文】 縱向變位時(shí)油位深度與油量的函數(shù)關(guān)系式計(jì)算儲(chǔ)油罐中油量的體積,我們將儲(chǔ)油罐分為三個(gè)部分，中間部分為圓柱體，兩端為球冠體。圓缺面積為，則： （）其中是圓缺的半徑，圓缺的高度。其中； （）記為：（4）當(dāng)時(shí)，極限容積為：的計(jì)算下面計(jì)算在右邊球冠中溶液的體積，以右端球冠的球心為坐標(biāo)原點(diǎn)，建立三維坐標(biāo)系，球的半徑，在直角三角形中，計(jì)算得。 （）為了便于后面數(shù)值解法書寫方便，我們記：所以（16）式變?yōu)椋涸谶M(jìn)行數(shù)值求解時(shí)，我們觀察出實(shí)驗(yàn)數(shù)據(jù)中且，于是我們放棄了計(jì)算中的第一種情況和第二種情況中的情況（4）。計(jì)算中采用本文改進(jìn)三點(diǎn)高斯公式，即：則有：二、建立疊加了橫向變位的油位深度和縱向變位時(shí)油位深度的函數(shù)關(guān)系，其中與的位置關(guān)系（如圖9）所示，（圖9）可得函數(shù)關(guān)系為: （）綜合可得三、找出函數(shù)關(guān)系將帶入到第二步的中，得到油液容積與疊加了橫向變位的油位深度的函數(shù)（具體程序見(jiàn)附件程序四）：。 算法描述： 第一步 賦初始值，定義步長(zhǎng)； 第二步 。并與實(shí)際體積對(duì)比，根據(jù)公式：計(jì)算出相對(duì)誤差（具體程序見(jiàn)附件程序六）。六、實(shí)際儲(chǔ)油罐罐容表的標(biāo)定根據(jù)函數(shù)，計(jì)算不同時(shí)的標(biāo)定值（具體程序見(jiàn)附件程序程序七）。最后將該公式應(yīng)用到2010年數(shù)學(xué)建模A題的求解中，取得了比較好的結(jié)果。為了指導(dǎo)我們的畢業(yè)論文他放棄了自己的休息時(shí)間，他的這種無(wú)私奉獻(xiàn)的敬業(yè)精神令人欽佩，在此我向他表示我誠(chéng)摯的謝意。在本論文的寫作過(guò)程中，我的導(dǎo)師唐利明老師傾注了大量的心血，從選題到開(kāi)題報(bào)告，從寫作提綱，到一遍又一遍地指出每稿中的具體問(wèn)題，嚴(yán)格把關(guān)，循循善誘，在此我表示衷心感謝。重慶科技學(xué)院本科生畢業(yè)設(shè)計(jì) 附錄附錄程序一、計(jì)算柯特斯公式的系數(shù)程序：function x=NewtonCotesxishu(n)。for i=1:n+1 l=1。for i=1:8 N(i,:)=NewtonCotesxishu(i)。B=simplify(g)程序三、數(shù)值算例程序clear。a=0。c3=sqrt(15)*(ba)/10。a=pi*A/180。L2=6。h1=(z) sqrt(^2z^2)。f1=(z) R^2*acos((Rh0(z))/R)(Rh0(z))*sqrt(2*h0(z)*Rh0(z)^2)。g2=r(z)^2*acos((r(z)h1(z))/r(z))(r(z)h1(z))*sqrt(2*h1(z)*r(z)h1(z)^2)。 bb=L1+h/tan(a)。 aa=。 V2=(bbaa)/2*(5/9*f2(z1z2)+8/9*f2(z1)+5/9*f2(z1+z2))+(bbaa)^7*subs(g22,z,(aa+bb)*cc)/2016000。h2*RL1*tan(a) aa=0。 V1=(bbaa)/2*(5/9*f1(z1z2)+8/9*f1(z1)+5/9*f1(z1+z2))+(bbaa)^7*subs(g11,z,(aa+bb)*cc)/2016000。 z2=sqrt(15)*(bbaa)/10。 z1=(aa+bb)/2。 else if h=2*Ramp。 z1=(aa+bb)/2。 bb=。 aa=。 V3=(bbaa)/2*(5/9*f2(z1z2)+8/9*f2(z1)+5/9*f2(z1+z2))+(bbaa)^7*subs(g22,z,(aa+bb)*cc)/2016000。format longdata2=data2./1000。min=inf。 for k=1:n sum=sum+(test2(A,B,data2(k))test2(A,B,data2(k+1))V2(k))^2。yanzh=yanzh./1000。程序七、罐容表標(biāo)定的程序：012：load 。附表一：流水號(hào)油罐號(hào)進(jìn)油量/L出油量/L顯示油高/mm顯示油量容積/L采集時(shí)間201020 20100801 08:00:49202020 20100801 08:15:42203020 20100801 08:23:41204020 20100801 08:38:14205020 20100801 08:53:08206020 20100801 09:08:01207020 20100801 09:32:06208020 20100801 09:56:18209020 20100801 10:10:09210020 20100801 10:26:28211020 20100801 10:48:35212020 20100801 11:02:48213020 20100801 11:36:04214020 20100801 13:16:24215020 20100801 13:24:25216020 20100801 13:37:33217020 20100801 14:00:03218020 20100801 14:36:02219020 20100801 14:52:08220020 20100801 15:06:14221020 20100801 15:16:43222020 20100801 15:48:22223020 20100801 16:06:31224020 20100801 16:25:54225020 20100801 16:41:24226020 20100801 17:06:08227020 20100801 17:45:54228020 20100801 18:03:34229020 20100801 18:14:52230020 20100801 18:35:24231020 20100801 19:12:55232020 20100801 20:34:44233020 20100801 21:05:57234020 20100802 07:18:34235020 20100802 07:23:30236020 20100802 07:39:21237020 20100802 07:48:36238020 20100802 08:04:44239020 20100802 08:28:32240020 20100802 08:33:36241020 20100802 08:46:28242020 20100802 08:58:39243020 20100802 09:03:54244020 20100802 09:14:47245020 20100802 09:25:55246020 20100802 09:35:54247020 20100802 09:59:33248020 20100802 10:13:32249020 20100802 10:18:24250020 20100802 10:26:36251020 20100802 10:38:33252020 20100802 10:46:54253020 20100802 10:58:37254020 20100802 11:05:34255020 20100802 11:26:38256020 20100802 11:58:37257020 20100802 12:09:46258020 20100802 12:36:44259020 20100802 12:55:32260020 20100802 13:18:34261020 20100802 13:26:42262020 20100802 13:38:24263020 20100802 13:45:37264020 20100802 13:55:34265020 20100802 14:18:36266020 20100802 14:27:33267020 20100802 14:48:24
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1