第一篇：Simulink一些知識(shí)總結(jié)
一、set_param 函數(shù)的應(yīng)用
（一）函數(shù)調(diào)用格式：
set_param(object,param,value)，其中object為模塊句柄，param 為參數(shù)，value 為參數(shù)值。
（二）舉例（1）啟動(dòng)模型仿真
set_param(gcs,'SimulationCommand','start')（2）精確控制模型仿真步長(zhǎng)
set_param(gcs, 'SimulationCommand', 'step')，step 為仿真步長(zhǎng)數(shù)。（3）設(shè)置模塊在模型中的位置
set_param('mymodel/Gain','Position',[50 250 125 275])（4）設(shè)置模塊的回調(diào)函數(shù)
set_param('mymodel/Gain','OpenFcn','my_open_cb')(5)設(shè)置模塊的朝向
set_param('mymodel/Gain','Orientation','left')
二、實(shí)時(shí)改變模塊的參數(shù)
單步仿真一次，改變一次參數(shù)的值，代碼如下：
set_param(gcs,'SimulationCommand','start')while 1 set_param(gcs,'SimulationCommand','pause')set_param(你的模塊名，'Resistance', 你想改變的值);%這一句你替換一下參數(shù)即可
set_param(gcs,'SimulationCommand','step')end
三、上次末狀態(tài)作為本次初始狀態(tài)進(jìn)行仿真
問：已經(jīng)建立好一個(gè)數(shù)據(jù)仿真模型 現(xiàn)在想單獨(dú)改變一個(gè)參數(shù)去觀察參數(shù)擾動(dòng)對(duì)結(jié)果的影響，例如參數(shù)變化范圍是1000 到1500，每7s 的時(shí)間讓參數(shù)變化100，用什么方法可以實(shí)現(xiàn) ？
答：分次仿真，每次仿真保存末狀態(tài)，下一次仿真前更改參數(shù)，再倒入上次仿真的末狀態(tài)作為初始狀態(tài),設(shè)置如下圖:
四、模型加上標(biāo)題
問：怎么給模型加上標(biāo)題？
答：雙擊模型空白位置，即可輸入text 文本，并可以設(shè)置字體大小顏色等。
五、模型仿真時(shí)間與實(shí)際運(yùn)行時(shí)間
Simulink 在Normal 模式下的仿真時(shí)間是非實(shí)時(shí)的，如果希望得知自己的模型轉(zhuǎn)換為實(shí)際產(chǎn)品時(shí)的運(yùn)行時(shí)間，需要將模型生成C 代碼，下載到硬件目標(biāo)中去運(yùn)行并使用Profiler Timer 進(jìn)行時(shí)間的記錄。如處理器在環(huán)仿真Processor in the Loop 可以測(cè)出模型實(shí)際的運(yùn)行時(shí)間。
六、關(guān)于模塊參數(shù)和模型解算配置的問題
（一）模塊參數(shù)的配置：根據(jù)用戶算法要求來(lái)配置，具體參數(shù)意義與使用方法可以雙擊模塊打開對(duì)話框后在help 按鈕獲取詳細(xì)幫助。
（二）模型結(jié)算方法配置：運(yùn)行仿真的算法一般由解算器提供，分為變步長(zhǎng)解算器和固定步長(zhǎng)解算器。默認(rèn)模型選擇變步長(zhǎng)解算器的ode45 算法，能夠適應(yīng)大多數(shù)模型應(yīng)用場(chǎng)景。如果用戶求解剛性模型，需要使用ode23s 之類算法。如果用戶是使用simulink 進(jìn)行控制算法建模，那么用戶一般使用固定步長(zhǎng)解算器以模擬與硬件芯片相同的時(shí)鐘頻率。
七、激活stateflow 狀態(tài)圖方法
（一）方法一：通過模型的調(diào)度時(shí)間周期性激活；
（二）方法二：通過trigger 信號(hào)去激活。
八、simulink 中過零問題討論
問：simulink 中有哪些過零檢測(cè)方法？是否有專門的過零檢測(cè)模塊？ 答：只有自適應(yīng)與非自適應(yīng)兩種選擇，沒有專門模塊，但是很多模塊里有是否使用過零檢測(cè)的checkbox，如Relational Operator，switch 模塊等。
九、模塊采樣時(shí)間與系統(tǒng)solver 最大步長(zhǎng)關(guān)系
當(dāng)模塊的Sample time 設(shè)置為-1 時(shí)表示繼承父層的采樣時(shí)間，如果模塊在子系統(tǒng)里就繼承子系統(tǒng)的采樣時(shí)間，如果模塊在最上層模型，就繼承系統(tǒng)solver 的步長(zhǎng)作為采樣時(shí)間間隔了。另外，如果模塊的采樣時(shí)間設(shè)置為-1 以外的值，如0 表示連續(xù)采樣，這種時(shí)刻需要solver 解算器類型支持連續(xù)狀態(tài)解算；當(dāng)模塊采樣時(shí)間為0.5 這樣的數(shù)字表示離散采樣時(shí)間，每隔0.5 計(jì)算一次，這時(shí)要求模塊的采樣時(shí)間必須是solver 計(jì)算步長(zhǎng)的整數(shù)倍。
十、solver 步長(zhǎng)間隔對(duì)仿真的影響
在變步長(zhǎng)解算器中，采樣時(shí)間間隔(即步長(zhǎng))會(huì)根據(jù)前后兩次計(jì)算值的誤差動(dòng)態(tài)改變。如果設(shè)置了最大步長(zhǎng)，那么每次的步長(zhǎng)就不會(huì)超出設(shè)置的這個(gè)上限。
十一、關(guān)于arduino 連接測(cè)溫模塊驅(qū)動(dòng)的問題
問：在arduino 上連接一個(gè)測(cè)溫模塊（如DS18B20），查看DS18B20 的技術(shù)手冊(cè)，說(shuō)是利用單總線通訊協(xié)議就可以讀取溫度，我想問問在有相關(guān)的庫(kù)函數(shù)及讀取溫度C 程序的情況下，怎么用simulink去驅(qū)動(dòng)測(cè)溫模塊工作，并通過串口將溫度傳回PC？
答：有相關(guān)的庫(kù)函數(shù)和C 例程的前提下，就把DS18B20 初始化的API 在model_initialize()中調(diào)用，在model_step()中使用讀取溫度的API、當(dāng)然，在Model.c 文件中包含DS18B20 庫(kù)函數(shù)所需要的頭文件。
十二、使用S-Function 創(chuàng)建自定義模塊如何避免重復(fù)初始化
問：使用S-Function 創(chuàng)建了一個(gè)自定義模塊，編寫了有關(guān)的模塊TLC 文件，并成功添加入Simulink 函數(shù)庫(kù)中，但是在使用時(shí)發(fā)現(xiàn)若同一模型文件中包含兩個(gè)以上該模塊，代碼生成時(shí)會(huì)重復(fù)進(jìn)行初始化，即TLC 文件中%functionStart(block, system)Output 中的內(nèi)容，如何解決？ 答：1.將這部分初始化放到BLockType 的tlc 函數(shù)中；
2.使用全局變量對(duì)Start 中的TLC 語(yǔ)句進(jìn)行判斷，僅第一次調(diào)用時(shí)生成，第二次以上屏蔽。
問：必須要在Start 中進(jìn)行，全局變量怎么定義？在模塊TLC 中進(jìn)行定義嗎？ 答：是的，當(dāng)這個(gè)變量為TLC_TRUE 時(shí)則不再執(zhí)行Start 函數(shù)，直接跳出。否則執(zhí)行。以保證此函數(shù)在整個(gè)編譯過程中只執(zhí)行一次。
十三、關(guān)于代碼空間的討論 問：在生成代碼后，如何查看代碼大概會(huì)占用多少RAM ROM，還有需要多少堆棧空間？
答：C 代碼看不出來(lái)的，在編譯完之后生成了map文件，就可以看到了。
十四、Embedded Coder、Simulink Coder 的關(guān)系 問：Embedded Coder 和Simulink Coder 有什么關(guān)系？
答：如果你想將Simulink 模型生成C 代碼，那么你需要使用Simulink Coder，而Simulink Coder的運(yùn)行需要依賴MATLAB Coder，所以，完整的講，如果你需要將Simulink 模型生成C 代碼，需要MATLAB Coder + Simulink Coder；如果你想將Simulink 模型生成為嵌入式C 代碼，那么，除了上述兩個(gè)Coder 之外，你還需要Embedded Coder。
十五、關(guān)于測(cè)量算法在處理器上執(zhí)行時(shí)間問題
問：最近在看論文，看到作者利用示波器探針測(cè)量電機(jī)MTPA 算法在DSP 上在線計(jì)算時(shí)間。我現(xiàn)在做了一個(gè)標(biāo)幺化定點(diǎn)后的電流指令表，怎么用示波器測(cè)量算法在處理器上執(zhí)行時(shí)間？
答：算法在目標(biāo)處理器上的執(zhí)行時(shí)間是一個(gè)非常重要的數(shù)據(jù)，我們更多時(shí)候關(guān)心這個(gè)數(shù)據(jù)的最大值，示波器測(cè)量的方式，我覺得不一定能采集到最長(zhǎng)時(shí)間的情形。在使用MBD 開發(fā)軟件的時(shí)候，建議使用PIL 獲取這個(gè)時(shí)間值，可以通過大量的測(cè)試用例，獲得這個(gè)最大時(shí)間值。
答：是的，matlab 的help 關(guān)于PIL 的用處也講到其中之一就是測(cè)量 execution time，把所需要的測(cè)量模塊封裝起來(lái)，給定一個(gè)輸入向量，下載到目標(biāo)板上，最后在matlab 窗口上輸入命令，能夠以圖像形式給出關(guān)于在此輸入向量下 average time maximum time。
十六、關(guān)于Simulink 模型的比較
問：在MBD 項(xiàng)目的工程管理中，模型的比較和差分是一個(gè)常見的問題。國(guó)外一些公司設(shè)計(jì)了專門用來(lái)對(duì)Simulink 模型進(jìn)行比較和充分的工具，比如德國(guó)的medini unite。我想請(qǐng)教一下，各位在進(jìn)行MBD 設(shè)計(jì)時(shí)，模型比較和差分的問題怎樣解決，都有哪些好用的工具，是否有開放代碼的工具。
答： MATLAB 產(chǎn)品提供了Simulink Report Generator，這個(gè)工具具有模型比較功能，也能實(shí)現(xiàn)一些模型合并。
十七、S-function builder 參數(shù)傳遞問題
問：怎么把s-function builder 中定義的參數(shù)傳遞到所集成的C 語(yǔ)言源文件中？ 答:在Data Properties 的頁(yè)面中定義參數(shù)的變量名何參數(shù)，如pvar。然后在Outputs, Discrete Update 等子頁(yè)面中都可以使用變量名和下標(biāo)訪問參數(shù)值，如pvar[0]。使用這個(gè)變量名編寫C代碼，完成所有配置之后點(diǎn)擊build，就可以將你寫的C 代碼打包集成到生成的C 文件中去了。
十八、快速原型和硬件在環(huán)的區(qū)別 問：快速原型和硬件在環(huán)的區(qū)別？
答：快速原型，也就是控制器快速原型的簡(jiǎn)稱，顧名思義，就是對(duì)Simulink 算法進(jìn)行控制器功能的快速實(shí)現(xiàn)。一般來(lái)講，你可以通過工控機(jī)，快速實(shí)現(xiàn)控制算法，而不必對(duì)控制算法進(jìn)行代碼生成相關(guān)的配置和設(shè)置。然后，你可以用這個(gè)工控機(jī)實(shí)現(xiàn)的控制器和被控對(duì)象相連，形成閉環(huán)，實(shí)現(xiàn)在真實(shí)被控對(duì)象上對(duì)控制算法的驗(yàn)證或者參數(shù)調(diào)節(jié)。硬件在環(huán)，是對(duì)控制器進(jìn)行測(cè)試的一個(gè)環(huán)節(jié)，在實(shí)現(xiàn)了控制器軟件和硬件的集成之后，這時(shí)候我們并不能確?？刂破饕欢ㄊ钦_的，如果這個(gè)時(shí)候直接把控制器和一些被控對(duì)象連接，可能造成財(cái)產(chǎn)或者人身傷害，在這種情況下，我們一般會(huì)通過工控機(jī)去模擬一個(gè)被控對(duì)象，然后把控制器和這個(gè)工控機(jī)相連，進(jìn)行控制器相關(guān)測(cè)試。
十九、tlc 文件該怎么編寫 問：tlc 文件該怎么編寫？
答：在Simulink Browser 中搜索S-Function,從example 模塊中找到帶有tlc 的例子，參考編寫即可。另外可以閱讀Embedded Coder 的User Guide 里面詳細(xì)介紹了tlc 的編寫方式。
二十、C MEX S-Function 延時(shí)功能
問：現(xiàn)在需要做一個(gè)PWM 波形輸出的S 函數(shù)，想利用延時(shí)功能進(jìn)行高低位的切換，但是不知道C-MEX 的S 函數(shù)怎么寫延時(shí)。
答：
1、利用狀態(tài)變量記錄輸入，在下幾個(gè)采樣時(shí)間后再輸出；
2、直接輸出，使用unit delay 實(shí)現(xiàn)延時(shí)。
二十一、S-Function 生成代碼如何添加自定義的頭文件
問：我用S 函數(shù)做了一個(gè)簡(jiǎn)單的s12xdp512 的Digital Output 的驅(qū)動(dòng)模塊，生成代碼時(shí)一直報(bào)錯(cuò)Error Test.c: 28 undeclared identifier PORTA' Error Test.c: 41 undeclared identifier `DDRA'。顯示的“PORTA”和“DDRA”都是在單片機(jī)的頭文件“mc9s12xdp512.h”中完成定義的，如何在tlc 文件中定義該頭文件，以便生成代碼時(shí)不會(huì)報(bào)錯(cuò)，代碼生成后不用手動(dòng)添加？
答：PORTA”和“DDRA”都是在單片機(jī)的頭文件“mc9s12xdp512.h”中完成定義的,你的Test.c 中有必要使用它們嗎？當(dāng)生成代碼時(shí)S 函數(shù)主要完成將GUI 用戶設(shè)定參數(shù)通過一系列中介文件傳遞到tlc 文件中去，最后生成的C 文件跟已有頭文件Merge 之后再使用IDE 編譯生成目標(biāo)文件下載到硬件中去。二
十二、PIL 測(cè)試時(shí)主機(jī)和DSP 的連接問題
問：PIL 測(cè)試時(shí)，主機(jī)和DSP 的連接，到底是用串口，USB 口，還是下載程序的JTAG 口？
答：串口和USB 口都是可以用的，串口的話需要
提供DSP 的串口驅(qū)動(dòng)代碼，USB 的話，PC 端和DSP端的都需要編寫，JTAG 就不確定了。
二十三、simulink 封裝的模塊Initialization有什么作用 問：simulink 封裝模塊Initialization 有什么作用，如何寫這個(gè)Initialization？
答：模塊的 Initialization Commands 執(zhí)行的時(shí)刻包括以下幾個(gè)情況：
1、在 Icon draw Commands 或 Initialization Commands 里更改封裝的參數(shù)時(shí)；
2、當(dāng) Icon draw Commands 有內(nèi)容時(shí)，翻轉(zhuǎn)或旋轉(zhuǎn)模塊；
3、雙擊打開模塊的參數(shù)對(duì)話框或點(diǎn)擊參數(shù)對(duì)話框上的 Apply/OK/Cancel 關(guān)閉對(duì)話框時(shí)。
二十四、simulink 生成C 代碼的問題
問：如果在simulink 中通過將模型生成C 代碼，但又想通過其他軟件如LabVIEW 調(diào)用此C 代碼，那么在simulink 中生成的C 代碼需要進(jìn)行一些手動(dòng)修改，還是說(shuō)自動(dòng)生成的C 代碼就可以進(jìn)行使用？
答：在選擇系統(tǒng)TLC 文件時(shí)要根據(jù)你的需要進(jìn)行選擇。如果是為了生成嵌入式C，則選擇ert.tlc文件； 如果是為了VC++ 環(huán)境使用，則選擇grt.tlc(for VC++)，Labview 里面你要調(diào)用此代碼，是指調(diào)用其中某個(gè)函數(shù)吧，如model_step()，將生成的代碼放到可搜索路徑下直接調(diào)用應(yīng)該就可以了。二
十五、simulink 生成C 代碼再轉(zhuǎn)化為DLL 文件的問題
問：通過simulink 中編寫仿真模型，再將模型生成為C 代碼，在VC 中生成DLL 文件，最后通過LabVIEW 調(diào)用DLL 文件，以實(shí)現(xiàn)混合編程的目的，應(yīng)該怎樣實(shí)現(xiàn)？ 答：系統(tǒng)目標(biāo)選擇ert_shrlib.tlc，可以在生成代碼的同時(shí)生成dll 文件。二
十五、Simulink 多個(gè)模塊生成代碼順序問題
問：在自動(dòng)代碼生成的時(shí)候，把一個(gè)模塊的代碼放到最前面，想知道模塊生成代碼的放置順序是通過什么來(lái)配置的？
答：如果有信號(hào)線連接的話，執(zhí)行順序按照信號(hào)流向（前提是DirectFeedThrough 設(shè)置正確）。沒有信號(hào)連接的模塊，郵件屬性可以設(shè)置Priority。數(shù)字大表示優(yōu)先級(jí)低。
二十六、如何封裝手寫底層代碼
問：有沒有一種方法能夠?qū)⑹謱懙牡讓觕 代碼封裝在matlab 的庫(kù)中，以后可以與simulink 自動(dòng)生成的代碼做接口的？
答：TSP 開發(fā)比較花費(fèi)時(shí)間。Legacy Code Tool對(duì)既有驅(qū)動(dòng)代碼的耦合性有要求，太過復(fù)雜的結(jié)構(gòu)無(wú)法直接應(yīng)用LCT。但是可以手動(dòng)開發(fā)。一般包括幾個(gè)步驟： 1.模塊和參數(shù)對(duì)話框制作 2.模塊S 函數(shù)開發(fā) 3.模塊TLC 文件開發(fā)
除了驅(qū)動(dòng)模塊之外，還有系統(tǒng)目標(biāo)文件和一系列回調(diào)函數(shù)的開發(fā)。共同構(gòu)成支持某一個(gè)目標(biāo)芯片的驅(qū)動(dòng)庫(kù)。二
十七、生成代碼該注意什么問題 問：生成代碼該注意什么問題
答：要生成代碼必須使用支持生成代碼的模塊，連續(xù)模塊一般不支持代碼生成，示波器等顯示模塊也不行。要生成代碼，需要確定生成代碼的類型，是嵌入式代碼還是一般Windows 應(yīng)用程序代碼或其他類型代碼，可以選擇對(duì)應(yīng)的System Target File，嵌入式通用代碼選擇ert.tlc, 快速原型代碼使用grt.tlc.如果是為了使用各個(gè)芯片廠商提供的芯片，最好還要安裝它們提供的TSP 工具箱和工具鏈，那里面會(huì)有配套的tlc 文件。二
十八、編寫S 函數(shù)時(shí)，采樣時(shí)間的設(shè)計(jì)
問：在編寫S 函數(shù)時(shí)，初始化模塊有關(guān)采樣時(shí)間ts=[x1 x2]，連續(xù)系統(tǒng)時(shí)設(shè)置為ts=[0 0];如果是離散系統(tǒng)，采樣時(shí)間為h，那么可以寫為ts=[h 0]嗎？例如20us 采樣一次，該怎么設(shè)計(jì)呢？
答：系統(tǒng)solver 設(shè)置為固定步長(zhǎng)，設(shè)置為20us，然后在S 函數(shù)里設(shè)置采樣ts = [-1, 0]來(lái)繼承系統(tǒng)采樣時(shí)間。二
十九、連續(xù)狀態(tài)是否支持自動(dòng)代碼生成 問：連續(xù)狀態(tài)是否支持自動(dòng)代碼生成？ 答1：含有連續(xù)狀態(tài)的模塊不支持代碼生成。因?yàn)橐纱a解算器必須是固定步長(zhǎng)的，但是此種解算器不支持連續(xù)狀態(tài)解算。答2：定步長(zhǎng)連續(xù)1 階積分好像可以生成代碼，剛開始用dSPACE 的時(shí)候全選的是這個(gè)，后來(lái)才開始用離散。
答3：是的，ode3 等定步長(zhǎng)算法也支持代碼生成?？磥?lái)連續(xù)模塊不能生成代碼的原因還需要討論。
三
十、sine wave 模塊中Sample based 和 Time based 的區(qū)別 問：sine wave 模塊中Sample based Time based的區(qū)別在哪里？
答：Sample based 基于采樣，Time based 基于時(shí)間，兩種計(jì)算正弦的數(shù)學(xué)公式，根據(jù)你使用情況需要選擇，如果你只有輸入時(shí)間t 采用基于時(shí)間。如果是離散系統(tǒng)，一般使用采樣點(diǎn)的序號(hào)來(lái)計(jì)算正弦。三
十一、simulink 中如何輸入帶參數(shù)的傳遞函數(shù)
問：如何輸入分子為s+a，分母為s^2+2bs+b^2的傳遞函數(shù)？a,b 其實(shí)不是變量，只是由其他模塊計(jì)算而來(lái)的。
答：transfer fcn 模塊支持輸入變量，變量只要定義在Workspace 里即可。三
十二、rtw 文件與tlc 文件的生成順序 問：rtw 文件與tlc 文件的生成順序是怎樣的？
答：先寫C 文件，編譯成mex 之后才可以被模型調(diào)用，S 函數(shù)模塊里的參數(shù)配置要正確，之后build 才能生產(chǎn)rtw。tlc 再?gòu)膔tw 里面獲取信息,進(jìn)行代碼生成。三
十三、M 語(yǔ)言編寫的S 函數(shù)是否支持代碼生成 問：M 語(yǔ)言編寫的S 函數(shù)是否支持代碼生成？
答：M 語(yǔ)言編寫的S 函數(shù)也是有兩種的，Level 1和Level2，Level1 僅支持仿真，不支持代碼生成。Level2 支持代碼生成。三
十四、如何在模型里面嵌入C 代碼
問：某部分內(nèi)容是c 代碼封裝的，如何在代碼生成時(shí)把該部分代碼直接連接添加？
答：可以使用custom code 模塊，直接添加C 代碼到模型中并生成到指定的位置，不進(jìn)行仿真也不確保代碼正確性。三
十五、Stateflow 如何代碼生成共用體
問1：stateflow 的RTW 自動(dòng)代碼生成，怎么設(shè)置才能生成c 語(yǔ)言的共用體？ 答1：樓主可以在Custom code 里面定義共同體，將此文件作為頭文件包含，然后在Stateflow 模型里使用此共同體變量或其成員，生成代碼時(shí)會(huì)自動(dòng)Merge 在一起。
問2：Stateflow 使用的時(shí)候，signal 的屬性怎么弄？比如：data type。答2：Stateflow 的運(yùn)行都是要先生成C 代碼，再編譯成動(dòng)態(tài)鏈接庫(kù)文件（mex），然后才能運(yùn)行的。所以，在custom code 設(shè)置通過C 代碼或者C 文件定義出來(lái)的函數(shù)或者變量，只要符合C 語(yǔ)言語(yǔ)法規(guī)則，就可以被用到Stateflow 模型的各種action 里。
三
十六、如何將simulink 自動(dòng)生成的代碼中如何使用外部RAM 問：因?yàn)槌绦蛱?，想在SIMULINK 自動(dòng)生成代碼時(shí)直接將PAGE1 定位到外部RAM，請(qǐng)問如何設(shè)置？
答：通過數(shù)據(jù)對(duì)象的方式定義各類數(shù)據(jù)，使其定義到期望的RAM 區(qū)域。在手工編碼的時(shí)候，遇到你說(shuō)的這種情況，我們一般會(huì)通過在link 文件中定義RAM 段，然后在C 文件中定義變量的時(shí)候，通過pragma 關(guān)鍵字，把這些變量定義到預(yù)先設(shè)定的地址段里。所以，你做自動(dòng)代碼生成的時(shí)候，也需要生成出來(lái)通過pragma 定義的變量，這需要你通過數(shù)據(jù)對(duì)象的方式實(shí)現(xiàn)，并且，這種情況下，需要自定義Package，以及package 里的類，在類的屬性storage class 里設(shè)置pragma。然后生成的代碼里，就可以有類似于以下的代碼： #pragma section_A begin int8_T var1;int16_T var2;............#pragma section_A end 這樣var1、var2 等就被定義到section_A 存儲(chǔ)區(qū)域了。
三
十七、如何拆分生成的代碼中model_data.c里面定義的parameter 結(jié)構(gòu)體？ 問：在使用simulink embedded coder 生成的代碼中，會(huì)將模型中的lookup table、constant 常數(shù)等統(tǒng)一放在了一個(gè)叫parameter 的結(jié)構(gòu)體內(nèi)(也就是在model_data.c 里面定義的數(shù)據(jù))。由于整個(gè)模型中的lookup table 很多，導(dǎo)致這個(gè)結(jié)構(gòu)體及其龐大，初步估計(jì)有7K byte 大小。請(qǐng)問有沒有辦法拆分這個(gè)結(jié)構(gòu)體的體積？
答：使用數(shù)據(jù)對(duì)象定義這些參數(shù)，就可以不用放到model_data.c 文件了。三
十八、目標(biāo)文件ert.tlc 與grt.tlc 的區(qū)別 問：grt 與ert 的區(qū)別在哪里？各自作用如何？
答：grt 中的g 表示general，ert 中的e 表示embedded，兩個(gè)文件都是系統(tǒng)目標(biāo)文件，使用grt.tlc 作為系統(tǒng)目標(biāo)文件，生成的代碼可以用作快速原型，而ert.tlc 生成的代碼可以用于嵌入式系統(tǒng)。ert.tlc 生成的是更高效的代碼。三
十九、關(guān)于觸發(fā)子系統(tǒng)的總結(jié) 當(dāng)系統(tǒng)內(nèi)部不存在積分環(huán)節(jié)時(shí)：
1.觸發(fā)子系統(tǒng)內(nèi)外不能使用goto 和from 模塊傳遞信號(hào)，這是由atomic 子系統(tǒng)屬性決定的
2.觸發(fā)子系統(tǒng)內(nèi)的模塊采樣時(shí)間只能是inf 或者-1，即繼承采樣時(shí)間；當(dāng)子系統(tǒng)內(nèi)部有信號(hào)源時(shí)需注意該點(diǎn)
3.觸發(fā)子系統(tǒng)的工作特性是在觸發(fā)的一瞬間輸出系統(tǒng)內(nèi)該時(shí)刻的值，然后保持；觸發(fā)子系統(tǒng)只可選擇held 信號(hào)，而不能reset 信號(hào)。
4.觸發(fā)子系統(tǒng)有信號(hào)源與外界接入信號(hào)源效果是一致的，但是內(nèi)部信號(hào)源要修改采樣時(shí)間為-1 5.子系統(tǒng)內(nèi)的信號(hào)顯示很有意思，是第一次觸發(fā)開始的值到最后一次觸發(fā)結(jié)束的值；中間數(shù)值都是保持狀態(tài)。當(dāng)系統(tǒng)內(nèi)部存在積分環(huán)節(jié)時(shí)：觸發(fā)子系統(tǒng)內(nèi)部有積分環(huán)節(jié)并且需要改變積分初值，這個(gè)問題沒有解決例如testforsubsystems中triggered subsystem。四
十、關(guān)于使能子系統(tǒng)的總結(jié)
1.當(dāng)系統(tǒng)內(nèi)部不存在積分環(huán)節(jié)時(shí)：使能觸發(fā)子系統(tǒng)相當(dāng)于一個(gè)觸發(fā)子系統(tǒng)，使能環(huán)節(jié)的作用應(yīng)包含在了觸發(fā)子系統(tǒng)之內(nèi)。
2.當(dāng)系統(tǒng)內(nèi)部存在積分環(huán)節(jié)時(shí)：使能觸發(fā)子系統(tǒng)解決了觸發(fā)子系統(tǒng)的有積分環(huán)節(jié)時(shí)不能解決的問題。四
十一、switch 的使用說(shuō)明
該模塊共有三個(gè)輸入，其中第二個(gè)（即u2）為控制信號(hào)的輸入，將其與閥值或者0 進(jìn)行判斷，判斷類型有：大于等于閥值、大于閥值、不等于0 三種類型，默認(rèn)判斷為真時(shí)輸出第一個(gè)信號(hào)，否則輸出第三個(gè)信號(hào)。
第二篇：通信原理課設(shè)文檔simulink
《通信原理》課程設(shè)計(jì)
開機(jī)，啟動(dòng)Matlab，在工作窗口Simulink。
模塊庫(kù)——>各個(gè)模塊 Simulink：基本模塊庫(kù)
Commuinications Blockset：通信模塊庫(kù) Signal Processing Blockset：信號(hào)處理模塊庫(kù)
基本原理：
1． 信源
通信工具箱：Communications Blockset——Comm Sources/Rondom Data Sources 選擇：Bernoulli Binary Generator二進(jìn)制伯努利序列產(chǎn)生器
Random Integer Generator 隨機(jī)整數(shù)產(chǎn)生器 Signal Processing Blockset——signal Processing Sources
Sine Wave正弦波 2． 噪聲源
Communications Blockset——Comm Sources/Noise Generators 選擇：Gaussian Noise Generator高斯白噪聲產(chǎn)生器
Rayleigh Noise Generator瑞利噪聲產(chǎn)生器
Uniform Noise Generator均勻噪聲產(chǎn)生器 3． 示波器
Simulink——〉Sinks——〉Scope 4． 相加器、增益、相乘器 Simulink——〉Math Operations——〉A(chǔ)dd、Gain、Produc 5． 信道
高斯白噪聲加性信道
Communications Blockset—〉Comm Sources—〉Noise Generators—〉Gaussian Noise Generator 6． 濾波器
Signal Processing Blockset——〉Filtering——〉Filter Designs——〉Digital Filter Design 低通濾波器的帶寬：基帶信號(hào)的帶寬
帶通濾波器的帶寬：載波頻率±基帶信號(hào)的帶寬 7． Relay（抽樣判決器）
Simulink——〉Discontinuities——〉Relay 這是一個(gè)滯環(huán)比較器，一般用來(lái)作為調(diào)節(jié)有波動(dòng)的系統(tǒng)時(shí)的緩沖。
parameters里的設(shè)置為：switch on point-閾值上限，switch off point-閾值下限，output when on-閾值上限輸出值，output when off-閾值下限輸出值。
例如某個(gè)信號(hào)的閾值上限為1.5，下限為0.5，上限輸出為1，下限輸出為0.當(dāng)該信號(hào)上升到大于1.5時(shí)，relay的輸出為1。若信號(hào)下一個(gè)周期小于1.5，但仍大于0.5時(shí)，其輸出仍保持為1。只有當(dāng)信號(hào)下降到小于0.5時(shí)，才會(huì)輸出0。
它的意義是使信號(hào)調(diào)節(jié)有一個(gè)足夠的范圍寬度，而不至于因?yàn)槊總€(gè)周期都進(jìn)行調(diào)節(jié)。8． 眼圖、星座圖、誤碼率
Communications Blockset—〉Comm Sinks——〉Discrete-Time Eye Diagram Scope ? Samples per symbol：每個(gè)符號(hào)的采樣點(diǎn)數(shù)，自己計(jì)算1個(gè)二進(jìn)制碼元在系統(tǒng)中有多少個(gè)采樣點(diǎn)，也就是要保證眼圖窗口的時(shí)間寬度為整數(shù)個(gè)碼元寬度（最好1、2個(gè)）。? Offset：非負(fù)整數(shù)，小于Samples per symbol和Samples per trace的乘積，指定在畫第一個(gè)點(diǎn)時(shí)省略的采樣點(diǎn)數(shù)。可調(diào)。
? Samples per trace：正數(shù)，指定每一條軌跡的符號(hào)數(shù)目。如2，顯示兩個(gè)符號(hào)周期。? Traces displayed：疊加的軌跡的數(shù)目。
? New traces per display：正整數(shù)，小于Traces displayed，指定每一次顯示時(shí)新軌跡的數(shù)目。
? Rendering Properties：線型和顏色設(shè)計(jì) ? Axes Properties：x、y軸范圍設(shè)計(jì)
? Figure Properties：眼圖顯示的數(shù)據(jù)類型設(shè)計(jì)
Communications Blockset—〉Comm Sinks——〉Discrete-Time Scatter Plot Scope Communications Blockset—〉Comm Sinks——〉Error Rate Calculation
第三篇：實(shí)驗(yàn)報(bào)告 5 Simulink仿真
實(shí)驗(yàn)五 Simulink仿真
（一）一、實(shí)驗(yàn)?zāi)康?/p>
1、熟悉Simulink仿真環(huán)境
2、了解Simulink基本操作
3、了解Simulink系統(tǒng)建?；痉椒?/p>
3、熟悉Simulink仿真系統(tǒng)參數(shù)設(shè)置和子系統(tǒng)封裝的基本方法
二、實(shí)驗(yàn)內(nèi)容
1、在matlab命令窗口中輸入simulink,觀察其模塊庫(kù)的構(gòu)成；
2、了解模塊庫(kù)中常用模塊的使用方法；
3、已知單位負(fù)反饋系統(tǒng)的開環(huán)傳遞函數(shù)為
G?100s?2s(s?1)(s?20)
建立系統(tǒng)的模型，輸入信號(hào)為單位階躍信號(hào)，用示波器觀察輸出。
4、建立一個(gè)包含Gain、Transfer Fcn、Sum、Step、Sine Wave、Zero-Pole、Integrator、Derivative等模塊構(gòu)成的自定義模塊庫(kù)Library1；
5、建立如圖7-12所示的雙閉環(huán)調(diào)速系統(tǒng)的Simulink的動(dòng)態(tài)結(jié)構(gòu)圖，再把電流負(fù)反饋內(nèi)環(huán)封裝為子系統(tǒng)，建立動(dòng)態(tài)結(jié)構(gòu)圖。
三、實(shí)驗(yàn)結(jié)果及分析：
圖5-1
圖5-2
圖5-3
圖5-4雙閉環(huán)調(diào)速系統(tǒng)的Simulink的動(dòng)態(tài)結(jié)構(gòu)圖
圖5-5把電流負(fù)反饋內(nèi)環(huán)封裝為子系統(tǒng)的動(dòng)態(tài)結(jié)構(gòu)圖
雙擊Subsystem模塊，編輯反饋電流環(huán)Subsystem子系統(tǒng)，如圖5-6所示：
圖5-6
分析：Simulink是Mathworks開發(fā)的MATLAB中的工具之一，主要功能是實(shí)現(xiàn)動(dòng)態(tài)系統(tǒng)建模、仿真與分析?？梢栽趯?shí)際系統(tǒng)制作出來(lái)之前，預(yù)先對(duì)系統(tǒng)進(jìn)行仿真與分析，并可對(duì)系統(tǒng)做適當(dāng)?shù)倪m時(shí)修正或按照仿真的最佳效果來(lái)調(diào)試及整定控制系統(tǒng)的參數(shù)，達(dá)到提高系統(tǒng)性能。減少涉及系統(tǒng)過程中的反復(fù)修改的時(shí)間、實(shí)現(xiàn)高效率地開發(fā)系統(tǒng)的目標(biāo)。Simulink提供了建模、分析和仿真各種動(dòng)態(tài)系統(tǒng)的交互環(huán)境，包括連續(xù)系統(tǒng)、離散系統(tǒng)和混雜系統(tǒng)，還提供了采用鼠標(biāo)拖放的方法建立系統(tǒng)框圖模型的圖形交互界面。
第四篇：Matlab與Simulink仿真學(xué)習(xí)心得
Matlab與Simulink仿真學(xué)習(xí)心得
班級(jí)：07610 學(xué)號(hào)：072016 姓名：呂天雄
一 Matlab學(xué)習(xí)心得體會(huì)與編程實(shí)踐
<1>學(xué)習(xí)Matlab的心得體會(huì)
真正開始接觸Matlab是大二上就開始了，到現(xiàn)在已經(jīng)一年多了，在此之間，Matlab的確為我提供了很多便利。Matlab的確不愧成為是草稿紙上的語(yǔ)言。我們不必去為很簡(jiǎn)單的顯示效果圖形去找一些什么其他軟件或者研究比較復(fù)雜的計(jì)算機(jī)圖形學(xué)，一個(gè)plot或者別的函數(shù)往往就可以得到很滿意的效果。
其實(shí)最初開始學(xué)習(xí)matlab的時(shí)候感覺這個(gè)東西和C沒什么兩樣，但是后來(lái)具體到一些 東西，比如信號(hào)處理和數(shù)學(xué)建模上以后才感覺到使用matlab編寫程序去驗(yàn)證結(jié)果比C要節(jié) 省很多時(shí)間，而且matlab寫東西基本都是按照自己的思路平鋪直敘很少去考慮什么函數(shù)的 嵌套調(diào)用或者指針等等很頭疼的東西。
關(guān)于matlab的學(xué)習(xí)，我感覺其實(shí)百度和matlab自帶的help基本能夠解決絕大數(shù)問題，而且一些比較好的論壇比如004km.cn都會(huì)為你產(chǎn)生很大的幫助，關(guān)鍵是在于多動(dòng)手實(shí)踐，多思考。但是matlab畢竟只是一個(gè)工具，原理和一些基本的編程素質(zhì)還是必須有的，否則matlab最多也只能是驗(yàn)證一些別人的東西而已，根本幫不上什么忙的。
<2>遇到的一些問題的思考方式與解決辦法
最開始用matlab的時(shí)候是在大物實(shí)驗(yàn)，實(shí)驗(yàn)要求去根據(jù)測(cè)量得到的數(shù)據(jù)作出圖。但是 手動(dòng)用鉛筆去畫確實(shí)很麻煩，所以用matlab確實(shí)可以省去很大的麻煩。但是第一次遇到問 的時(shí)候是有關(guān)極化坐標(biāo)下的曲線擬合。
首先是一個(gè)物理實(shí)驗(yàn)的問題；在做一個(gè)關(guān)于光的偏振的實(shí)驗(yàn)的時(shí)候，最后的結(jié)果要在一個(gè)極化坐標(biāo)下顯示出來(lái)；因?yàn)閿?shù)據(jù)是離散的，所以顯示出來(lái)的圖像是一個(gè)折來(lái)折去的一個(gè)東東； 然后很自然的想法是對(duì)這個(gè)曲線進(jìn)行插值處理。
但是極化坐標(biāo)下MATLAB并未提供插值處理的函數(shù)，interp1這個(gè)函數(shù)只能在笛卡爾坐標(biāo)系，也就是直角坐標(biāo)系下使用。
然后就想到把極坐標(biāo)的數(shù)據(jù)轉(zhuǎn)換的直角坐標(biāo)系下，pol2cart可以實(shí)現(xiàn)這個(gè)想法，但是隨后而來(lái)，也就是最后導(dǎo)致整個(gè)問題失敗的關(guān)鍵也在這里。
pol2cart以后產(chǎn)生的一串?dāng)?shù)據(jù)中出現(xiàn)了重復(fù)的數(shù)據(jù)，那么interp1這個(gè)東西也就無(wú)能為力了，因?yàn)閕nterp1不能處理一串?dāng)?shù)據(jù)中有重復(fù)出現(xiàn)的情況。最后的處理辦法是把這些數(shù)據(jù)c os，sin這些東西變換一下后，使其大致規(guī)則，然后再用polar畫出極坐標(biāo)下的圖形。
接下來(lái)這個(gè)問題就有點(diǎn)超過我的范圍了，可能會(huì)有點(diǎn)敘述不清楚。問題可以概要為：人臉網(wǎng)格插值。
這是一個(gè)用三角形網(wǎng)格表示的一個(gè)人臉模型。需要通過插值使其變得光滑，當(dāng)然問題的復(fù)雜之處在于，插值會(huì)改變?cè)镜木W(wǎng)格結(jié)構(gòu)。
對(duì)于這個(gè)問題許多人都給出了解決的辦法，當(dāng)然是一些關(guān)于人臉識(shí)別技術(shù)和運(yùn)動(dòng)圖像處理的范疇之內(nèi)的。
主要有兩種，第一種是face—ls算法。這種算法是基于RBF（徑向插值函數(shù)）和loop細(xì)分原則的一種算法，當(dāng)然其精確程度比較差點(diǎn)。但是速度挺好。原理是：基于網(wǎng)格上的點(diǎn)，產(chǎn)生出頂點(diǎn)的迭代函數(shù)，從而產(chǎn)生出新的定點(diǎn)，進(jìn)行細(xì)化。然后進(jìn)行插值。也就是引進(jìn)新的頂點(diǎn)。分為兩個(gè)過程內(nèi)部插值和邊緣插值。具體可以去看蘭州大學(xué)信息工程學(xué)院的學(xué)報(bào)。
第二種也就是基本的老辦法RBF。首先解釋一下什么叫做徑向插值 徑向基函數(shù)
主要是考慮多維空間的數(shù)據(jù)插值問題，徑向基函數(shù)在三維圖形的變形中常常被采用，用其來(lái)變分三維人臉的時(shí)候還應(yīng)結(jié)合人臉的固有特征。確定該方法來(lái)實(shí)現(xiàn)三維人臉變形是可行的，但要構(gòu)造好的基函數(shù)，以及解這些高維方程的解仍是難點(diǎn)。插值方法：
假設(shè)大多數(shù)人臉的形狀都可以由一個(gè)拓?fù)湓妥兓脕?lái)，那么，通過調(diào)整一個(gè)一般模型的構(gòu)造參數(shù)可以建立不同的面部模型。但是，這種參數(shù)模型僅僅局限于那些構(gòu)造參數(shù)已知的情況，并且對(duì)特定人臉參數(shù)的調(diào)整非常困難。在離散數(shù)據(jù)的多變量插值問題方面，徑向基函數(shù)（radial basis function，RBF）插值方法是一個(gè)行之有效的工具，所以也適用于類似人臉這樣高維曲面的近似或平滑插值?，F(xiàn)有的許多方法使用了基于RBF的插值技術(shù)，將一般人臉網(wǎng)格變化到特定人臉的形狀。這種方法的優(yōu)點(diǎn)在于：（1）通過插值可以得到丟失的數(shù)據(jù)點(diǎn)，所以源網(wǎng)格和目標(biāo)網(wǎng)格不需要相同數(shù)目的結(jié)點(diǎn)；（2）如果選擇了合適的匹配點(diǎn)，數(shù)學(xué)上可以保證能夠?qū)⒃淳W(wǎng)格變形到目標(biāo)網(wǎng)格。
當(dāng)然我也不是很懂，只是理解了部分。這種方法的關(guān)鍵是找到一個(gè)很好的核函數(shù)來(lái)計(jì)算出新的頂點(diǎn)。
最后問題解決的方式是查著了一些圖書館的學(xué)報(bào)。然后而且在網(wǎng)上找到了部分代碼。然后做出來(lái)的。這是我校大四一個(gè)同學(xué)畢設(shè)中的一個(gè)東西。
其實(shí)這個(gè)事件最好的解決辦法是在3D-MAX中，不過因?yàn)闆]能找到MATLAB與3D-MAX是怎么接口的，所以作罷了。
第三次遇到的問題是我在數(shù)學(xué)建模課上遇到的一個(gè)關(guān)于矩陣LU分解的問題：
Matlab作出的結(jié)果和手算的結(jié)果竟然不一樣。
10??1????1???1?
U=?1? ?一個(gè)矩陣 A=? 經(jīng)過三角分解為L(zhǎng)=101????????11?因?yàn)?/?為一個(gè)極大數(shù) 則1—1/?可以看成—1/?；
設(shè)?=0.***000000001%10的負(fù)24次方
從而L*U=則最后一位1沒有了 變成了0
其中L=
U=Matlab 程序?yàn)?/p>
epusino=0.***000000001;U=[epusino 1;0-1/epusino];L=[1 0;1/epusino 1];L*U ans =
0.***
1.***
1.***
0 而直接用lu函數(shù)對(duì)矩陣A分解的結(jié)果為
A=[epusino 1;1 1] A =
0.0000
1.0000
1.0000
1.0000
>> [L,U]=lu(A)L =
0.0000
1.0000
1.0000
0
U =
0 我們老師當(dāng)時(shí)的解釋是這是由于我們所用matlab是破解版的緣故，破解版的計(jì)算精度 沒正版的高所以導(dǎo)致運(yùn)算結(jié)果的錯(cuò)誤，后來(lái)我看了一些LU這個(gè)函數(shù)的help文檔 psychologically lower triangular matrix"(i.e.a product of lower triangular and permutation matrices)in L。翻譯過來(lái)是L是一個(gè)心理上的下三角矩陣，其實(shí)是下三角矩陣和置換矩陣也就是最后的結(jié)論。P*A = L*U.。由此可見這并不是什么精度的問題導(dǎo)致的結(jié)果只是LU分解的矩陣意義和書本上的不同而已。
以下是我自己寫的一些程序的代碼
1：網(wǎng)絡(luò)隨機(jī)拓?fù)鋱D
目的是要生成一個(gè)度數(shù)隨機(jī)、權(quán)值隨機(jī)的拓?fù)鋱D。首先的理解是利用鄰接矩陣，先隨機(jī)生成一個(gè)矩陣，該矩陣為一個(gè)對(duì)稱矩陣。然后畫出這個(gè)矩陣就行了。
隨機(jī)拓?fù)鋱D10026516398191 ***0430 ***065
%產(chǎn)生數(shù)組A用來(lái)存放表示兩點(diǎn)之間權(quán)值的矩陣A，也就是臨接矩陣，那么兩點(diǎn)之間權(quán)值不為零元素的個(gè)數(shù)即為該點(diǎn)的度數(shù)
DEF=5;%設(shè)定一個(gè)東東 方便改變隨機(jī)點(diǎn)的個(gè)數(shù) A=rand(DEF,DEF);%產(chǎn)生DEF*DEF的隨機(jī)矩陣 for i=1:DEF A(i,i)=0%將對(duì)角線上的數(shù)置為0 end A=10*A;A=floor(A);%向下去整 for i=1:DEF for j=1:i A(j,i)=A(i,j)%將A矩陣變?yōu)橐粋€(gè)上三角或者下三角矩陣 end end x=100*rand(1,DEF);y=100*rand(1,DEF);%產(chǎn)生10個(gè)隨機(jī)的點(diǎn) plot(x,y,'r+');
for i=1:DEF a=find(A(i,:)>0)%將A矩陣每行大于0的數(shù)的在該行的地址找出來(lái)放在a中
for j=1:length(a)
c=num2str(A(i,j));%將A中的權(quán)值轉(zhuǎn)化為字符型
if c~='0'%不顯示為0的值 因?yàn)锳矩陣為零代表兩點(diǎn)不相連 text((x(i)+x(j))/2,(y(i)+y(j))/2,c,'Fontsize',18);%將權(quán)值顯示在兩點(diǎn)連線中間 end hold on;
line([x(i)x(a(j))],[y(i)y(a(j))]);%連線 end end
title('隨機(jī)拓?fù)鋱D');e=num2str(DEF);legend(e);%左上角顯示節(jié)點(diǎn)的個(gè)數(shù) for m=1:DEF A(m,m)=m;f=num2str(A(m,m));hold on;text((x(m)+x(m))/2,(y(m)+y(m))/2,f,'Fontsize',18);%將權(quán)值顯示在兩點(diǎn)連線中間 end 接下來(lái)是一個(gè)在信息安全課上寫的一個(gè)關(guān)于256色圖分層的程序，一副圖像有m*n個(gè)像素然后每個(gè)像素是一個(gè)8bit的二進(jìn)制數(shù)據(jù)換算為十進(jìn)制是0-255之間。要做的就是把每bit的信息提取出來(lái)。
I=imread('bupt副本.bmp');%讀入源圖像
I=double(I);%將圖像轉(zhuǎn)換為double類型便于MATLAB中的計(jì)算 %確定圖像的長(zhǎng)寬 M=size(I,1);%長(zhǎng) N=size(I,2);%寬 l=M*N;%圖像長(zhǎng)寬之積 for n=1:8;%剔除每層 for i=1:M
for j=1:N
B=numdec2bin(I(i,j),8);%先將每個(gè)像素轉(zhuǎn)換為8位2進(jìn)制序列
B(n)=[];%刪除相應(yīng)位置的元素
y(i,j)=numbin2dec(B);%將剩下的7個(gè)元素轉(zhuǎn)換為10進(jìn)制數(shù)放入y矩陣中
end end
y=uint8(y);%將10進(jìn)制的矩陣轉(zhuǎn)換為無(wú)符號(hào)整形
subplot(4,2,9-n);%顯示n個(gè)圖像，每層是相應(yīng)的剔除該層的結(jié)果
imshow(y);%顯示這個(gè)圖像
title(strcat('去除第',num2str(9-n),'層后的結(jié)果'));%顯示標(biāo)題 end
附帶的兩個(gè)函數(shù)：
function y=numdec2bin(x,n);%從函數(shù)將十進(jìn)制數(shù)x轉(zhuǎn)化為n位二進(jìn)制 y=zeros(1,n);a=x;i=n;while a>0 y(i)=mod(a,2);a=floor(a/2);i=i-1;end function y=numbin2dec(x)%將二進(jìn)制數(shù)轉(zhuǎn)換為十進(jìn)制數(shù) a=0;lx=length(x);for i=1:lx a=a+x(i)*2^(lx-i);end y=a;
去除第1層后的結(jié)果去除第2層后的結(jié)果去除第3層后的結(jié)果去除第4層后的結(jié)果去除第5層后的結(jié)果去除第6層后的結(jié)果去除第7層后的結(jié)果去除第8層后的結(jié)果
Matlab的RGB轉(zhuǎn)換為YCbCr之后轉(zhuǎn)換回來(lái)和原圖不符
這個(gè)是做JPEG圖像壓縮的時(shí)候遇到的一個(gè)問題，當(dāng)時(shí)的代碼是用C寫的但是結(jié)果不對(duì)，后來(lái)想到用matlab進(jìn)行驗(yàn)證結(jié)果發(fā)現(xiàn)通過公式轉(zhuǎn)換得到的圖是錯(cuò)的。
原理：BMP圖像壓縮為JPEG的第一步是將RGB色彩空間通過這個(gè)公式映射到Y(jié)CbCr空間上。
Y=0.299 R + 0.587 G + 0.114 B Cb =0.3313G + 0.5 B + 128 Cr = 0.5 R0.0813 B + 128 然后再進(jìn)行量化，DCT，編碼等等步驟。JPEG解壓時(shí)需要將YCbCr空間的圖像又轉(zhuǎn)化回來(lái)到RGB上。
R = Y + 1.402(Cr0.34414(Cb-128)128)
B = Y + 1.772(Cb-128)
但是通過書上給的公式和網(wǎng)上大部分公式卻發(fā)現(xiàn)根本轉(zhuǎn)化不會(huì)來(lái)。下面是寫的程序代碼以及測(cè)試得到的結(jié)果圖像。
clc,clear;Source=imread('hl.jpg');%讀入原始RGB圖像
figure(1);subplot(1,2,1);imshow(Source):title('original image');%顯示圖像 [r c d]=size(Source);%計(jì)算圖像大小 %------計(jì)算紅色分量并顯示分解圖------% R(:,:,1)=Source(:,:,1);R(:,:,2)=zeros(r,c);R(:,:,3)=zeros(r,c);R=uint8(R);whos;figure(2);subplot(1,3,1);imshow(R)title('Red Component');%-------計(jì)算綠色分量并顯示分解圖-------% G(:,:,2)=Source(:,:,2);G(:,:,1)=zeros(r,c);G(:,:,3)=zeros(r,c);G=uint8(G);figure(2);subplot(1,3,2);imshow(G)title('Green Component');%--------計(jì)算藍(lán)色分量并顯示分解圖-------% B(:,:,3)=Source(:,:,3);B(:,:,1)=zeros(r,c);B(:,:,2)=zeros(r,c);B=uint8(B);figure(2);subplot(1,3,3)imshow(B)title('Blue Component');%------------合成-------------% Comp(:,:,1)=R(:,:,1);Comp(:,:,2)=G(:,:,2);Comp(:,:,3)=B(:,:,3);figure(1);subplot(1,2,2);imshow(Comp):title('composition image');Y=0.229*R+0.587*G+0.114*B;Cb=0.5*B-0.1687*R-0.3313*G+128;Cr=0.5*R-0.4187*G-0.0813*B+128;red=Y+1.402*(Cr-128);green=Y-0.34414*(Cb-128)-0.71414*(Cr-128);blue=Y+1.772*(Cb-128);Comp2(:,:,1)=red(:,:,1);Comp2(:,:,2)=green(:,:,2);Comp2(:,:,3)=blue(:,:,3);figure(3);imshow(Comp2);title('RGB轉(zhuǎn)換為YCrCb后又轉(zhuǎn)換為RGB的圖像');RD=R(:,:,1)-red(:,:,1);GD=G(:,:,2)-green(:,:,2);BD=B(:,:,3)-blue(:,:,3);figure(4);subplot(1,3,1);imshow(RD);title('紅色分量差異');subplot(1,3,2);imshow(GD);title('綠色分量差異');subplot(1,3,3);imshow(BD);title('藍(lán)色分量差異');
然而可以通過figure(4)的圖像清楚看到R 和G分量轉(zhuǎn)換后恢復(fù) 與以前的差別很大。而且看了一下matlab自帶的rgb2ycbcr和ycbcr2rgb并且用了這兩個(gè)函數(shù)測(cè)試后圖像是一樣的。很是不解
實(shí)驗(yàn)得到的圖為: original imagecomposition image
將源圖像分為RGB三個(gè)部分。
Red ComponentGreen ComponentBlue Component
將源圖像的RGB分量轉(zhuǎn)換到Y(jié)CbCr上后又轉(zhuǎn)化為RGB得到的三個(gè)分量差異的圖。
紅色分量差異綠色分量差異藍(lán)色分量差異 公式轉(zhuǎn)換的y分量公式轉(zhuǎn)換的cb分量公式轉(zhuǎn)換的Cr分量matlab函數(shù)轉(zhuǎn)換的y分量matlab函數(shù)轉(zhuǎn)換的Cb分量matlab函數(shù)轉(zhuǎn)換的Cr分量
從最后的結(jié)果可以看出原圖轉(zhuǎn)換到Y(jié)CbCr空間之后又轉(zhuǎn)化回來(lái)得到的圖像與原圖不符….RGB轉(zhuǎn)換為YCrCb后又轉(zhuǎn)換為RGB的圖像
以上只是我自己用Matlab做過的部分問題，還有一些比如場(chǎng)聲源定位Music算法仿真，DES加密算法。神經(jīng)網(wǎng)絡(luò)濾波器等等，鑒于代碼量太大，就不再敖述。二 Simulink學(xué)習(xí)心得與編程實(shí)踐
SIMILINK模塊庫(kù)按功能進(jìn)行分類，包括以下8類子庫(kù)： Continuous（連續(xù)模塊）Discrete（離散模塊）
Function&Tables（函數(shù)和平臺(tái)模塊）Math（數(shù)學(xué)模塊）
Nonlinear（非線性模塊）
Signals&Systems（信號(hào)和系統(tǒng)模塊）Sinks（接收器模塊）Sources（輸入源模塊）
對(duì)于simulink來(lái)說(shuō)其實(shí)沒什么好說(shuō)的，需要什么就把什么拖到Model里面，連線之后設(shè)定參數(shù)。不是很麻煩的。但是重要的是對(duì)于系統(tǒng)的設(shè)計(jì)與仿真思想才是最重要的。而且里面的S-Function也是比較好的一個(gè)設(shè)計(jì)，用S-Function可以自定義一些組件，使自己的仿真看上去更加清晰與精簡(jiǎn)。
Simulink的工作方式：
（1）模塊內(nèi)的參數(shù)值首先會(huì)送到Matlab中進(jìn)行計(jì)算，得到的參數(shù)值會(huì)用來(lái)當(dāng)做以后需要調(diào)用的參數(shù)。
（2）模型系統(tǒng)中的各個(gè)層級(jí)將被平展開來(lái)，每一個(gè)子系統(tǒng)將被相應(yīng)的模塊所代替。（3）模塊按被處理的順序排列，此時(shí)代數(shù)回路結(jié)構(gòu)也被檢查出來(lái)，此種排列產(chǎn)生一個(gè)列表，以確保具有代數(shù)回路的模塊驅(qū)動(dòng)輸入的模塊被更新后才更新。
（4）檢查塊之間的鏈接，是否每一個(gè)塊的輸出端口與它所連接的模塊輸入端口有相同的信號(hào)寬度。
現(xiàn)在可以準(zhǔn)備執(zhí)行仿真操作，仿真時(shí)使用數(shù)值迭代求的的結(jié)果，每種數(shù)值積分模型提供的連續(xù)狀態(tài)的微分能力。
Simulink中的模型都是分級(jí)的，因此可以通過自上而下或者自下而上的方式建立模型。定義了一個(gè)模型以后，就可以通過Simulink的菜單或者在Matlab的Command中輸入命令進(jìn)行仿真。
關(guān)于學(xué)習(xí)的心得就寫到這里了，接下來(lái)是一些自己做的仿真。
1：觀察一個(gè)信號(hào)與積分之后的區(qū)別，目的是了解Scope的功能和用法。
仿真結(jié)果，左邊為原始信號(hào)，右邊為積分后的信號(hào)。
用XYGraph看到的對(duì)比
Rossler吸引子產(chǎn)生仿真：
Rossler 吸引子產(chǎn)生是服從下面這個(gè)方程的 1
212 331
仿真圖 ???(x?x)x??x?axx??b?x(x?c)xa?b?0.2c?5.7
仿真結(jié)果
LMS自適應(yīng)濾波器： 仿真圖：
遇到的問題
問題的解決辦法：
將Spectrum Scope中的buffer input打勾就行了。
仿真結(jié)果圖形：
1：Spectrum Scope的輸出圖形
2：Vector Scope的輸出結(jié)果
3：Scope的輸出結(jié)果圖形
第五篇：MATLAB實(shí)驗(yàn)六《SIMULINK交互式仿真》
《計(jì)算機(jī)仿真及應(yīng)用》實(shí)驗(yàn)教案
實(shí)驗(yàn)六 SIMULINK交互式仿真
一、實(shí)驗(yàn)?zāi)康?/p>
1、熟悉SIMULINK交互式仿真集成環(huán)境。
2、掌握連續(xù)時(shí)間系統(tǒng)建模的方法。
二、實(shí)驗(yàn)主要儀器與設(shè)備
裝配有MATLAB7.6軟件的計(jì)算機(jī)
三、預(yù)習(xí)要求
做實(shí)驗(yàn)前必須認(rèn)真復(fù)習(xí)第七章SIMULINK交互式仿真集成環(huán)境。
四、實(shí)驗(yàn)內(nèi)容及實(shí)驗(yàn)步驟
示的SIMULINK模塊庫(kù)瀏覽器。
《計(jì)算機(jī)仿真及應(yīng)用》實(shí)驗(yàn)教案
圖7.2 SIMULINK模塊庫(kù)瀏覽器
4、開啟空白（新建）模型窗
單機(jī)Simulink模塊庫(kù)瀏覽器工具條上的圖標(biāo)“去的）。
5、從模塊庫(kù)復(fù)制所需模塊到空白（新建）模型窗（見圖7.3）
把連續(xù)模塊子庫(kù)中的積分器
”；或在選擇下拉菜單項(xiàng){File：
New}，引出如圖7.3所示的空白模型窗（注：窗中的模塊是后一個(gè)步驟復(fù)制進(jìn)
《計(jì)算機(jī)仿真及應(yīng)用》實(shí)驗(yàn)教案
7.3 制進(jìn)庫(kù)模塊的新建模型窗
6、新建模型窗中的模型再?gòu)?fù)制
在該新建模型窗中，可以通過再?gòu)?fù)制，產(chǎn)生建模所需的2個(gè)積分模塊和2個(gè)增益 模塊。具體操作方法是：按住[Ctrl]鍵，用鼠標(biāo)“點(diǎn)亮并拖拉”積分模塊
見一個(gè)所需的增益模塊，可采用類似方法獲得。
7、模塊連接線的形成方法隨信號(hào)起始端位置不同而不同，具體如下。
起始端位于模塊輸出口的信號(hào)線生成法。使光標(biāo)靠近模塊輸出口；待光標(biāo)變 為“單線十字叉”時(shí)，按下鼠標(biāo)左鍵；移動(dòng)十字叉，拖出一根“虛連線”；光標(biāo) 與另一個(gè)模塊輸出口靠近到一定程度，單十字叉變?yōu)殡p十字叉；放開鼠標(biāo)左鍵，”虛連線”便變?yōu)閹Ъ^的信號(hào)連線。
起始端位于其他信號(hào)線上的信號(hào)線生成法。使光標(biāo)置于已經(jīng)存在的其他信號(hào)線 之上；按下鼠標(biāo)右鍵，光標(biāo)變?yōu)椤皢尉€十字叉”；運(yùn)動(dòng)鼠標(biāo)，引出”虛連線”；當(dāng) 鼠標(biāo)與待連接模塊輸入口靠得足夠近，單十字叉變?yōu)殡p十字叉；放開鼠標(biāo)右鍵，“虛連線”便變?yōu)閹Ъ^的信號(hào)連線。
雙擊相應(yīng)的信號(hào)線，給信號(hào)線作標(biāo)注，如x’’,x’,x等。
《計(jì)算機(jī)仿真及應(yīng)用》實(shí)驗(yàn)教案
圖7.4
已構(gòu)建完成的新建模型窗
8、根據(jù)理論數(shù)學(xué)模型設(shè)置模塊參數(shù)
為使構(gòu)造的Simulink模型與理論數(shù)學(xué)模型一致，須對(duì)模塊進(jìn)行如下參數(shù)設(shè)置。
設(shè)置增益模塊
《計(jì)算機(jī)仿真及應(yīng)用》實(shí)驗(yàn)教案
修改求和模塊輸入口的代數(shù)符號(hào)：雙擊求和模塊，引出如圖7.1-6所示的參數(shù)設(shè)置窗； 把符號(hào)列表欄中的默認(rèn)符號(hào)（++）修改成代數(shù)符號(hào)（--）；單擊【OK】鍵，完成設(shè)置；
圖7.6 改變輸出入口符號(hào)的求和模塊參數(shù)設(shè)置窗
據(jù)初始位移x(0)=0.05m對(duì)積分模塊
《計(jì)算機(jī)仿真及應(yīng)用》實(shí)驗(yàn)教案
圖7.7 實(shí)現(xiàn)初始化位移0.05設(shè)置的Integrator1設(shè)置窗
9仿真運(yùn)行參數(shù)采用默認(rèn)解算器“ode45”、默認(rèn)“變步長(zhǎng)”和默認(rèn)仿真終止時(shí)間10.10把新建模型保存為exm070101.mdl.11試運(yùn)行，以便發(fā)現(xiàn)問題加以改善。
雙擊
”仿真啟動(dòng)鍵，使該模型運(yùn)行；在示波器上呈現(xiàn)
圖7.8 坐標(biāo)范圍設(shè)置不當(dāng)時(shí)的信號(hào)
《計(jì)算機(jī)仿真及應(yīng)用》實(shí)驗(yàn)教案
單擊Scope顯示窗上的縱坐標(biāo)范圍自動(dòng)設(shè)置圖標(biāo)“”，示波器顯示窗改變?yōu)閳D7.9所
示。在顯示窗中，可以看到位移x（t）的變化曲線。同時(shí)可以發(fā)現(xiàn)：縱坐標(biāo)的適當(dāng)范 圍大致在【-0.06,0.06】；仿真時(shí)間取【0,5】即可；顯示的曲線不夠光滑。
圖7.9 采用軸自動(dòng)設(shè)置功能的信號(hào)顯示
12、據(jù)試運(yùn)行結(jié)果，進(jìn)行仿真參數(shù)的再設(shè)置
示波器縱坐標(biāo)設(shè)置：用鼠標(biāo)單擊示波器的黑色顯示屏，在彈出菜單中選擇Axes Properties，引出縱坐標(biāo)設(shè)置對(duì)話窗7.10所示；把縱坐標(biāo)的下限、上限分別設(shè)置為（-0.06）和(+0.06);單 擊【OK】，完成設(shè)置。
示波器時(shí)間顯示范圍的修改:單擊示波器的參數(shù)設(shè)置圖標(biāo)“
”，引出示波器參數(shù)設(shè)置窗；
在General卡片的Axes區(qū)的Time range欄中，填寫5或auto；單擊【OK】，完成設(shè)置。
圖7.10 對(duì)顯示屏的縱坐標(biāo)范圍進(jìn)行設(shè)置
7.11 對(duì)示波器時(shí)間顯示范圍的設(shè)置
·仿真終止時(shí)間最簡(jiǎn)捷的修改方法：在exam070101模型窗“仿真終止時(shí) 間”欄“
”中的默認(rèn)值10改變?yōu)?。
《計(jì)算機(jī)仿真及應(yīng)用》實(shí)驗(yàn)教案
·顯示曲線的光滑化設(shè)置：選中exam070101模型窗的下拉菜單項(xiàng)SimulationConfiguration Parameters，引出仿真參數(shù)配置窗，如圖7.12所示；再在該窗左側(cè)的選擇欄中，選中Data Import/Export項(xiàng)，與之相應(yīng)的參數(shù)設(shè)置欄便出現(xiàn)在窗口的右側(cè)；把右半窗下方Save options 區(qū)中Refine factor欄中的默認(rèn)值1改為5；單擊OK，完成設(shè)置。
7.12 通過仿真參數(shù)配置窗設(shè)置輸出光滑因子
《計(jì)算機(jī)仿真及應(yīng)用》實(shí)驗(yàn)教案
完成以上修改后的模型窗如圖7.13所示。再運(yùn)行exam070101,可得比較滿意的位移變化曲線（見圖7.14所示）。順便指出：模型運(yùn)行后，在模塊輸出口的信號(hào)線上會(huì)出現(xiàn)double字樣。這 表示該信號(hào)是采用“雙精度”類型數(shù)據(jù)運(yùn)算。如果用戶不希望這類標(biāo)識(shí)出現(xiàn)，那么應(yīng)取消對(duì)下 拉菜單項(xiàng)FormatPort signal DisplaysPort Data Types的勾選。
圖7.13
仿真參數(shù)調(diào)整運(yùn)行后的exam070101模型
在模型窗中，x’’,x’,x等信號(hào)名稱是模型創(chuàng)建者根據(jù)需要寫入的。標(biāo)識(shí)信號(hào)名稱的操作方法 是：用鼠標(biāo)雙擊信號(hào)附近的適當(dāng)位置后，就會(huì)出現(xiàn)一個(gè)虛線框，該虛線框中允許輸入任何 文字。
13、仿真結(jié)果顯示
原本比較稀疏的解計(jì)算器數(shù)據(jù)，經(jīng)設(shè)定的“光滑因子”下的插值，增加了用于描繪曲線的 數(shù)據(jù)點(diǎn)數(shù)，因此使示波器顯示出更為光滑的曲線，如圖7.14所示。
《計(jì)算機(jī)仿真及應(yīng)用》實(shí)驗(yàn)教案
圖7.14 適當(dāng)?shù)仫@示仿真所得的位移變化曲線
五、實(shí)驗(yàn)小結(jié)
點(diǎn)擊咨詢 