簡(jiǎn)介：機(jī)電系統(tǒng)動(dòng)力學(xué)仿真實(shí)驗(yàn)報(bào)告機(jī)電系統(tǒng)動(dòng)力學(xué)仿真實(shí)驗(yàn)報(bào)告機(jī)工機(jī)工1004班201003130150宋康宋康習(xí)題一習(xí)題一13ABCWU_2004兩個(gè)是合法的。14（1）（2）（3）CLEARX8058YXXONESSIZEYXYYONESSIZEX11110505101050510050051XYZ習(xí)題二習(xí)題二21

簡(jiǎn)介：_畢業(yè)設(shè)計(jì)畢業(yè)論文題目畢業(yè)論文題目基于基于MATLABMATLAB的指紋識(shí)指紋識(shí)別院系班級(jí)學(xué)生姓名學(xué)生姓名學(xué)號(hào)指導(dǎo)教師指導(dǎo)教師時(shí)間時(shí)間_442WALK函數(shù)32443LAST1函數(shù)33第五章特征點(diǎn)匹配3551紋線長(zhǎng)度匹配3552三角形邊長(zhǎng)匹配3553點(diǎn)類型匹配36第六章實(shí)驗(yàn)過(guò)程38第七章總結(jié)45參考文獻(xiàn)45致謝47

簡(jiǎn)介：1擴(kuò)頻通信原理課程設(shè)計(jì)報(bào)告擴(kuò)頻通信原理課程設(shè)計(jì)報(bào)告題目直接擴(kuò)頻系統(tǒng)仿真題目直接擴(kuò)頻系統(tǒng)仿真班級(jí)班級(jí)0110910和0110911姓名姓名詹曉丹（詹曉丹（2009210432）姜微（微（2009210503）張建華（張建華（2009210336）指導(dǎo)老師李兆玉指導(dǎo)老師李兆玉34模塊功能分析（1）直擴(kuò)系統(tǒng)的調(diào)制功能模塊（都包含模塊框圖和不同調(diào)制、解調(diào)方式介紹、分析）（A）擴(kuò)頻調(diào)制模塊用擴(kuò)頻碼發(fā)生器產(chǎn)生一個(gè)偽隨機(jī)碼PN（這里用的是M序列），與信源信息碼序列XT相乘，實(shí)現(xiàn)頻譜的展寬（B）BPSK調(diào)制模塊

簡(jiǎn)介：MATLABMATLAB繪制三維圖形的方法繪制三維圖形的方法PLOT3函數(shù)與PLOT函數(shù)用法十分相似，其調(diào)用格式為PLOT3X1Y1Z1選項(xiàng)1X2Y2Z2選項(xiàng)2XNYNZN選項(xiàng)N，其中每一組XYZ組成一組曲線的坐標(biāo)參數(shù)，選項(xiàng)的定義和PLOT函數(shù)相同。當(dāng)XYZ是同維向量時(shí)，則XYZ對(duì)應(yīng)元素構(gòu)成一條三維曲線。當(dāng)XYZ是同維矩陣時(shí)，則以XYZ對(duì)應(yīng)列元素繪制三維曲線，曲線條數(shù)等于矩陣列數(shù)。例繪制三維曲線。程序如下T0PI10020PIXSINTYCOSTZTSINTCOSTPLOT3XYZGRIDTITLELINEIN3DSPACEXLABELXYLABELYZLABELZ如下圖10500511050051402002040XLINEIN3DSPACEYZ三維曲面三維曲面程序如下XYMESHGRID8058ZSINSQRTX2Y2SQRTX2Y2EPSSUBPLOT221MESHXYZTITLEMESHXYZSUBPLOT222MESHCXYZTITLEMESHCXYZSUBPLOT223MESHZXYZTITLEMESHZXYZSUBPLOT224SURFXYZTITLESURFXYZ如下圖1001010010101MESHXYZ1001010010101MESHCXYZ1001010010101MESHZXYZ1001010010101SURFXYZ3標(biāo)準(zhǔn)三維曲面SPHERE函數(shù)的調(diào)用格式為XYZSPHERENGENERATESTHREEN1BYN1MATRICESSOTHATSURFXYZPRODUCESAUNITSPHERECYLINDER函數(shù)的調(diào)用格式為XYZCYLINDERRN，其中R為圓周半徑，N為組成圓周的點(diǎn)數(shù)。MATLAB還有一個(gè)PEAKS函數(shù)，稱為多峰函數(shù)，常用于三維曲面的演示。例繪制標(biāo)準(zhǔn)三維曲面圖形。

簡(jiǎn)介：MATLAB特殊圖形和高維可視化20091020010674特殊圖形和高維可視化741特殊圖形指令例示7411面域圖AREA【例74111】面域圖指令A(yù)REA。該指令的特點(diǎn)是在圖上繪制多條曲線時(shí)，每條曲線（除第一條外）都是把“前”條曲線作基線，再取值繪制而成。因此，該指令所畫的圖形，能醒目地反映各因素對(duì)最終結(jié)果的貢獻(xiàn)份額。注意（1）AREA的第一輸入宗量是單調(diào)變化的自變量。第二輸入宗量是“各因素”的函數(shù)值矩陣，且每個(gè)“因素”的數(shù)據(jù)取列向量形式排放。第三輸入宗量是繪圖的基準(zhǔn)線值，只能取標(biāo)量。當(dāng)基準(zhǔn)值為0（即以X軸為基準(zhǔn)線）時(shí)，第三輸入宗量可以缺省。（2）本例第條指令書寫格式XY，強(qiáng)調(diào)沿列方向畫各條曲線的事實(shí)。CLFX22注意自變量要單調(diào)變化Y352413452154325各因素的相對(duì)貢獻(xiàn)份額CUM_SUMCUMSUMY各曲線在圖上的絕對(duì)坐標(biāo)AREAXY0LEGEND因素A因素B因素CGRIDONCOLMAPSPRINGX21012Y352413452154325CUM_SUM352416976211131087圖74122三維直方圖7413餅圖PIEPIE3【例74131】餅圖指令PIEPIEPIE3PIE3用來(lái)表示各元素占總和的百分?jǐn)?shù)。該指令第二輸入宗量為與第一宗量同長(zhǎng)的01向量，1使對(duì)應(yīng)扇塊突出。A116120821SUBPLOT121PIEA10100LEGEND12345SUBPLOT122PIE3AAMINACOLMAPCOOL圖74131餅形統(tǒng)計(jì)圖7414填色圖FILLFILL3【例74141】讀者試驗(yàn)本例時(shí)，注意三點(diǎn)MATLAB畫任意多邊形的一種方法；保證繪圖數(shù)據(jù)首尾重合，使勾畫多邊形封閉；使用圖柄對(duì)圖形的屬性進(jìn)

簡(jiǎn)介：1MATLAB及應(yīng)用上機(jī)作業(yè)及應(yīng)用上機(jī)作業(yè)學(xué)院名稱機(jī)械工程學(xué)院專業(yè)班級(jí)測(cè)控1201學(xué)生姓名學(xué)生學(xué)號(hào)201年4月3作業(yè)作業(yè)11、用MATLAB可以識(shí)別的格式輸入下面兩個(gè)矩陣，12342357135732391894A144367723355422675342189543IIBI再求出它們的乘積矩陣，并將矩陣的右下角子矩陣賦給矩陣。賦值完成CC23D后，調(diào)用相應(yīng)的命令查看MATLAB工作空間的占有情況。解A12342357135732391894B14I436772335542I267I5342189543B14I436772335542I267I5342189543CABDC45462、設(shè)矩陣，求，，，，，并求矩陣16231351110897612414152A1A3A12AA13AAA的特征值和特征向量。解A16231351110897612414152DETAINVA

