簡介：計量地理學（徐建華，高等教育出版社，2005）配套實習指導13919利用MATLAB編程計算最短路徑及中位點選址1、最短路問題兩個指定頂點之間的最短路徑。例如，給出了一個連接若干個城鎮(zhèn)的鐵路網(wǎng)絡，在這個網(wǎng)絡的兩個指定城鎮(zhèn)間，找一條最短鐵路線。以各城鎮(zhèn)為圖的頂點，兩城鎮(zhèn)間的直通鐵路為圖相應兩頂點間的邊，GG得圖。對的每一邊，賦以一個實數(shù)直通鐵路的長度，稱為的權，GGEEWE得到賦權圖。的子圖的權是指子圖的各邊的權和。問題就是求賦權圖中GGG指定的兩個頂點間的具最小權的軌。這條軌叫做間的最短路，它的00VU00VU權叫做間的距離，亦記作。00VU00VUD求最短路已有成熟的算法迪克斯特拉（DIJKSTRA）算法，其基本思想是按距從近到遠為順序，依次求得到的各頂點的最短路和距離，直至（或0U0UG0V直至的所有頂點），算法結束。為避免重復并保留每一步的計算信息，采用了G標號算法。下面是該算法。I令，對，令，，。00UL0UVVL00US0III對每個（），用ISVIISVSMINUVWULVLISU代替。計算，把達到這個最小值的一個頂點記為，令VLMINVLISV1IU計量地理學（徐建華，高等教育出版社，2005）配套實習指導141存放由始點到第點最短通路的值。IDI求第一個城市到其它城市的最短路徑的MATLAB程序如下CLEARCLCM10000A1050M402510A2ZEROS121520M25A3ZEROS131020MA4ZEROS141025A5ZEROS1555A6ZEROS16AAAPB1LENGTHA0PB11D1LENGTHAMD10TEMP1WHILESUMPBLENGTHATBFINDPB0DTBMINDTBDTEMPATEMPTBTMPBFINDDTBMINDTBTEMPTBTMPB1PBTEMP1ENDD運行輸出，第一個城市到其它城市的最短路徑長度，即D035453525102、選址問題－以中位點選址為例中位點選址問題的質量判據(jù)為使最佳選址為止所在的定點到網(wǎng)絡圖中其他頂點的最短路徑距離的總和（或者以各個頂點的載荷加權求和）達到最小。例2某縣下屬七個鄉(xiāng)鎮(zhèn)，各鄉(xiāng)鎮(zhèn)所擁有的人口數(shù)AVI（I1，2，，7），以及各鄉(xiāng)鎮(zhèn)之間的距離WIJ（I，J1，2，，7）如圖所示?，F(xiàn)在需要設立一個中心郵局，為全縣所轄的七個鄉(xiāng)鎮(zhèn)共同服務。試問該中心郵局應該設在哪一個

簡介：VISUALC是當前主流的應用程序開發(fā)環(huán)境之一，開發(fā)環(huán)境強大，開發(fā)的程序執(zhí)行速度快。但在科學計算方面函數(shù)庫顯得不夠豐富、讀取、顯示數(shù)據(jù)圖形不方便。MATLAB是一款將數(shù)值分析、矩陣計算、信號處理和圖形顯示結合在一起，包含大量高度集成的函數(shù)可供調用，適合科學研究、工程設計等眾多學科領域使用的一種簡潔、高效的編程工具。不過由于MATLAB使用的是解釋性語言，大大限制了它的執(zhí)行速度和應用場合?；赩C和MATLAB混合編程是很多熟悉VC編程而又需要進行科學計算、數(shù)據(jù)仿真的科研人員常用的一種方式，其中最簡單也最直接的方法就是調用MATLAB引擎。本文以下部分將詳細介紹通過VC60調用MATLAB65引擎來達到VC與MATLAB數(shù)據(jù)共享編程的方法。1什么是MATLAB引擎所謂MATLAB引擎（ENGINE），是指一組MATLAB提供的接口函數(shù)，支持CC、FTRAN等語言，通過這些接口函數(shù)，用戶可以在其它編程環(huán)境中實現(xiàn)對MATLAB的控制?？梢灾饕δ苡小锎蜷_關閉一個MATLAB對話；★向MATLAB環(huán)境發(fā)送命令字符串；★從MATLAB環(huán)境中讀取數(shù)據(jù)；★向MATLAB環(huán)境中寫入數(shù)據(jù)。與其它各種接口相比，引擎所提供的MATLAB功能支持是最全面的。通過引擎方式，應用程序會打開一個新的MATLAB進程，可以控制它完成任何計算和繪圖操作。對所有的數(shù)據(jù)結構提供100的支持。同時，引擎方式打開的MATLAB進程會在任務欄顯示自己的圖標，打開該窗口，可以觀察主程序通過ENGINE方式控制MATLAB運行的流程，并可在其中輸入任何MATLAB命令。實際上，通過引擎方式建立的對話，是將MATLAB以ACTIVEX控件方式啟動的。在MATLAB初次安裝時，會自動執(zhí)行一次MATLABREGSERVER將自己在系統(tǒng)的控件庫中注冊。如果因為特殊原因，無法打開MATLAB引擎，可以在DOS命令提示符后執(zhí)行上述命令，重新注冊。2配置編譯器要在VC中成功編譯MATLAB引擎程序，必須包含引擎頭文件ENGINEH并引入MATLAB對應的庫文件LIBMXLIB、LIBMATLIB、LIBENGLIB。具體的說，打開一個工程后，做如下設置（以VC6為例）1通過菜單工程選項，打開設置屬性頁，進入DIRECTIES頁面，在目錄下拉列表框中選擇INCLUDEFILES，添加路徑“CMATLABEXTERNINCLUDE“（假定MATLAB安裝在CMATLAB目錄）。MESSAGEBOX“CANTSTARTMATLABENGINE“EXIT1ENGCLOSEEP關閉MATLAB引擎。32向MATLAB發(fā)送命令字符串ENGEVALSTRING－發(fā)送命令讓MATLAB執(zhí)行。函數(shù)聲明INTENGEVALSTRINGENGINEEPCONSTSTRING參數(shù)EP為函數(shù)ENGOPEN返回的引擎指針，字符串STRING為要MATLAB執(zhí)行的命令。函數(shù)返回值為0表示成功執(zhí)行，返回1說明執(zhí)行失敗（如命令不能被MATLAB正確解釋或MATLAB引擎已經(jīng)關閉了）。33獲取MATLAB命令窗口的輸出要在VC中獲得函數(shù)ENGEVALSTRING發(fā)送的命令字符串被MATLAB執(zhí)行后在MATLAB窗口中的輸出，可以調用ENGOUTPUTBUFFER函數(shù)。函數(shù)聲明INTENGOUTPUTBUFFERENGINEEPPINTN參數(shù)EP為MATLAB引擎指針，P為用來保存輸出結構的緩沖區(qū)，N為最大保存的字符個數(shù)，通常就是緩沖區(qū)P的大小。該函數(shù)執(zhí)行后，接下來的ENGEVALSTRING函數(shù)所引起的命令行輸出結果會在緩沖區(qū)P中保存。如果要停止保存，只需調用代碼ENGOUTPUTBUFFEREPNULL0。34讀寫MATLAB數(shù)據(jù)341從MATLAB引擎工作空間中獲取變量。MXARRAYENGGETVARIABLEENGINEEPCONSTNAME參數(shù)EP為打開的MATLAB引擎指針，NAME為以字符串形式指定的數(shù)組名。函數(shù)返回值是指向NAME數(shù)組的指針，類型為MXARRAY（MXARRAY數(shù)據(jù)類型在本文第4節(jié)詳細簡介）。

