1、數學實驗2013年2月實驗1 matlab基本特性與基本運算【實驗目的】了解Matlab基本特性與基本運算【實驗要求】1、熟悉MATLAB語言編程環境；2、熟悉MATLAB語言命令；3、熟悉Matlab基本運算命令【實驗原理】MATLAB 是以復雜矩陣作為基本編程單元的一種程序設計語言。它提供了各種矩陣的運算與操作，并有較強的繪圖功能。1.1 基本規則(1) 一般MATLAB命令格式為 輸出參數1，輸出參數2，=（命令名）（輸入參數1，輸入參數2，） 輸出參數用方括號，輸入參數用圓括號如果輸出參數只有一個可不使用 括號。 (2) %后面的任意內容都將被忽略，而不作為命令執行，一般用于為代碼加注

2、釋。(3) 可用、鍵來重現已輸入的數據或命令。用、鍵來移動光標進行修改。(4) 所有MATLAB命令都用小寫字母。大寫字母和小寫字母分別表示不同的變量。(5) 常用預定義變量，如pi 、Inf、NaN、ans(6) 矩陣的輸入要一行一行的進行，每行各元素用空格或“，”分開，每行用“；”分開。如 MATLAB書寫格式為A=1 2 3 ；4 5 6 ；7 8 9 在MATLABZ中運行如下程序可得到A矩陣 a=1 2 3;4 5 6;7 8 9 a = 1 2 3 4 5 6 7 8 9(7) 需要顯示命令的計算結果時，則語句后面不加“；”號，否則要加“；”號。 運行下面兩種格式可以看出它們的區別

3、：a=1 2 3;4 5 6;7 8 9 a=1 2 3;4 5 6;7 8 9;a= 1 2 3 不顯示結果 4 5 6 7 8 9(8) 當輸入語句過長需要換行時，應加上“”后再回車，則可續行輸入。1.2 文件管理常用命令(1) 幫助（HELP）命令 MATLAB有很多命令，因此很不容易記憶。使用HELP命令可以得到有關命令的屏幕幫助信息。如 在MATLAB環境下直接運行HELP命令就會在屏幕上給出MATLAB的基本命令以及相關的工具箱軟件命令（請自行驗證）。 而下面的格式可以得到指定命令的幫助信息，和詳細說明。 help 命令名 (2) save 命令該命令存儲定義的變量或演算結果，也可

4、以用來存儲指定的變量。 命令格式為 save 文件名 變量名表 -append -ascii (3) what 命令 該命令可以在當目錄下顯示MATLAB文件和MAT數據文件(4) dir 命令 顯示當前目錄下的所有文件.(5) clear命令 如果輸入 calear a b c ，則表示清除工作空間中指定變量a,b,c；如果僅僅輸入calear命令，則清除整個工作空間。 1.3 矩陣運算 MATLAB具有強大的矩陣運算功能,這里我們只作簡單介紹.(1) 在MATLAB中表示一個矢量要用方括號， 而列矢量的輸入只需在行矢量輸入格式基礎上加轉置符“”即可。 如 x=1 2 3;4 5 6 x =

5、 1 2 3 4 5 6 而 x=1 2 3;4 5 6' (加轉置符) x = 1 4 2 5 3 6(2) 下面三條命令可以產生一個行矢量1 a=linspace(x,y,n) 2 a=logspace(x,y,n) 3 a=x:n:y第一條命令可以在線性空間產生一個值在x至y之間間隔點數為n的行矢量(一組數據)。其行矢量的起始值是x, 終值為y, 點數為n。第二條命令可以在對數空間產生一個值在10x至 10y之間等間隔的行矢量(一組數據)。第三條命令產生x至y步長為n的行矢量。但是, 三個命令之間存在差別，下面的例子可以說明這一點。 例1 x=logspace(0,5,6) x

6、= 1 10 100 1000 10000 100000 例2 x=linspace(0,10,11) x = 0 1 2 3 4 5 6 7 8 9 10 例3 x=0:1:10 x = 0 1 2 3 4 5 6 7 8 9 10 通過上面三個例子可以看出例一, 例二中n代表選取的點數。而在例三 中n則表示步長. 注意它們的區別。(3) 矩陣的加, 減, 乘, 除等, 和其它語言書寫一樣。 但要注意的是在運算 符前面加有“.”則表示是數組運算，為對應元素的相應操作。 例4 請注意下面的不同情況. A=1 2;3 4; b=2 2;2 2;c=a*b c=a.*b c = c = 6 6 2