簡(jiǎn)介：控制系統(tǒng)仿真實(shí)驗(yàn)MATLAB部分實(shí)驗(yàn)結(jié)果1實(shí)驗(yàn)一實(shí)驗(yàn)一MATLAB基本操作基本操作1用MATLAB可以識(shí)別的格式輸入下面兩個(gè)矩陣12332357135732391894A144367823355422675342189543IIBI再求出它們的乘積矩陣C，并將C矩陣的右下角23子矩陣賦給D矩陣。賦值完成后，調(diào)用相應(yīng)的命令查看MATLAB工作空間的占用情況。A12332357135732391894B14I436782335542I267I5342189543CABDC4546WHOSNAMESIZEBYTESCLASSATTRIBUTESA5X4160DOUBLEB4X6384DOUBLECOMPLEXC5X6480DOUBLECOMPLEXD2X396DOUBLECOMPLEX2選擇合適的步距繪制出下面的圖形，其中1SINT11TT1011YSIN1TPLOTTY

簡(jiǎn)介：課程設(shè)計(jì)實(shí)驗(yàn)體會(huì)課程設(shè)計(jì)實(shí)驗(yàn)體會(huì)學(xué)生姓名學(xué)生姓名李祥勝李祥勝學(xué)生學(xué)號(hào)學(xué)生學(xué)號(hào)20120704專業(yè)班級(jí)專業(yè)班級(jí)光信息科學(xué)與技術(shù)光信息科學(xué)與技術(shù)指導(dǎo)老師指導(dǎo)老師MISSCHEN學(xué)院院信息工程學(xué)院信息工程學(xué)院題目MATLAB學(xué)期實(shí)驗(yàn)總結(jié)學(xué)期實(shí)驗(yàn)總結(jié)操作系統(tǒng)上運(yùn)行。（6）MATLAB的圖形功能強(qiáng)大。在FTRAN和C語(yǔ)言里，繪圖都很不容易，但在MATLAB里，數(shù)據(jù)的可視化非常簡(jiǎn)單。MATLAB還具有較強(qiáng)的編輯圖形界面的能力。（7MATLAB的缺點(diǎn)是，它和其他高級(jí)程序相比，程序的執(zhí)行速度較慢。由于MATLAB的程序不用編譯等預(yù)處理，也不生成可執(zhí)行文件，程序?yàn)榻忉寛?zhí)行，所以速度較慢。8功能強(qiáng)勁的工具箱是MATLAB的另一重大特色。MATLAB包含兩個(gè)部分核心部分和各種可選的工具箱。核心部分中有數(shù)百個(gè)核心內(nèi)部函數(shù)。其工具箱又可分為兩類功能性工具箱和學(xué)科性工具箱。功能性工具箱主要用來(lái)擴(kuò)充其符號(hào)計(jì)算功能、圖示建模仿真功能、文字處理功能以及與硬件實(shí)時(shí)交互功能。功能性工具箱能用于多種學(xué)科。而學(xué)科性工具箱是專業(yè)性比較強(qiáng)的，如CONTROL、TOOLBOX、SIGNALPROCESSINGTOOLBOX、COMMUNICATIONTOOLBOX等。這些工具箱都是由該領(lǐng)域內(nèi)的學(xué)術(shù)水平很高的專家編寫的，所以用戶無(wú)需編寫自己學(xué)科范圍內(nèi)的基礎(chǔ)程序，而直接進(jìn)行高、精、尖的研究。下表列出了MATLAB的核心部分及其工具箱等產(chǎn)品系列的主要應(yīng)用領(lǐng)域。9源程序的開放性。開放性也許是MATLAB最受人們歡迎的特點(diǎn)。除內(nèi)部函數(shù)以外，所有MATLAB的核心文件和工具箱文件都是可讀可改的源文件，用戶可通過(guò)對(duì)源文件的修改以及加入自己的文件構(gòu)成新的工具箱。112112優(yōu)勢(shì)優(yōu)勢(shì)（1）友好的工作平臺(tái)和編程環(huán)境MATLAB由一系列工具組成。這些工具方便用戶使用MATLAB的函數(shù)和文件，其中許多工具采用的是圖形用戶界面。包括MATLAB桌面和命令窗口、歷史命令窗口、編輯器和調(diào)試器、路徑搜索和用于用戶瀏覽幫助、工作空間、文件的瀏覽器。隨著MATLAB的商業(yè)化以及軟件本身的不斷升級(jí)，MATLAB的用戶界面也越來(lái)越精致，更加接近WINDOWS的標(biāo)準(zhǔn)界面，人機(jī)交互性更強(qiáng)，操作更簡(jiǎn)單。而且新版本的MATLAB提供了完整的聯(lián)機(jī)查詢、幫助系統(tǒng)，極大的方便了用戶的使用。簡(jiǎn)單的編程環(huán)境提供了比較完備的調(diào)試系統(tǒng)，程序不必經(jīng)過(guò)編譯就可以直接運(yùn)行，而且能夠及時(shí)地報(bào)告出現(xiàn)的錯(cuò)誤及進(jìn)行出錯(cuò)原因分析。（2）簡(jiǎn)單易用的程序語(yǔ)言MATLAB一個(gè)高級(jí)的矩陣陣列語(yǔ)言，它包含控制語(yǔ)句、函數(shù)、數(shù)據(jù)結(jié)構(gòu)、輸入和輸出和面向?qū)ο缶幊烫攸c(diǎn)。用戶可以在命令窗口中將輸入語(yǔ)句與執(zhí)行命令同步，也可以先編寫好一個(gè)較大的復(fù)雜的應(yīng)用程序（M文件）后再一起運(yùn)行。新版本的MATLAB語(yǔ)言是基于最為流行的C＋＋語(yǔ)言基礎(chǔ)上的，因此語(yǔ)法特征與C＋＋語(yǔ)言極為相似，而且更加簡(jiǎn)單，更加符合科技人員對(duì)數(shù)學(xué)表達(dá)式的書寫格式。使之更利于非計(jì)算機(jī)專業(yè)的科技人員使用。而且這種語(yǔ)言可移植性好、可拓展性極強(qiáng)，這也是MATLAB能夠深入到科學(xué)研究及工程計(jì)算各個(gè)領(lǐng)域的重要原因。（3）強(qiáng)大的科學(xué)計(jì)算機(jī)數(shù)據(jù)處理能力MATLAB是一個(gè)包含大量計(jì)算算法的集合。其擁有600多個(gè)工程中要用到的數(shù)學(xué)運(yùn)算函數(shù)，可以方便的實(shí)現(xiàn)用戶所需的各種計(jì)算功能。函數(shù)中所使用的算法都是科研和工程計(jì)算中的最新研究成果，而前經(jīng)過(guò)了各種優(yōu)化和容錯(cuò)處理。