簡介：一、實驗目的一、實驗目的1自行編寫功能函數(shù)實現(xiàn)DTFT的變化，并運用該函數(shù)實現(xiàn)之后的功能。2利用編寫的DTFT函數(shù)，觀察函數(shù)的離散傅里葉變換，并在三維空間中觀察函數(shù)的頻響特性。3利用MATLAB設計奇對稱系統(tǒng)觀察群延時對輸入信號的影響。二、實驗條件二、實驗條件PC機，MATLAB703、實驗內容實驗內容1）編寫功能函數(shù)MYDTFT。根據(jù)公式，將HN用樣值向量H和位置向量N描述，NSNJFFENHH2則DTFT可以用一個向量和一個矩陣的相乘來實現(xiàn)。函數(shù)如下FUNCTIONHMYDTFTHOMGALLENGTHHN0L1MEXP1INOMGAHHM2）做出DTFT，由于MATLAB無法計算連續(xù)變量，只能在某個區(qū)間范圍內將OMGA賦值為很密的、長度很長的向量來近似為連續(xù)變量。由可知，當FS夠高時，HN可近似于HT。FETHHTJ2FS50T2A5P2PIFS40OMGAPI8PPI8SZEROS1FSAONES1FLOTFS1構造門函數(shù)SMYDTFTSOMGAPLOTOMGAABSSGRIDONFIGUREPLOT3OMGAREALSIMAGSGRIDONFIGUREP2PIFS40OMGA02PIPPI22SMYDTFTSOMGAPLOTOMGAABSSFIGURESDSSSSMYDTFTSDOMGASTEMOMGAABSSGRIDONFIGURE通過三維空間將幅頻特性和相頻特性同時放在實軸和虛軸上觀察其變化。XY軸的圖像實部頻譜函數(shù)圖

簡介：2012大學生物理實驗研究論文利用利用MATLAB和EXCEL進行光電效應測普朗克常量實驗中的數(shù)據(jù)處理進行光電效應測普朗克常量實驗中的數(shù)據(jù)處理東南大學東南大學學院學院南京南京211189摘要摘要本文回顧了光電效應中用反向遏制電壓法測量普朗克常數(shù)的實驗原理，討論了可能影響實驗精度的原因，提出了用EXCEL和MATLAB軟件繪制實驗曲線確定普朗克常量的方法。關鍵詞關鍵詞光電效應普朗克常量MATLABEXCELUSINGMATLABEXCELTORESEARCHPHOTOELECTRICEFFECTEXPERIMENTOFTHEDETERMINATIONOFPLANCKCONSTANTABSTRACTTHISPAPERINTRODUCESTHEEXPERIMENTALPRINCIPLEOFPHOTOELECTRICEFFECTEXPERIMENTOFTHEDETERMINATIONOFPLANCKCONSTANTTHENLEADSTOTHEEXPERIMENTALDATAUSINGEXCELSOFTWAREMATLABSOFTWAREONANEWTREATMENTMETHODSTOIMPROVEACCURACYKEYWDSPRINCIPLEOFPHOTOELECTRICEFFECTPLANCKCONSTANTEXCELMATLAB引言引言一般我們是用制圖紙繪制零電流法測普朗克常量的線性直線，由于人為的誤差，筆的粗細，取點時候的偏差都會直接影響到精度的測定，所以本文比較使用EXCEL和MATLAB進行這些數(shù)據(jù)的處理，和繪圖分析，可以完全減少人為的誤差，提高實驗的準確精度，同時也可以讓我更加熟悉EXCEL和MATLAB的使用。正文正文實驗原理實驗原理1光電效應當一定頻率的光照射到某些金屬表面上時，可以使電子從金屬表面逸出，這種現(xiàn)象稱為光電效應。所產(chǎn)生的電子，稱為光電子。光電效應是光的經(jīng)典電磁理論所不能解釋的。當金屬中的電子吸收一個頻率為的光子時，便獲得這光子的全部能量，如果這VHV能量大于電子擺脫金屬表面的約束所需要的脫出功W，電子就會從金屬中逸出。按照能量守恒原理有（1）Ν1222012大學生物理實驗研究論文圖3是利用光電管進行光電效應實驗的原理圖。頻率為Ν、強度為P的光線照射到光電管陰極上，即有光電子從陰極逸出。如在陰極K和陽極A之間加正向電壓，它使K、A之間建立起的電場對從光電管陰極逸出的光電子起加速作用，隨著電壓的增加，到達陽極的光電子將逐漸增多。當正向電壓增加到時，光電流達到最大，不再增加，此時即稱為飽和狀態(tài)，對應的光電流即稱為飽和光電流。圖3光電效應原理圖由于光電子從陰極表面逸出時具有一定的初速度，所以當兩極間電位差為零時，仍有光電流I存在，若在兩極間施加一反向電壓，光電流隨之減少；當反向電壓達到截止電壓時，光電流為零。愛因斯坦方程是在同種金屬做陰極和陽極，且陽極很小的理想狀態(tài)下導出的。實際上做陰極的金屬逸出功比作陽極的金屬逸出功小，所以實驗中存在著如下問題1）暗電流和本底電流。當光電管陰極沒有受到光線照射時也會產(chǎn)生電子流，稱為暗電流。它是由電子的熱運動和光電管管殼漏電等原因造成的。室內各種漫反射光射入光電管造成的光電流稱為本底電流。暗電流和本底電流隨著K、A之間電壓大小變化而變化。（2）陽極電流。制作光電管陰極時，陽極上也會被濺射有陰極材料，所以光入射到陽極上或由陰極反射到陽極上，陽極上也有光電子發(fā)射，就形成陽極電流。由于它們圖4入射光頻率不同的I－U曲線圖5入射光強度不同的I－U曲線

