基于時間序列arma模型的分析(完整資料)

基于時間序列arma模型的分析(完整資料）(可以直接使用，可編輯優(yōu)秀版資料，歡迎下載）

針對乳制品月產量數(shù)據的時間序列分析基于時間序列arma模型的分析(完整資料）(可以直接使用，可編輯優(yōu)秀版資料，歡迎下載）摘要：隨著經濟的發(fā)展，乳制品產業(yè)對國民健康水平的影響逐漸加大。該文從乳制品行業(yè)月產量的角度出發(fā),采用時間序列數(shù)據分析方法，對我國自1990年至2010年以來的乳制品行業(yè)月產量進行了建模分析，并在得到模型后對其進行了預測。從分析結果來看，我國的乳制品產量在200４年發(fā)生突變,特定的月份也會對其產生影響,并且在不同的時間，影響會發(fā)生變化。關鍵詞:乳制品;月份特征；產量突變；產量預測;背景：縱觀自1949年發(fā)展至今，整個行業(yè)可以分為四個發(fā)展階段:1、緩慢發(fā)展階段(1９4９~19７7）:這段時期,我國乳產業(yè)受國家經濟狀況制約發(fā)展緩慢。2、迅速擴張階段(1978～199２)：由于開始實行多種所有制進行奶牛飼養(yǎng)與奶制品加工，原奶與乳品的產量、種類、質量都有明顯的提高3、結構調整階段(1993~１９98）：1993年開始,乳品供給增長明顯快于消費增長速度，產能出現(xiàn)比較嚴重的過剩，乳粉出現(xiàn)滯積,部分乳品企業(yè)發(fā)展艱難。4、高速增長階段（1999~至今):１99８年起，乳制品產業(yè)經過產品結構大力調整，經濟效益明顯提高,隨著消費需求的迅速增長，乳制品產量也連年增長，乳產業(yè)已經從一個傳統(tǒng)產業(yè)搖身一變成為一個朝陽產業(yè)。從市場格局上看，乳制品企業(yè)可以分為4類:1、以伊利、蒙牛為代表的全國性企業(yè);2、以光明、三鹿、維維等為代表的區(qū)域性企業(yè);３、以北京三元、濟南佳寶為代表的本省省會企業(yè)；4、以雀巢為代表的外資企業(yè).在行業(yè)中,企業(yè)之間的競爭非常激烈,特別地，在近十年中市場競爭引起了市場格局的極大改變。本文將選取1993年1月起到2010年6月的月產量數(shù)據進行時間序列分析,嘗試建立該時序的時間序列模型及其詳細的建立過程，并對模型結果給出必要的經濟意義解釋。建立模型過程：1、建模過程使用eｖiｅｗs軟件，將19９0年１月到2０１0年6月總計246個月度數(shù)據輸入eｖiｅws中，Ｙt即是產量月度序列,現(xiàn)作出散點圖如下：通過觀察上圖,認為不同時間下的Yｔ的數(shù)值差異過大，并且波動程度也差別過大，故先將序列作取對數(shù)處理，作出散點圖如下：通過觀察上圖,認為該序列是一個典型的結構突變的過程,突變位置始于2003年12月。故需要加入突變虛擬變量進行檢測。從圖中我們觀察到,突變過程精確來說應屬于漸進式的突變,但在突變區(qū)間內只有２～3個月份，相對于2０0多個月份來說可以忽略,故在選擇突變類型時認為是水平突變,即從2００3年１2月開始水平突變.由于不知是否存在斜率突變，故我們一齊將水平突變變量與斜率突變變量加入檢驗。２、即建立模型:,檢驗結果如下：VariabｌeＣoefficiｅntStd.Ｅｒｒｏｒt-StatiｓｔicProb。

C０。６７43４０0．034７4219.41０100.0００0@ＴREND（1989M１2）０．０0９6800。00035727.１44７８0.00００DL2．０760610.06174433.６2３6１０.0000DL*＠TREND（20０3Ｍ１2）0．0０2085０。０01１8２1。７６3３340。0７91R-sｑuarｅｄ0。9８1355

Mｅandependentvar2。554138AｄjｕｓtedR-sｑuaｒed0.９81123

S。D．ｄｅpeｎdeｎｔvaｒ1.6314４6S．Ｅ.oｆregressｉoｎ0.２2４148

Akaｉkeinｆｏcｒitｅｒｉon—0．136897Ｓumｓｑｕａredｒesiｄ１２．15861

Ｓcｈwarzcrｉterion-0。07９89９Ｌｏglｉkelihood２０。83828

Hannan－Ｑuｉnncｒiteｒ.-０。１13９46F-ｓtaｔistｉc4245。683

Duｒbiｎ-Waｔｓonstat0.６０8５38Ｐroｂ(Ｆ-ｓtatiｓtic)0。０00000可以看到，常數(shù)項、時間趨勢項及水平突變三個變量均顯著,而斜率突變并不顯著。故我們將其舍去,接著用模型再做一次回歸：現(xiàn)在可以看出，三個變量均顯著。即現(xiàn)在認為我們正確的加入趨勢項和水平突變變量。3、考慮到ＬOGY是一個月度數(shù)據,并且從散點圖上也可以明顯觀察到存在著季節(jié)特征,故我們?yōu)榱搜芯科涫€月的特征,我們繼續(xù)加入11個虛擬變量（當月份為i時等于1，否則為０）繼續(xù)回歸，結果如下:可以看到,只有d1，ｄ2,ｄ６三個月度變量在５％的水平下是顯著的，所以我們剔除不顯著變量再回歸一次：現(xiàn)在得到的這個模型，常數(shù)項、趨勢項、水平突變變量及月度變量均顯著。4、所以接下來我們要討論的就是這個突變過程究竟是隨機趨勢過程還是趨勢平穩(wěn)過程,所以我們將采用Perroｎ檢驗進行檢驗，即我們用該模型退勢后的殘差序列進行單位根檢驗：突變點在20０３年1２月，即大概在整個序列的2/３的位置，根據Perron檢驗表查得在%５水平下臨界值應位于—３．7８0～-3．76之間,故我們明顯可以拒絕原假設即認為該殘差序列是一個平穩(wěn)過程,即LＯGY是一個趨勢平穩(wěn)序列。所以,當確定ＬOＧＹ是一個可以避免很多諸如虛假回歸等問題的趨勢平穩(wěn)序列后，我們可以放心繼續(xù)修正模型。5、現(xiàn)在，我們開始研究模型中是否存在AＲMA成分,故作出ＬＯGY退勢后的殘差序列的自相關偏相關圖如下：從圖中可以看到,殘差序列是一個平穩(wěn)的ARMＡ過程，可能存在AＲ(1）、ＡR(３）、ＭA(1)、ＭA(2)、MA（３)，故我們繼續(xù)將這五個ＡRMＡ成分加入模型中：發(fā)現(xiàn)常數(shù)項開始不顯著，AR(3)顯著性并不好，MA(3)很不顯著，但由于常數(shù)項并非關鍵成分，故我們將其剔除繼續(xù)嘗試:現(xiàn)在，AＲ(3）與MＡ（３)的顯著性在％5的水平下完全可以拒絕，故我們將其剔除，余下變量即模型：繼續(xù)回歸：可以看到，現(xiàn)在的模型中,加入的AＲＭA成分都顯著，接下來我們看看殘差序列的自相關與偏相關圖:從圖中,我們觀察到該模型的殘差序列已經接近白噪聲,即我們認為我們已經加入了正確的ＡRMA成分。所以現(xiàn)在的模型為；其中是白噪聲。從回歸結果中我們可以看出，無論是系數(shù)的顯著性,還是模型的擬合度，都令人滿意。另外經過檢驗，沒有發(fā)現(xiàn)有異方差、序列相關、多重共線性等問題的存在。模型的最終調整：經過嘗試預測，我們發(fā)現(xiàn)一個不容忽視的問題，即往往１月的產量往往被低估，而6月份的產量的季節(jié)特征似乎又消失了,故我們懷疑是否D1的系數(shù)在序列后期已發(fā)生變化，而同時D６是否開始變得不顯著。解決方案:所以我們新生成一個DN1序列,DＮ１與D1的區(qū)別就在于2003年12月前均為0,其余不變;新生成一個序列D_6,在2003年１２月前,D_６與Ｄ6一樣，后期其余均為０.我們在原來模型中加入DN1，用D_6代替D６進行回歸:從結果發(fā)現(xiàn),月度特征果然發(fā)生了變化，1月虛擬變量的系數(shù)果然需要被調高０。102,而6月的季節(jié)特征在２00４年后就不再顯著。從而驗證了我們的猜想。經過嘗試預測檢測，沒有出現(xiàn)比較異常的擬合值,故認為已經得到目前最好的模型,即：模型解釋：經過層層分析,我們得到最后的模型,即使一個帶有結構突變、季節(jié)虛擬變量、ARMA成分的模型：首先對于時間趨勢項t，其系數(shù)為0．０17，證明logＹt是有一個穩(wěn)定的時間增長趨勢，對于其它條件保持不變,下一月的產量相對于本月將增加0。017個百分點（可看做一個增長模型),說明我國的乳產業(yè)多年來一直保持向上增長的勢頭.并且由于logＹｔ的線性增長，Yt將呈現(xiàn)指數(shù)增長。季節(jié)虛擬變量Ｄ１、D2、Ｄ６的系數(shù)分別為—0。258、-0.154、0.０95，經過查驗資料表明，1月、2月是產奶淡季，這時候奶牛的產量會因為氣候因素而降低，但是降低的程度在２０0３年1２月以后會減弱;而6月是一個產奶旺季,而且此時社會對乳制品的訂單會迅速增加，這反過來也會刺激產量提升，但這個特征會在20０3年12月以后消失。即在２０03年以前,保持其他條件不變,1月會降低產量的0.258個百分點，２月會降低產量的0．１54個百分點，6月會增加0.095個百分點；在20０3年12月以后，1月會發(fā)生改變，只會降低產量的0。19１個百分點，6月份的季節(jié)特征會消失。對水平突變變量DL,說明從0４年左右開始，乳制品產量出現(xiàn)了一個短時間內的飛躍，這與乳產業(yè)經歷數(shù)年滯積與低迷后發(fā)生迅速的行業(yè)調整有著緊密的關系。蒙牛企業(yè)就是在04年占有的市場份額迅速增加。至于模型預測：1、樣本內預測:我們采用靜態(tài)方法進行樣本內預測,下圖為預測值和真實值的點線圖及誤差結果：其中,logY為實際值,ｌogYF為預測值。我們可以看到，采用樣本內靜態(tài)預測得到的預測值十分接近實際值,誤差比較理想。下圖是取了反函數(shù)得到Y的預測值與真實Y的值：2、接下來我們采用靜態(tài)預測方法進行樣本外一期預測(左序列為預測值，右序列為真實值）：預測值為18４ｔ.通過近一年的真實數(shù)據和預測值的比較,我們發(fā)現(xiàn)預測值的吻合度還是相當高的,預測誤差基本控制在5%以內.通過網上查閱資料,發(fā)現(xiàn)2０1０年7月的產量17９噸,預測誤差為（184—17９)／1７9=0。０２7，即２。7％的預測誤差。故可以認為我們得到一個預測效果很好的模型。課程設計任務書學生姓名:專業(yè)班級：電子科學與技術070１班指導教師：劉金根工作單位:信息工程學院題目：基于ＭATＬAB的連續(xù)時間系統(tǒng)的頻域分析初始條件:ＭＡTLＡB６．5微機要求完成的主要任務:深入研究連續(xù)時間信號和系統(tǒng)時域分析的理論知識.利用MＡTLAB強大的圖形處理功能、符號運算功能以及數(shù)值計算功能,實現(xiàn)連續(xù)時間系統(tǒng)頻域分析.1.利用ＭATLＡB分析系統(tǒng)的頻率特性；2．用ＭＡTLAＢ實現(xiàn)連續(xù)時間信號的采樣及重構；3.撰寫《MAＴLＡＢ應用實踐》課程設計說明書。時間安排：學習ＭATLAＢ語言的概況第１天學習MAＴLＡＢ語言的基本知識第2、3天學習MATLAB語言的應用環(huán)境,調試命令,繪圖能力第4、５天課程設計第６－9天答辯第1０天指導教師簽名：年月日系主任(或責任教師）簽名:年月日目錄摘要…………………Ⅰ１。緒論………………12.對課題內容的分析…………………2２．1連續(xù)時間信號概述………………22.2采樣定理…………22。3總體思路…………23．MATLAB的仿真實現(xiàn)……………３3．1利用MATLAB分析系統(tǒng)的頻率特性…………33。1。1低通濾波器3。1.２高通濾波器3。1．3全通濾波器３.1。4帯通濾波器3.2用MAＴＬAB實現(xiàn)連續(xù)時間信號的采樣及重構…………………6３.２。13.2.23．２．34.心得體會…………1０5。參考文獻…………11附錄…………12摘要本文介紹了基于MＡTLＡB的連續(xù)時間系統(tǒng)的頻域分析。首先利用MＡTLAB分析了系統(tǒng)的頻率特性，分別分析了基于連續(xù)時間系統(tǒng)的低通濾波器、高通濾波器、帶通濾波器、全通濾波器的頻域特性，并依次做出了它們的時域沖激響應波形、頻域內幅頻特性波形、相頻特性波形。在編程過程中分別用到了y=abs()、ｙ＝angle（)、h=fｒeqs(b,a，ｗ）等函數(shù)。然后用MAＴLＡＢ實現(xiàn)了連續(xù)時間信號的采樣及重構，并以f(t）=Sａ（t)為例，分別以過采樣、等采樣、欠采樣三種情況，繪出原信號、采樣信號、重構信號的時域波形圖.關鍵詞：連續(xù)時間系統(tǒng)；頻特性；采樣；重構ＡbstｒａctＴｈｉsarticlｅintroｄucedbａsｅdｏntｈeMＡTＬABrun-onｔiｍesystｅmfrｅqｕencyｒangｅanalysｉs．FirsthasaｎalyzedthesysｔemｆreｑuencycharａcteｒisticusｉｎgMAＴLＡＢ，anａlyzeｄｓeparaｔｅｌybａｓｅdonｔherun—ontｉｍesｙsｔｅmlowpassfilｔer，haspａssｅｄthｅfｉｌｔer，ｔｈebａndpassfiｌteｒｈｉgh,allpaｓsｅsｔhefiltｅrthefreｑｕenｃｙｒａｎgechaｒacteristic，anｄhａsｉntuｒnmadeiｎtheirtiｍｅdomaiｎimpulseresponsepｒｏｆiｌe，theｆｒeｑuencyrangetheａmｐlitude-ｆrequeｎcyｃhａｒａcteｒisticprofｉlｅ，thｅｆrｅqｕenｃycharacterisｔｉcprofiｌｅ.Ｕsedy=abssepａrａt(yī)elyｉntheprogrammiｎgpｒoceｓs(),y＝anglｅ(),h=freqｓ（b，a,w）ａnｄsoｏntｈｅfunctｉons.Ｔhｅnｈａsreａlizedtherun-ontimｅsignａlsamplｉｎgａndresｔructｕｒingwiｔhＭATLAB，aｎdｔａkeｆ（t)=Sa（t）ａstheｅxamｐle，rｅspｅctivｅlybytｈｅsamｐling,ａndsoｏｎthesａｍpｌing,haｓowedｔhesamｐlingｔhｒeekｉｎdoｆsituationｓ，dｒawsthｅorigｉnalsｉｇnal,tｈｅsａｍpｌingｓｉgｎaｌ,ｔhereｓtructuriｎgsignａｌｔimedomａinosｃiｌlogｒam。Kｅｙwｏrd：Run—ontiｍesｙｓtem;Frｅqｕencychａｒacteristｉc;Saｍpliｎg;Restruｃturing１緒論ＭATＬAB是國際上公認的優(yōu)秀科技應用軟件，它的基本功能是數(shù)值計算、符號運算、圖形控制,它的出現(xiàn)給“信號與系統(tǒng)”課程的計算機輔助教學帶來了福音，使利用計算機輔助學生完成“信號與系統(tǒng)”課程的數(shù)值計算、信號與系統(tǒng)分析的可視化建模及仿真調試成為可能。該軟件由公司于是1984年推出,經過十幾年的發(fā)展與完善,目前已成為科技界最流行的應用軟件。它的主要特點是:（１）高效的數(shù)值計算及符號計算功能,能使用戶從繁雜的數(shù)學運算分析中解脫出來.（2)完備的圖形處理功能，實現(xiàn)計算結果和編程的可視化。（3)友好的用戶界面及接近數(shù)學表達式的自然化語言，使學者易于學習和掌握。（４)功能豐富的應用工具箱（如信號處理工具箱），為用戶提供了大量方便實用的處理工具。運用MATLAB對信號與線性系統(tǒng)進行分析與實現(xiàn)的具體方法和過程,其目的在于:（1)讓學生在學習“信號與系統(tǒng)”課程的同時,掌握MATＬAB的應用，對MATLＡＢ語言在低年級學生中的推廣應用起到促進作用。(２)學會應用MAＴLAＢ的數(shù)值計算功能，將學生從繁瑣的數(shù)學運算中解脫出來,從而將便多的時間留于對信號與系統(tǒng)的基本分析方法和應用的理解與思考。（３)讓學生將課程中的重點、難點及部分課后練習用MＡＴLAＢ進行形象、直觀的可視化計算機模擬與仿真實現(xiàn),從而加深對信號與系統(tǒng)基本原理、方法及應用的理解，以培養(yǎng)學生主動獲取知識和獨立解決問題的能力,為學習后繼專業(yè)課打下堅實的基礎.本文將以MATLＡB為工具,對信號與系統(tǒng)在聯(lián)系時間系統(tǒng)的頻域進行分析。２對課題內容的分析２。1連續(xù)時間信號概述在某一時間區(qū)間內,除若干個不連續(xù)點外，如果任意識可都可給出確定的函數(shù)值，則稱該信號為連續(xù)時間信號，簡稱為連續(xù)信號。從嚴格意義上講，ＭATLAＢ數(shù)值計算的方法并不能處理連續(xù)時間信號，然而，可利用連續(xù)信號在等時間間隔的取樣值來近似表示連續(xù)信號，即當取樣時間間隔足夠小時,這些離散樣值能夠被MAＴLＡB處理，并且能較好地近似表示連續(xù)信號。２。2采樣定理對于一個有限頻寬信號進行理想采樣，當采樣頻率時，采樣值唯一確定；當此采樣信號通過截止頻率的理想低通濾波器后，原始信號可以完全重建。通常把最低允許的采樣頻率2稱為奈奎斯特頻率，把最大允許的采樣間隔稱為奈奎斯特間隔。2。3總體思路利用MATLAB強大的圖形處理功能、符號運算功能以及數(shù)值計算功能，實現(xiàn)連續(xù)時間系統(tǒng)頻域分析。其中分析系統(tǒng)的頻率特性時用到了ｙ＝ａbs（）、y=anｇｌe(）、h=freｑｓ(b，a,w）等函數(shù)。實現(xiàn)連續(xù)時間信號的采樣及重構時把采樣分為了過采樣、等采樣、欠采樣三種情況，分別做它們的原信號、采樣信號、重構信號的時域波形圖，其中用到Sa(t）=sinc（t/ｐi）、fa=ｆ＊Tｓ*wc/pi*ｓinｃ(（ｗc／ｐi)＊（oneｓ（lengtｈ(nTｓ）,1)*t－nTs＇*onｅs(１，ｌenｇtｈ(t）))）等函數(shù)。３設計內容３。1利用MＡTＬAB分析系統(tǒng)的頻率特性３。1.1低通濾波器的頻率特性低通濾波器模型:取k＝70,b=10,ｃ=９0可得圖３.1。１圖3．１。１低通濾波器時域沖激響應波形、頻域內幅頻特性、相頻特性。對于低通濾波器，可以從圖3。1。1看出,當w〈10時1．2＞Ｈ(s）〉０．8，而其相頻特性為隨著w的增加而減小并逐漸趨近于—3．1.2高通濾波器高通濾波器模型：以下解法只供參考,方法不唯一:若令１)當,，可以得出：2）考慮一種情況：當，由以上兩條件可得＝—９０－157ｊ，=-9０+1５7j，ｋ=1，并由此可通過matlａb得圖3．1。2圖3。１.２高通濾波器時域沖激響應波形、頻域內幅頻特性、相頻特性。對于高通濾波器，可以從圖３.1.2看出,當ｗ〉50時1.２〉Ｈ（s）>0。8,而其相頻特性為隨著ｗ的增加而減小并逐漸趨近于０。3。1.3全通通濾波器的頻率特性全通濾波器模型：可以取=4,通過matlａb可得圖３。1．3圖3。1。3全通通濾波器時域沖激響應波形、頻域內幅頻特性、相頻特性。當w從0增大時,H（jw）的幅頻特性是一條數(shù)值為１的水平線，即對輸入信號的各頻率分量都進行等值傳輸；而ｊw從0開始下降，最終趨于—.這種網絡稱為全通網絡，在傳輸系統(tǒng)中常用來進行相位校正，如作相位均衡器或移相器。3.1。4帶通濾波器的頻率特性帶通濾波器模型：可?。?5,b=40,=1８,=306由matｌaｂ可得圖3.１。4圖3．1。４帶通濾波器時域沖激響應波形、頻域內幅頻特性、相頻特性。對于帶通濾波器，可以從圖3。1。4看出,當４0＞w〉10時1。2〉H（s)>0．8，而其相頻特性為隨著w的增加先增加而后減小并逐漸趨近于-/23.2用MATLＡB實現(xiàn)連續(xù)時間信號的采樣及重構我們選取信號f(t)Sａ(t)作為被采樣的信號，是因為：第一,f（ｔ）是一個帶限信號，其1;第二，它是一個典型的信號,是分析其他信號的基礎，因此完全有必要對它信號特征詳加了解。此外，應該指出的是,實際信號中,絕大多數(shù)都不是嚴格意義上的帶限信號,這時根據實際精度要求來確定信號的帶寬。對于一個有限頻寬信號進行理想采樣，當采樣頻率時,采樣值唯一確定；當此采樣信號通過截止頻率的理想低通濾波器后,原始信號可以完全重建。通常把最低允許的采樣頻率2稱為奈奎斯特頻率,把最大允許的采樣間隔稱為奈奎斯特間隔.3.2．1過采樣過采樣是使用遠大于奈奎斯特采樣頻率的頻率對輸入信號進行采樣。對于f（t）＝Ｓａ(ｔ),它的帶寬=1,為了由ｆ(t）的采樣信號不失真的重構ｆ（t)，由時域采樣定理知采樣間隔,取,即為過采樣。利用MATＬＡB中的抽樣函數(shù)Ｓinｃ（t)＝sｉｎ(）/來表示Sａ（ｔ），有Sa（t）=Sｉｎc(t／）.圖３。２.1過采樣原信號、采樣信號、重構信號的時域波形圖3.２.２等采樣由采樣定理知f（t)=Sa（ｔ)的采樣間隔,取=,即為等采樣圖3。２。2等采樣3。2。3欠采樣欠采樣的信號頻率(大于fｓ／２),由采樣定理知ｆ（t）=Sa（t）的采樣間隔，取=１．5，即為欠采樣。圖3。2．3欠采樣４.心得體會MATLAB對自己以后的工作和學習很有幫助。在剛開始設計時,有點茫然,連設計的題意都不太明白,在逐步的深入學習、了解之后慢慢有所了解.開始寫程序時也是無從下手，是通過查閱一些圖書資料和網上資料后開始進入正題。在這次課程設計中，我通過多方面地搜集資料，成功地用ＭATＬAB編寫出低通、高通、全通、帯通濾波器時域沖激響應波形、頻域內幅頻特性、相頻特性的仿真實現(xiàn)，及過采樣、等采樣、欠采樣的原信號、采樣信號、重構信號的時域波形圖。通過這次的實踐，我明白了要將理論與實際相結合的道理，盡管這個過程會有一些辛苦,但通過努力實現(xiàn)后，就能大大深化我對知識的理解程度，增長實踐經驗。這表現(xiàn)在我對連續(xù)時間系統(tǒng)的頻域的理解的加深、ＭATLＡB的功能特性都有了進一步的認識?？傊?我在本次課程設計中學到了很多關于ＭATLAB的知識,獲益良多。5.參考文獻［1］孫祥,徐流美，吳清。MAＴLAB７。０基礎教程.北京:清華大學出版社，20０６[2］劉泉，姜雪梅.信號與系統(tǒng).北京：高等教育出版社,2006[3］唐向宏，岳恒立，鄧雪峰．ＭＡTLAB及在電子信息類課程中的應用.北京：電子工業(yè)出版社，２0０6[4］趙靜，張瑾，高新科．基于MATＬＡＢ的通信系統(tǒng)仿真。北京：北京航空航天大學出版社。2００７[5]梁虹．信號與線性系統(tǒng)分析基于MATLＡB的方法與實現(xiàn).北京：電子工業(yè)出版社，20０6.5附錄低通b=90;ａ＝［11０90]；sys＝tf（ｂ,a）;t=0:0.1：5；h=imｐulse(sys，t）;subpｌot（２21）；plot(h）；ｇrｉdxlabeｌ（’ｔ＇）；ylａbel(’h（t)’）;titlｅ（’ｈ（t）’）；［ｈ,w］＝frｅｑs（b，a，100）;h1=abs（h）；h2=anglｅ(ｈ)；ｓuｂploｔ（222）；pｌoｔ（w，h1）；ｇｒidxlａbel（＇角頻率（w）'）;ylabel(’幅度');titｌe(’H（ｊw)的幅頻特性’);suｂpｌot（２23)；plot(ｗ，ｈ2*1８0/pi)；grｉdxlabｅl（'角頻率（w)’）;yｌabel（＇相位（度)'）;tiｔle('Ｈ（jｗ)的相頻特性’)；高通k=1；ｂ=[１00]＊k；a=cｏnv（[1９0+157j],［190-15７j］)；sys=ｔf(b,ａ）；t=0：0.１:10;h=impulse(sys,ｔ）;sｕｂplot（22１);ｐlot（ｈ);gｒidxlabｅl(＇t’）;ylabel（＇h(t）'）;tｉｔle(’h(ｔ）’);［h，ｗ］=ｆｒeｑs（b，a,100)；ｈ1=abs(h);h2＝angle(h)；ｓuｂplot（2２２）;ｐlot（w，h1)；ｇridxlabel(’角頻率(ｗ)'）；ｙｌabel(＇幅度’）;tｉtlｅ（’Ｈ（jw)的幅頻特性’）；subｐlｏt（2２3）；plot(w,h2＊180/ｐｉ)；gridxlabel('角頻率(w）＇）;yｌａｂｅl（＇相位（度）');tiｔle（＇H(jｗ)的相頻特性＇);全通b＝[-14］；a=［１4］；sｙs=ｔf（ｂ，a）;t=０:0.1：10；h＝impulsｅ(syｓ,t）;subplot（２21）;pｌot（h)；ｇridｘlａbel（'t');ｙlabeｌ(’ｈ（t）')；title(’h(ｔ)'）；[ｈ,w］＝freｑs（b，a,100）；h1=abs（h);h２=aｎgｌe（ｈ）;ｓubplot（２22）;ｐlot（w,h1)；griｄaｘis([0，１00,０,1.5]）ｘlabｅl(＇角頻率(w)’）；ｙlabeｌ(＇幅度＇）；ｔiｔｌe（’Ｈ（jｗ)的幅頻特性’）;subｐlot（223）；plot(w,h2*180/pi）；grｉdxlaｂel（'角頻率（w)’）；ｙlａbel(＇相位（度)')；ｔiｔlｅ(’H（jw）的相頻特性’）；帯通k=１；b=［02540]＊k；a=coｎv(［１９+15j］，［１9-15j］）；ｓｙs=ｔｆ(b,a)；t＝0：0.1：１0;h=iｍpulse（ｓys,ｔ）;subpｌｏt（221);plot(h）；ｇridｘlaｂｅｌ(’ｔ'）；ylabeｌ（’h(t)’)；tｉtle(’h(t）’)；［h,ｗ］＝freqs（b,ａ，10０）；h1＝abs(h）；h2=ａｎgle（ｈ)；sｕｂplot（22２)；plｏt(w，h１）;griｄxlａbel(’角頻率(w）')；yｌａbel（＇幅度')；titlｅ（'H（jw)的幅頻特性＇);ｓuｂpｌｏt(2２３）；ｐloｔ(w，h２＊180/pｉ）;ｇridxlabel(＇角頻率（w）’);ylａbel（＇相位(度）＇）;titlｅ(＇H（jw）的相頻特性');過采樣ｔ=—１5：0.０1：15;ｆ=ｓｉn（t）./t;suｂｐloｔ(221）;ｐlｏt(ｔ，f)；xｌａbｅｌ(＇t’);ylaｂel('ｆ(t)’）;titｌe('f（t）=ｓin(t)/t的過采樣原信號')；griｄwm=１；ｗc＝wｍ；Ts=0.02＊pi./ｗm;ws=2＊ｐi．/Ts；n＝－100：100;nTs=ｎ＊Tsｆ=sｉｎc（nTs/pi）;Dt=0。0０5；t＝-１５：Dt:１5；fａ＝f＊Ts＊wc/pi＊ｓｉｎc((wc/pi）＊(ones（lｅｎgｔh(nTs），1)＊t-ｎTｓ’*oｎｅs（1,lenｇｔh（t））））;t1=－1５:0.5：15;f1=sinｃ（t1/pｉ)；sｕｂplｏt（222）；ｓtem(t1，ｆ1)；ｘlaｂeｌ（’kTｓ');ylabel(’ｆ(kTs)’）；title（’sａ（t)=sinc(t/pi）的過采樣采樣信號’）;subplot（2２3)；plot(ｔ，fａ)xlabel(’ｔ’)；ｙlabel（＇ｆa（t）’);tiｔlｅ(＇由sa（ｔ）=sinｃ（t／pi)的過采樣信號重構ｓa(t）’)；gｒｉd；等采樣t=-15:０．0１:15；f＝siｎ（t）．/t;suｂplot（２21)；plｏt（t,f)；xlabel（’t’）;ylａbｅｌ(＇f(t）'）；titｌe（’ｆ(t）＝sin(t）／t的等采樣信號＇）;gridwｍ=1;wc=ｗｍ;Ts＝0.2*ｐi。/wm；wｓ=2＊pｉ。/Tｓ;n=—100:10０；nTｓ＝n＊Tsf=sｉnc（ｎTs／pi)；Dt＝０.０0５；t=-１5：Ｄt：１５;ｆa＝f*Ts＊wc/pi*ｓinc(（wc/pi）*(ｏｎｅs（lｅngｔh（nTs),１)＊t－ｎTs’＊onｅｓ（1,lｅｎgｔh（t））)）；t1＝—１5:0.5：15;f１=sｉｎc（t1／pi）;ｓuｂplot(222)；sｔｅｍ(t1,ｆ1)；xlabｅｌ（'ｋTs＇）；ylａbｅl（'f（ｋTs）’）；titｌe(’ｓa（t)=sｉnc(t/ｐi）的等采樣信號＇）；suｂplｏt（223)；plｏｔ（t，fa)ｘｌabｅl（＇t＇）；ylabｅｌ('fa(t）');title(’由ｓa（t)＝sinc（ｔ/pｉ）的等采樣信號重構sa(ｔ）');grｉd;欠采樣t＝—1５：0.01：15;f=ｓｉn(t）。/t;subpｌｏt（22１）;ｐｌot（t，f）；ｘlabel（＇t’）；ylａbel(’f(t)'）；title(’f(t）＝sin（t)／t的欠采樣原信號'）；ｇｒｉdwm=1；wc＝wm;Ｔs=1。5＊pi./wm;wｓ=2*pｉ．／Ｔs；n=－１00:100；ｎTs=n*Tｓf＝sinc(ｎＴs／pi);Dt=０.００5;ｔ=－15：Ｄt：１5;fａ＝f＊Ts＊wｃ/ｐi＊ｓiｎc(（wｃ/pi）*(ｏnｅs（length（nTs)，１)＊t-nTs'＊oｎes（1,leｎｇｔh(ｔ）)）)；ｔ1＝-１５：0.5:15;f1=sｉnc(ｔ1／pi);ｓubplot（2２2)；ｓtem（t１，f1）；xlabeｌ（＇ｋTs'）;ylaｂeｌ（＇ｆ（ｋＴｓ）’）；title（’sa(ｔ）=ｓiｎｃ（t/ｐi)的欠采樣采樣信號'）;ｓｕbplot（２23)；plｏt（ｔ，fa）ｘlaｂel（’t’）;ylabel（’ｆa（t）’)；titｌe（'由sa（ｔ）=sinc（ｔ／ｐi）的欠采樣信號重構sa(t)'）;ｇriｄ；本科生課程設計成績評定表姓名陳講重性別男專業(yè)、班級電子科學與技術07０1課程設計題目:基于MATLAB的連續(xù)時間系統(tǒng)的頻域分析課程設計答辯或質疑記錄：成績評定依據：最終評定成績（以優(yōu)、良、中、及格、不及格評定)指導教師簽字：年月日《擴頻通信原理》課程設計報告題目：直接擴頻系統(tǒng)仿真班級：01１09１0和0110911姓名：詹曉丹（2009２１0432)姜微(2009210503）張建華(2０09２1033６）指導老師：李兆玉課程設計目的了解、掌握直接擴頻通信系統(tǒng)的組成、工作原理;了解、熟悉擴頻調制、解調、解擴方法,并分析其性能；學習、掌握Ｍａｔlab相關編程知識并用其實現(xiàn)仿真的直接擴頻通信系統(tǒng)；課程設計實驗原理直接擴頻通信系統(tǒng)工作原理：直接序列擴頻,就是直接用高碼率的擴頻碼序列在發(fā)端去擴展信號的頻譜，在收端用相同的擴頻碼去解擴，把展寬的擴頻信號還原成原始的基帶信號。在發(fā)端輸入的信息與擴頻碼發(fā)生器產生的偽隨機碼序列(這里使用的是m序列）進行波形相乘,得到復合信號，實現(xiàn)信號頻譜的展寬，展寬后的信號再調制射頻載波發(fā)送出去。由于采用平衡調制可以提高系統(tǒng)抗偵波的能力，所以直接序列擴頻調制一般都采用二相平衡調制方式。一般擴頻調制時一個信息碼包含一個周期的偽碼，用擴頻后的復合信號對載波進行二相相移監(jiān)控（BＰＳK)調制，當ｇｔ從“0"變成“1"或從“1”變到“0”時，載波相位發(fā)生1８0度相移。接收端的本振信號與發(fā)射端射頻載波相差一個中頻,接收端收到的寬帶射頻信號與本振信號混頻、低頻濾波后得到中頻信號,然后與本地產生的與發(fā)端相同并且同步的擴頻碼序列進行波形相乘，實現(xiàn)相關解擴，再經信息解調，恢復出原始信號。建立模型描述直接擴頻通信系統(tǒng)組成框圖:直接擴頻通信系統(tǒng)波形圖:模塊功能分析直擴系統(tǒng)的調制功能模塊:(都包含模塊框圖和不同調制、解調方式介紹、分析）擴頻調制模塊用擴頻碼發(fā)生器產生一個偽隨機碼ｐn(這里用的是m序列），與信源信息碼序列xt相乘，實現(xiàn)頻譜的展寬ＢＰSＫ調制模塊調制的方式可以有二相相移監(jiān)控BPSK、四相相移鍵控QPＳK、偏移四相相移監(jiān)控OＱPSＫ、最小頻移監(jiān)控ＭSK。QPＳＫ調制的目的是節(jié)省頻譜，但在擴頻系統(tǒng)中有時候帶寬的利用并不是最重要的;OQPSK的優(yōu)點就是調制信號的相位改變沒有倒π現(xiàn)象；MSＫ調制信號時可以避免相位突變，由于以上調制方式實現(xiàn)比較復雜，所以我們選用擴頻系統(tǒng)中最常用的BＰSK調制方式。直擴系統(tǒng)的解調功能模塊:BPＳＫ解調模塊在常規(guī)數(shù)字通信中,解調可以用鎖相環(huán)解調器、平方環(huán)解調器、科斯塔斯環(huán)解調器。在直擴系統(tǒng)中，一般擴頻調制方式是用抑制載波雙平衡調制來產生BPＳＫ信號的，而對于BPSK信號,不管是絕對相移還是差分相移。其載波分量都被抑制了幾十分貝,并且直擴信號譜密度都很低,而大氣噪聲及接收機內部噪聲又很大，有用信號常淹沒在噪聲中，所以用一般的鎖相環(huán)難以提取載波。而平方環(huán)雖然便于載波提取，但環(huán)路工作在二倍頻后的頻率上，工作頻率高，環(huán)路穩(wěn)定性較差。我們選用的是科斯塔斯環(huán)，因為它的突出優(yōu)點是能夠解調移相鍵控信號和抑制了載波的信號，且環(huán)路的工作頻率與輸入信號載波頻率完全相同。擴頻解調模塊解擴方式有相關解擴、直接式相關器解擴、外差式相關器解擴、序列匹配濾波器解擴。直接式相關器的優(yōu)點是結構簡單,缺點是對干擾信號有直通和碼速率泄露現(xiàn)象外差式相關器的抗干擾能力較低；由于相關解擴在性能上很好，在接收端產生的本地ｐｎ’碼,可以用科斯塔斯環(huán)實現(xiàn)與發(fā)端的pn碼精確的同步。模塊源代碼及調試過程（1）直擴系統(tǒng)的調制模塊（a)信息碼生成模塊cｏde_lｅngth=20;％信息碼元個數(shù)Ｎ=1:codｅ_ｌｅngtｈ；rand(＇ｓeｅｄ’,0）；ｘ=ｓign(ｒand（1，coｄe_leｎgth)-０.5）;%信息碼從0、1序列變成—1、1序列x1=ｒeｃtpuｌse（ｘ，8０0）;％每個碼元內采樣800個點ploｔ（x１)；ａxｉs(［０１60０0－１．5１.5］)；ｔiｔle('信源信息碼序列'）;ｇrｉｄon；生成信息碼的波形圖(b）偽隨機碼生成模塊ｆunｃtiｏnｙ=mgeｎ(g,ｓtatｅ,N)%輸入g：m序列生成多項式（1０進制輸入）％state:寄存器初始狀態(tài)（１０進制輸入)%N:輸出序列長度figuｒｅ（1）g＝19;ｓtａt(yī)e=8；N=2000；ｇen＝dｅc2bin（g）-48；Ｍ=lenｇth（gｅｎ）；cuｒStａt(yī)e=dｅc2bin(ｓtatｅ，Ｍ—1)－48;foｒk=1:Nｙ（k）＝cuｒState（Ｍ—１）；a=rｅm（sum（gｅｎ（2:enｄ）。＊cｕrStaｔｅ),2);curStａt(yī)e=［acuｒSｔａｔe(1：M—２)];ｅnｄx＿ｃｏｄe=sign(y—0．５);pn=rｅctｐulse（ｘ_code，8）；%每個偽碼元內采樣8個點ｐｌot(ｐn）;ａｘis(［060０—1。5１。5]）；ｔitle（'偽隨機碼序列’）;gridon；生成偽隨機碼的波形圖（c）擴頻調制模塊gt=ｘ1。*pn;plot(gt);axiｓ（［０１０00—1.51．5］）；ｔiｔｌe（’復合碼序列’);ｇｒidoｎ；生成的復合碼波形圖（d）BＰSK調制模塊%用BPSK調制ｆs=2０e6;ｆ０＝30ｅ6；foｒｉ=１:2000ＡI＝2；dt=fｓ/f0;n＝0:dｔ/７：ｄt;%一個載波周期內采樣8個點cI＝AI*cｏs(2*pｉ＊f0*n/fｓ);x_bｐsk（(1+（i－１）＊８)：ｉ＊8）＝ｇｔ(（１+（ｉ-1）*8)：i＊8）.＊cIｅndplot(x_ｂpsｋ)；axｉs([02０0—2.52。5]);titlｅ('BPSK調制后的波形’）;gｒidon;生成BＰSK調制后的波形圖(２)加噪模塊ｓigma＝0。１;nt=sigma*randn（1,20）；nｔ1=reｃtpｕlｓｅ(nt，800）；ｇt1=gt＋ｎt１；ｆs=20e６；f0＝30e6;fｏrｉ=1：200０AI=２；dt=ｆs/f0;n=0：ｄt/7:dt;％一個載波周期內采樣8個點cＩ＝AI＊cｏs(2＊pi＊f0*n/fｓ）;x_ｂｐsk1(（1+(i-1）*8）：i*８)＝gt1(（１+（ｉ—1)＊８）：ｉ*8）.＊ｃIendｐloｔ（x＿bpsｋ１）;axis（［０2０0—２.５2．5］)；ｔitlｅ(＇加噪后已調波的波形'）;gridon;生成加噪后已調波的波形圖(3）直擴系統(tǒng)的解調模塊（a）ＢPSＫ解調模塊AI=1；dt=fｓ/f0；n=０：dt/7：ｄt；％一個載波周期內采樣八個點ｃI=AＩ*cｏs(2＊pi*f０＊n/fs);fｏｒi=1：２00０s（(1+（ｉ—1）＊8)：ｉ＊8）=x_bpsk1（（1+（ｉ—1）＊8）：i*8)．＊ｃI；endplot(s）；axｉs([0200—2．52.5]）；title（＇解調后的波形＇）;gridon;生成BPSK解調后的波形圖（b）解擴模塊％相關解擴jｋ＿code=s．*pn；%低通濾波ｗｎ＝５/1000000;%截止頻率wｎ=fn/(fｓ/2），這里fn為擴頻碼的帶寬５Mb=fir1(16，wｎ)；H=freqz（b，1，１６000）；ｘx=filter（b，１，jk_ｃｏde）；ｐlot（ｘx）；axis（［０1６0０0-1。51。5]）；title（'解擴并濾波后的波形＇);gｒiｄon;生成解擴后的波形圖調試分析及結論調試分析(a）信源信息碼與收端恢復出的波形的比較(b)頻譜分析調制過程中出現(xiàn)的問題及解決方式：在調試程序中出現(xiàn)“undefinedfunctｉｏｎorｖａrｉａblｅ‘ｍgen’”,后來檢測是子函數(shù)調用格式出現(xiàn)錯誤,經改正后得以修改;復合碼產生程序中，出現(xiàn)“veｃｔｏrｍusｔbｅthesamelｅnｇｔhｓ",后檢測得出原因是信息碼和偽隨機碼長度不相同時不能直接相乘；由于所加噪聲為高斯白噪聲,是用raｎｄｎ函數(shù)產生的，具有隨機性，所以不同時刻加噪了的已調波具有不同波形；結論通過調試分析及研究，我們發(fā)現(xiàn)直擴系統(tǒng)具有較多的特點,將其運用于通信系統(tǒng)能發(fā)揮其很大的優(yōu)勢。首先,直擴系統(tǒng)具有很高的處理增益，因此直擴系統(tǒng)具有抗寬帶干擾、抗多頻干擾及單頻干擾的能力;其次，直擴信號的功率譜密度低,具有隱蔽性和低的截獲概率,從而抗截獲的能力強，另外功率污染少;第三，直擴偽隨機序列的偽隨機性和密鑰量具有保密性，即系統(tǒng)本身具有加密的能力；第四,利用直擴偽隨機序列碼型的正交性,可構成直擴碼分多址系統(tǒng)；第五，利用直擴信號的相關接收,具有抗多徑干擾的能力;第六，利用直擴信號可實現(xiàn)精確的測距定位。心得體會通過這次直接擴頻系統(tǒng)的仿真實驗，我們加深了對擴頻理論知識的理解。懂得了直接擴頻系統(tǒng)的原理所在，并且經過分析，自己應用MＡTLAＢ實現(xiàn)了直接擴頻系統(tǒng)的仿真,掌握了一定的軟件仿真能力.通過頻譜分析,更使我們了解到了直接擴頻的抗干擾性能。在實踐過程中,我們小組也遇到了一些問題，但是我們通過討論以及查閱資料，解決了難題.同時,也激發(fā)了我們對擴頻這門課程學習的積極性。在此次直接擴頻系統(tǒng)的仿真實驗中，我們收獲頗豐,記憶深刻，并全面培養(yǎng)了我們的能力。參考文獻劉煥淋.擴展頻譜通信［M].北京：北京郵電大學出版社，２008田日才.擴頻通信[M］.北京：清華大學出版社，2００7曾一凡，李暉.擴頻通信原理[Ｍ]。北京：機械工業(yè)出版社，2０05附件一：程序源代碼ｆuｎctionｙ=mｇｅn（g，stａｔe，N）％輸入g:m序列生成多項式（10進制輸入)％sｔate:寄存器初始狀態(tài)(１０進制輸入）%N:輸出序列長度figｕrｅ（１）ｇ=19；ｓｔａt(yī)e＝８；Ｎ＝2000；gｅn=dec２bin(g）-4８；Ｍ=length（gen)；ｃurState＝dec2bin（stａt(yī)e，Ｍ—１）－48;ｆｏrk=1:Ｎｙ（ｋ)＝ｃurStatｅ(M—１)；a=rem（suｍ（gen（2:end）。*curStaｔe）,2）；curState=［acuｒSｔａｔe（1：Ｍ—２）］；eｎdx_code=ｓign（y－0.5)；pｎ=rｅctpulｓe(x_cｏｄe，8)；%每個偽碼元內采樣８個點ploｔ(pn）;axis([06００－１.51.５])；tiｔlｅ（'偽隨機碼序列＇）；gridon;figurｅ（2)code_ｌｅngtｈ=20；%信息碼元個數(shù)N＝1:code_leｎgth；raｎd(’seeｄ'，0）；ｘ＝ｓiｇn（ｒand(1,ｃode_lｅngｔh)—０。5)；％信息碼從0、1序列變成-1、1序列ｘ1＝rｅctpulse（x，８00)；%每個碼元內采樣800個點plｏｔ(ｘ1)；axis([016０00—1。51。5]);tｉtｌe（'信源信息碼序列'）；ｇridon；figure（3)gt=x1.＊pn；pｌot（gt）;aｘis（［01000-１。51.5］）;tiｔｌe（'復合碼序列'）;gridｏn;ｆiguｒe（４）;%用BPSK調制ｆｓ=２0ｅ６；f０=30e6；fｏrｉ=1:2000ＡI＝2；dt=fs/f0；n=0：dt/7：dt;％一個載波周期內采樣8個點ｃＩ=ＡI＊ｃoｓ（2*pi＊ｆ0*n/fｓ）；x_bpｓｋ(（１+（i-1）＊８）:i＊8）=gｔ(（1＋（i—1)*8）:i＊8）.＊cIendplot(ｘ_bｐｓk）；axis(［0２00—2．52.5］);ｔitle（＇BPSK調制后的波形'）；griｄoｎ;ｆiｇure（5）sｉgma=0.1；nｔ=sigma＊rａｎdn（1，２0）；nt1=recｔpulsｅ(ｎt，８00);gt１=ｇt＋ｎｔ1;fs=20e6；f０=３0e6;fori＝1:２００0AI=2；dt＝fｓ/f0；n=０:ｄt/7：dt；％一個載波周期內采樣８個點cI=AI＊cｏs（2*pi＊f0*n/fs);ｘ＿bpsk１((1＋（i－1)＊8）：i＊8)=gt1（（1＋（i-１)＊8）：ｉ*8)。＊cＩeｎdploｔ(x_bpｓk1）；ａxｉs（［02００-2。52.5］）；tｉtle（’加噪后已調波的波形'）；gｒidon；fiｇｕre（６）％解調ＡI=1；dt=fs/f0；n=0：dt/7：dt;％一個載波周期內采樣八個點cＩ=ＡI*cos(2＊pi＊ｆ0＊n/fs）；forｉ=１：2000ｓ（（1+（i—1)＊8）:i＊8）=x_bpsk１（（１+（i-1）*８)：i＊８）.*cI;endploｔ(s）;ａxiｓ（[0200－2.52.５]）；title（'解調后的波形＇）;ｇridon；figuｒe（7）％相關解擴jk＿ｃodｅ＝s.＊pｎ;％低通濾波wn＝５/1000000；%截止頻率ｗn=fｎ/(fs／２)，這里ｆｎ為擴頻碼的帶寬5Ｍb＝fir１(16,ｗｎ）；H=fｒeｑz(b，1,16000）;xx=ｆilter(b,1，jk_codｅ)；plot（ｘx）;ａxis（［016０0０-1。５１。5］）;tｉtle（'解擴并濾波后的波形＇)；grｉｄｏn；%信源信息碼與收端恢復出的波形的比較figuｒe（8）subplot（2，１，1）;plｏt（x1）;aｘiｓ（[01６00０—1.５1．５]）；tiｔle（＇信源信息碼序列’）;gridon；subplｏt(2，１，2)；ｐlot(xｘ)；axｉs（[016000-1．51．5])；ｔitｌｅ(’收端接收到的波形＇)；gridoｎ；％頻譜圖figｕre（9)；Ｔ=1；N＝２０；Ｎ_saｍpｌｅ=800；df=１/（N＊Ｔ)；f=－（Ｎ＊N_saｍple/2）＊df+df:df：（N＊N_sａmｐlｅ/2）＊df；ｓubplｏt(２，2，１)；pｌoｔ(f,ａbs(fftshift(fｆt(x1)）));ａxiｓ([—1０01000２０00］）;title（'信源信息碼序列的頻譜圖');griｄｏｎ;T＝1;Ｎ=2000;N_sample=8；df=１／(N＊T）；f=－（N*Ｎ_saｍｐlｅ/2)*df+df：df：（N＊Ｎ_saｍple/2)＊df；sｕbｐlｏt(2，２,2）；ｐlot（ｆ，abｓ(ｆftshｉft（ffｔ（pn）））);aｘis（［-4４０4００0］);ｔitlｅ（’偽隨機碼序列的頻譜圖’)；ｇｒidoｎ；T＝1;N=20０0;N_sampｌe=8；ｄｆ=1/(N＊T）;f＝—(N＊N＿ｓample/2)＊ｄf＋df:df:（Ｎ*Ｎ_sａmple/2）*ｄf；subpｌｏｔ（２，２,３);ploｔ（f,abｓ（fｆtｓhifｔ（ｆft(ｇt））)）;axis([—4402000］）；ｔitlｅ（'擴頻碼序列的頻譜圖’）；gｒidoｎ；T=１；N＝２０00;Ｎ_samｐle=８;df=１／(N*Ｔ）;f=-(N＊N_sａｍpｌe/2)＊df＋df:df：（N＊Ｎ_ｓample／2）*df；subplot（2，２，4）;plot(f，aｂｓ(ｆftsｈift(fft（x_bpsk））））；axis(［—440２000］)；titｌe（’bpsｋ已調波的頻譜圖’);gｒidｏn;figｕrｅ（1０）；T＝１;N=2000;N_sａmpｌe=8;dｆ＝１/(N*Ｔ）；f＝-(N＊N_samplｅ／2)＊dｆ＋dｆ：df：(N＊N_sａｍplｅ/２)＊df；subplot（2,２，１)；ｐlot（f,abs(fｆtｓhｉft（ffｔ（x_ｂpsｋ１））））；axis([-4402000]);ｔitle(’加噪已調波的頻譜圖＇）；gridon;T=1；N=20０0;N_sample=8；df=1/（N＊T)；f＝—（N＊N_sａmpｌe/2）*df＋dｆ：ｄf：(N＊N_sａｍple/2）＊dｆ；sｕｂplot(2,2，2）；ｐlot（f,abs(ｆftｓhift（ffｔ(s））））;axis(［-４4020００]）；title（’解調后的頻譜圖’)；grｉdon;T=１；N=２0；N_samplｅ=８00;df=１/(N*Ｔ）；ｆ=—(Ｎ*N＿sample/２）＊df+df:ｄf：（N＊N＿saｍple/２）＊df;subplot(２,2,3）;plｏt(f,abs（fftｓhift（ｆfｔ(xx））））；aｘis（[－100１００020０0］）；tiｔle(’解擴后收端恢復出的信號的頻譜圖＇）；gridon;成都信息工程學院畢業(yè)論文基于關鍵鏈的項目時間管理方法研究論文作者姓名：XＸX專業(yè)：XＸXＸXXXＸＸ（講師）論文提交日期:2０11年0５月2５日基于關鍵鏈的項目時間管理方法研究摘要“按時、保質地完成項目”必定是每一位項目經理最希望做到的，但工期拖延的情況卻時常發(fā)生。因而合理地安排項目時間是項目管理中一項關鍵內容，它的目的是保證按時完成項目、合理分配資源、發(fā)揮最佳工作效率。本文主要是基于關鍵鏈的項目時間管理方法的研究,關鍵鏈項目管理方法是約束理論在項目管理中的應用.本文主要分析的是在對關鍵鏈項目管理方法的基本思想和理論基礎分析介紹的基礎上，針對如何利用關鍵鏈將項目計劃與項目進度控制更好的結合進行了研究和評價，利用項目管理領域相關知識,補充了緩沖區(qū)管理的方法和實施步驟,并運用了實例進行了求證。關鍵詞：項目時間管理；項目進度制定方法；關鍵路徑；關鍵鏈;案例分析BａseｄoｎｋｅｙｃhainsｏｆthepｒoｊectｔiｍemaｎagemeｎtmeｔhoｄＡｂｓtract”Coｎditionandontｉmｅcompletｉonｏfpｒojｅｃts"evｅryprｏｊｅcｔｍaｎagerｉsprｏbablythｅmｏｓthope。However,thｅsｉtuaｔｉonｉsscｈeduledelaｙswerefrequenｔ.Thusａreasonableｔiｍetoarrangeｔhepｒojectisakeyiｎprojｅctmanagemｅｎt，itspurposeｉstoenｓureｔimelyｃｏmpｌetioｎ，ｒatｉoｎalａllｏｃatｉoｎｏｆrｅsouｒces，ｔｏachｉeｖetheｂestｅfｆｉｃiencｙ．Thｉspaperisｍａiｎｌｙbaｓｅdonkeｙchaiｎprojｅｃttiｍemanageｍenｔｍethodsofｒeｓeａrcｈ，kｅychainprｏｊecｔmaｎageｍentmetｈodiｓcoｎstraiｎｔtｈeｏryaｐplｉcａt(yī)ｉoｎinpｒojectmanagｅment．Thｉspａpeｒmaiｎlyanaｌｙsisoｆkeycｈaiｎｉnthｅprｏjｅcｔｍanａgemｅｎtmethｏｄofbasｉcideｏlｏgｉｃalandｔheoretｉcalbasisofａnalysｉsｉntroｄucｔｉon，basedonhowtoｕｓetｈｅkeychaiｎwiｌｌprojecｔplanninｇａnｄｐｒｏjｅｃtsｃheduleconｔrｏlbetterisstｕdｉedandthecombinａt(yī)ionofevaluatｉoｎｏｆprojectmａnageｍentdｏmainknowledge，ｔocｏｍpｌemｅntｔｈebuｆfermａｎagｅmeｎｔmｅｔｈodanｄtheiｍｐlemeｎtationsteps，aｎdusingtheexampleｖeｒifｉcatｉon.Kｅｙwords：ProｊｅctTimeＭａnagemｅｎｔ;projectscheduledｅvelopmeｎｔmethod；criｔiｃａlｐatｈ;kｅychaｉn；Casｅstudy目錄論文總頁數(shù)：1６頁TOC＼o"1-3"\h\ｚ\uＨYPEＲLIＮK是對人們平時所形成的一種學習工作行為和習慣的通俗描述。面對一項任務，人們在做計劃安排時一般習慣盡量讓時間留有余地（通?？紤]了不確定性和風險因素)，極力爭取更多的安全時間。然而實際在計劃執(zhí)行時這些影響工作的風險因素未必真正發(fā)生，而因為人有惰性，知道自己有充足的時間,開始時便覺得不必著急,沒必要全力以赴,等到時間變得緊迫時，才開始努力去做,這樣使原本富裕的時間被浪費掉，出現(xiàn)前松后緊的現(xiàn)象,一旦有意外發(fā)生，往往造成工期的延誤.（二）多任務效應，反應了共享資源地項目執(zhí)行的約束作用。假設某項目成員（共享資源)同一時間內需要完成A、B和Ｃ三個任務,每個任務周期都為10天.但在多項目環(huán)境、情緒壓力或人為干預等因素影響下，上訴任務通常以輪換形式執(zhí)行,那么每個任務周期將延長一倍.可見項目的多任務效應將直接導致項目延誤。（三)工期變動效應的不對稱性,在項目實施進程中,任務的延遲通常會由于工序的邏輯關系產生級聯(lián)效應，但項目卻很難從某些任務的提前完成中受益。例如:A、Ｂ、C是三個可并發(fā)執(zhí)行的任務，Ｄ為Ａ、B、Ｃ的緊后任務,A、B、C的任務周期為５天,D的任務周期為１0天。通常該項目需要15天的時間就可以完成.但是任務Ｂ在實施過程中出現(xiàn)了意外推遲了3天完成,從而導致整個項目周期延期至1８天。如果現(xiàn)在反過來思考，如果A、B、C能夠提前2天完成,那么整個項目周期是否就變成了13天呢?事實上答案是否定的。分析如下:首先如果項目組成員提前完成了任務,他不會因此而得到獎勵，反而卻導致項目經理對類似任務的預估完成時間將減少，顯然這對他沒有好處，因此他不會愿意提前完工.另外即使A、Ｂ、Ｃ能夠提前完成，但D也許因為資源尚未就位等原因而無法開始，所以項目就無法提前完成［１１].關鍵鏈技術就是在項目的執(zhí)行過程中，把每道工序節(jié)省下來的安全時間綜合利用起來,通過設置緩沖（ PB），輸入緩沖(ＦＢ）和資源緩沖(RB）來降低風險，保障項目的順利運行。①項目緩沖區(qū)的設置在傳統(tǒng)的項目管理中，認為一個項目能否如期完成與每一個活動完成日期有關，因此都會在每一個工序上加上很多的安全時間，來確保項目能準時完成,但人們因為惰性,知道自己又充足的時間，開始不全力以赴，等到時間變得緊迫時,才開始努力去做，這樣使原本富裕的時間被浪費，而且不會累積到下一個活動，因此刪減每個活動的預估時間,便會釋放出足夠的時間來設置項目緩沖（PB）。項目緩沖一般設置在關鍵鏈的尾端,通過設置項目緩沖區(qū)將延誤控制在預定的范圍內,保證項目如期交付，以關鍵鏈上各工序節(jié)省時間總和的一半作為其大小。②輸入緩沖區(qū)的設置雖然關鍵路徑有了項目緩沖，但當有非關鍵鏈路徑與關鍵路徑會合時,就必須在非關鍵路徑有可能延遲關鍵路徑上的活動，因此必須在非關鍵路徑與關鍵路徑會合的地方插入輸入緩沖ＦB．ＦB主要是保護關鍵路徑，可以消化非關鍵鏈工序帶來的工期延遲。輸入緩沖一般設置在非關鍵鏈與關鍵鏈的匯合處，主要是為了保證關鍵鏈上的工序如期開始而不受非關鍵鏈上工序延誤的影響,以非關鍵鏈上節(jié)省時間總和的一般作為其大小.③資源緩沖區(qū)的設置由于關鍵鏈上任何活動的延誤都會引起整個項目的延誤,關鍵鏈上的工序一般被賦予最高優(yōu)先級。資源緩沖ＲB主要放置于有可能發(fā)生資源爭用的活動上。資源緩沖是一種時間緩沖,本質上是一種預警機制,當關鍵鏈上的活動準備使用某瓶頸資源時，通過預告機制即在工序來臨前發(fā)出預先警告提醒人們。項目管理者可以通過對緩沖區(qū)的監(jiān)控了解緩沖的剩余情況和項目的進展情況,進而及時采取措施進行調控。資源緩沖一般設置在關鍵鏈上，通過制定一定的前置時間來進行設置，以保證關鍵鏈上所需要資源能按時到位[9]。④關鍵鏈技術的實施步驟根據關鍵鏈技術的基本思想，對于項目進度計劃可以按照以下步驟編制:繪制網絡圖，找出關鍵工序，運用工作分