7、 4 14 14 6 8 c=a/b c=a./b c = c = 1.0e+015 * 0.5000 1.0000 -2.2518 2.2518 1.5000 2.0000 -2.2518 2.25181.4 常見數學函數函 數 名數 學 計 算 功 能函 數 名數 學 計 算 功 能abs（x）實數的絕對值或復數的幅值floor（x）對x朝-方向取整acos（x）反余弦arcsingcd（m，n）求正整數m和n的最大公約數acosh（x）反雙曲余弦arccoshimag（x）求復數x的虛部angle（x）在四象限內求復數 x 的相角lcm（m，n）求正整數m和n的最小公倍數asin（x）反

8、正弦arcsinlog（x）自然對數（以為底數）asinh（x）反雙曲正弦arcsinhlog10（x）常用對數（以10為底數）atan（x）反正切arctanreal（x）求復數x的實部atan2（x,y）在四象限內求反正切rem（m，n）求正整數m和n的m/n之余數atanh（x）反雙曲正切arctanhround（x）對x四舍五入到最接近的整數ceil（x）對x朝+方向取整sign（x）符號函數：求出x的符號conj（x）求復數x的共軛復數sin（x）正弦sincos（x）余弦cossinh（x）反雙曲正弦sinhcosh（x）雙曲余弦coshsqrt（x）求實數x的平方根：exp（x）

9、指數函數 tan（x）正切tanfix（x）對x朝原點方向取整tanh（x）雙曲正切tanh例5：輸入 x=-4.85 -2.3 -0.2 1.3 4.56 6.75，則： ceil(x)= -4 -2 0 2 5 7 fix(x) = -4 -2 0 1 4 6 floor(x) = -5 -3 -1 1 4 6 round(x) = -5 -2 0 1 5 71.4 關系及邏輯運算 (1) 關系運算允許常量(或矩陣中的元素與元素) 之間的比較. 如果比較結果為真,則答案為1, 否則為0. 常用的關系運算符如下:< 小于 , > 大于 , = 等于 <= 小于等于, >

10、;= 大于等于, = 不等于對于復數”=” 和”=” 既比較實部也比較虛部, 而其它運算僅比較實部。例6 比較 a=5:1:15 與 b=1 2 8 8 7 10 12 11 13 14 15 中的對應元素 是否相同,則運行下面程序a=5:1:15;b=1 2 8 8 7 10 12 11 13 14 15;y=a=b y = 0 0 0 1 0 1 0 0 1 1 1 (2) MATLAB有三個邏輯運算符. & 邏輯與 | 邏輯或 邏輯非在矩陣運算中,以上命令是將兩個矩陣中的對應元素進行運算, 得到的結 果是具有同樣元素的矩陣。例7 分別對下面兩個矩陣進行邏輯運算. X=0 1;1

11、0; Y=0 0;1 0; X&Y ans = 0 0 1 0 X|Y ans = 0 1 1 0 X ans = 1 0 0 1 MATLAB在進行混合運算時, 其順序為: 算術運算 關系比較 邏輯運算例8 計算 x&y+x>y 時, 程序首先執行x+y運算, 然后將相加結果與y比較, 最后將比較結果和x作與運算，運行結果如下。 x=0 1;1 0; y=0 0;1 0; x&y+x>y ans = 0 1 1 0【實驗內容】熟悉MATLAB語言編程環境。在這一環境中，系統提供了許多編寫，調試和執行MATLAB程序的便利工具。熟悉變量、MATLAB基本數學

12、函數使用，矩陣操作練習，主要包括矩陣的建立，算數運算，邏輯運算和關系運算。（一）熟悉Matlab的運行環境以及HELP命令，運行demo命令，觀看matlab所具有的功能。（二）驗證例題中所有操作。（三）完成以下各題，將步驟命令寫入實驗報告，并在機器上運行結果驗證。1輸入A=7 1 5；2 5 6；3 1 5，B=1 1 1; 2 2 2; 3 3 3，在命令窗口中執行下列表達式，掌握其含義：A(2, 3) A(:,2) A(3,:) A(:,1:2:3) A(:,3).*B(:,2) A(:,3)*B(2,:) A*B A.*B A2 A.2 B/A B./A2輸入C=1:2:20，則C（i