簡介：安裝好MATLAB2012后再安裝目錄下點擊SETUPEXE會出現(xiàn)“查找安裝程序類時出錯，查找類時出現(xiàn)異常“的錯誤提示。該錯誤的解決方法是進入安裝目錄下的BIN文件夾雙擊MATLABEXE對安裝程序進行激活。這是可以對該MATLABEXE創(chuàng)建桌面快捷方式，以后運行程序是直接雙擊該快捷方式即可。信號運算信號運算1、信號加信號加MATLABMATLAB實現(xiàn)實現(xiàn)XX1X2XX1X22、信號延遲信號延遲YNXNKYNXNK3、信號乘信號乘XX1X2XX1X24、信號變化幅度信號變化幅度YKXYKX5、信號翻轉信號翻轉YFLIPLRXYFLIPLRX6、信號采樣和信號采樣和數(shù)學描述數(shù)學描述YYMATLABMATLAB實現(xiàn)實現(xiàn)YSUMXN1N2YSUMXN1N221NNNNX7、信號采樣積信號采樣積數(shù)學描述數(shù)學描述MATLABMATLAB實現(xiàn)實現(xiàn)YPRODXN1N2YPRODXN1N221NNNNXY8、信號能量信號能量數(shù)學描述數(shù)學描述NXNXE2||MATLABMATLAB實現(xiàn)實現(xiàn)EXSUMABSX2EXSUMABSX2數(shù)字濾波器的特性分析數(shù)字濾波器的特性分析1、脈沖響應、脈沖響應IMPZ函數(shù)函數(shù)調用方式調用方式（1）HTIMPZBA返回參數(shù)返回參數(shù)H是脈沖響應的數(shù)據(jù)，是脈沖響應的數(shù)據(jù)，T是脈沖響應是脈沖響應的時間間隔。的時間間隔。（2）HTIMPZBANN用來指定脈沖信號的長度。用來指定脈沖信號的長度。（3）HTIMPZBANFSFS用來指定脈沖信號的采樣頻率用來指定脈沖信號的采樣頻率4HTIMPZBAFS不再指定指定脈沖信號的長度。不再指定指定脈沖信號的長度。例例BABUTTER4005IMPZBA1002、頻率響應（幅頻響應和相頻響應）、頻率響應（幅頻響應和相頻響應）（1）數(shù)字濾波器頻率響應）數(shù)字濾波器頻率響應FREQZ函數(shù)函數(shù)調用方式調用方式HWFREQZBAN返回濾波器的返回濾波器的N點復頻率響應，點復頻率響應，BA分別是濾波器系數(shù)的分子和分別是濾波器系數(shù)的分子和分母向量。分母向量。H是復頻率響應，是復頻率響應，W是頻率點。是頻率點。HWFREQZBAN’WHOLE’采用單位圓上的采用單位圓上的N個點。個點。