簡(jiǎn)介：MATLAB程式設(shè)計(jì)入門篇程式碼與記憶體之最佳化,張智星JANGMIRLABORGHTTP//MIRLABORG/JANG清大資工系多媒體檢索實(shí)驗(yàn)室,重點(diǎn)說(shuō)明65版之前，MATLAB執(zhí)行效率非常倚賴程式碼的向量化65（含）版後，MATLAB引進(jìn)了JIT加速器，程式碼向量化的重要性相對(duì)降低。同學(xué)們?nèi)皂毑t解程式碼向量化的實(shí)作方法、優(yōu)缺點(diǎn)。本章範(fàn)例之執(zhí)行環(huán)境OSWINDOWS7CPUINTELCOREI72670QM,220GHZRAM8GBMATLAB7120635R2011A,21程式碼的向量化,65版之前，MATLAB程式環(huán)境是一個(gè)傳統(tǒng)的解譯器（INTERPRETER），在執(zhí)行MATLAB的程式碼時(shí)，會(huì)進(jìn)行下列動(dòng)作逐列對(duì)程式碼轉(zhuǎn)換為PCODE，這是MATLAB可以解讀的格式。對(duì)產(chǎn)生的PCODE進(jìn)行逐列執(zhí)行。,MATLAB程式運(yùn)作流程I,在執(zhí)行每一列PCODE時(shí)，都還包含一些「經(jīng)常開銷」（OVERHEAD）若該列執(zhí)行需要大量運(yùn)算，經(jīng)常開銷就會(huì)顯得很渺小而不會(huì)拖累程式執(zhí)行速度。若該列只是簡(jiǎn)單的運(yùn)算，那經(jīng)常開銷的比例就會(huì)相對(duì)提高。若該列是放在迴圈內(nèi)，那麼程式碼的執(zhí)行速度就會(huì)被這些經(jīng)常開銷大幅拖慢。,MATLAB程式運(yùn)作流程II,若要加快執(zhí)行速度MATLAB65版之前，就要盡量不使用迴圈，而盡量改用向量化的運(yùn)算，以降低經(jīng)常開銷的比例。MATLAB65版（含）之後，已經(jīng)加入了JITACCELERATOR，有效地降低了經(jīng)常開銷的比例，同時(shí)也使得向量化程式碼的重要性越來(lái)越低。JITJUSTINTIME對(duì)向量化運(yùn)算來(lái)進(jìn)行說(shuō)明，主要著眼點(diǎn)在於很多舊的程式碼還是包含向量化的運(yùn)算。MATLAB65之前的版本，向量化的運(yùn)算還是提高執(zhí)行效率的主要關(guān)鍵。某些情況下，JITACCELERATOR的效能無(wú)法完全發(fā)揮，還是需要配合向量化的運(yùn)算才可以讓程式碼更有效率。,MATLAB程式碼的向量化,傳統(tǒng)範(fàn)例I,計(jì)算N項(xiàng)調(diào)和數(shù)列的總和，若使用FOR迴圈，當(dāng)N100,000時(shí)，可計(jì)算其執(zhí)行時(shí)間如下範(fàn)例21HSUM01M,TICN100000TOTAL0FORI1NTOTALTOTAL1/IENDTOC,ELAPSED_TIME263000,由上可知，當(dāng)N100,000時(shí)，使用FOR迴圈的程式碼約需263秒才能執(zhí)行完畢。若改用向量化的運(yùn)算，可計(jì)時(shí)如下範(fàn)例22HSUM02M,TICN100000SEQUENCE1NTOTALSUM1/SEQUENCETOC,ELAPSED_TIME00600,傳統(tǒng)範(fàn)例II,提示,上述的測(cè)試時(shí)間是根據(jù)PENTIUM450,256MBRAM在MATLAB61所得到的結(jié)果。如果你的電腦更快，計(jì)算時(shí)間的差異性就不會(huì)那麼顯著。若使用我目前環(huán)境進(jìn)行測(cè)試HSUM01?0009773SECHSUM02?0008938SEC在MATLAB65版之後，採(cǎi)用了JITJUSTINTIME的編譯技術(shù)，因此「向量化」和「非向量化」的程式執(zhí)行速度差異較小。,要能夠熟練地運(yùn)用向量化的運(yùn)算，有下列三個(gè)訣竅對(duì)矩陣的索引（INDEXING）非常熟悉。對(duì)MATLAB可用的內(nèi)建（BUILTIN）指令非常熟悉。對(duì)問(wèn)題本身的解法非常瞭解。,程式碼向量化的要訣,若要對(duì)一個(gè)矩陣X的每一個(gè)直行乘上向量A的每一個(gè)元素，我們可使用內(nèi)建的DIAG指令來(lái)達(dá)成此功能，例如範(fàn)例23COLMULTIPLY01M此為「以空間換取時(shí)間」的範(fàn)例。,X123111A321YXDIAGA,Y343321,矩陣直行乘法,若要對(duì)一個(gè)矩陣X的每一個(gè)橫列乘上向量A的每一個(gè)元素，也可以使用DIAG指令，例如範(fàn)例24ROWMULTIPLY01M此為「以空間換取時(shí)間」的範(fàn)例。,X123456A31YDIAGAX,Y369456,矩陣橫列乘法,將一個(gè)行向量拷貝N份，可執(zhí)行如下範(fàn)例25COPYCOLMAONES1,N也可以產(chǎn)生相同的效果。,N10A137REPA,ONES1,N,REP111111111133333333337777777777,行向量拷貝,將一個(gè)列向量拷貝N份，可執(zhí)行如下範(fàn)例26COPYROWMONESN,1A也可以產(chǎn)生相同的效果。,N5A148REPAONES1,N,,REP148148148148148,列向量拷貝,將一個(gè)MN矩陣A拷貝MN次,並安排成維度大小為MMNN的大矩陣，可輸入如下範(fàn)例27COPYMATRIXM亦可使用內(nèi)建指令REPMAT來(lái)達(dá)成同樣的功能。,A123456MSIZEA,1NSIZEA,2M3N2MINDEX1MONES1,MNINDEX1NONES1,NBAMINDEX,NINDEX,B123123456456123123456456123123456456,矩陣拷貝,22JIT加速器,65版後，MATLAB引進(jìn)了JITJUSTINTIME的編譯技術(shù)，簡(jiǎn)稱JIT加速器（JITACCELERATOR）大幅地提高了迴圈的執(zhí)行效率。使得一般程式碼的撰寫，可以使用程式設(shè)計(jì)師最能代表演算法精神的方式來(lái)呈現(xiàn)。不必為了提高效率、勉強(qiáng)使用向量化運(yùn)算，造成程式碼維護(hù)的困難。,在JIT加速器的運(yùn)作下，向量化運(yùn)算的光環(huán)不如往日耀眼，測(cè)試範(fàn)例如下範(fàn)例28HSUM03M,JIT加速器可以加速迴圈運(yùn)算,NS10001100FORI1LENGTHNS第一種方法FORLOOPOPERATIONTICTOTAL0FORJ1NSITOTALTOTALSIN1/JENDTIME1TOC,由此圖可以看出，TIME1/TIME2的值都大約在4左右，顯示JIT加速器的效能已經(jīng)大幅提升。（若無(wú)JIT，此值約在40左右。）,第二種方法VECTORIZEDOPERATIONTICTOTALSUMSIN1/1NTIME2TOC計(jì)算倍數(shù)RATIOITIME1/TIME2ENDPLOTNS,RATIO,GRIDONXLABELNYLABELTIME1/TIME2,22JIT加速器,JIT加速器的主要功能,JIT加速器的主要功能，列出如下JIT加速器可將每一列MATLAB程式碼直接轉(zhuǎn)成「原機(jī)指令」，而不再直接使用PCODE，因此省去大部分PCODE會(huì)遇到的經(jīng)常開銷。JIT加速器也會(huì)對(duì)INTELX86為主的WINDOWS及LINUX進(jìn)行程式碼最佳化。JIT加速器可以分成兩部分JUSTINTIMECODEGENERATION將PCODE轉(zhuǎn)成原機(jī)指令。RUNTIMETYPEANALYSIS對(duì)資料型態(tài)的分析，以便用於JUSTINTIMECODEGENERATION,有了JIT加速器，是否還要使用向量化運(yùn)算呢基本的判斷原則如下所用的向量化運(yùn)算是根基於「以空間換取時(shí)間」，就必須經(jīng)過(guò)測(cè)試，才能決定是否值得使用。所用的向量化運(yùn)算並非根基於「以空間換取時(shí)間」，那還是要盡量使用（尤其是要盡量使用內(nèi)建指令），因?yàn)檫@種向量化運(yùn)算不需要額外的空間來(lái)暫存資訊，因此效率會(huì)比JIT加速器更好。碰到較複雜的資料結(jié)構(gòu)（例如稀疏矩陣及維度超過(guò)3維的陣列），JIT加速器並不會(huì)轉(zhuǎn)成原機(jī)指令，因此在這種情況下，也要盡量使用向量化運(yùn)算。,JIT加速器與向量化運(yùn)算,使用「不需要額外空間」的向量化運(yùn)算來(lái)提昇效率範(fàn)例29MATMULTIPLY01M,「不需要額外空間」的向量化運(yùn)算I,NS50120FORJ1LENGTHNSNNSJARANDNBRANDN第一種方法FORLOOPOPERATIONCZEROSNTICFORP1NFORQ1NCP,QAP,QBP,QEND,ENDTIME1JTOC第二種方法VECTORIZEDOPERATIONTICCABTIME2JTOCENDSUBPLOT2,1,1PLOTNS,TIME1,V,NS,TIME2,GRIDONLEGENDTIME1FORFORLOOPCODE,TIME2FORVECTORIZEDCODEXLABELNYLABELSECONDSUBPLOT2,1,2PLOTNS,TIME1/TIME2,GRIDONXLABELNYLABELTIME1/TIME2,「不需要額外空間」的向量化運(yùn)算II,「不需要額外空間」的向量化運(yùn)算III,,第一個(gè)圖畫出迴圈運(yùn)算和向量化運(yùn)算所花的時(shí)間。第二個(gè)圖則是時(shí)間的比值，有此可以看出加速的倍數(shù)大約都在10～20之間。,為了使JIT加速器發(fā)揮最大效能，我們的程式碼也要盡量配合，以下是幾點(diǎn)注意事項(xiàng)迴圈中的索引變數(shù)盡量使用純量。迴圈中的變數(shù)盡量是簡(jiǎn)單的資料型態(tài)，並以不超過(guò)二維為主。迴圈內(nèi)所呼叫的函式僅限於MATLAB內(nèi)建函式。,如何發(fā)揮JIT加速器的最大功能,MATLAB可隨時(shí)改變矩陣的維度，以容納新的資料，但是每一次增大矩陣的容量時(shí)，可能造成下列現(xiàn)象MATLAB必須要向電腦的作業(yè)系統(tǒng)索取記憶體，因此會(huì)造成程式執(zhí)行效率的降低造成記憶體使用的分散現(xiàn)象，可能導(dǎo)致「實(shí)際雖仍有充足的記憶體，但卻未有連續(xù)空間以處理較大矩陣」的現(xiàn)象。基於以上種種原因，若預(yù)先知道使用矩陣的維度大小，則最好實(shí)施矩陣的預(yù)先配置,23矩陣的預(yù)先配置,計(jì)算亂數(shù)矩陣的各個(gè)元素平方範(fàn)例210PREALLOCATE01M,N10000000XRAND1,NCLEARY清除變數(shù)YTICFORI1NYIXI2ENDTIME1TOC,矩陣預(yù)先配置的範(fàn)例I,YZEROS1,N預(yù)先配置所須矩陣YTICFORI1NYIXI2ENDTIME2TOCFPRINTFTIME1G,TIME2G,TIME1/TIME2G\N,TIME1,TIME2,TIME1/TIME2,矩陣預(yù)先配置的範(fàn)例II,TIME1286088,TIME20152491,TIME1/TIME218761,上述範(fàn)例中，使用YZEROS1,N來(lái)事先宣告Y所佔(zhàn)用的空間，因此程式碼的執(zhí)行效率快了將近20倍。,24記憶體管理,MATLAB提供4個(gè)和記憶體管理直接相關(guān)的指令,執(zhí)行MATLAB程式碼時(shí)，發(fā)生「OUTOFMEMORY」的警告訊息，即表示已無(wú)足夠空間可容納新變數(shù)，可採(cǎi)取的解決方案為隨時(shí)使用CLEAR指令清除不再須要的變數(shù)。執(zhí)行PACK指令，減小記憶體分散使用的現(xiàn)象。增加虛擬記憶體的大小。選用適合的資料型態(tài)，例如DOUBLE?UINT8。如果上述方法仍無(wú)法解決，那麼只有減少所用資料量的大小。,「記憶體不足」之對(duì)策,

