廣告費用與銷售價格調(diào)整程序設(shè)計

1、數(shù)學(xué)建模實踐數(shù)學(xué)建模課程設(shè)計（程序設(shè)計和論文）題目 1 對函數(shù)進行麥克勞林展開及誤差分析 2 無變位油罐中油量確定及誤差分析 3 評卷成績調(diào)整程序設(shè)計 4 廣告費用與銷售價格調(diào)整程序設(shè)計 班 級 學(xué) 號 學(xué) 生 姓 名 指 導(dǎo) 教 師 41沈陽航空航天大學(xué)課 程 設(shè) 計 任 務(wù) 書課 程 名 稱 數(shù)學(xué)建模實踐 院（系） 理學(xué)院 專業(yè) 信息與計算科學(xué) 班級 學(xué)號 姓名課程設(shè)計題目 1 對函數(shù)進行麥克勞林展開及誤差分析 2 無變位油罐中油量確定及誤差分析 3 評卷成績調(diào)整程序設(shè)計 4 廣告費用與銷售價格的調(diào)整程序設(shè)計 課程設(shè)計時間: 2021 年 6 月 27 日至 2021 年 7 月 15 日

2、課程設(shè)計的內(nèi)容及要求：內(nèi)容1（1）求函數(shù)（2）編寫對任意固定的n計算多項式函數(shù)值的函數(shù)M文件（3）任取n，在同一平面內(nèi)畫出函數(shù)的圖形，并進行比較。2無變位油罐中油量確定設(shè)油罐中油量V與高度h的關(guān)系是其中，（1）編寫計算體積V(h)的函數(shù)M文件fv；（2）根據(jù)“無變位實驗采集數(shù)據(jù)表”中的無變位進油表中的數(shù)據(jù)計算公式V(h)與實驗數(shù)據(jù)之間的誤差WC(h)，并用多項式擬合確定函數(shù)WC(h)表達式。（3）用誤差WC(h)調(diào)整V(h)，并用“無變位實驗采集數(shù)據(jù)表”中的無變位出油表中的數(shù)據(jù)檢驗調(diào)整結(jié)果。3評卷成績調(diào)整程序設(shè)計 設(shè)個專家分別對名學(xué)生的試卷進行評閱，設(shè)表示教師對學(xué)生的試卷所給定的成績，這樣形成

3、成績矩陣。由于各專家的評分標(biāo)準(zhǔn)不一致，因此需要對成績進行一致性調(diào)整，具體方法如下： 設(shè)分別表示整體成績，是教師j的平均成績和標(biāo)準(zhǔn)差，即第j列數(shù)據(jù)的平均值和標(biāo)準(zhǔn)差。調(diào)整后的成績?yōu)?形成調(diào)整后的成績矩陣，則的平均值就是第i個學(xué)生的最后綜合成績。而是綜合成績向量，依此確定學(xué)生獲獎情況。 （1）編寫函數(shù)M文件，收入成績矩陣，輸出是綜合成績向量。 （2）根據(jù)下表是成績數(shù)據(jù) 學(xué)生編號專家1專家2專家3專家419090809729590809038585909048590809059580907069080907578090708588585708099570707510906060601175657060

4、127570606513607060501465605050156560505016606050501760505050用上述方法計算綜合成績向量，并由此確定1個特等獎，1個一等獎，2個二等獎；3個三等獎。4. 廣告模型某裝飾材料公司欲以每桶2元的價錢購進一批彩漆以供日后銷售。為了盡快收回資金并獲得較多的贏利，公司經(jīng)理李先生打算做廣告，于是便找到廣告公司的王經(jīng)理進行咨詢。李經(jīng)理認(rèn)為，隨彩漆售價的提高，預(yù)期銷售量將減少，并對此進行了估算（見表2）。他問王經(jīng)理廣告有多大的效應(yīng)。王經(jīng)理說“投入一定的廣告費后，銷售量將有一個增長，這由銷售增長因子來表示。例如，投資3萬元的廣告費，銷售增長因子為1.85