簡介：基于基于MATLAB的滯后的滯后超前校正器的設計超前校正器的設計摘要摘要對控制系統(tǒng)的校正設計方法進行了簡單的介紹；介紹了基于MATLAB的滯后超前校正器的設計過程，并用仿真實例驗證了該方法比傳統(tǒng)的方法節(jié)省了相當大的工作量，實現(xiàn)起來非常的方便。利用MATLAB軟件中的控制系統(tǒng)工具箱和SIMULINK工具箱可以很方便的對控制系統(tǒng)進行建模、分析和設計。關鍵詞關鍵詞MATLAB；滯后超前校正器；設計1引言引言MATLAB（MATRIXLABATY即“矩陣實驗室”）是集數(shù)值分析、矩陣運算、信號處理和圖形顯示于一體，構成了一個方便的、界面友好的用戶環(huán)境，其強大的科學計算與可視化功能，簡單易用的開放式可擴展環(huán)境，使得MATLAB成為控制領域內被廣泛采用的控制系統(tǒng)計算與仿真軟件?！白詣涌刂圃怼笔枪た祁悓I(yè)一門重要的課程，其所需數(shù)學基礎寬而深、控制原理抽象、計算復雜且繁瑣以及繪圖困難等原因，使學生學習感覺枯燥并有畏難情緒。將MATLAB軟件應用到該門課程教學中，可以解決深奧繁瑣的計算，簡單、方便又精確的繪圖，并可以用豐富多彩的圖形來說明抽象的控制原理，可以提高學生的學習興趣。早期的校正器設計利用試湊法，其計算量非常大，而且還要手工繪制系統(tǒng)的頻率特性圖，很難達到滿意的結果。將MATLAB軟件應用到校正器設計中，則大大提高了設計的效率，并能很方便的達到滿意的效果。本文介紹在MATLAB環(huán)境下進行滯后超前校正器的設計方法。2控制系統(tǒng)校正設計概述控制系統(tǒng)校正設計概述在經(jīng)典控制理論中，系統(tǒng)校正設計，就是在給定的性能指標下，對于給定的對象模型，確定一個能夠完成系統(tǒng)滿足的靜態(tài)與動態(tài)性能指標要求的控制器（常稱為校正器或補償控制器），即確定校正器的結構與參數(shù)?？刂葡到y(tǒng)經(jīng)典校正設計方法有基于根軌跡校正設計法、基于頻率特性的BODE圖校正設計法及PID校正器設計法。按照校正器與給定被控對象的連接方式，控制系統(tǒng)校正可分為串聯(lián)校正、反饋校正、前饋校正和復合校正四種。串聯(lián)校正控制器的頻域設計方法中，使用的校正器有超前校正器、滯后校正器、滯后超前校正器等。超前校正設計方法的特點是校正后系統(tǒng)的截止頻率比校正前的大，系統(tǒng)的快速性能得到提高，這種校正設計方法對于要求穩(wěn)定性好、超調量小以及動態(tài)過程響應快的系統(tǒng)被經(jīng)常采用。滯后校正設計方法的特點是校正后系統(tǒng)的截止頻率比校正前的小，系統(tǒng)的快速性能變差，但系統(tǒng)的穩(wěn)定性能卻得到提高，因此，在系統(tǒng)快速性要求不是很高，而穩(wěn)定性與穩(wěn)態(tài)精度要求很高的場合，滯后校正設計方法比較適合。滯后超前校正設計是指既有滯后校正作用又有超前校正作用的校正器設計。它既具有了滯后校正高穩(wěn)定性能、高精確度的好處，又具有超前校正響應快、超調小的優(yōu)點，這種設計方法在要求較高的場合經(jīng)常被采用。3滯后滯后超前校正器的設計超前校正器的設計（6）用MATLAB軟件繪制校正后系統(tǒng)的BODE圖，并檢驗系統(tǒng)的各項性能指標是否滿足設計要求，若不滿足則可適當增大ΩC，重新執(zhí)行步驟（4）、（5）、（6），直至滿足要求。3333仿真實例仿真實例下面通過一實例來說明基于MATLAB軟件的滯后超前校正器的設計過程以及利用MATLAB軟件對控制系統(tǒng)進行時域分析和頻域分析。已知一單位反饋系統(tǒng)的開環(huán)傳遞函數(shù)為要求系統(tǒng)校正后滿足下列指標要求靜態(tài)速度誤差系數(shù)KV1801S，截止頻率ΩC≥32RADS，相角裕度?！?5，幅值裕度H≥10DB，試設計一個串聯(lián)滯后超前校正裝置。設計過程（1）求解K。KKV180（2）用MATLAB軟件繪制校正前原系統(tǒng)的BODE圖與階躍響應曲線，檢查是否滿足題目的設計要求。MATLAB程序如下NUM1180DEN1CONV10CONV011051SYS1TFNUM1DEN1MARGINSYS1FIGURE2SYS2FEEDBACKSYS11STEPSYS2程序運行后，可得未校正系統(tǒng)的BODE圖如圖1所示，還有未校正系統(tǒng)的階躍響應曲線如圖2所示。由圖1可讀出數(shù)據(jù)幅值裕度H235DB（即GM）、穿越頻率ΩX447RADS、相角裕度Γ471（即PM）及截止頻率ΩC143RADS，性能指標不滿足題目的要求。由圖2可看出原系統(tǒng)是不穩(wěn)定的，因此必須對原系統(tǒng)進行校正。

簡介：第九章對應分析第九章對應分析安慶師范學院安慶師范學院胡云峰胡云峰習題習題用對應分析法，對我國16個地區(qū)農民1982年支出情況的6項指標進行分類。16個地區(qū)農民生活水平的調查數(shù)據(jù)（單位元）食品衣著燃料住房生活用品及其他文化生活服務支出地區(qū)X1X2X3X4X5X6北京19033437797360544901904天津1352364104744161349394河北95212283932244228128山西104782511649891817325內蒙12841276389412582399327遼寧145683283177927293909347吉林159373338183711812529522黑龍江116222957132413762175604上海2211138641253115655082589江蘇1449829121167426273574浙江1699232751272471234355安徽153112309156223541818639福建144922126169619522175673江西14054215176419191597494山東11584302612233613377385河南10118232684620220543解第一步求出協(xié)方差矩陣AZ’Z程序X19033437797360544901904135203640104744161349394952122839302244228128010478251164098918173251284127638941258239932714568328317792729390934715937333818371181252952211622295713241376217560422111386412531156550825891449829121167426027305741699232751272471234355001531123091562235418186391449221261696195221756731405421517641919159749411584302612203361337738510118232684620220543NPSIZEXTSUMSUMXXLIEHEZEROS1PXHANGHEZEROS1NFK1PXLIEHEKSUMXKEND程序FZEROSP2FT11PFT10035705XT6FT20006905XT5ENDFGZF輸出結果F005620026000232002890058200198016740019100161006740018300089G0007900021000540002400001000130007500011000800001200035000240012800001000830000700246000120004400006000420000000062000270007700018000800002800030000250001800007最后我們在因子軸平面上左變量點和樣品點程序H1F1G1F2H2G1G2G2PLOTH1G1OHOLDONPLOTH2G2

