【正文】 縱向變位時油位深度與油量的函數(shù)關(guān)系式計算儲油罐中油量的體積,我們將儲油罐分為三個部分，中間部分為圓柱體，兩端為球冠體。圓缺面積為，則： （）其中是圓缺的半徑，圓缺的高度。其中； （）記為：（4）當時，極限容積為：的計算下面計算在右邊球冠中溶液的體積，以右端球冠的球心為坐標原點，建立三維坐標系，球的半徑，在直角三角形中，計算得。 （）為了便于后面數(shù)值解法書寫方便，我們記：所以（16）式變?yōu)椋涸谶M行數(shù)值求解時，我們觀察出實驗數(shù)據(jù)中且，于是我們放棄了計算中的第一種情況和第二種情況中的情況（4）。計算中采用本文改進三點高斯公式，即：則有：二、建立疊加了橫向變位的油位深度和縱向變位時油位深度的函數(shù)關(guān)系，其中與的位置關(guān)系（如圖9）所示，（圖9）可得函數(shù)關(guān)系為: （）綜合可得三、找出函數(shù)關(guān)系將帶入到第二步的中，得到油液容積與疊加了橫向變位的油位深度的函數(shù)（具體程序見附件程序四）：。 算法描述： 第一步 賦初始值，定義步長； 第二步 。并與實際體積對比，根據(jù)公式：計算出相對誤差（具體程序見附件程序六）。六、實際儲油罐罐容表的標定根據(jù)函數(shù)，計算不同時的標定值（具體程序見附件程序程序七）。最后將該公式應(yīng)用到2010年數(shù)學建模A題的求解中，取得了比較好的結(jié)果。為了指導我們的畢業(yè)論文他放棄了自己的休息時間，他的這種無私奉獻的敬業(yè)精神令人欽佩，在此我向他表示我誠摯的謝意。在本論文的寫作過程中，我的導師唐利明老師傾注了大量的心血，從選題到開題報告，從寫作提綱，到一遍又一遍地指出每稿中的具體問題，嚴格把關(guān)，循循善誘，在此我表示衷心感謝。重慶科技學院本科生畢業(yè)設(shè)計 附錄附錄程序一、計算柯特斯公式的系數(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。程序七、罐容表標定的程序：012：load 。附表一：流水號油罐號進油量/L出油量/L顯示油高/mm顯示油量容積/L采集時間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