簡(jiǎn)介：本科畢業(yè)設(shè)計(jì)本科畢業(yè)設(shè)計(jì)論文論文論文題目論文題目基于MATLAB的圖像處理與分析英文IMAGEPROCESSINGANALYSISBASEDONMATLAB所在學(xué)院電子信息學(xué)院專業(yè)班級(jí)電氣工程及其自動(dòng)化學(xué)生姓名學(xué)號(hào)指導(dǎo)教師職稱完成日期年月日ABSTRACTABSTRACTALONGWITHTHERAPIDDEVELOPMENTOFSCIENTIFICTECHNOLOGYITNEEDSHIGHERIMAGEPROCESSINGCONTROLPRECISIONOFTHESYSTEMHIGHERRESPONSESPEEDMEFUNCTIONTHEABILITYTOADJUSTTHEDEMFMOSTCOMPLEXPROCESSINGOBJECTTHEUSEOFTHETRADITIONALPROCESSINGMETHODISOFTENDIFFICULTTORECEIVESATISFACTYRESULTITWILLBEMECONVENIENTQUICKTOUSECOMPUTERTOHLETHESEIMAGESINTHISPAPERWEDOSOMEANALYSISRESEARCHFIMAGEPROCESSINGSYSTEMBASEDONTHEMATLABLANGUAGEWEMAINLYINTRODUCETHEBASICCONCEPTTHEDEVELOPMENTAPPLICATIONACTERISTICSOFIMAGEPROCESSINGWEDOSOMERESEARCHONIMAGEPROCESSINGSYSTEMOFBASICSTRUCTUREPRINCIPLEIMAGEPROCESSINGTHEMANMACHINEINTERFACEDESIGNMETHODOFEXCHANGEWEDISCUSSONIMAGEPROCESSINGSYSTEMINTHEMODERNSOCIETYWEDESIGNTHEIMAGEPROCESSINGSYSTEMTOREALIZETHEFUNCTIONSREADPRESERVETHEIMAGEJOINNOISETOIMAGEDENOISEIMAGEZOOMINREDUCETHEIMAGESOONFINALLYWESIMULATETHEFUNCTIONSTHROUGHMATLABTHESIMULATIONRESULTSSHOWTHATTHEIMAGEPROCESSINGTHROUGHMATLABPRACTICALFEASIBLEKEYKEYWDSWDSMATLABIMAGEPROCESSINGFUNCTIONGUI