簡介：1PRINCOMP功能主成分分析格式PCPRINCOMPXPCSCELATENTTSQUAREPRINCOMPX說明PCSCELATENTTSQUAREPRINCOMPX對數(shù)據(jù)矩陣X進行主成分分析，給出各主成分PC、所謂的Z得分SCE、X的方差矩陣的特征值LATENT和每個數(shù)據(jù)點的HOTELLINGT2統(tǒng)計量TSQUARE。2PCACOV功能運用協(xié)方差矩陣進行主成分分析格式PCPCACOVXPCLATENTEXPLAINEDPCACOVX說明PCLATENTEXPLAINEDPCACOVX通過協(xié)方差矩陣X進行主成分分析，返回主成分PC、協(xié)方差矩陣X的特征值LATENT和每個特征向量表征在觀測量總方差中所占的百分數(shù)EXPLAINED。3PCARES功能主成分分析的殘差格式RESIDUALSPCARESXNDIM說明PCARESXNDIM返回保留X的NDIM個主成分所獲的殘差。注意，NDIM是一個標量，必須小于X的列數(shù)。而且，X是數(shù)據(jù)矩陣，而不是協(xié)方差矩陣。4BARTTEST功能主成分的巴特力特檢驗格式NDIMBARTTESTXALPHANDIMPROBCHISQUAREBARTTESTXALPHA說明巴特力特檢驗是一種等方差性檢驗。NDIMBARTTESTXALPHA是在顯著性水平ALPHA下，給出滿足數(shù)據(jù)矩陣X的非隨機變量的N維模型，NDIM即模型維數(shù)，它由一系列假設檢驗所確定，NDIM1表明數(shù)據(jù)X對應于每個主成分的方差是相同的；NDIM2表明數(shù)據(jù)X對應于第二成分及其余成分的方差是相同的。第一種方法用MATLAB的各個函數(shù)組合得到的結果CLCCLEARX281110050521200123、PCSCELATENTTSQUAREPRINCOMP結果中PC是主成分系數(shù)，LATENT是特征值，應該和第一種方法中得到的特征向量和特征值相同；SCE是主成分得分，應該和Y相同。4、SPSS中得到的結果不同是因為其計算方法和MATLAB稍有差異，兩種結果應該成某個倍數(shù)關系（特征值取法不一樣），但是結論是一致的，不影響分析5、MATLAB中對原始數(shù)據(jù)標準化可以直接用ZSCE函數(shù)主成分分析和因子分析的區(qū)別1、因子分析中是把變量表示成各因子的線性組合，而主成分分析中則是把主成分表示成個變量的線性組合。2、主成分分析的重點在于解釋個變量的總方差，而因子分析則把重點放在解釋各變量之間的協(xié)方差。3、主成分分析中不需要有假設（ASSUMPTIONS），因子分析則需要一些假設。因子分析的假設包括各個共同因子之間不相關，特殊因子（SPECIFICFACT）之間也不相關，共同因子和特殊因子之間也不相關。4、主成分分析中，當給定的協(xié)方差矩陣或者相關矩陣的特征值是唯一的時候，的主成分一般是獨特的；而因子分析中因子不是獨特的，可以旋轉得到不到的因子。5、在因子分析中，因子個數(shù)需要分析者指定（SPSS根據(jù)一定的條件自動設定，只要是特征值大于1的因子進入分析），而指定的因子數(shù)量不同而結果不同。在主成分分析中，成分的數(shù)量是一定的，一般有幾個變量就有幾個主成分。和主成分分析相比，由于因子分析可以使用旋轉技術幫助解釋因子，在解釋方面更加有優(yōu)勢。大致說來，當需要尋找潛在的因子，并對這些因子進行解釋的時候，更加傾向于使用因子分析，并且借助旋轉技術幫助更好解釋。而如果想把現(xiàn)有的變量變成少數(shù)幾個新的變量（新的變量幾乎帶有原來所有變量的信息）來進入后續(xù)的分析，則可以使用主成分分析。當然，這中情況也可以使用因子得分做到。所以這中區(qū)分不是絕對的。總的來說，主成分分析主要是作為一種探索性的技術，在分析者進行多元數(shù)據(jù)分析之前，用主成分分析來分析數(shù)據(jù)，讓自己對數(shù)據(jù)有一個大致的了解是非常重要的。主成分分析一般很少單獨使用，通常用做以下用途A了解數(shù)據(jù)（SCREENINGTHEDATA）；B和CLUSTERANALYSIS一起使用；C和判別分析一起使用，比如當變量很多，個案數(shù)不多，直接使用判別分析可能無解，這時候可以使用主成份發(fā)對變量簡化（REDUCEDIMENSIONALITY）；

簡介：練習61設計增益K，使系統(tǒng)的相位裕度為60。TRANSFERFUNCTION417S0417S315S2467S0417GMINFPM600610WCGINFWCP25711K041605010SSK110SS校正后分析有階躍響應圖可知，校正之后系統(tǒng)更加穩(wěn)定，調節(jié)時間更短。練習63設被控對象開環(huán)傳遞函數(shù)為1010110SSSKSGG系統(tǒng)的靜態(tài)速度誤差系數(shù)，相位穩(wěn)定裕度，截止角頻率≥451100SKV030C。1S要求要求（1）設計滿足上述要求的超前控制器；（2）畫出校正前后的BODE圖；（3）用SIMULINK對原系統(tǒng)和校正后的系統(tǒng)進行仿真，畫出校正前后的階躍響應圖；（4）分析設計結果。