13、）表示什么？其中i=1,2,3,10；3查找已創建變量的信息，刪除無用的變量；4輸入以下程序，察看結果，如果出錯給出處理方法并驗證。（1）a=magic(3);b=243; （2）a=magic(3);b=243;save myvar a b; save myvar a b ascii;clear; pause;clear;load myvar; load myvar ascii;whoa,b5求下列表達式的值。（1）（2），其中a=3.5,b=5,c=-9.8。6欲通過系統做一平面圖，請查找相關的命令與函數，獲取函數的幫助信息。實驗2 matlab繪圖【實驗目的】學會用MATLAB語言繪圖【

14、實驗要求】1、熟悉matlab繪制二維圖形2、熟悉matlab繪制三維圖形3、了解特殊圖形的繪制4、熟悉圖形的標識與修飾方法【實驗原理】1、 plot 該命令在線性坐標系下繪制y對應于x的軌跡。若其中之一為矩陣則該命令將對應于矩陣的行或者列繪制一簇曲線。該命令也可以在同一坐標軸下繪制多條曲線。例1 運行程序 t=0:1:10; x1=t; plot(x1)可在圖形窗口得到下面圖形: 例2 x=0:0.5:7; y=sin(x); y1=2*x; y2=3*x; plot(y,x,y1,x,y2,x)從上面的例子可以看出, 如果要在同一坐標軸下繪制多條曲線, 可用如下格式： plot(x1,y1

15、,x2,y2.xn,yn) % 對應于不同矢量時 其中x1,y1.xn,yn是對應的關系。2、 log log 在兩個對數（標準對數）坐標間繪制圖形。 3、 semilogx 或semilogy 在x軸或y軸為對數坐標, 而另一軸為線性軸的空間繪制圖形。4、polar 在極坐標空間繪制圖形。5、 grid 在圖形中加柵格。6、 subplot (r,c,p) 該命令將屏幕分成r*c個窗口,p表示在第幾個窗口. 同時繪制多個圖形窗口時,繪圖語句間應用逗號分開。 例3 下面的程序產生四個窗口,且分別繪制相應圖形. t=0:0.01:2*pi; y=sin(t); x=cos(t); subplot

16、(221),plot(t,y),grid,subplot(222),plot(t,x),grid,subplot(223),polar(t,y),subplot(224),polar(t,x)從上面的程序可以看出,subplot語句中的r,c,p可以不用逗號分開。 上面程序執行后得到如下圖形：7、 gtext 該命令用于鼠標定位的文字注釋. 執行完繪圖命令后再執行gtext(說明文字) 命令, 就可在屏幕上得到一個光標, 然后,用鼠標選擇說明文字的位置。例4 t=0:1:10; x1=2*t; x2=0.5*t; plot(x1,t,x2,t); gtext('x=2t');

17、gtext('x=0.5t') 8、 如果需要同時在不同窗口繪制圖形, 可用figure語句來完成。 例5 t=0:1:10; x1=2*t; x2=3*t; figure,plot(x1),figure,plot(x2) 9、有關繪圖方面的還有： ² Title 給當前圖像加標題。² xlabel 或 ylabel 給x或y軸加標注。² text 在當前圖像上指定位置加注釋。10、在同一座標系中繪制多個曲線時, 曲線的顏色可以自動確定，也可以自己選擇。其使用格式如下:plot(x1,y1,選項1,x2,y2,選項2,xn,yn,選項n)11、三維

18、圖形繪制函數plot3函數與plot函數的用法類似。還可以用mesh函數繪制三維網格圖；用surf函數繪制三維曲面圖。MATLAB提供了豐富的三維繪圖函數，可以通過help命令查看。 例6 繪制三維曲線：t=0:pi/100:20*pi;x=sin(t) y=cos(t);z=t.*sin(t).*cos(t);plot3(x,y,z) title(line in 3-D space) xlabel(X) ;ylabel(Y);zlabel(z);grid on運行結果如下：【實驗內容】（一）驗證例題中的操作。（二）完成以下各題，將答案寫入實驗報告，并在機器上運行