簡(jiǎn)介：摘要圖像增強(qiáng)是指依據(jù)圖像所存在的問(wèn)題，按特定的需要突出一幅圖像中的某些信息，同時(shí)，削弱或去除某些冗余信息的處理方法。其主要目的是使得處理后的圖像對(duì)給定的應(yīng)用比原來(lái)的圖像更加有效同時(shí)可以有效的改善圖像質(zhì)量。圖像增強(qiáng)技術(shù)主要包含直方圖修改處理、圖像平滑化處理、圖像尖銳化處理和彩色處理技術(shù)等。圖像增強(qiáng)單純從技術(shù)上可分成兩大類一類是頻域處理法；一類是空域處理法。頻域處理法的采用的是卷積定理，它利用修改圖像傅立葉變換的方法實(shí)現(xiàn)對(duì)圖像的增強(qiáng)處理；空域處理法則是對(duì)圖像中的像素進(jìn)行直接的處理，大多數(shù)是以灰度映射變換為基礎(chǔ)的，所用的映射變換取決于增強(qiáng)的目的。在本論文中，利用MATLAB提供的若干函數(shù)，用于圖像類型的轉(zhuǎn)換。MATLAB支持五種圖像類型，即二值圖像、索引圖像、灰度圖像、RGB圖像和多幀圖像陣列。在MATLAB中要查詢一個(gè)圖像文件的信息，只要用IMFINFO指令加上文件及其完整路徑名即可。MATLAB提供了兩個(gè)重要的用于圖像文件的讀寫的指令，分別是從圖像文件中讀取數(shù)據(jù)的IMREAD，以及將數(shù)據(jù)寫入到圖像文件中的IMWRITE。至于多個(gè)圖像的顯示，則可分為兩個(gè)方面在不同的圖形窗口顯示不同的圖像，可以用FIGURE指令來(lái)實(shí)現(xiàn)；在同一個(gè)圖形窗口顯示多圖，可以用SUBPLOT來(lái)實(shí)現(xiàn)。關(guān)鍵字MATLAB圖像增強(qiáng)灰度濾波器畢業(yè)設(shè)計(jì)（論文）原創(chuàng)性聲明和使用授權(quán)說(shuō)明原創(chuàng)性聲明本人鄭重承諾所呈交的畢業(yè)設(shè)計(jì)（論文），是我個(gè)人在指導(dǎo)教師的指導(dǎo)下進(jìn)行的研究工作及取得的成果。盡我所知，除文中特別加以標(biāo)注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過(guò)的研究成果，也不包含我為獲得及其它教育機(jī)構(gòu)的學(xué)位或?qū)W歷而使用過(guò)的材料。對(duì)本研究提供過(guò)幫助和做出過(guò)貢獻(xiàn)的個(gè)人或集體，均已在文中作了明確的說(shuō)明并表示了謝意。作者簽名日期指導(dǎo)教師簽名日期使用授權(quán)說(shuō)明本人完全了解大學(xué)關(guān)于收集、保存、使用畢業(yè)設(shè)計(jì)（論文）的規(guī)定，即按照學(xué)校要求提交畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版本；學(xué)校有權(quán)保存畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版，并提供目錄檢索與閱覽服務(wù)；學(xué)校可以采用影印、縮印、數(shù)字化或其它復(fù)制手段保存論文；在不以贏利為目的前提下，學(xué)?？梢怨颊撐牡牟糠只蛉?jī)?nèi)容。作者簽名日期

簡(jiǎn)介：基于MATLAB的PID參數(shù)調(diào)整方法的仿真研究院系系物理與電氣工程學(xué)院物理與電氣工程學(xué)院專業(yè)業(yè)自動(dòng)化自動(dòng)化小組成員小組成員指導(dǎo)老師指導(dǎo)老師完成時(shí)間完成時(shí)間工程整定法。它不需要得到過(guò)程模型，主要依賴工程經(jīng)驗(yàn)，在控制系統(tǒng)的試驗(yàn)中直接進(jìn)行參數(shù)整定。方法簡(jiǎn)單實(shí)用，計(jì)算簡(jiǎn)便且易于掌握，可以解決一般的實(shí)際問(wèn)題，在工程實(shí)際中被廣泛采用。PID控制器參數(shù)的工程整定法，主要有臨界比例度法（又稱穩(wěn)定邊界法）、反應(yīng)曲線法和41衰減法。其共同點(diǎn)都是通過(guò)實(shí)驗(yàn)，然后按照工程經(jīng)驗(yàn)公式對(duì)控制器參數(shù)進(jìn)行整定。然而，無(wú)論采用哪一種方法整定所得到的控制參數(shù)，都需要在實(shí)際運(yùn)行中進(jìn)行最后的調(diào)整與完善。理論和實(shí)踐證明，即便是整定得很好的PID參數(shù)值，系統(tǒng)響應(yīng)的快速性與超調(diào)量之間也存在矛盾，二者不可能同時(shí)達(dá)到最優(yōu)，且系統(tǒng)在跟蹤設(shè)定值與抑制擾動(dòng)方面對(duì)控制參數(shù)的要求也是矛盾的。下面從系統(tǒng)穩(wěn)定性、響應(yīng)速度、超調(diào)量和控制精度等各方面特性來(lái)分析PID三參數(shù)對(duì)PID控制品質(zhì)的影響。比例系數(shù)的作用在于加快系統(tǒng)的響應(yīng)速度，提高系統(tǒng)調(diào)節(jié)精度。越PKPK大，系統(tǒng)的響應(yīng)速度越快，但將產(chǎn)生超調(diào)和振蕩，甚至導(dǎo)致系統(tǒng)不穩(wěn)定，因此值不能取的過(guò)大；如果值較小，則會(huì)降低調(diào)節(jié)精度，使響應(yīng)速度變慢，從PKPK而延長(zhǎng)調(diào)節(jié)時(shí)間，使系統(tǒng)動(dòng)、靜態(tài)特性變壞。積分環(huán)節(jié)作用系數(shù)的作用在于消除系統(tǒng)的穩(wěn)態(tài)誤差。越大，積分速度越ITIT快，系統(tǒng)靜差消除越快。但過(guò)大，在響應(yīng)過(guò)程的初期以及系統(tǒng)在過(guò)渡過(guò)程中，IT會(huì)產(chǎn)生積分飽和現(xiàn)象，從而引起響應(yīng)過(guò)程出現(xiàn)較大的超調(diào)，使動(dòng)態(tài)性能變差。若過(guò)小，積分作用變?nèi)?，則系統(tǒng)的靜差難以消除，過(guò)渡過(guò)程時(shí)間加長(zhǎng)，不能較快IT的達(dá)到穩(wěn)定狀態(tài)，影響系統(tǒng)的調(diào)節(jié)精度和動(dòng)態(tài)特性。微分環(huán)節(jié)作用系數(shù)的作用在于改善系統(tǒng)的動(dòng)態(tài)特性。因?yàn)镻ID控制器的微DT分環(huán)節(jié)只影響系統(tǒng)偏差的變化率，其作用主要是在響應(yīng)過(guò)程中抑制偏差向任何方向的變化，對(duì)偏差變化進(jìn)行提前制動(dòng)，降低超調(diào)，增加系統(tǒng)的穩(wěn)定性。但過(guò)DT大，則會(huì)使響應(yīng)過(guò)程過(guò)分提前制動(dòng)，從而拖長(zhǎng)調(diào)節(jié)時(shí)間，而且系統(tǒng)的抗干擾性也會(huì)變差。