簡介：05050511第五講第五講第五講第五講MATLABMATLABMATLABMATLAB可視化可視化可視化可視化三三圖像顯示圖像顯示圖像顯示圖像顯示MATLABMATLABMATLAB進進行行圖圖像像處處理的步理的步理的步理的步驟驟如下如下如下如下讀讀取取圖圖像到像到像到像到MATLABMATLABMATLAB工作空工作空工作空工作空間間↓↓對圖對圖對圖對圖像矩像矩像矩像矩陣進陣進陣進陣進行格式行格式行格式行格式轉換轉換轉換轉換↓↓對圖對圖對圖對圖像矩像矩像矩像矩陣陣數(shù)據(jù)數(shù)據(jù)數(shù)據(jù)數(shù)據(jù)進進行行處處理理↓↓↓↓顯顯示示圖圖像像保存保存保存保存圖圖像像【【目錄目錄目錄目錄】】一、圖像文件格式一、圖像文件格式一、圖像文件格式一、圖像文件格式2211、調色板、調色板、調色板2222、圖像類型、圖像類型、圖像類型2233、圖像文件格式、圖像文件格式、圖像文件格式33二、讀圖像和圖像信息二、讀圖像和圖像信息二、讀圖像和圖像信息二、讀圖像和圖像信息3311、讀取圖像、讀取圖像、讀取圖像3322、讀取圖像信息、讀取圖像信息、讀取圖像信息44三、圖像類別與數(shù)據(jù)格式三、圖像類別與數(shù)據(jù)格式三、圖像類別與數(shù)據(jù)格式三、圖像類別與數(shù)據(jù)格式55四、圖像顯示四、圖像顯示四、圖像顯示四、圖像顯示7711、、IMSHOWINIMSHOWINIMSHOWIN7722、、IMSHOWILOWHIGHIMSHOWILOWHIGHIMSHOWILOWHIGH8833、、IMSHOWBWIMSHOWBWIMSHOWBW9944、、IMSHOWXMAPIMSHOWXMAPIMSHOWXMAP12121255、、IMSHOWRGBIMSHOWRGBIMSHOWRGB13131366、、IMSHOWIMSHOWIMSHOWFILENAMEFILENAMEFILENAME14141477、、SUBIMAGESUBIMAGESUBIMAGE141414五、保存圖像五、保存圖像五、保存圖像五、保存圖像1515151511、、IMWRITEIMWRITEIMWRITE函數(shù)函數(shù)函數(shù)151515六、圖像數(shù)據(jù)格式轉換六、圖像數(shù)據(jù)格式轉換六、圖像數(shù)據(jù)格式轉換六、圖像數(shù)據(jù)格式轉換1616161611、索引圖像、索引圖像、索引圖像16161622、灰度圖像、灰度圖像、灰度圖像1616160505053333、圖像文件格式、圖像文件格式、圖像文件格式、圖像文件格式文件文件文件文件類型類型類型類型說明說明說明說明11位位單色單色單色單色88位位灰度灰度灰度灰度88位位彩色彩色彩色彩色16161616位位灰度灰度灰度灰度24242424位位彩色彩色彩色彩色48484848位位彩色彩色彩色彩色BMPBMPBMP微微軟視軟視軟視軟視窗系窗系窗系窗系統(tǒng)圖統(tǒng)圖統(tǒng)圖統(tǒng)圖像文件格式像文件格式像文件格式像文件格式√√√√√√√√GIFGIFGIF用于網(wǎng)用于網(wǎng)用于網(wǎng)用于網(wǎng)頁頁上的上的上的上的動動畫、透明畫、透明畫、透明畫、透明√√√√√√TIFFTIFFTIFF可可儲儲存多幅存多幅存多幅存多幅圖圖像像√√√√√√√√√√√√PCXPCXPCXPAINTBRUSHPAINTBRUSHPAINTBRUSH軟軟件件圖圖像格式像格式像格式像格式√√√√√√√√JPGJPGJPG連續(xù)連續(xù)連續(xù)連續(xù)色色調調靜靜態(tài)圖態(tài)圖態(tài)圖態(tài)圖像數(shù)字像數(shù)字像數(shù)字像數(shù)字壓縮壓縮壓縮壓縮√√√√PSDPSDPSDPHOTOSHOPPHOTOSHOPPHOTOSHOP軟軟件件圖圖像格式像格式像格式像格式√√√√√√√√√√√√二、讀圖像和圖像信息二、讀圖像和圖像信息二、讀圖像和圖像信息二、讀圖像和圖像信息11、讀取圖像、讀取圖像、讀取圖像、讀取圖像函數(shù)函數(shù)函數(shù)函數(shù)IMREADIMREADIMREAD可以從任何可以從任何可以從任何可以從任何MATLABMATLABMATLAB支持的支持的支持的支持的圖圖像文件格式中，以任意位深度像文件格式中，以任意位深度像文件格式中，以任意位深度像文件格式中，以任意位深度讀讀取一幅取一幅取一幅取一幅圖圖像。格式像。格式像。格式像。格式為為XMAPIMREADFILENAMEFMTXMAPIMREADFILENAMEFMTXMAPIMREADFILENAMEFMT，其中，其中，其中，其中FILENAMEFILENAMEFILENAME－－為為需要需要需要需要讀讀入的入的入的入的圖圖像文件名稱，像文件名稱，像文件名稱，像文件名稱，F(xiàn)MTFMTFMT－－為圖為圖為圖為圖像格式。像格式。像格式。像格式?！尽纠俊繄D像讀取演示圖像讀取演示圖像讀取演示圖像讀取演示X1MAP1IMREADX1MAP1IMREADX1MAP1IMREAD演示演示演示演示圖圖像－像－像－像－11位黑白位黑白位黑白位黑白TIFTIFTIFX2MAP2IMREADX2MAP2IMREADX2MAP2IMREAD演示演示演示演示圖圖像－像－像－像－88位灰度位灰度位灰度位灰度TIFTIFTIFX3MAP3IMREADX3MAP3IMREADX3MAP3IMREAD演示演示演示演示圖圖像－像－像－像－256256256色色TIFTIFTIF

