1、數(shù)值分析課程設(shè)計(jì)報(bào)告第 15頁，共 15頁一、課程設(shè)計(jì)目的：1 .訓(xùn)練學(xué)生靈活應(yīng)用所學(xué)數(shù)值分析知識，獨(dú)立完成問題分析， 結(jié)合數(shù)值分析理論知識，編寫程序求解指定問題。2 .初步掌握解決實(shí)際問題過程中的對問題的分析、系統(tǒng)設(shè)計(jì)、程序編碼、測試等基本 方法和技能；3 .提高綜合運(yùn)用所學(xué)的理論知識和方法獨(dú)立分析和解決問題的能力；4 .訓(xùn)練用數(shù)值分析的思想方法和編程應(yīng)用技能模擬解決實(shí)際問題，鞏固、深化學(xué)生的理論知識，提高學(xué)生對數(shù)值分析的認(rèn)知水平和編程水平，并在此過程中培養(yǎng)他們嚴(yán)謹(jǐn)?shù)目茖W(xué)態(tài)度和良好的工作作風(fēng)二、課程設(shè)計(jì)任務(wù)與要求：課程設(shè)計(jì)題目：計(jì)算水塔的水流量【問題描述】某居民區(qū)的民用自來水是由一個(gè)圓柱形的

2、水塔提供，水塔高12.2米，直徑17.4米。水塔是由水泵根據(jù)水塔內(nèi)水位高低自動(dòng)加水，一般每天水泵工作兩次， 現(xiàn)在需要了解該居民區(qū)用水規(guī)律與水泵的工作功率。按照設(shè)計(jì)，當(dāng)水塔的水位降至最低水位，約8.2米時(shí)，水泵自動(dòng)啟動(dòng)加水；當(dāng)水位升高到一個(gè)最高水位，約 10.8米時(shí)，水泵停止工作?？梢钥紤]采用用水率 （單位時(shí)間的用水量） 來反映用水規(guī)律，并通過間隔一段時(shí)間測量 水塔里的水位來估算用水率，原始數(shù)據(jù)表是某一天的測量記錄數(shù)據(jù)，測量了28個(gè)時(shí)刻，但是由于其中有3個(gè)時(shí)刻遇到水泵正在向水塔供水，而無水位記錄。試建立合適的數(shù)學(xué)模型，推算任意時(shí)刻的用水率、一天的總用水量。進(jìn)一步：可自己增加一些新的計(jì)算功能?！締?/p>

3、題假設(shè)】1 .水塔中水流量是時(shí)間的連續(xù)光滑函數(shù)，與水泵工作與否無關(guān)， 并忽略水位高度對水流速度的影響。2 .水泵工作與否完全取決于水塔內(nèi)水位的高度。3 .水塔為標(biāo)準(zhǔn)的圓柱體。體積 V=PI*D*D*h/4 其中D為底面直徑，h為水位高。4 . 水泵第一次供水時(shí)間段為 8.967,10.954，第二次供水時(shí)間段為20.839,22.958。 【實(shí)驗(yàn)數(shù)據(jù)】原始數(shù)據(jù)（單位：時(shí)刻（小時(shí)），水塔中水位（米）時(shí)刻t00.9211.8432.9493.8714.9785.900水位h9.6779.4799.3089.1258.9828.8148.686時(shí)刻t7.0067.9288.9679.981110.9

4、2510.95412.032水位h8.5258.3888.220泵水泵水10.82010.500時(shí)刻t12.95413.87514.98215.90316.82617.93119.037水位h10.2109.9369.6539.4099.1808.9218.662時(shí)刻t19.95920.83922.01522.95823.88024.98625.908水位h8.4338.220泵水10.82010.59110.35410.180【實(shí)現(xiàn)提示】由問題的要求，關(guān)鍵在于確定用水率函數(shù)，即單位時(shí)間內(nèi)用水體積，記為f，又稱水流速度。如果能夠通過測量數(shù)據(jù)，產(chǎn)生若干個(gè)時(shí)刻的用水率，也就是f(t)在若干個(gè)點(diǎn)的函