19、結果驗證。1. 在同一坐標系中繪制余弦曲線y=cos(t-0.25)和正弦曲線y=sin(t-0.5)， t0，2 2. 繪制0，4區間上的x=10sint曲線，并要求： （1）線形為點劃線、顏色為紅色、數據點標記為加號； （2）坐標軸控制：顯示范圍、刻度線、比例、網絡線 （3）標注控制：坐標軸名稱、標題、相應文本。 3. 繪制不同插值算法（nearest、linear、cubic、spline）對peaks函數的插值效果圖，要求使用subplot命令在一個圖形窗口顯示四幅圖。4. 繪制三維圓柱螺線 ： ，要求給相應的坐標軸和標題附加標注，螺線為蘭色虛線；并使用stem函數繪制三維火柴桿圖。5

20、.利用圖形顯示命令作出下列函數的圖形：（1），其中（2），其中(3) ，其中(4) ，其中(5) 6.制作平面振動動畫（利用函數）7.我們選用球體和圓柱體，其中實驗3 matlab程序設計【實驗目的】了解Matlab編程原理和編程方法【實驗要求】1、熟悉Matlab for-end、while-end循環結構2、數學Matkab if-end、if-else-end選擇判斷結構【實驗原理】MATLAB提供三種決策或控制流結構： for循環、while循環、if-else-end結構. 這些結構經常包含大量的MATLAB命令，故經常出現在MATLAB程序中，而不是直接加在MATLAB提示符下.

21、1、for循環：允許一組命令以固定的和預定的次數重復 for x=array commands end 在for和end語句之間的命令串commands按數組（array）中的每一列執行一次. 在每一次迭代中，x被指定為數組的下一列，即在第n次循環中，x=array(：，n) 例 對n=1,2,10,求xn= 的值for n=1:10 x(n)=sin(n*pi/10);endx2、While循環與for循環以固定次數求一組命令相反，while循環以不定的次數求一組語句的值. while expression commands end 只要在表達式(expression)里的所有元素為真，就執

22、行while和end語句之間的命令串commands. 例 設銀行年利率為11.25%。將10000元錢存入銀行，問多長時間會連本帶利翻一番？money=10000years=0while money<20000 years=years+1 money=money*(1+11.25/100)end3、If-Else-End結構（1）有一個選擇的一般形式是： if expression commands end 如果在表達式(expression)里的所有元素為真，就執行if和end語句之間的命令串commands. 先建立M文件fun1.m定義函數f（x），再在Matlab命令窗口輸入f

23、un1(2),fun1(-1)即可。function f=fun1(x)if x>1 f=x2+1endif x<=1 f=2*xend【實驗內容】對以下問題,編寫M文件:(1)用起泡法對10個數由小到大排序. 即將相鄰兩個數比較,將小的調到前頭.(2)有一個矩陣,編程求出其最大值及其所處的位置.(3)編程求(4)一球從100米高度自由落下,每次落地后反跳回原高度的一半,再落下. 求它在第10次落地時,共經過多少米?第10次反彈有多高?(5)有一函數 ,寫一程序,輸入自變量的值,輸出函數值.實驗4 matlab基礎數學實驗高等數學極限【實驗目的】1、熟悉函數極限的概念2、掌握各種求