簡介：Ⅰ1統(tǒng)計工具箱函數(shù)統(tǒng)計工具箱函數(shù)表Ⅰ1概率密度函數(shù)函數(shù)名對應分布的概率密度函數(shù)BETAPDF貝塔分布的概率密度函數(shù)BINOPDF二項分布的概率密度函數(shù)CHI2PDF卡方分布的概率密度函數(shù)EXPPDF指數(shù)分布的概率密度函數(shù)FPDFF分布的概率密度函數(shù)GAMPDF伽瑪分布的概率密度函數(shù)GEOPDF幾何分布的概率密度函數(shù)HYGEPDF超幾何分布的概率密度函數(shù)NMPDF正態(tài)（高斯）分布的概率密度函數(shù)LOGNPDF對數(shù)正態(tài)分布的概率密度函數(shù)NBINPDF負二項分布的概率密度函數(shù)NCFPDF非中心F分布的概率密度函數(shù)NCTPDF非中心T分布的概率密度函數(shù)NCX2PDF非中心卡方分布的概率密度函數(shù)POISSPDF泊松分布的概率密度函數(shù)RAYLPDF雷利分布的概率密度函數(shù)TPDF學生氏T分布的概率密度函數(shù)UNIDPDF離散均勻分布的概率密度函數(shù)UNIFPDF連續(xù)均勻分布的概率密度函數(shù)WEIBPDF威布爾分布的概率密度函數(shù)表Ⅰ2累加分布函數(shù)函數(shù)名對應分布的累加函數(shù)BETACDF貝塔分布的累加函數(shù)BINOCDF二項分布的累加函數(shù)CHI2CDF卡方分布的累加函數(shù)EXPCDF指數(shù)分布的累加函數(shù)FCDFF分布的累加函數(shù)GAMCDF伽瑪分布的累加函數(shù)GEOCDF幾何分布的累加函數(shù)HYGECDF超幾何分布的累加函數(shù)LOGNCDF對數(shù)正態(tài)分布的累加函數(shù)NBINCDF負二項分布的累加函數(shù)NCFCDF非中心F分布的累加函數(shù)NCTCDF非中心T分布的累加函數(shù)表Ⅰ4隨機數(shù)生成器函數(shù)函數(shù)對應分布的隨機數(shù)生成器TARND貝塔分布的隨機數(shù)生成器BINND二項分布的隨機數(shù)生成器CHI2RND卡方分布的隨機數(shù)生成器EXPRND指數(shù)分布的隨機數(shù)生成器FRNDF分布的隨機數(shù)生成器GAMRND伽瑪分布的隨機數(shù)生成器GEND幾何分布的隨機數(shù)生成器HYGERND超幾何分布的隨機數(shù)生成器LOGNRND對數(shù)正態(tài)分布的隨機數(shù)生成器NBINRND負二項分布的隨機數(shù)生成器NCFRND非中心F分布的隨機數(shù)生成器NCTRND非中心T分布的隨機數(shù)生成器NCX2RND非中心卡方分布的隨機數(shù)生成器NMRND正態(tài)（高斯）分布的隨機數(shù)生成器POISSRND泊松分布的隨機數(shù)生成器RAYLRND瑞利分布的隨機數(shù)生成器TRND學生氏T分布的隨機數(shù)生成器UNIDRND離散均勻分布的隨機數(shù)生成器UNIFRND連續(xù)均勻分布的隨機數(shù)生成器WEIBRND威布爾分布的隨機數(shù)生成器表Ⅰ5分布函數(shù)的統(tǒng)計量函數(shù)函數(shù)名對應分布的統(tǒng)計量BETASTAT貝塔分布函數(shù)的統(tǒng)計量BINOSTAT二項分布函數(shù)的統(tǒng)計量CHI2STAT卡方分布函數(shù)的統(tǒng)計量EXPSTAT指數(shù)分布函數(shù)的統(tǒng)計量FSTATF分布函數(shù)的統(tǒng)計量GAMSTAT伽瑪分布函數(shù)的統(tǒng)計量GEOSTAT幾何分布函數(shù)的統(tǒng)計量HYGESTAT超幾何分布函數(shù)的統(tǒng)計量LOGNSTAT對數(shù)正態(tài)分布函數(shù)的統(tǒng)計量NBINSTAT負二項分布函數(shù)的統(tǒng)計量NCFSTAT非中心F分布函數(shù)的統(tǒng)計量NCTSTAT非中心T分布函數(shù)的統(tǒng)計量NCX2STAT非中心卡方分布函數(shù)的統(tǒng)計量

簡介：MATLABMATLAB語言與控制系統(tǒng)仿真實訓教程語言與控制系統(tǒng)仿真實訓教程參考答案參考答案第7章173控制系統(tǒng)根軌跡分析控制系統(tǒng)根軌跡分析MATLAB仿真實訓仿真實訓731實訓目的實訓目的1掌握運用MATLAB繪制180度、0度根軌跡圖的編程方法；2學會通過根軌跡圖獲取相關信息的方法；3利用仿真結果對系統(tǒng)根軌跡進行分析；4學會通過根軌跡圖分析和解決一些實際問題。732實訓內容實訓內容1單位反饋系統(tǒng)的開環(huán)傳遞函數(shù)為321SSSSKSGG試編程繪制閉環(huán)系統(tǒng)的根軌跡。并回答（1）閉環(huán)系統(tǒng)穩(wěn)定的的取值范圍；GK（2）系統(tǒng)的階躍響應有超調的的取值范圍；GK（3）分離點的坐標。GGJ01MZ1P023K1SYSZPKZPKRLOCUSSYS系統(tǒng)穩(wěn)定的的取值范圍；GK0GK系統(tǒng)的階躍響應有超調的的取值范圍GK4190GK分離點的坐標472D2設單位反饋控制系統(tǒng)的開環(huán)傳遞函數(shù)為1365312SSSSSKSG試編程繪制閉環(huán)系統(tǒng)的根軌跡。并回答（1）閉環(huán)系統(tǒng)穩(wěn)定的的取值范圍；KMATLABMATLAB語言與控制系統(tǒng)仿真實訓教程語言與控制系統(tǒng)仿真實訓教程參考答案參考答案第7章3由圖上數(shù)據(jù)可知閉環(huán)系統(tǒng)穩(wěn)定的的取值范圍；K2620K根軌跡與虛軸的交點坐標；183J分離點的坐標，。21D452J232，D4設單位反饋控制系統(tǒng)的開環(huán)傳遞函數(shù)為21SSSKSG的變化范圍為，試編程繪制閉環(huán)系統(tǒng)的根軌跡。并回答K0（1）閉環(huán)系統(tǒng)穩(wěn)定的的取值范圍；K（2）根軌跡與虛軸的交點坐標；（3）系統(tǒng)的單位階躍響應無超調的的取值范圍。K解系統(tǒng)的根軌跡方程121SSSK根軌跡仿真實訓第4題N11DCONV1012RLOCUSND