簡(jiǎn)介：,基于CMATLAB的圖像識(shí)別技術(shù)仿真研究,姓名周易德專業(yè)電子信息工程學(xué)號(hào)2012214709班級(jí)0121208,,一個(gè)完整的圖像識(shí)別系統(tǒng)包含三個(gè)部分圖像預(yù)處理圖像到圖像特征提取圖像到特征向量模式識(shí)別特征向量到抽象含義,圖像識(shí)別是針對(duì)圖像的模式識(shí)別。圖像識(shí)別的目的是使得計(jì)算機(jī)擁有像人類一樣能從圖片中認(rèn)識(shí)事物分析事物的能力。,圖像預(yù)處理的目的就是將輸入圖像處理成輸出圖像，從而達(dá)到圖像消噪、圖像增強(qiáng)和圖像恢復(fù)等目的。需要注意的是我在這里將圖像分割也歸為預(yù)處理的一部分，圖像分割輸出的是原始圖片的子集合。,原始圖片,按行分割,分割結(jié)果,特征提取是對(duì)圖像中的目標(biāo)采取合適的數(shù)據(jù)結(jié)構(gòu)來(lái)進(jìn)行表達(dá)。圖像識(shí)別中的特征提取就是一種圖像表達(dá)技術(shù)。一般需要通過(guò)對(duì)識(shí)別對(duì)象的深入研究經(jīng)驗(yàn)，選取特征；有時(shí)也可以通過(guò)圖像壓縮技術(shù)對(duì)圖像進(jìn)行特征提取。,模式識(shí)別是圖像識(shí)別的核心技術(shù)。模式廣義的來(lái)講就是一種完美的概念，是從某類事物中歸納總結(jié)出來(lái)的一種認(rèn)識(shí)。模式識(shí)別是讓計(jì)算機(jī)像人腦一樣，從千變?nèi)f化的圖像中識(shí)別出抽象的模式。模式識(shí)別的核心是分類器設(shè)計(jì)。,基于最小錯(cuò)誤的貝葉斯決策,則,概率神經(jīng)網(wǎng)絡(luò)（PNN）是一種徑向基神經(jīng)元實(shí)現(xiàn)最小錯(cuò)誤的貝葉斯決策的方法。,概率神經(jīng)網(wǎng)絡(luò),概率神經(jīng)網(wǎng)絡(luò),概率神經(jīng)網(wǎng)絡(luò),,,D3使用了MATLAB2015B引入的NEURALNETWORKTOOLBOX。使用了兩層AUTOENCODER來(lái)壓縮數(shù)據(jù)，一層SOFTMAX來(lái)計(jì)算后驗(yàn)概率。,D3結(jié)果統(tǒng)計(jì),完整圖像識(shí)別系統(tǒng)的流圖,謝謝答辯組老師,

簡(jiǎn)介：46異步電動(dòng)機(jī)變頻調(diào)速系統(tǒng)與MATLAB實(shí)現(xiàn),461異步電動(dòng)機(jī)變頻調(diào)速系統(tǒng)介紹,交流異步電動(dòng)機(jī)所采用的變頻控制方式有V/F控制、矢量控制以及直接轉(zhuǎn)矩控制等實(shí)現(xiàn)方式。,462變頻調(diào)速控制方式,三相異步電動(dòng)機(jī)定子每相電動(dòng)勢(shì)的有效值為EG444F1N1RN1ΦM（420）式中EG氣隙磁通在定子每相中感應(yīng)電動(dòng)勢(shì)的有效值，單位為V；F1定子頻率（HZ）；N1定子每相繞組串聯(lián)匝數(shù)；RN1與繞組結(jié)構(gòu)有關(guān)的常數(shù)；ΦM每極氣隙磁通量WB。,,1．基頻以下變頻調(diào)速控制基本方式,2．基頻以下變頻調(diào)速控制基本方式,,,463矢量控制變頻調(diào)速系統(tǒng),464交流異步電動(dòng)機(jī)變頻矢量控制系統(tǒng),,465矢量控制變頻調(diào)速仿真,,1,CALCULATION子模塊,（2）DQ－ABC子模塊,,（3）TEATCALCULATION子模塊,（4）ＡＢＣＤＱ子模塊,,（5）,CALCULATION子模塊,（6）CURRENTREGULATOR子模塊,,（7）ASR子模塊,異步電機(jī)參數(shù)線電壓380；額定頻率5０ＨＺ；定子內(nèi)阻0087Ω定子漏感08ＭＨ；轉(zhuǎn)子內(nèi)阻0228Ω；轉(zhuǎn)子漏感08ＭＨ；定、轉(zhuǎn)子漏感347；極對(duì)數(shù)為4。為逆變器直流電源VDC設(shè)為780V。仿真方法選擇為固定步長(zhǎng)（FIXEDSTEP）。仿真時(shí)間設(shè)為0－315秒。,,電機(jī)轉(zhuǎn)速,,電磁轉(zhuǎn)矩,,三相定子電流,

