數學建模實驗報告試驗名稱：生日問題問題背景描述：在100個人的團體中，如果不考慮年齡的差異，研究是否有兩個以上的人生日相同。假設每人的生日在一年365天中的任意一天是等可能的，那么隨機找n個人（不超過365人）。求這n個人生日各不相同的概率是多少？從而求這n個人中至少有兩個人生日相同這一隨機事件發生的概率是多少？實驗目的：用計算機求解概率計算問題；當冪方次數較大時用多項式擬合方法確定求概率的近似計算公式；了解隨機現象的計算機模擬技術。實驗原理與數學模型：這是一個古典概率問題，n個人中每一人的生日都可能在365天中任何一天，樣本空間中樣本點總數為365，考慮n個人的生日兩兩不同，第一個人的生日可能在365天中任一天，第二個人的生日不能與第一個人生日相同，第二個人生日可能在364天中任何一天，類推可得，n個人生日兩兩不同的這一事件的總共有365*364*(365-n+1). 故這n個人的生日各不相同的概率（可能性）以下面公式計算：P （1）因而，n個人中至少有兩人生日相同這一隨機事件發生的概率為：P(n)= （2）但是在利用公式進行計算時,所用的乘法次數和除法次數較多,可以考慮用多項式做近似計算。這需要解決多項式擬合問題。主要內容（要點）：1、 求出n個人中至少有兩個人生日相同的概率P（n）的近似公式；2、 根據P（n）的近似公式，用計算機分別計算出當團體人數取n=1，2，100時的概率值：P（1），P（2），P（100）。在Matlab環境下用指令plot(p)繪制圖形，描述概率值隨團體人數變化的規律；3、 特殊概率值的計算。在有40個學生的班上，至少有2個同學生日相同的概率是多少？60個人的團體中，至少有兩個人生日在同一天的概率又是多少？在80個人的團體中，情況又如何？4、 用5次多項式擬合方法尋找一個近似計算概率的公式；5、 考慮團體總人數對概率值的影響；計算機仿真（數值模擬）。實驗過程記錄（含：基本步驟、主要程序清單及異常情況記錄等）：1、 利用（2），用計算機分別計算出當團體人數取n=1,2,100時的概率值：P（1），P（2），P（100），并繪制圖形。Matlab程序具體如下：for k=1:100p(k)=1-prod(365-k+1:365)/365k;endplot(p)并以shengriyi.m為文件名保存，然后在Matlab工作環境下輸入如下指令： shengriyi結果所得圖形如下：2、 特殊概率值的計算。由于前面已經計算了概率值P(k)(k=1,2,100),所以只需鍵入P（40），P（60），P（80）即可。如輸入如下指令： p(40)ans =0.8912一個40個同學的學生班上，至少有兩個同學生日相同的概率是P（40）=0.8912；同理可求出60個人的團體中，至少有兩個人生日相同的概率是P（60）= 0.9941；在80個人的團體中，至少有兩人生日相同的概率是P（80）=0. 0.9999。3、 參考上圖，用五次多項式擬合方法尋找近似計算概率的公式。在Matlab環境下鍵入下列指令（該指令為求五次多項式擬合的多項式系數）： n=1:100; c5=polyfit(n,p,5)c5 =-0.0000 0.0000 -0.0001 0.0023 -0.0046 -0.0020該多項式即為：在Matlab環境下繼續鍵入下列指令： p5=polyval(c5,n); /用多項式近似計算100個概率值 plot(n,p,n,p5,.) /畫出擬合多項式的圖象與概率曲線作比較結果所得的圖象如下所示：用五次多項式作近似計算P（30）、P（50）和P（70），指令和結果如下： p5(40)ans =0.8895 p5(60)ans = 0.9985 p5(80)ans =0.99434、在某團體中，要保證“至少有兩人生日相同”的概率大于99%，可以利用第一個步驟以算出的100個概率值，鍵入如下指令： find(p0.99)，可得結果為：ans = Columns 1 through 27 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 Columns 28 through 44 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100從結果可看出，該團體總人數若超過57人，則這個團體中至少有兩人生日相同的概率將大于99%。5、 算機仿真（數值模擬）。隨機產生30個正整數，介于1到365之間（用這30個數代表一個學生班的30個同學的生日），然后統計數據，觀察是否有兩人以上的人生日相同。當30人中有兩人生日相同時，計算機輸出為“1”，否則輸出為“0”。如此重復觀察100次，可得頻率。下面是做計算機模擬的Matlab源程序： n=0;for m=1:100 y=0; x=1+fix(365*rand(1,30); for i=1:29 for j=i+1:30 if (x(i)=x(j),y=1; break,end end end n=n+y;endf=n/mshengrier.m為文件名保存在Matlab工作空間中，并在Matlab環境下鍵入j，回車，可輸出結果：=0.6500實驗結果報告與實驗總結：通過本試驗的學習，對一般較簡單的Matlab語句有了更深得了解,對一些循環語句也有了一定的認識,但對于語句與語句之間在循環判斷條件下如何進行連接,以及如何寫出正確的語句還存在著一定的困難。然而從這個實驗中也有了不少的收獲，在Matlab環境下計算概率值，但當冪方很大的時候，就較難用乘冪直接求出，其已超出計算機的最大數，最終只能作近似計