簡介：單鉆頭退火算法單鉆頭退火算法MATLAB編程編程CLEARCLCA0999溫度衰減函數(shù)的參數(shù)T097TF3TT0MARKOV_LENGTH2800MARKOV鏈長度CODINATESCODINATES1AMOUNTSIZECODINATES1城市的數(shù)目通過向量化的方法計算距離矩陣DIST_MATRIXZEROSAMOUNTAMOUNTCO_X_TMP1CODINATES1ONES1AMOUNTCO_X_TMP2CO_X_TMP1CO_Y_TMP1CODINATES2ONES1AMOUNTCO_Y_TMP2CO_Y_TMP1DIST_MATRIXSQRTCO_X_TMP1CO_X_TMP22CO_Y_TMP1CO_Y_TMP22SOL_NEW1AMOUNT產(chǎn)生初始解SOL_NEW是每次產(chǎn)生的新解；SOL_CURRENT是當前解；SOL_BEST是冷卻中的最好解；E_CURRENTINFE_BESTINFE_CURRENT是當前解對應的回路距離；E_NEW是新解的回路距離；E_BEST是最優(yōu)解的SOL_CURRENTSOL_NEWSOL_BESTSOL_NEWP1WHILETTFFR1MARKOV_LENGTHMARKOV鏈長度產(chǎn)生隨機擾動IFR05隨機決定是進行兩交換還是三交換兩交換IND10IND20WHILEIND1IND2IND1CEILRAMOUNTIND2CEILRAMOUNTENDTMP1SOL_NEWIND1DIST_MATRIXSOL_NEWAMOUNTSOL_NEW1IFE_NEWE_CURRENTE_CURRENTE_NEWSOL_CURRENTSOL_NEWIFE_NEWE_BEST把冷卻過程中最好的解保存下來E_BESTE_NEWSOL_BESTSOL_NEWENDELSE若新解的目標函數(shù)值小于當前解的，則僅以一定概率接受新解IFREXPE_NEWE_CURRENTTE_CURRENTE_NEWSOL_CURRENTSOL_NEWELSESOL_NEWSOL_CURRENTENDENDENDTTA控制參數(shù)T（溫度）減少為原來的A倍ENDDISP最優(yōu)解為DISPSOL_BESTDISP最短距離DISPE_BESTFIGURESETGCFNAMEMONITUIHUOSOL_BESTCOLRNLENGTHSOL_BESTSCATTERCODINATES1CODINATES250FILLEDHOLDONPLOTCODINATESSOL_BEST11CODINATESSOL_BESTN1CODINATESSOL_BEST12CODINATESSOL_BESTN2SETGCACOLGHOLDONFI2NPLOTCODINATESSOL_BESTI11CODINATESSOL_BESTI1CODINATESSOL_BESTI12CODINATESSOL_BESTI2HOLDONEND

簡介：第二講第二講導數(shù)與微分方程導數(shù)與微分方程一、實驗內容1、實際引例（牛頓冷卻模型）警察上午9點鐘發(fā)現(xiàn)一被謀殺者，并測得尸體溫度為324℃，一小時以后，尸體的溫度變?yōu)?17℃，尸體所在房間的溫度是20℃。如果人的正常體溫為365℃，并知道熱物體冷卻速度與自身溫度與外界溫度之差成正比，試推斷謀殺時間。解設T為尸體溫度，T為時間（十進制，如105代表10點30分）溫度的變化率（）比例常數(shù)為KDTDT解題關鍵熱物體冷卻速度與自身溫度與外界溫度之差成正比解題關鍵熱物體冷卻速度與自身溫度與外界溫度之差成正比。即KT20DTDT該方程就為微分方程，那么如何求解喃以下為MATLAB求解過程TDSOLVEDTKT20T9324T（9）表示在9點時的溫度為324℃ANS62EXPKT5EXP9K20由于由于1010點室溫為點室溫為317℃317℃，帶入后可求，帶入后可求KK命令如下命令如下先輸入SOLVE然后把上面得到的表達式復制進去62EXPKT5EXP9K20KSOLVE62EXPK105EXP9K20317（帶入T10T317K005810763080728074591965065204834T62EXPKT5EXP9K20T2091915298056906312642973833139EXP005810763080728074591965065204834T20TVPAT6保留幾位數(shù)T209192EXP00581076T200TSOLVE209192EXP00581076T200365ANS4083923903950647782688258218643400839236050354即謀殺時間即謀殺時間凌晨凌晨4點5分左右分左右ANSTANPI4T4、導數(shù)、積分的求法。導數(shù)、積分的求法。1、導數(shù)先要定義要用的字母SYMSXYADIFFY’X’N其中Y為待求導運算的函數(shù)，‘X’為求導變量，N為求導階次例1求YASINX對X求一階導二階導，以及對A求一階導。SYMSXAY定義變量DIFFASINXX1ANSACOSX對X求一階導數(shù)DIFFASINXX2ANSASINX對X求二階導數(shù)DIFFASINXA1ANSSINX對A求導，把SINX當做常數(shù)2、積分INTY’X’’A’’B’其中Y為待求函數(shù)，‘X’為積分變量，’A’’B’為積分區(qū)域，若區(qū)間為用字母INF，即INFINF詳見練習（7）例2DXCOSSYMSXAY定義變量INTCOSXANSCOSX過手練習（1）求Y’Y”2YXXY2331221X（大家要習慣MATLAB表示法和熟悉寫法的轉化）

簡介：134、兩路信號相加，并進行BPSK調制395、觀察調制信號頻譜，并經(jīng)AWGN信道加高斯白噪和窄帶強干擾406、接收機信號乘以恢復載波，觀察時域和頻域427、信號經(jīng)凱薩爾窗低通濾波器438、對濾波后信號分別用M1和M2進行解擴449、對兩路信號分別采樣，判決45八、產(chǎn)生隨機序列八、產(chǎn)生隨機序列GOLD碼和正交碼和正交GOLD碼471、產(chǎn)生GOLD碼并仿真其自相關函數(shù)482、產(chǎn)生正交GOLD碼并仿真其互相關函數(shù)50九、實驗心得體會九、實驗心得體會51