簡(jiǎn)介：信號(hào)與系統(tǒng),SIGNALSANDSYSTEMS,國(guó)家精品課程主教材、北京市精品教材信號(hào)與系統(tǒng)第2版陳后金，胡健，薛健清華大學(xué)出版社，2005年,利用MATLAB進(jìn)行信號(hào)與系統(tǒng)分析,,MATLAB簡(jiǎn)介信號(hào)的MATLAB表示利用MATLAB進(jìn)行系統(tǒng)的時(shí)域分析利用MATLAB進(jìn)行信號(hào)的頻域分析利用MATLAB分析系統(tǒng)的頻率特性利用MATLAB進(jìn)行連續(xù)系統(tǒng)的S域分析利用MATLAB進(jìn)行離散系統(tǒng)的Z域分析利用MATLAB進(jìn)行系統(tǒng)的狀態(tài)變量分析,MATLAB簡(jiǎn)介,,MATRIXLABORATORY,MATLAB的工作方式如何獲取幫助表達(dá)式變量、數(shù)值、算數(shù)運(yùn)算符、關(guān)系運(yùn)算符、邏輯運(yùn)算符、冒號(hào)運(yùn)算符數(shù)組及其運(yùn)算函數(shù)文件循環(huán)（FOR、WHILE循環(huán)）基本繪圖語(yǔ)句,一、MATLAB的工作方式,（1）窗口命令方式（2）運(yùn)行以M為擴(kuò)展名磁盤文件,工作方式舉例,用PLOT函數(shù)畫一個(gè)方波T100113F001100PLOTT,FXLABELTYLABELFTAXIS1302,直接在命令窗口輸入以上命令,建一個(gè)名字為MY_FILEM的文件，然后在命令窗口輸入文件名及回車。,二、獲取幫助,命令窗口輸入HELP函數(shù)名例如HELPPLOT,三、表達(dá)式,不需要變量的類型說(shuō)明變量名的第一個(gè)字符必須是字母變量名長(zhǎng)度不超過(guò)31個(gè)字符大寫和小寫的字母視為不同的字符例如NUM_STUDENTS25特殊變量PI表示圓周率，INF表示無(wú)窮大，NANNOTANUMBER表示不定量，如0/0。,變量,三、表達(dá)式,數(shù)值,MATLAB用常規(guī)的十進(jìn)制表示數(shù)值用I或J作為后綴來(lái)表示復(fù)數(shù)的虛部例1235E5表示1235?105，X23JABSX求復(fù)數(shù)X的模ANGLEX求復(fù)數(shù)X的相角弧度REALX求復(fù)數(shù)X的實(shí)部IMAGX求復(fù)數(shù)X的虛部CONJX求復(fù)數(shù)X的共軛,三、表達(dá)式,運(yùn)算符號(hào),算數(shù)運(yùn)算符,加減乘/除乘方‘矩陣的復(fù)共軛轉(zhuǎn)置,三、表達(dá)式,運(yùn)算符號(hào),邏輯運(yùn)算符,A34。,用LINSPACE產(chǎn)生數(shù)組,四、數(shù)組,1數(shù)組的構(gòu)造,MATLAB提供了一些產(chǎn)生基本矩陣的函數(shù)ZEROS產(chǎn)生矩陣元素全為0的矩陣ONES產(chǎn)生矩陣元素全為1的矩陣RAND產(chǎn)生0,1均勻分布隨機(jī)數(shù)矩陣RANDN產(chǎn)生正態(tài)分布隨機(jī)數(shù)矩陣,四、數(shù)組,2數(shù)組的運(yùn)算,數(shù)組和一個(gè)標(biāo)量相加或相乘例YX1Z3X2個(gè)數(shù)組的對(duì)應(yīng)元素相乘除/例ZXY確定數(shù)組大小的函數(shù)SIZEA返回值數(shù)組A的行數(shù)和列數(shù)（二維）LENGTHB確定數(shù)組B的元素個(gè)數(shù)（一維）,五、函數(shù)文件,M文件的第一行包含F(xiàn)UNCTION功能建立一個(gè)函數(shù)，可以同MATLAB的庫(kù)函數(shù)一樣使用。,五、函數(shù)文件,例編一個(gè)繪制圖示波形的函數(shù)。,FUNCTIONYTRITYABSTEPSSS1/NNNN1END,例計(jì)算的值，且誤差小于106,FPRINTFS5F\N,S,八、PLOT函數(shù)繪圖函數(shù)CONTINUOUS,TLINSPACE0,4PI,512PLOTT,SINT,T,COST,TITLEMYFIGUREXLABELTLEGENDSINT,COST,,八、PLOT函數(shù)繪圖函數(shù)CONTINUOUS,九、STEM函數(shù)繪圖函數(shù)DISCRETE,K039STEMK,COS09PIKTITLECOS09\PIK,,九、STEM函數(shù)繪圖函數(shù)DISCRETE,COS09?K波形,信號(hào)的MATLAB表示,,基本信號(hào)的MATLAB表示指數(shù)信號(hào)AEAT、指數(shù)序列AK、抽樣函數(shù)SAT、正弦型信號(hào)、矩形脈沖信號(hào)、三角脈沖信號(hào)信號(hào)基本運(yùn)算的MATLAB實(shí)現(xiàn)尺度變換、翻轉(zhuǎn)、時(shí)移、相加、相乘、差分與求和、微分與積分,一、基本信號(hào)的MATLAB表示,指數(shù)信號(hào)AEATYAEXPAT指數(shù)序列AK冪運(yùn)算AK實(shí)現(xiàn)正弦型信號(hào)內(nèi)部函數(shù)COS和SIN抽樣函數(shù)SATSINCT矩形脈沖信號(hào)YRECTPULST,WIDTH三角波脈沖信號(hào)YTRIPULST,WIDTH,SKEW,一、基本信號(hào)的MATLAB表示,,?CAYINGEXPONENTIAL,T000110A1A04FTAEXPATPLOTT,FT,T00110A1A04FTAEXPATSTEMT,FT,一、基本信號(hào)的MATLAB表示,,RECTPULS,T000014T1FTRECTPULST2T,TPLOTT,FTAXIS0,4,05,15,一、基本信號(hào)的MATLAB表示,,TRIPULS,T300013FTTRIPULST,4,05PLOTT,FT,FTTRIPULST,4,1,一、基本信號(hào)的MATLAB表示,,UNITIMPULSSEQUENCE,K5050DELTAZEROS1,50,1,ZEROS1,50STEMK,DELTA,FUNCTIONF,KIMPSEQK0,K1,K2產(chǎn)生FKDELTAKK0；K11,WLINSPACE6PI,6PI,512NLENGTHWFZEROS1,NFORK1NFKQUAD8SF1,1,1,,,WKENDFIGURE1PLOTW,REALFTITLEXLABEL\OMEGAYLABELFJ\OMEGAFIGURE2PLOTW,REALFSINCW/2/PI2XLABEL\OMEGATITLE計(jì)算誤差,例3試用數(shù)值方法近似計(jì)算三角波信號(hào)的頻譜,運(yùn)行結(jié)果,利用MATLAB進(jìn)行系統(tǒng)頻域分析,,連續(xù)系統(tǒng)頻響特性的計(jì)算周期信號(hào)通過(guò)系統(tǒng)的響應(yīng)離散系統(tǒng)頻響特性的計(jì)算,一、連續(xù)系統(tǒng)頻響特性的計(jì)算,B分子多項(xiàng)式系數(shù)A分母多項(xiàng)式系數(shù)W需計(jì)算的HJW的抽樣點(diǎn)數(shù)組W中少需包含兩個(gè)W的抽樣點(diǎn)。,計(jì)算頻響的MATLAB函數(shù),HFREQSB,A,W,,一、連續(xù)系統(tǒng)頻響特性的計(jì)算,例1三階歸一化的BUTTERWORTH低通濾波器的系統(tǒng)函數(shù)為,WLINSPACE0,5,200B1A1221HFREQSB,A,WSUBPLOT2,1,1PLOTW,ABSHSUBPLOT2,1,2PLOTW,ANGLEH,試畫出|HJW|和?W。,,一、連續(xù)系統(tǒng)頻響特性的計(jì)算,三階BUTTERWORTH低通濾波器的幅度響應(yīng)和相位響應(yīng),,二、周期信號(hào)通過(guò)系統(tǒng)的響應(yīng),例2周期方波通過(guò)RC系統(tǒng)的響應(yīng)。,,二、周期信號(hào)通過(guò)系統(tǒng)的響應(yīng),例2周期方波通過(guò)RC系統(tǒng)的響應(yīng)。,P5_2PERIODICSIGNALPASSLTISYSTEM,T4W02PI/TRC01T60016N51C005XNC0ONES1,LENGTHT?FORN12NEVENHARMONICSAREZEROHABS1/1JRCW0NPHIANGLE1/1JRCW0NXNXNHCOSW0NTPHISINCN05ENDPLOTT,XNXLABELTIMERC,NUM2STRRCGRIDSETGCA,XTICK,5310135,,二、周期信號(hào)通過(guò)系統(tǒng)的響應(yīng),例2周期方波通過(guò)RC系統(tǒng)的響應(yīng)。,三、離散系統(tǒng)頻率響應(yīng)的計(jì)算,計(jì)算頻率響應(yīng)的MATLAB函數(shù),B分子的系數(shù)A分母系數(shù),W抽樣的頻率點(diǎn)至少2點(diǎn)，W在02P之間,幅頻特性ABS,相頻特性ANGLE,HFREQZB,A,W,三、離散系統(tǒng)頻響特性的計(jì)算,,B1,A1109A2109,WLINSPACE0,2PI,512,H1FREQZB,A1,W,H2FREQZB,A2,W,PLOTW/PI,ABSH1,W/PI,ABSH2,,LEGEND\ALPHA09,\ALPHA09,三、離散系統(tǒng)頻響特性的計(jì)算,,利用MATLAB進(jìn)行連續(xù)系統(tǒng)的復(fù)頻域分析,,部分分式展開的MATLAB實(shí)現(xiàn)HS的零極點(diǎn)與系統(tǒng)特性的MATLAB計(jì)算,一、部分分式展開的MATLAB實(shí)現(xiàn),R,P,KRESIDUENUM,DEN,NUM,DEN分別為FS分子多項(xiàng)式和分母多項(xiàng)式的系數(shù)向量。,R為部分分式的系數(shù)，P為極點(diǎn)，K為多項(xiàng)式的系數(shù)。若為真分式，則K為零。,二、HS的零極點(diǎn)與系統(tǒng)特性的MATLAB計(jì)算,計(jì)算多項(xiàng)式根ROOTS的函數(shù)可用于計(jì)算HS的零極點(diǎn)。,RROOTSN計(jì)算多項(xiàng)式N的根,HS零極點(diǎn)分布圖可用PZMAP函數(shù)畫出，調(diào)用形式為,PZMAPSYS,表示畫出SYS所描述系統(tǒng)的零極點(diǎn)圖。,例1用部分分式展開法求FS的反變換。,PROGRAM7_1FORMATRAT將結(jié)果數(shù)據(jù)以分?jǐn)?shù)的形式輸出NUM12DEN1430R,PRESIDUENUM,DEN,運(yùn)行結(jié)果為R1/6,1/2,2/3P3,1,0,故FS可展開為,例2用部分分式展開法求FS的反變換。,PROGRAM7_2NUM2305DENCONV11,112將因子相乘的形式轉(zhuǎn)換成多項(xiàng)式的形式R,P,KRESIDUENUM,DENMAGRABSR求R的模ANGRANGLER求R的相角,例2用部分分式展開法求FS的反變換。,運(yùn)行結(jié)果為,R2000011339I,2000011339I,30000P0500013229I,0500013229I,10000K2MAGR2299,22991,30000ANGR26258,26258,0,故FS可展開為,例3試畫出系統(tǒng)的零極點(diǎn)分布圖，求其單位沖激響應(yīng)HT和頻率響應(yīng)HJ?，并判斷系統(tǒng)是否穩(wěn)定。,NUM1DEN1221SYSTFNUM,DENPOLESROOTSDENFIGURE1PZMAPSYST000210HIMPULSENUM,DEN,TFIGURE2PLOTT,HTITLEIMPULSERESPONEH,WFREQSNUM,DENFIGURE3PLOTW,ABSHXLABEL\OMEGATITLEMAGNITUDERESPONE,運(yùn)行結(jié)果,利用MATLAB進(jìn)行離散系統(tǒng)的Z域分析,,部分分式展開的MATLAB實(shí)現(xiàn)HZ的零極點(diǎn)與系統(tǒng)特性的MATLAB計(jì)算,一、部分分式展開的MATLAB實(shí)現(xiàn),R,P,KRESIDUEZNUM,DEN,NUM,DEN分別為FZ分子多項(xiàng)式和分母多項(xiàng)式的系數(shù)向量。,R為部分分式的系數(shù)，P為極點(diǎn)，K為多項(xiàng)式的系數(shù)。若為真分式，則K為零。,二、HZ的零極點(diǎn)與系統(tǒng)特性的MATLAB計(jì)算,利用TF2ZP函數(shù)計(jì)算HZ的零極點(diǎn)，調(diào)用形式為,Z,P,KTF2ZPB,A,HZ零極點(diǎn)分布圖可用ZPLANE函數(shù)畫出，調(diào)用形式為,ZPLANEB,A,B和A分別為HZ分子多項(xiàng)式和分母多項(xiàng)式的系數(shù)向量。,返回值Z為零點(diǎn)、P為極點(diǎn)、K為增益常數(shù)。,例1將FZ用部分分式展開。,PROGRAM8_1NUM18DEN18341R,P,KRESIDUEZNUM,DEN,運(yùn)行結(jié)果為R03600,02400,04000P05000,03333,03333K,故FZ可展開為,例2試畫出系統(tǒng)的零極點(diǎn)分布圖，求其單位沖激響應(yīng)HK和頻率響應(yīng)HEJW。,PROGRAM8_2B121A105000503FIGURE1ZPLANEB,ANUM0121DEN105000503HIMPZNUM,DENFIGURE2STEMHXLABELKTITLEIMPULSERESPONEH,WFREQZNUM,DENFIGURE3PLOTW/PI,ABSHXLABELFREQUENCY\OMEGATITLEMAGNITUDERESPONE,運(yùn)行結(jié)果,利用MATLAB進(jìn)行系統(tǒng)狀態(tài)變量分析,,微分方程到狀態(tài)方程的轉(zhuǎn)換狀態(tài)方程系統(tǒng)函數(shù)矩陣HS的計(jì)算MATLAB求解連續(xù)系統(tǒng)的狀態(tài)方程MATLAB求解離散系統(tǒng)的狀態(tài)方程,一、微分方程到狀態(tài)方程的轉(zhuǎn)換,NUM,DEN分別表示系統(tǒng)函數(shù)HS的分子和分母多項(xiàng)式；A,B,C,D分別為狀態(tài)方程的矩陣。,A,B,C,DTF2SSNUM,DEN,二、狀態(tài)方程系統(tǒng)函數(shù)矩陣HS的計(jì)算,A,B,C,D分別表示狀態(tài)方程的矩陣。K表示函數(shù)SS2TF計(jì)算的與第K個(gè)輸入相關(guān)的系統(tǒng)函數(shù)，即HS的第K列。NUM表示HS第K列的M個(gè)元素的分子多項(xiàng)式DEN表示HS公共的分母多項(xiàng)式。,NUM,DENSS2TFA,B,C,D,K,,三、MATLAB求解連續(xù)系統(tǒng)的狀態(tài)方程,SYS由函數(shù)SS構(gòu)造的狀態(tài)方程模型T需計(jì)算的輸出樣本點(diǎn),T0DTTFINALF,K系統(tǒng)第K個(gè)輸入在T上的抽樣值X0系統(tǒng)的初始狀態(tài)可缺省Y,K系統(tǒng)的第K個(gè)輸出TO實(shí)際計(jì)算時(shí)所用的樣本點(diǎn)；X系統(tǒng)的狀態(tài),獲得連續(xù)系統(tǒng)狀態(tài)方程的計(jì)算機(jī)表示模型,SYSSSA,B,C,D,求解狀態(tài)方程,Y,TO,XLSIMSYS,F,T,X0,,四、MATLAB求解離散系統(tǒng)的狀態(tài)方程,SYS由函數(shù)SS構(gòu)造的狀態(tài)方程模型F,K系統(tǒng)第K個(gè)輸入序列X0系統(tǒng)的初始狀態(tài)可缺省Y,K系統(tǒng)的第K個(gè)輸出N序列的下標(biāo)X系統(tǒng)的狀態(tài),獲得離散狀態(tài)方程的計(jì)算機(jī)表示模型,SYSSSA,B,C,D,,求解狀態(tài)方程,Y,N,XLSIMSYS,F,,X0,例1寫出系統(tǒng)的狀態(tài)方程。,由A,B,C,DTF2SS1,1510可得,所以系統(tǒng)的狀態(tài)方程為,C0D0,例2已知某連續(xù)系統(tǒng)的狀態(tài)方程和輸出方程為,其初始狀態(tài)和輸入分別為,求該系統(tǒng)的系統(tǒng)函數(shù)矩陣HS和輸出。,計(jì)算系統(tǒng)函數(shù)矩陣HS,A2301B0110C1101D1010NUM1,DEN1SS2TFA,B,C,D,1NUM2,DEN2SS2TFA,B,C,D,2運(yùn)行結(jié)果NUM1101120DEN1112NUM2011000DEN2112所以系統(tǒng)函數(shù)矩陣HS為,計(jì)算輸出,PROGRAM9_1A2301B0110C1101D1010X021DT001T0DT2F,1ONESLENGTHT,1F,2EXP3TSYSSSA,B,C,DYLSIMSYS,F,T,X0SUBPLOT2,1,1PLOTT,Y,1,RYLABELY1TXLABELTSUBPLOT2,1,2PLOTT,Y,2YLABELY2TXLABELT,運(yùn)行結(jié)果,例3已知某離散系統(tǒng)的狀態(tài)方程和輸出方程為,其初始狀態(tài)和輸入分別為,求該系統(tǒng)的輸出。,計(jì)算輸出,PROGRAM9_2A0123B01C1121DZEROS2,1X011N10FONES1,NSYSSSA,B,C,D,YLSIMSYS,F,,X0SUBPLOT2,1,1Y1Y,1STEM0N1,Y1XLABELKYLABELY1SUBPLOT2,1,2Y2Y,2STEM0N1,Y2XLABELKYLABELY2,運(yùn)行結(jié)果,