5、，即銷售量將是預(yù)期銷售量的1.85倍。根據(jù)經(jīng)驗，廣告費與銷售增長因子的關(guān)系有表3?！?售價2.002.503.003.504.004.505.005.506.00預(yù)期銷售量（千桶）413834322928252220表2 售價與預(yù)期銷售量廣告費（元）0100002021 03000040000500006000070000銷售增長因子1.001.401.701.851.952.001.951.80表3 廣告費與銷售增長因子 問李經(jīng)理如何確定彩漆的售價和廣告費，才能使公司獲得的利潤最大？ 要求1、學(xué)習(xí)態(tài)度要認(rèn)真，要積極參與課程設(shè)計，鍛煉獨立思考能力；2、嚴(yán)格遵守上機時間安排；3、按照MATLAB

6、編程訓(xùn)練的任務(wù)要求來編寫程序；4、根據(jù)任務(wù)來完成數(shù)學(xué)建模論文；5、報告書寫格式要求按照沈陽航空航天大學(xué)“課程設(shè)計報告撰寫規(guī)范”；7、報告上交時間：課程設(shè)計結(jié)時上交報告。8、嚴(yán)謹(jǐn)抄襲行為。指導(dǎo)教師 年 月 日負(fù)責(zé)教師 年 月 日學(xué)生簽字 年 月 日沈陽航空航天大學(xué)課 程 設(shè) 計 成 績 評 定 單課 程 名 稱 數(shù)學(xué)建模實踐 院（系） 理學(xué)院 專業(yè) 信息與計算科學(xué) 課程設(shè)計題目1 對函數(shù)進行麥克勞林展開及誤差分析 2 無變位油罐中油量確定及誤差分析 3 評卷成績調(diào)整程序設(shè)計 4 廣告費用與銷售價格的調(diào)整程序設(shè)計 學(xué)號 2021041401002 姓名 郭 婧 指導(dǎo)教師評語：課程設(shè)計成績 指導(dǎo)教師

7、簽字 年 月 日目 錄目 錄V摘 要1正文31 題目一31.1 問題重述31.2 問題求解31.3 題目結(jié)果42 題目二52.1 問題重述52.2 問題求解62.3 題目結(jié)果103 題目三163.1 問題重述163.2 問題求解173.3 題目結(jié)果204 題目四214.1 問題重述214.2 問題求解224.3 題目結(jié)果24參考文獻26源程序27摘 要在本次課程設(shè)計中，我的課程設(shè)計題目是四道題。第一道題目里的第一個問號是用Matlab編寫函數(shù)，根據(jù)人為設(shè)定的n，函數(shù)可以任意展開，并且在Matlab運行界面顯示的是展開的多項式。第二個問號里要求在任意設(shè)定的n階下，帶入自變量的值，然后求出的麥克勞

8、林展式的函數(shù)值。第三個問號里要求我們通過畫圖對的本來的式子、麥克勞林展開的式子、以及作比較。在題中已經(jīng)給定畫圖區(qū)間，在這個區(qū)間內(nèi)畫出圖形，進行比較。第二道題目中，給出了一個Excel表格，里面有無變位進油量表和無變量出油量表。我們知道在一個油罐中罐中油的高度和體積是有一定的關(guān)系的，題中就把這種關(guān)系式給了我們，式中的一些參量已經(jīng)給出，編寫這個式子的程序即可。帶入不同的高度可以輸出不同的體積。接下來根據(jù)“無變位實驗采集數(shù)據(jù)表”中的無變位進油表中的數(shù)據(jù)，計算公式V(h)與實驗數(shù)據(jù)之間的誤差WC(h)。我們可以先把體積數(shù)據(jù)保存在Matlab中，然后用表中已經(jīng)給了的高度帶入V(h)，這時可以求出一系列的