5、數(shù)值，則f(t)的計(jì)算問題就可以轉(zhuǎn)化為插值或擬合問題?！緦?shí)現(xiàn)要求】1、在處理每個(gè)題目時(shí)，要求分別從數(shù)據(jù)處理階段和程序設(shè)計(jì)階段兩個(gè)主要階段實(shí)現(xiàn)課程設(shè)計(jì)，詳細(xì)的通過文字以及插圖等形式，按需求分析、數(shù)據(jù)處理、算法設(shè)計(jì)、代碼、計(jì)算結(jié)果和程序執(zhí)行的截圖 等若干步驟完成題目，最終寫出完整的分析報(bào)告。前期準(zhǔn)備工作 完備與否直接影響到后序上機(jī)調(diào)試工作的效率。在程序設(shè)計(jì)階段應(yīng)盡量利用已有的標(biāo)準(zhǔn)函 數(shù)，加大代碼的重用率。2、設(shè)計(jì)的題目要求達(dá)到一定工作量，并具有一定的深度和難度。3、程序設(shè)計(jì)語言推薦使用 C/C+,程序書寫規(guī)范，源程序需加必要的注釋；4、每位同學(xué)需提交可獨(dú)立運(yùn)行的程序；5、每位同學(xué)需獨(dú)立提交設(shè)計(jì)報(bào)告

6、書(每人一份)，要求編排格式統(tǒng)一、規(guī)范、內(nèi)容充實(shí)； 6、課程設(shè)計(jì)實(shí)踐作為培養(yǎng)學(xué)生動(dòng)手能力的一種手段，單獨(dú)考核。三、課程設(shè)計(jì)說明書【需求分析】根據(jù)一組關(guān)于若干個(gè)時(shí)刻的水塔內(nèi)水位高度的數(shù)據(jù)，經(jīng)過必要的處理，得到一組若干個(gè)時(shí)刻的用水率的數(shù)據(jù)，然后根據(jù)這一組數(shù)據(jù)提出假設(shè)，建立合適的數(shù)學(xué)模型，通過C語言程序?qū)⒃撃Ｐ偷囊恍﹨?shù)求出來，最后得到一條能夠擬合這一組若干個(gè)時(shí)刻的用水率的光滑的函數(shù)，進(jìn)而推算出任意一個(gè)時(shí)刻的用水率，并計(jì)算出一天的總用水量。如果還有余力可以增加更多的實(shí)際計(jì)算要求?！緮?shù)據(jù)處理概要】1 .關(guān)于原始數(shù)據(jù)中有三個(gè)時(shí)刻水泵是在泵水的處理：由問題假設(shè)可以知道水塔中水 流量是時(shí)間的連續(xù)光滑函數(shù)，與

7、水泵工作與否無關(guān)，如果將這三個(gè)時(shí)刻的數(shù)據(jù)忽 略不計(jì)，是不影響設(shè)計(jì)結(jié)果的。所以設(shè)計(jì)中忽略這三個(gè)時(shí)刻的數(shù)據(jù)。2 .關(guān)于確定求解哪些時(shí)刻的用水率問題：我采用每兩個(gè)時(shí)刻的中間時(shí)刻和這兩個(gè)時(shí)刻之間的時(shí)間段內(nèi)的平均用水率作為最后的若干時(shí)刻的用水率。3 .關(guān)于數(shù)據(jù)處理過程中三個(gè)無數(shù)據(jù)時(shí)刻的處理：按照處理需要，適當(dāng)?shù)纳釛夁吘壍臄?shù)據(jù)和錯(cuò)誤的結(jié)尾數(shù)據(jù) (比如要計(jì)算兩個(gè)值的差， 最后一項(xiàng)后面已經(jīng)沒有數(shù)據(jù)了， 系統(tǒng)會(huì)在該數(shù)據(jù)后產(chǎn)生一個(gè)零，所以是錯(cuò)誤的數(shù)據(jù))?！驹敿?xì)設(shè)計(jì)】本設(shè)計(jì)報(bào)告按照我的設(shè)計(jì)步驟完成:1、將所給的數(shù)據(jù)中的時(shí)刻輸入到Excel表格中，如下圖所示:E I . P G |K| 小庫豐水的體和外限出間用水體耙

8、中間時(shí)利用水率l85$5T«glcllla1314lElGLLllJ13lgM22233現(xiàn)SS6ZTE829小期A B _ C | 上麗1常劃差普解函市間時(shí)到京位00. 9211.0452. 9493. ST14. 9T87. OOG 7. 9SS S. 967 亂物10. 925 10, 054 12, 032 12. 95-1 13. 8T5 U.9RS 15. 9Q3 16. 82& 37.931.13 03?19.湖 20. 339 22. 01E.22, 95823. BE 2% 4脂 國婚，h- Sluegtl / Sfaeeta / Sh»t3/2、計(jì)