24、極限的方法3、會用Matlab命令求極限【實驗要求】熟悉Matlab中求極限的命令limit【實驗原理】例1：（1），（2）實驗過程：（1）>> syms x a b>> limit(sin(a*x)/sin(b*x), x,0)運行結果：ans =a/b（2）>>syms x>> limit(1-cos(x)/(x*sin(x),x,0)運行結果：ans =1/2【實驗內容】1、計算以下極限（1），（2），（3），（4）2、設計實驗復利問題復利，即利滾利。不僅是一個經濟問題，而且是一個古老又現代的經濟社會問題。隨著商品經濟的發展，復利計算將日益

25、普遍，同時復利的期限將日益變短，即不僅用年息、月息，而且用旬息、日息、半日息表示利息率。現在我們已進入電子商務時代，允許儲戶隨時存款或取款，如果一個儲戶連續不斷存款和取款，結算本息的頻率趨于無窮大，每次結算后將本息全部存入銀行，這意味著銀行不斷地向儲戶支付利息，稱為連續復利問題。 若銀行一年活期年利率為0.06，那么儲戶存10萬元的人民幣，如果銀行允許儲戶在一年內可任意次結算，在不計利息稅的情況下，由于復利，顯然這比一年結算一次要多，因為多次結算增加了復利。結算越頻繁，獲利越大。連續復利會造成總結算額無限增大嗎？隨著結算次數的無限增加，一年后該儲戶是否會成為百萬富翁？實驗5matlab基礎數學

26、實驗高等數學微分【實驗目的】1.熟悉基本求導公式，掌握初等函數的求導方法.2.會求函數在給定點處的導數值.3.掌握隱函數求導的方法和步驟.4.掌握參量函數求一階導數和二階導數的方法和公式.【實驗要求】1.熟悉，Matlab中的求導命令diff2.熟悉Matlab中解方程的命令solve和求導命令diff【實驗原理】例1：求下列函數的導數（1），（2）例2.已知函數 ，求 ； 例3：設 ，求例4.已知 ，求實驗步驟：1.（1）>> syms x>> y=exp(x)*(sin(x)+cos(x);>> diff(y)運行結果：ans =exp(x)*(sin(

27、x)+cos(x)+exp(x)*(cos(x)-sin(x) 即函數的導數為 （2）>> syms x>>y=log(x3+1)/(x2+1);>>diff(y)運行結果：ans =(3*x2/(x2+1)-2*(x3+1)/(x2+1)2*x)/(x3+1)*(x2+1)即函數的導數化簡得 2、（1）>>syms x;>>f=1/x;>>f1=diff(f,x);>>ff=inline(f1);>>x=1;>>ff(1)運行結果：ans =-1>> x=-2;>&g

28、t; ff(-2)運行結果：ans = -0.25003、（1）解法一：>> syms x y;>> f=solve('x2+y2-R2=0',y);>> diff(f,x)運行結果：ans = -1/(-x2+R2)(1/2)*x 1/(-x2+R2)(1/2)*x即 或說明：對于能從方程中求出函數顯示形式的題可以采用這種做法。解法二：>> syms x y R;>> f=x2+y2-R2;>> f1=diff(f,x);>> f2=diff(f,y);>> -f1/f2運行結果：

29、ans =-x/y 即 說明：對于不能從方程中解出函數顯示形式的題要采用這種做法。4.(1)>> syms t;>> x=t2;>> y=4*t;>> f=diff(y,t);f1=diff(x,t);>> f2=f/f1運行結果：f2 =2/t 即 【實驗內容】1、計算下列導數（1）（2）（3）已知函數，求，（4）設，求（5）已知，求2、設計實驗最優價格問題某房地產公司擁有100套公寓當每套公寓的月租金為1000元時，公寓全部租出。當月租金每增加25元時，公寓就會少租出一套。1.請你為公司的月租金定價，使得公司的收益最大，并檢驗結論

30、2.若租出去的公寓每月每套平均花費20元維護費，又應該如何定價出租，才能使公司收益最大實驗6 matlab基礎數學實驗高等數學積分【實驗目的】1.掌握求函數的原函數的方法.2.熟悉基本積分公式和積分方法.3.掌握求函數定積分的方法.4.會求變上限函數的導數和帶有變上限函數的極限.【實驗要求】1.掌握Matlab中積分命令int.2.熟悉Matlab中求定積分的命令.【實驗原理】例1、 求下列函數的一個原函數（1）（2）例2、 1.求下列定積分例3、 求變上限函數的導數例4、 求下列極限實驗步驟：1.（1）>> syms x;>> f=1/x4;>> int(

31、f,x)運行結果：ans =-1/3/x3 即函數 的一個原函數為 .（2）>> syms x;>> f=exp(x)/(1+exp(x);>> int(f,x)運行結果：ans =log(1+exp(x) 即函數 的一個原函數為 .2、（1）>> syms x;>> f=sqrt(1-x2);>> int(f,x,0,1)運行結果：ans =1/4*pi 3、2.（1）>> syms t x;>> y=sin(t)/t;>> diff(int(y,t,0,x),x)運行結果：ans =

32、sin(x)/x 即 4、3.（1）>> syms x t;>> f=cos(t2);>> int(f,t,sin(x),0);>> f1=diff(int(f,t,sin(x),0),x)>> f2=f1/1>> limit(f2)運行結果：ans =-1【實驗內容】1、求下列函數的一個原函數(1)、 (2) 2、求下列定積分（1）、（2）、3、求變上限函數的導數（1）、（2）、4、求極限5、現購買一棟別墅價值300萬元，若首付50萬元，以后分期付款，每年付款數目相同。10年付清，年利率為6%，按連續復利計算，問每年應付

33、款多少？（e -0.60.5448）實驗7 matlab基礎數學實驗概率論與數理統計隨機變量及其分布【實驗目的】1會利用 Matlab 軟件計算離散型隨機變量的概率、連續型隨機變量概率密度值, 以及產生離散型隨機變量的概率分布(即分布律).2掌握常見分布的隨機數產生的有關命令.【實驗要求】1掌握常見分布的分布律和概率密度的產生命令，如binopdf，normpdf等2掌握常見分布的分布分布函數命令，如binocdf，normcdf等掌握常見分布的隨機數產生命令，如binornd，normrnd等【實驗原理】1.1 常見分布1. 事件 A 在每次試驗中發生的概率是 0.3，計算（1）在 10 次

34、試驗中 A 恰 好發生 6 次的概率；（2）生成事件A發生次數的概率分布；（3）在 10 次試驗中 A 至少發生 6 次的概率；（4）設事件A發生次數為X，且X的分布函數為F(x)，求F(6.1)；又已知F（x）=0.345，求x。1（1）>> binopdf(6,10,0.3)運行結果為：ans = 0.0368（2）>>binopdf(0:10,10,0.3)運行結果為：ans = Columns 1 through 8 0.0282 0.1211 0.2335 0.2668 0.2001 0.1029 0.0368 0.0090 Columns 9 through

35、 11 0.0014 0.0001 0.0000 （3）>>binocdf(6,10,0.3)運行結果為：ans = 0.98944（4）>> binocdf(6.1,10,0.3)運行結果為：ans = 0.9894>> binoinv(0.345,10,0.3)運行結果為：ans = 2或者：>> icdf('bino',0.345,10,0.3)運行結果為：ans = 22 設隨機變量 X服從參數是 3 的泊松分布， 求（1）概率 PX=6；（2）X的分布律前七項；（3）設X的分布函數為F(x)，求F(6.1)；又已知F（x

36、）=0.345，求x。（1）>> poisspdf(6,3)運行結果為：ans = 0.0504（2）>> poisspdf(0:6,3)運行結果為：ans = 0.0498 0.1494 0.2240 0.2240 0.1680 0.1008 0.0504（3）>> poisscdf(6.1,3)運行結果為：ans = 0.9665>> poissinv(0.345,3)運行結果為：ans = 2或者：>> icdf('poiss',0.345,3)運行結果為：ans = 23設隨機變量 X服從區間2, 6上的均勻分

37、布，求（1）X=4 時的概率密度值； （2） PX5；（3）若PXx=0.345，求x。(1)>> unifpdf(4,2,6)運行結果為：ans = 0.2500 （2）>> unifcdf(5,2,6)運行結果為：ans =0.7500（3）>> unifinv(0.345,2,6)運行結果為：ans = 3.3800或者：>> icdf('unif',0.345,2,6)運行結果為：ans = 3.38004設隨機變量 X 服從參數是 6 的指數分布， 求（1）X=0，1，2，3，4，5，6 時的概率密度值；（2）PX5；

38、（3）若PXx=0.345，求x。（1）>> exppdf(0:6,6)運行結果為：ans =0.1667 0.1411 0.1194 0.1011 0.0856 0.0724 0.0613(2) >> expcdf(5,6)運行結果為：ans = 0.5654（3）>> expinv(0.345,6)運行結果為：ans = 2.5387或者：>> icdf('exp',0.345,6)運行結果為：ans = 2.53875設隨機變量 X 服從均值是 6，標準差是2 的正態分布，求：（1）X=3，4，5，6，7，8，9 時的概率密

39、度值；（2）X=3，4，5，6，7，8，9 時的分布函數值；（3）若PXx=0.345，求x；（4）求標準正態分布的上0.05分為點。（1）>> normpdf(3:9,6,2)運行結果為：ans =0.0648 0.1210 0.1760 0.1995 0.1760 0.1210 0.0648（2）>> normcdf(3:9,6,2)運行結果為：ans = 0.0668 0.1587 0.3085 0.5000 0.6915 0.8413 0.9332（3）>> norminv(0.345,6,2)運行結果為：ans = 5.2023或者：>>

40、; icdf('norm',0.345,6,2)運行結果為：ans = 5.2023（4）>> norminv(0.95,0,1)運行結果為：ans = 1.6449或者：>> icdf('norm',0.95,0,1)運行結果為：ans = 1.64491.2 隨機數的生成1產生參數為 20, 概率為 0.25 的二項分布的隨機數。(1) 產生 1 個隨機數； (2) 產生 7個隨機數； (3) 產生 21(要求 3 行7 列)個隨機數。(1) >> binornd(20,0.25)運行結果為：ans = 4（2）>&

41、gt; binornd(20,0.25,1,7)運行結果為：ans = 4 7 3 7 7 5 7（3）>> binornd(20,0.25,3,7)運行結果為：ans = 4 9 2 0 6 5 7 5 6 3 4 4 5 2 2 5 2 0 5 6 5 或者：>> random('bino',20,0.25,3,7)運行結果為：ans = 2 9 5 1 3 4 4 3 4 2 5 3 7 1 5 4 9 6 7 6 42產生7個服從參數為6的泊松分布的隨機數。>> poissrnd(6,1,7)運行結果為：ans = 3 10 5 9

42、3 4 6或者：>> random('poiss',6,1,7)運行結果為：ans = 6 8 8 6 7 13 73產生區間(1, 3)上的連續型均勻分布的隨機數。 (1) 產生 6×6個隨機數； (2) 產生 21(要求 3 行7列)個隨機數。（1）>> unifrnd(1,3,6)運行結果為：ans = 2.0649 1.8077 1.4863 1.3516 2.2505 2.5248 2.4330 2.0971 1.3083 1.7207 2.0861 2.1521 1.3586 1.0975 2.9128 1.3776 1.8781

43、2.4953 1.6731 2.1055 2.8713 1.0024 1.5749 2.2911 1.3754 1.5496 2.6374 1.6328 2.0033 1.24641.6439 1.4830 2.4565 2.3992 2.5231 2.0088（2）>> unifrnd(1,3,3,7)運行結果為：ans = 1.6945 1.3963 2.3888 2.0646 2.8129 2.3429 1.1139 1.1843 2.3445 1.5136 1.5588 1.7854 2.6743 1.9006 1.2957 1.8630 1.0195 2.8925 1.0

44、497 2.9430 2.1649或者：>> random('unif',1,3,3,7)運行結果為：ans = 2.9571 1.9324 1.4606 2.1998 2.0276 1.9198 2.6108 2.6972 1.6513 2.1598 1.8969 1.8155 1.9018 2.4017 1.1013 2.2604 2.2063 1.0708 1.2161 2.1023 2.74454產生服從均值為4，均方差為1的正態分布的隨機數。（1）產生7個隨機數；(2) 產生 21(要求 3 行7列)個隨機數。（1）>> normrnd(4,1

45、,1,7)運行結果為：ans =3.5674 2.3344 4.1253 4.2877 2.8535 5.1909 5.1892（2）>> normrnd(4,1,3,7)運行結果為：ans = 3.9624 3.8133 6.1832 5.0668 3.1677 4.7143 4.8580 4.3273 4.7258 3.8636 4.0593 4.2944 5.6236 5.2540 4.1746 3.4117 4.1139 3.9044 2.6638 3.3082 2.4063或者：>> random('norm',4,1,3,7)運行結果為：ans = 2.5590 4.6900 5.2902 2.7975 2.3959 5.4151 4.2193 4.5711 4.8156 4.6686 3.9802 4.2573 3.1949 3.0781 3.6001 4.7119 5.1908 3.8433 2.9435 4.5287 1.8293【實驗內容】1、驗證以上例題。2、請將教材上的部分題目用計算機實現。3、下圖中每一個黑點表示釘在板上的一顆釘子，每排釘子等距排列，下一排的每個釘子恰在上一排兩相鄰釘子之間。假設有