9、體積，然后與真實值進行作差，得到的數(shù)據(jù)即為誤差。再次用多項式擬合確定誤差函數(shù)WC(h)表達式。最后，用誤差函數(shù)WC(h)調(diào)整V(h)，并用“無變位實驗采集數(shù)據(jù)表”中的無變位出油表中的數(shù)據(jù)檢驗調(diào)整結(jié)果。第三道題目是評卷成績調(diào)整程序設(shè)計，題目中給出了一些學(xué)生的由不同專家給出的閱卷成績。要求先求出每一個學(xué)生由不同專家給出的成績的平均值，然后求出標(biāo)準(zhǔn)差。再求出第j個專家給出成績的平均值，然后求出標(biāo)準(zhǔn)差。調(diào)整后的成績?yōu)椤Ｐ纬烧{(diào)整后的成績矩陣，則的平均值就是第i個學(xué)生的最后綜合成績。而是綜合成績向量，依此確定學(xué)生獲獎情況。由此最后求出的列向量確定1名特等獎，1名一等獎，2名二等獎；3名三等獎。第四題是一個

10、廣告模型，某裝飾材料公司欲以每桶2元的價錢購進一批彩漆以供日后銷售。為了可以很快的收益并且收回大量的資金，公司經(jīng)理李先生打算做廣告，于是便找到廣告公司的王經(jīng)理進行咨詢。李經(jīng)理認(rèn)為，隨彩漆售價的提高，預(yù)期銷售量將減少，并對此進行了估算，見表格2，隨著銷售價格的增加，銷售量下降。他問王經(jīng)理廣告有多大的效應(yīng)。王經(jīng)理說：“投入一定的廣告費后，銷售量將有一個增長，這由銷售增長因子來表示。通過表格3可知，隨著廣告費用的升高，銷售因子先上升后下降。問李經(jīng)理如何確定彩漆的售價和廣告費，才能使公司獲得的利潤最大。關(guān)鍵詞：擬合函數(shù)；誤差分析；調(diào)整矩陣；利益最大化正文1 題目一1.1 問題重述第一個問號是用Matl

11、ab編寫函數(shù)，根據(jù)人為設(shè)定的n，函數(shù)可以任意展開，并且在Matlab運行界面顯示的是展開的多項式。第二個問號里要求在任意設(shè)定的n階下，帶入自變量的值，然后求出的麥克勞林展式的函數(shù)值。第三個問號里要求我們通過畫圖對的本來的式子、麥克勞林展開的式子、以及作比較。在題中已經(jīng)給定畫圖區(qū)間，在這個區(qū)間內(nèi)畫出圖形，進行比較。1.2 問題求解（1）根據(jù)數(shù)學(xué)分析課程中學(xué)到的麥克勞林展開的定義，可知= （） （1） （） （2） （） （3）有了公式（1）（2）（3）就可以對編程進行麥克勞林展開，把、用Matlab語言進行編輯，然后作差即得得麥克勞林展式。（2）再上一個問號中已經(jīng)把麥克勞林展式求出來，在第二個程

12、序中只需把任意自變量值代入求函數(shù)值。把Matlab中M文件的函數(shù)名由function Tx=myfun1(n)改為function Tx=myfun2(x,n)，輸入任意的x和n就可以求出任意階展式的任意函數(shù)值。（3）第三個問號是畫圖比較，x的區(qū)間已經(jīng)給出，用plot命令可以直接畫出圖形調(diào)用格式為a(k)=log(1-x(k)/(1+x(k); plot(x,a,'*')。畫用麥克勞林展開的式子調(diào)用格式為Tn(k)=myfun2(x(k),n); plot(x,Tn,'*')。畫圖形的調(diào)用格式為y(k)=a(k)-Tn(k); plot(x,y,'*&#

13、39;)。為了更直觀的觀察圖形之間的關(guān)系和差距，最后把三個圖形畫到同一個圖形中，可以用subplot(m,n,p)命令把四個畫到同一圖中，分為四個小子圖，m，n為畫幾乘幾的子圖，p為第幾個圖。1.3 題目結(jié)果（1）第一個問號的運行結(jié)果，n=10時的的麥克勞林展式。圖1.3.1（2）第一個問號運行結(jié)果，x=10，n=2時的麥克勞林展式的函數(shù)值。圖1.3.2（3）第三個問號運行結(jié)果，取n=1時的、麥克勞林展開式、以及同時在一個圖形時的圖形。圖1.3.32 題目二2.1 問題重述第二道題目中，給出了一個Excel表格，里面有無變位進油量表和無變量出油量表。我們知道在一個油罐中罐中油的高度和體積是有一

14、定的關(guān)系的，題中就把這種V與h的大致關(guān)系式給了我們：，式中的一些參量已經(jīng)給出，a=17.8/2、b=12/2、L1=0.4、L2=2.05，編寫這個式子的程序即可。帶入不同的高度可以輸出不同的體積。接下來根據(jù)“無變位實驗采集數(shù)據(jù)表”中的無變位進油表中的數(shù)據(jù)，計算公式V(h)與實驗數(shù)據(jù)之間的誤差WC(h)。我們可以先把體積數(shù)據(jù)保存在Matlab中，然后用表中已經(jīng)給了的高度帶入V(h)，這時可以求出一系列的體積，然后與真實值進行作差，得到的數(shù)據(jù)即為誤差。再次用多項式擬合確定誤差函數(shù)WC(h)表達式。最后，用誤差函數(shù)WC(h)調(diào)整V(h)，并用“無變位實驗采集數(shù)據(jù)表”中的無變位出油表中的數(shù)據(jù)檢驗調(diào)整

15、結(jié)果。2.2 問題求解（1）第一個問號中要求編輯計算體積的公式的函數(shù)M文件。根據(jù)題目中給出的a、b、L1、L2的值以及公式V(h)進行編輯。程序為：function Vh=myfun4(h)a=17.8/2;b=12/2;L1=0.4;L2=2.05;Vh=a*b*(L1+L2)*asin(h-b)/b)+(h-b)/b)*sqrt(1-(h-b)2/b2)+pi/2*10;（2）要求根據(jù)“無變位實驗采集數(shù)據(jù)表”中的無變位進油表中的數(shù)據(jù)計算公式V(h)與實驗數(shù)據(jù)之間的誤差WC(h)，并用多項式進行擬合。誤差值=|真實值-公式求解的函數(shù)值| （4）所以我用以下M文件進行求解誤差，此M文件可以求解

16、每一個進油高度所對應(yīng)的誤差。function myfun5(X1,X2)for l=1:length(X1) a=17.8/2; b=12/2; L1=0.4; L2=2.05;Vh(l)=a*b*(L1+L2)*asin(X1(l)-b)/b)+(X1(l)-b)/b)*sqrt(1-(X1(l)-b)2/b2)+pi/2*10; a(l)=abs(X2(l)+262-Vh(l); disp(a(l)end根據(jù)高度和誤差進行曲線擬合，擬合命令為：x=X1;y=Y;plot(x,y,*)polyfit(x,y,n)hold onfplot(fx,0,12,r)我分別對曲線進行了二次，三次，四次

17、，五次擬合，得到以下擬合曲線：二次擬合曲線：圖2.2.1三次擬合曲線：圖2.2.2四次擬合曲線：圖2.2.3五次擬合曲線：圖2.2.4由以上曲線擬合可知：進行三次和四次多項式擬合的曲線較好，更貼合原圖。（3）用誤差WC(h)調(diào)整V(h)，并用“無變位實驗采集數(shù)據(jù)表”中的無變位出油表中的數(shù)據(jù)檢驗調(diào)整結(jié)果。把編輯的M文件的求誤差的語句中的絕對值去掉后，誤差都為負(fù)值，所以應(yīng)在V(h)后減去WC(h)，所以，調(diào)整后的V(h)= V(h)-WC(h)。分別用三次多項和四次多項式進行求解誤差。2.3 題目結(jié)果（1）h取10分米時的結(jié)果：圖2.3.1（2）求解的誤差值：圖2.3.2圖2.3.3三次擬合出來的

18、曲線為：WC(h)= -0.084* +1.5065* +5.8216* -1.7108四次擬合出來的曲線為：WC(h)= -0.0025*-0.0167* +0.8876*+8.0826*-4.3828（3）用三次多項式擬合出的WC(h)調(diào)整的結(jié)果：圖2.3.4用三次多項式擬合出的WC(h)圖形：圖2.3.5用四次多項式擬合出的WC(h)調(diào)整的結(jié)果：圖2.3.6用四次多項式擬合出的WC(h)圖形：圖2.3.73 題目三3.1 問題重述第三道題目是評卷成績調(diào)整程序設(shè)計，題目中給出了一些學(xué)生的由不同專家給出的閱卷成績。要求先求出每一個學(xué)生由不同專家給出的成績的平均值，然后求出標(biāo)準(zhǔn)差。再求出第j個

19、專家給出成績的平均值，然后求出標(biāo)準(zhǔn)差。調(diào)整后的成績?yōu)?。形成調(diào)整后的成績矩陣，則的平均值就是第i個學(xué)生的最后綜合成績。而是綜合成績向量，依此確定學(xué)生獲獎情況。由此確定1個特等獎，1個一等獎，2個二等獎；3個三等獎。3.2 問題求解這道題目就是對矩陣進行變換，根據(jù)矩陣求出要求解出的數(shù)值。如每個學(xué)生由不同的專家評卷得出的成績的平均值，每個學(xué)生成績的標(biāo)準(zhǔn)差。每個專家對不同學(xué)生評卷成績的平均值，每個專家評卷成績的標(biāo)準(zhǔn)差。求出這些必要的數(shù)值后，根據(jù)題中所給的公式求出Y矩陣，把矩陣的每一行在作和求平均值，最后得到一個列向量，列向量的每一行就為最后評定學(xué)生成績的標(biāo)準(zhǔn)，再對這個列向量的每一行的數(shù)值進行有大到小的

20、排列，由此最后求出的列向量就可以得出得特等獎、一等獎、二等獎、三等獎的人。標(biāo)準(zhǔn)差的求解公式為： （5）平均值的求解公式為： （6）求解每個學(xué)生成績的平均值：B=A(:,1);for i=2:4 B=B+A(:,i);endB=B/4;求解每個學(xué)生成績的標(biāo)準(zhǔn)差：for j=1:4 D(:,j)=(A(:,j)-B).2;endF=0;for l=1:4 F=F+D(:,l);endF=sqrt(F/4);求解每位專家評卷成績的平均值：C=A(1,:);for j=2:17 C=C+A(j,:);endC=C/17;求解每位專家評卷成績的標(biāo)準(zhǔn)差：for k=1:17 E(k,:)=(A(k,:)-

21、C).2./17;endG=0;for q=1:17 G=G+E(q,:);endG=sqrt(G/17);求解最后的公式：for p=1:4 H(:,p)=(A(:,p)-C(p)./G(p).*F+B;endI=0;for n=1:4 I=I+H(:,n);endI=I./4;對最后求出的列向量進行排列并輸出學(xué)生號：a=size(I);a=a(1);for h=1:a J(h,1)=h;endfor w=1:16 for s=w:17 if I(w)<I(s) t=I(w); I(w)=I(s); I(s)=t; z=J(w); J(w)=J(s); J(s)=z; end ende

22、nd3.3 題目結(jié)果（1）每個學(xué)生由不同的專家評卷得出的成績的平均值，每個學(xué)生成績的標(biāo)準(zhǔn)差。每個專家對不同學(xué)生評卷成績的平均值，每個專家評卷成績的標(biāo)準(zhǔn)差。學(xué)生1234567889.2588.7587.586.5083.7583.7581.2580.009101112131415161777.5067.5067.5067.5060.0056.0056.2555.0052.50表3.3.1學(xué)生123456786.05705.44862.50004.14589.60146.49527.39516.12379101112131415161710.307812.99045.59025.59027.701

23、16.49526.49525.00004.3301表3.3.2專家123479.411873.823568.823569.8235專家12343.08143.12733.41813.9585表3.3.3（2）最后求出經(jīng)過排序的綜合成績向量。學(xué)6332115.6118114.0837104.7201102.162598.890397.924992.71039111210131617141591.642358.569958.479847.861233.286728.456626.051224.403924.4039表3.3.4所以得到特等獎的學(xué)生為學(xué)生1，得到一等獎的學(xué)生

24、為學(xué)生5，得到二等獎的學(xué)生為學(xué)生2和學(xué)生6，得到三等獎的學(xué)生為學(xué)生4，學(xué)生7，學(xué)生3。4 題目四4.1 問題重述第四題是一個廣告模型，某裝飾材料公司欲以每桶2元的價錢購進一批彩漆以供日后銷售。為了可以很快的收益并且收回大量的資金，公司經(jīng)理李先生打算做廣告，于是便找到廣告公司的王經(jīng)理進行咨詢。李經(jīng)理認(rèn)為，隨彩漆售價的提高，預(yù)期銷售量將隨之減少，并對此進行了估算，銷售價格與銷售量之間的關(guān)系見表格2：售價2.002.503.003.504.004.505.005.506.00預(yù)期銷售量（千桶）413834322928252220表2 售價與預(yù)期銷售量通過表格2知道，銷售量隨著銷售價格的增加而減少。他

25、問王經(jīng)理廣告有多大的效應(yīng)。王經(jīng)理說：“投入一定的廣告費后，銷售量將有一個增長，這有銷售因子決定，如投入的廣告費用為2021 0元時，銷售量將為原銷售量的1.7倍。廣告費用與銷售增長因子的關(guān)系見表格3：廣告費（元）0100002021 03000040000500006000070000銷售增長因子1.001.401.701.851.952.001.951.80表3 廣告費與銷售增長因子通過表格3可知，隨著廣告費用的升高，銷售因子先上升后下降。問李經(jīng)理如何確定彩漆的售價和廣告費，才能使公司獲得的利潤最大。4.2 問題求解根據(jù)題中所給要求可知，此題為一個應(yīng)用線性規(guī)劃求解的題目，所以應(yīng)使用Lingo

26、軟件進行解題。題中問如何確定彩漆的售價和使用的廣告費用，才能使公司獲得的利益最大，根據(jù)這個要求可以寫出目標(biāo)函數(shù)。設(shè)x為廣告費用，則銷售增長因子為，y為彩漆售價，則為預(yù)期銷售量。根據(jù)表格中所給的數(shù)據(jù)可以擬合出和的曲線。擬合出來的=-0.0426*x2+0.4092*x+1.0288（這里的x以萬元為單位），=-5.1333*x+50.4222（這里的以千桶為單位）。銷售利潤=銷售量*銷售因子*銷售價格-廣告費用 (7)目標(biāo)函數(shù)為：max=(-5.1333*y+50.4222)*1000*(-0.0426*x2+0.4096*x+0.0188)*y-x*10000;因此只要求出這個目標(biāo)函數(shù)的最大值

27、就是利益最大組合。下面是擬合的和的圖像，根據(jù)圖像可知，擬合出的函數(shù)與實際數(shù)值有很好的相近關(guān)系，所以可以用這兩個函數(shù)擬合表格中的數(shù)據(jù)。圖4.2.1圖4.2.2應(yīng)用Lingo求解：圖4.2.3在此模型求解中，我假設(shè)廣告費用x可以在0元到70000元中任意取值，銷售價格可以在2元到6元中任意取值，在由線性規(guī)劃求解后的值向整數(shù)靠近，此時求得的目標(biāo)函數(shù)的值為線性規(guī)劃中與現(xiàn)實情況比較接近的解。4.3 題目結(jié)果由Lingo軟件解得的結(jié)果為：圖4.3.1由線性規(guī)劃后的結(jié)果可知：在廣告費用x=3.859588萬元，銷售價格y=4.911285元時，此時的目標(biāo)函數(shù)值最大：max=(-5.1333*y+50.422

28、2)*1000*(-0.0426*x2+0.4096*x+0.0188)*y-x*10000=80901.71元所以我們?nèi)V告費用x=4萬元，銷售價格y=5元，則此時公司獲得的利潤最大。最大的利潤為80758元。參考文獻1王正東,數(shù)學(xué)軟件與數(shù)學(xué)實驗.北京:科學(xué)出版社,2021 2劉玉璉,傅沛仁等,數(shù)學(xué)分析講義.北京.高等教育出版社,20213吳建國,數(shù)學(xué)建模案例精編.北京:中國水利水電出版社,2021 源程序1 題目一（1）function Tx=myfun1(n)syms xfx=1;for k=1:n fx=fx+(-1)(2*k-1)*xk/k;%x存在的區(qū)間為endgx=1;for l

29、=1:n gx=gx+(-1)(l-1)*xl/l; %x存在的區(qū)間為endTx=fx-gx;(2) function Tx=myfun2(x,n)fx=1;for k=1:n fx=fx+(-1)(2*k-1)*xk/k; %x存在的區(qū)間為endgx=1;for l=1:n gx=gx+(-1)(l-1)*xl/l; %x存在的區(qū)間為endTx=fx-gx;(3) function myfun3(n)x=linspace(-2/3,2/3,50);for k=1:length(x) Tn(k)=myfun2(x(k),n); a(k)=log(1-x(k)/(1+x(k); y(k)=a(k

30、)-Tn(k);endsubplot(2,2,1);plot(x,a,'*')title('ln(1-x)/(1+x)')subplot(2,2,2);plot(x,Tn,'*') title('麥克勞林展式')subplot(2,2,3);plot(x,y,'*')title('作差圖')subplot(2,2,4);plot(x,Tn,'r',x,a,'b',x,y,'g')2 題目二（1）function Vh=myfun4(h)a=17.8/2

31、;b=12/2;L1=0.4;L2=2.05;Vh=a*b*(L1+L2)*asin(h-b)/b)+(h-b)/b)*sqrt(1-(h-b)2/b2)+pi/2*10;（2）function myfun5(X1,X2)for l=1:length(X1) a=17.8/2; b=12/2; L1=0.4; L2=2.05; Vh(l)=a*b*(L1+L2)*asin(X1(l)-b)/b)+(X1(l)-b)/b)*sqrt(1-(X1(l)-b)2/b2)+pi/2*10; a(l)=abs(X2(l)+262-Vh(l); disp(a(l)end數(shù)據(jù)1：X1=159.02176.1

32、4192.59208.50223.93238.97253.66268.04282.16296.03309.69323.15336.44349.57362.56375.42388.16400.79413.32425.76438.12450.40462.62474.78486.89498.95510.97522.95534.90546.82558.72570.61582.48594.35606.22618.09629.96641.85653.75665.67677.63678.54690.53690.82702.85714.91727.03739.19751.42763.70764.16776.5

33、3788.99801.54814.19826.95839.83852.84866.00879.32892.82892.84906.53920.45934.61949.05963.80978.91994.431010.431026.991044.251062.371081.591102.331125.321152.361193.49./100數(shù)據(jù)2：X2=50100150200250300350400450500550600650700750800850900950100010501100115012001250130013501400145015001550160016501700175018

34、001850190019502021 20502053.832103.832105.062155.062205.062255.062305.062355.062404.982406.832456.832506.832556.832606.832656.832706.832756.832806.832856.832906.832906.912956.913006.913056.913106.913156.913206.913256.913306.913356.913406.913456.913506.913556.913606.913656.913706.91數(shù)據(jù)3：Y=10.8826 12.6

35、330 14.3649 16.1318 17.8519 19.6058 21.3520 23.0810 24.8468 26.5777 28.3288 30.0561 31.8023 33.5423 35.2971 37.0548 38.8081 40.5541 42.2939 44.0317 45.7749 47.4914 49.2354 50.9789 52.7376 54.4860 56.2425 57.9836 59.7297 61.4589 63.1929 64.9543 66.6790 68.4335 70.1978 71.9514 73.6741 75.4326 77.1627

36、78.8869 80.6703 80.7751 82.5410 82.5612 84.3353 86.0549 87.8238 89.5311 91.2775 93.0301 93.1098 94.8345 96.5893 98.3393 100.0883 101.8369 103.5815 105.3144 107.0620 108.8044 110.5536 110.5498 112.2694 114.0274 115.7643 117.5153 119.2572 121.0106 122.7597 124.5085 126.2281 127.9808 129.7393 131.4753

37、133.2285 134.9532 136.7038 138.4521（3）function myfun7(X3,X4)for k=1:length(X3) a=17.8/2; b=12/2; L1=0.4; L2=2.05; Vh(k)=a*b*(L1+L2)*asin(X3(k)-b)/b)+(X3(k)-b)/b)*sqrt(1-(X3(k)-b)2/b2)+pi/2*10; Wc(k)= -0.084*X3(k)3+1.5065*X3(k)2+5.8216*X3(k)-1.7108; b(k)=abs(X4(k)-Vh(k)+Wc(k)+262); disp(b(k);endfunct

38、ion myfun9(X3,X4)for m=1:length(X3) a=17.8/2; b=12/2; L1=0.4; L2=2.05; Vh(m)=a*b*(L1+L2)*asin(X3(m)-b)/b)+(X3(m)-b)/b)*sqrt(1-(X3(m)-b)2/b2)+pi/2*10; Wc(m)= -0.0025*X3(m)4-0.0167*X3(m)3+0.8876*X3(m)2+8.0826*X3(m)-4.3828; c(m)=abs(X4(m)-Vh(m)+Wc(m)+262); disp(c(m);end數(shù)據(jù)1：X3=1150.721123.991101.151080.

39、511061.361043.291026.081009.54993.57978.08962.99948.26933.84919.69905.78892.10878.61865.30852.15839.14826.27813.52800.87788.33775.88763.51751.21738.98726.81714.70702.64690.61678.63666.68654.75642.84630.96619.08607.21595.35583.48571.61559.72547.82535.90523.95511.97499.96487.90475.80463.65451.43439.15

40、426.80414.36401.84389.22376.49363.64350.67337.55324.27310.82297.18283.33269.24254.88240.21225.21209.81193.94177.54160.48142.62./100數(shù)據(jù)2：X4=3706.91-52.72102.72152.72202.72252.72302.72352.72402.72452.72502.72552.72602.72652.72702.72752.72802.72852.72902.72952.721002.721052.721102.721152.721202.721252.7

41、21302.721352.721402.721452.721502.721552.721602.721652.721702.721752.721802.721852.721902.721952.722021 .722052.722102.722152.722202.722252.722302.722352.722402.722452.722502.722552.722602.722652.722702.722752.722802.722852.722902.722952.723002.723052.723102.723152.723202.723252.723302.723352.723402

42、.723452.723502.723552.723602.723652.723702.723 題目三function I,J=myfun8(A)B=A(:,1);for i=2:4 B=B+A(:,i);endB=B/4;for j=1:4 D(:,j)=(A(:,j)-B).2;endF=0;for l=1:4 F=F+D(:,l);endF=sqrt(F/4);C=A(1,:);for j=2:17 C=C+A(j,:);endC=C/17;for k=1:17 E(k,:)=(A(k,:)-C).2./17;endG=0;for q=1:17 G=G+E(q,:);endG=sqrt(G/17);for p=1:4 H(:,p)=(A(:,p)-C(p)./G(p).*F+B;endI=0;for n=1:4 I=I+H(:,n);endI=I./4;a=size(I);a=a(1);for h=1:a J(h,1)=h