9、算出每兩個(gè)相鄰的時(shí)刻之間的時(shí)間差(表格列B(n尸A(n+1)- A(n),然后除以2 (表格列C(n尸B(n)/2 )得到相鄰兩個(gè)時(shí)刻的中間時(shí)刻:3、將原始數(shù)據(jù)中的水位數(shù)據(jù)輸入到該表格中:4、根據(jù)公式：V=PI*D*D*h/4 計(jì)算各個(gè)時(shí)刻水塔中水的體積(表格列E(n尸PI* D* D*表格列D(n)/4),然后求出各個(gè)時(shí)間段內(nèi)用水的體積(表格列F(n尸D(n+1)-D(n)：號1=J窩件# >M：i根型界熱相域制工田功近陰百口爐現(xiàn)助卿片的0澗胃P 4;Jj U 3町剪w口，此必蚓氣，修嘉章單一 必LJ11干,ft水庫中不班ABCbE|FHI占時(shí)多.時(shí)型差普段時(shí)間口間的到水位水陣中不由體

10、積各裳li阿用乩體斗中間時(shí)刻同也中2口 Q. 9Z1(.46059. 6T72301. 06623147. 081 S552930.奠 10, 9221.3S2-47勺22s 1始必耨40-&61G023_J_1.8431,1 好2. 3轉(zhuǎn)9.3閽2215.33277343.515048072169. EO77Z534. Q035G2161 5Z. 9dy 0. 92214L'?. L253. STI 1.1074我第 S. 32213&.圜如的胸.9朝剖加9工腐礴22匆7旃T54融74. 9相 0, 932瓦4融 g, 81435.91.1C66. 4K3 S. 6E

11、SZC65.dL916T33. 2B3730E2mT. 0060. 922孔4利 S. 525加卻.霓32. 5T戰(zhàn)油卸199d- 55S9T39. 9蛇£如1E6LD7. 93S1, 059比 4475 S. 3E311a 9&T L 014當(dāng)474皮a1韭1610356% 9B10,綱司10.例13.IO,?250.02910. 93KL410.9541.0T31L4弟1也峻2572.85641T6.的 13股副IE12. 0320. 922】£4驕10.5網(wǎng)9隊(duì)傳向53毋,95耽設(shè)機(jī)工427. EO67S165.1536T3M1612. 9540. '

12、?2113LU14510.21L；13. ST5 1.1071.4, 428523&?. 65310267. 2937«2ff714 . 9830, 921飄/ y, 6轉(zhuǎn)2295> 36中融5&國口嗎ML¥IS. 9的0.必16. 365'?. 4a3修圻.33耍襠M.45黃S6閱,2C-16. S2§ I.IDS17.37ES9.1W2132- 83M1B&L5&6ST15221793.11.105】£. 4網(wǎng)E. 92L2121.29934761.536E71321221 女 0J70. 92219,

13、4鮑 8, 66-2即5曳T122將Z. 45地防圜_2S.1軋959。,聆品口, 33S. 4處凰口區(qū)盤馴198L 5粥6呢切1954.610354；2420. 3391. 17621.427B.222E-22. 0150,綱322,堿更曜脂醇鼠因r 482轉(zhuǎn)腳比22, 9隔0, 9舞25.4J91爪混STZ3,BS 1.1C624. 43310. 591251E. '10365»56. 35555406j 1324.另% 0, 3225, 447 10. 3抬M&2,Q的134L匿為的M熊90S_ 1<Ql 12MW. H Sbeetl /ShflatS /

14、5、求中間時(shí)刻的用水率，即求該時(shí)間段(時(shí)亥I差)內(nèi)的平均用水量(表格列G(n尸F(xiàn)(n)/B(n):<| Merosofl Ltd型環(huán)吉MF小1 - a.文件舊 劇15 ffiEd）防風(fēng)工具B 范優(yōu)面 £«1 WlOlBi,串一fl£ 口E S 3 jd5雅，nsF ：F_， a ,aG1i'巾枷閩用水率ABcI?EFHL.12時(shí)則0時(shí)到差0. 921&段時(shí)間中間時(shí)豺 水位木即中禾的律相.各段時(shí)間用水體和也的05 9.67723m.崎鼠乳47.0S13552'9中間片刻艮永率5L 1203050. 9210. 9221.3824T922

15、53. &3437S鐘,兩1狗的4110152053上5： 7L0L1L2L5 ,U_ _LL.L£ LT L£ 及21 '比 23 別1,即3 gqg 3S714. 9T85,97, OK 7. 923SL967 務(wù)死110.泗 1 口涮32, OK 12. 954J 13 S7S14. 9B215.迎.1/106 Q. 922 LI婚 0.迤 1.106 0. 922 1. 039 1.014 0.甄4 0.婚 1.QT8,。陽 Q. 9ZI 1.1婚 0. 921 0.例23我9,303 我1工笈2m1.419.125Z169. B077Z5工養(yǎng)方g.

16、9822135.的區(qū) 4割8.8142C95. 85的整6. 4盼亂8蹣 凰6氏n91肝7. 4例B.525前室3M圜比 44T5E.3E31994. E5ES9T9. 4741B. 221 犯k 8103法10. 4E-3 位 939611. 4310.82372.85641用,4卷0,5 國為6岡E13. Jld5ID. ZL2-027.E067E111 式吃9. -93665310215u 44259.6532295.知羽匏1粒誨死S, 41?32237.處第T5好,515O4SOT 風(fēng)口。脂團(tuán)5 料.馴船q 口的 圜.隴T同制 38, 283730 52. 57楣2費(fèi) 39. q附工小

17、。隨慍 Q91837N65.15367S&9 67.23762875al 020D64 L 國YR256吃39, SM5需】 3E. ES021-925 球.尚的38M 33.01166479 3蟲同強(qiáng)掰83 55. 3327沂4733, «8網(wǎng)Ta 5gM的判 刊.門四同事| TA.T4Z3Z197 網(wǎng).789£3OTO262.的圈12方EB.995鋼非 2I1I0瞄 17 911 1 且 017 1SL 的 20.咽1,1051.1D6 09220.081.1T61T.W 脂&,13 282-83Mlfl】2 4M B. 9ZI 2121. 299347

18、 19.拈S S. 662 就E1曳71愛花 20. 3的 S. 433 2C0S. 25翼 19 21.427B.22 1玄5立也中將1，5前混 &1.5SE3T1323.止，二:上6n電 S4B&&2&1的,T34?251g 生.的43321】 EI皚的的M57. 55s2粥乳HE22, 0150,綱3浣.招的就 2T _2L ¥22 953Z3,B824.雌25. 聆0,9221.1D60.132223.4391Q.葭Z572. E5694124. 43310.591Z51B. J036BU25. 4710. 3H2462 04813Uh IS加

19、2優(yōu)即31%3.5325603訪.35555406如4星L咤汨別網(wǎng) 51. 9543B3S5 4L &7&23183H 4 就討* X5hagtX 5h""t2 Z 5h?(it3/l< L束.t曰詡niR研國£_T6、將中間時(shí)刻和中間時(shí)刻的用水率復(fù)制到另一張表中，以便處理（這時(shí)共產(chǎn)生22組數(shù)據(jù)）：到這里所有的數(shù)據(jù)處理已經(jīng)結(jié)束。接下來是研究數(shù)據(jù)所包含的規(guī)律，根據(jù)數(shù)據(jù)提出適當(dāng)?shù)臄?shù)學(xué)模型，擬合一條光滑的連續(xù)的曲線。7、在Excel中有擬合曲線的工具，所以我先用Excel中的工具對這組數(shù)據(jù)進(jìn)行了擬合*畫出散點(diǎn)圖：圖表向?qū)?4步鴕1 -圖表賽型標(biāo)段類

20、壬自京義類用圖表芟型位)二子圖麥類型d)按T不放可萱看示例四觸消 ' K-QU占咸g*確定橫縱坐標(biāo)的意義:圖吉向橋 4步景之3 -圖看姻標(biāo)題坐標(biāo)軸網(wǎng)格維I囪例F塔丁制總cst制女看三水塔水繞里數(shù)值軸時(shí)更I 小時(shí)3數(shù)值on油5：用水率t立方米/小暇) 次分類的軸(x).次數(shù)值軸yy取消 <LT®I 丁一步心L完成但)*形成散點(diǎn)圖:水塔水流量2BW/米松摑票旺*通過添加趨勢線，來擬合光滑曲線:*得到擬合的函數(shù)曲線，以及表達(dá)式用水率函數(shù)OOOOOOOOO 87654321 £w至后)si*系列1多項(xiàng)式（系列Dy = -0. OOOlx5 + 0. 0104x4 -

21、0.32541 + $ 35弱1-21.169x + 64. 71= 0. 78371015202530時(shí)刻（小時(shí)）在Excel中就可以很方便的擬合出比較符合要求的曲線，并且可以得到函數(shù)表達(dá)式。然后開始嘗試?yán)镁幊探鉀Q這個(gè)問題，也就是說用編程的方法得到相同的函數(shù)表達(dá)式！最小二乘法通常用于曲線擬合，所謂擬合是指已知某函數(shù)的若干離散函數(shù)值f1,f2,fn通過調(diào)整該函數(shù)中若干待定系數(shù)f（入1, 2,，入版導(dǎo)該函數(shù)與已知點(diǎn)集的差別（最小二乘意義）最小。如果待定函數(shù)是線性，就叫線性擬合或者線性回歸 （主要在統(tǒng)計(jì)中）,否則叫作非線性擬合或者非線性回歸。表達(dá)式也可以是分段函數(shù)，這種情況下叫作樣條擬合。通過求

22、解該函數(shù)中待而插值是指已知某函數(shù)的在若干離散點(diǎn)上的函數(shù)值或者導(dǎo)數(shù)信息,定形式的插值函數(shù)以及待定系數(shù)，使得該函數(shù)在給定離散點(diǎn)上滿足約束。插值函數(shù)又叫作基否則叫作分域基。如果約束條件中函數(shù)，如果該基函數(shù)定義在整個(gè)定義域上，叫作全域基,只有函數(shù)值的約束，叫作Lagrange插值，否則叫作Hermite插值。從幾何意義上將，擬合是給定了空間中的一些點(diǎn)，找到一個(gè)已知形式未知參數(shù)的連續(xù)曲面來最大限度地逼近這些點(diǎn)而插值是找到一個(gè)（或幾個(gè)分片光滑的）連續(xù)曲面來穿過這些 點(diǎn)。可以看出利用最小二乘法擬合比較方便，我采用的函數(shù)形式為多項(xiàng)式函數(shù)，最高次為5。擬合部分代碼執(zhí)行結(jié)果如下:經(jīng)過比較，可以發(fā)現(xiàn)在誤差允許的范圍

23、內(nèi)通過Excel擬合出來的函數(shù)和通過編程而擬合出來的多項(xiàng)式是相統(tǒng)一的。然后再將該函數(shù)應(yīng)用到Excel中，利用程序代碼所求出來的多項(xiàng)式，畫出一條光滑連續(xù)的曲線，如果這條曲線和最上面在Excel中，純粹用Excel擬合出來的曲線相同的話表明擬合成功！1、按照遞增的順序，將從 0-24的數(shù)據(jù)輸入到Excel中，如下圖所示。并在第二列中輸入公 式：=-0.00015*A2*A2*A2*A2*A2+0.012193*A2*A2*A2*A2-0.366317*A2*A2*A2+4.760591*A2*A2-22.665987*A2+66.02523如下圖所示。將數(shù)據(jù)填充完畢:利用和上面的相同的方法，直接選

24、擇光滑的曲線圖，而不是選擇散點(diǎn)圖:第向?qū)?d步騾之1-圖表美莖標(biāo)準(zhǔn)類型目定義類型是圖圖 理胱腕 舉柱條面XT三期肖一步上|越g然后就畫出了一條曲線:代碼擬合用水率系列1OOOOOOOO 7 6 5 4 3 2 1 (后W/米歸揭)助考旺51015202530時(shí)刻（小時(shí)）經(jīng)過比較，可以知道在誤差允許的范圍內(nèi)數(shù)學(xué)模型已經(jīng)成功建立！到這里該問題的數(shù)學(xué)模型已經(jīng)建立好了(采用編程所得的結(jié)果進(jìn)行建模和進(jìn)一步的計(jì)算)：f(t)=-0.00015t5+0.012193f-0.366317t3+4.760591t2-22.665987t+66.025230最后再通過編程計(jì)算任意時(shí)刻的用水率、一天的總用水量和水泵

25、的工作功率：任意時(shí)刻的用水率的算法簡單，只要給出任意時(shí)間，根據(jù)上面的式子就可以求出；一天的總用水量為上式在 0-24內(nèi)的積分V= / 024f(t)dt=-0.000025t6+0.0024386t5-0.09157925t4+1.58686367t3-11.3329935f+66.02523t+C024水泵的工作功率的計(jì)算方法是用一次的總泵水量除以泵水時(shí)間水泵泵一次水泵水量為：V=2*PI*D*D*(h2-h1)/4=2*3.1415926*17.4*17.4*(10.82-8.22)/4泵水時(shí)間為：T=t2-t1=(10.954-8.967)+(22.958-20.839)所以水泵的工作功

26、率為：P= p Vg(h2-h1)/2T其中p為水的密度(kg/m3), V為上面的體積(m3), g為重力單位(9.8N/kg)。然后將程序加以修改，將上述的三個(gè)求解功能加進(jìn)去就完成了本次設(shè)計(jì)。最終源代碼：#include<stdio.h>#include<math.h>#include<iostream.h>void yongShuiLv();void yongShuiLiang();void gongLv();void Approx(float,float,int,int,float);void main()int i;float a6; 多項(xiàng)式的待定

27、系數(shù)個(gè)數(shù)float x22=0.4605,1.382,2.396341,4.4245,5.439,6.4537467, 8.4475,11.493,12.493,13.4145,14.4285,15.4425,16.3645, 17.3785,18.484,19.498,20.399,23.419,24.433,25.447; 22 個(gè)時(shí)亥U的值float y22=51.12036405,44.10152093,39.3445281,36.88021926,36.08693844, 33.01166479,34.61458483,35.33279747,38.44874,70.58616637

28、,74.79205304, 70.74232197,60.78930702,62.99681227,58.99594462,55.73472518,55.68433211, 59.05993154,57.55529831,59.05993154,50.95438885,44.87523183; 1122 個(gè)時(shí)刻所對應(yīng)的用水率Approx(x,y,22,5,a);for(i=0;i<=5;i+)printf("a%d=%fn",i,ai);printf("擬合多項(xiàng)式為nf(t)=(%f)*t*t*t*t*t+(%f)*t*t*t*t+(%f)*t*t*t+(%

29、f)*t*t+(%f)*t+(%f)n",a5,a4,a3,a2,a1,a 0);yongShuiLiang();gongLv();char s;while(1)/通過判斷用戶輸入的字符是Y/N,來控制是否繼續(xù)求某一時(shí)刻的用水率cout<<"您是否還要查詢某日刻的用水率( Y/N ) ： "<<endl;cin>>s;switch(s) case 'Y':yongShuiLv();break;case 'N':break;if(s='N')break;void Approx(flo

30、at x,float y,int m,int n,float a) (int i,j,t;float *c=new float(n+1)*(n+2);float power(int,float);void ColPivot(float *,int,float);for(i=0;i<=n;i+)(for(j=0;j<=n;j+)(*(c+i*(n+2)+j)=0;for(t=0;t<=m-1;t+)*(c+i*(n+2)+j)+=power(i+j,xt);*(c+i*(n+2)+n+1)=0;for(j=0;j<=m-1;j+)*(c+i*(n+2)+n+1)+=yj*

31、power(i,xj);ColPivot(c,n+1,a);delete c;void ColPivot(float *c,int n,float x)(int i,j,t,k;float p;for(i=0;i<=n-2;i+)(k=i;for(j=i+1;j<=n-1;j+)if(fabs(*(c+j*(n+1)+i)>(fabs(*(c+k*(n+1)+i)k=j;if(k!=i)for(j=i;j<=n;j+)(p=*(c+i*(n+1)+j);*(c+i*(n+1)+j)=*(c+k*(n+1)+j);*(c+k*(n+1)+j)=p;for(j=i+1;j&

32、lt;=n-1;j+)(p=(*(c+j*(n+1)+i)/(*(c+i*(n+1)+i);for(t=i;t<=n;t+)*(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t);for(i=n-1;i>=0;i-)(for(j=n-1;j>=i+1;j-)(*(c+i*(n+1)+n)-=xj*(*(c+i*(n+1)+j);xi=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i);float power(int i,float v)(float a=1;while(i-)a*=v;return a;void yongShuiLv() 計(jì)算某時(shí)刻的用水率函數(shù)float t;cout<<"請輸入任意一個(gè)時(shí)刻，程序?qū)⒕退愠鲈摃r(shí)刻的用水率"<<endl;cin>>t;floatf=-0.00015*t*t*t*t*t+0.012193*t*t*t*t-0.