MATLAB基礎(chǔ)及應(yīng)用課件第78章

1、MATLAB基礎(chǔ)及應(yīng)用（第5版）MATLAB基礎(chǔ)及應(yīng)用（第5版）第7章 Simulink仿真 本章要點 Simulink系統(tǒng)的基本模塊 仿真模型的編輯 仿真參數(shù)的設(shè)置第7章 Simulink仿真 本章要點7.1 認(rèn)識Simulink Simulink 是MATLAB環(huán)境下對動態(tài)系統(tǒng)進(jìn)行建模、仿真和分析的一個軟件包。該系統(tǒng)的兩個主要功能就是Simu（仿真）和Link（連接）。 一、Simulink的啟動和退出 1Simulink的啟動在MATLAB操作桌面下，單擊工具欄中的“Simulink Library”按鈕 ；或在命令窗口鍵入命令“simulink”，單擊【Enter】鍵，都會彈出一個名為

2、“Simulink Start Page”的對話框。7.1 認(rèn)識Simulink Simulink 是MA圖7-1 “Simulink Start Page”對話框圖7-1 “Simulink Start Page”對話框2模型的創(chuàng)建 “Simulink Start Page”對話框有“New”（新建）和“Examples”（實例）兩個選項卡，其中“New”（新建）選項卡可以選擇新建“Blank Model”（空白模型）、“Blank Library”（空白庫）、“Blank Project”（空白項目）等；“Examples”（實例）選項卡是一些仿真項目實例，可供參考。 單擊“Simulin

3、k Start Page”對話框中的“Blank Model”（空白模型），會彈出一個名為Untitled（無標(biāo)題）的空白窗口，所有控制模塊都可以創(chuàng)建在這個窗口中，如圖7-2所示。 2模型的創(chuàng)建 “Simulink Star圖7-2 新建模型圖7-2 新建模型 Simulink模塊庫瀏覽器 單擊“新建模型”窗口工具欄中的按鈕 ，可以打開Simulink模塊庫瀏覽器。 Simulink模塊庫瀏覽器 單擊“新建模型1信號源模塊（Sources） 信號源模塊是輸入模塊，提供各種類型的信號或數(shù)據(jù)。1信號源模塊（Sources） 信號源模塊Signal Generator信號發(fā)生器Chirp Signa

4、l快速正弦掃描信號模塊功能模塊功能In1創(chuàng)建輸入端Ground接地From File從文件讀數(shù)據(jù)From Workspace從工作空間讀數(shù)據(jù) Constant常數(shù)Enumerated Constant枚舉類型常數(shù)Signal Builder信號生成器Ramp斜波信號Step階躍信號Sine Wave正弦波Signal Generator信號發(fā)生器Chirp Sig Random Number隨機(jī)信號Uniform Random Number均勻隨機(jī)信號Band-Limited White Noise帶限白噪聲Digital Clock數(shù)字時鐘Pulse Generator脈沖發(fā)生器Repeat

5、ing Sequence重復(fù)序列信號Repeating Sequence Stair重復(fù)階梯序列信號Repeating Sequence Interpolated重復(fù)曲線序列信號Clock當(dāng)前時間Digital Clock數(shù)字時間Counter Free-Running自動運行計數(shù)器（溢出時自動清零）Counter Limited有限計數(shù)器（可自定義計數(shù)上限） Random Number隨機(jī)信號Uniform Rand2輸出模塊（Sinks） 輸出模塊也可稱為接收模塊，用于顯示仿真結(jié)果或輸出仿真數(shù)據(jù)。2輸出模塊（Sinks） 輸出模塊也可稱為接2輸出模塊（Sinks）模塊功能模塊功能Outl創(chuàng)

6、建輸出端Terminator通用終端To File輸出到文件To Workspace輸出到工作空間Scope示波器Floating Scope浮點格式的示波器XY GraphXY關(guān)系圖Display實時數(shù)值顯示Stop Simulation輸出不為0時停止仿真2輸出模塊（Sinks）模塊功能模塊功能Outl創(chuàng)建輸出端3連續(xù)系統(tǒng)模塊（Continuous） 連續(xù)系統(tǒng)模塊提供積分、導(dǎo)數(shù)等連續(xù)系統(tǒng)仿真的常用子模塊。3連續(xù)系統(tǒng)模塊（Continuous） 連續(xù)系連續(xù)系統(tǒng)模塊的功能 模塊功能模塊功能Integrator不定積分Integrator Limited定積分Integrator Second-

7、Order二重不定積分Integrator Second-Order Limited二重定積分State-Space狀態(tài)方程Transfer Fcn傳遞函數(shù)Zero-Pole零極點PID ControllerPID控制器PID Controller（2D OF）雙自由度PID控制器Transport Delay傳輸延時Variable Time Delay可變時間延時Variable Transport Delay可變傳輸延時Derivative導(dǎo)數(shù)連續(xù)系統(tǒng)模塊的功能 模塊功能模塊功能Integrator不定8數(shù)學(xué)運算模塊（Math Operations） 數(shù)學(xué)運算模塊提供了基本數(shù)學(xué)運算函數(shù)、

8、三角函數(shù)、復(fù)數(shù)運算函數(shù)以及矩陣運算函數(shù)。 在Simulink模塊庫瀏覽器左側(cè)的資源列表選中“Math Operations”、或雙擊右側(cè)窗口中的Math Operations模塊圖標(biāo)，如圖7-7所示。8數(shù)學(xué)運算模塊（Math Operations） MATLAB基礎(chǔ)及應(yīng)用課件第7-8章表7-4 數(shù)學(xué)運算模塊的功能模塊功能模塊功能Sum求和Add加法Subtract減法Sum of Elements各元素的累加和Bias偏移量WeightedSampleTimeMath對信號經(jīng)過加權(quán)時間采樣的運算Gain常數(shù)增益Slider Gain可變增益Product乘法Divide除法Product of

9、 Elements各元素的累積Dot Product點乘表7-4 數(shù)學(xué)運算模塊的功能模塊功能模塊功能Sum求和Ad表7-4 數(shù)學(xué)運算模塊的功能模塊功能模塊功能Sign符號函數(shù)Abs求絕對值Unary Minus單一元素的負(fù)數(shù)Math Function數(shù)學(xué)運算函數(shù)Rounding Function取整函數(shù)Polynomial多項式MinMax求最大值或最小值MinMaxRunningResettable求最大值或最小值，帶復(fù)位功能Trigonometric function三角函數(shù)Sine Wave Function正弦波形函數(shù)Algebraic Constraint代數(shù)上的約束常量Sqrt平方

10、根Signal Sqrt信號的的平方根Reciprocal Sqrt平方根的倒數(shù)表7-4 數(shù)學(xué)運算模塊的功能模塊功能模塊功能Sign符號函7.2 Simulink模塊操作 一、模塊的編輯處理 1模塊的操作（1）添加模塊 當(dāng)要把一個模塊添加到模型中，先在Simulink模塊庫中找到它，然后直接將這個模塊拖入模型窗口中即可。（2）選取模塊 當(dāng)模塊已經(jīng)位于模型窗口中時，只要用鼠標(biāo)在模塊上單擊就可以選中該模塊，這時模塊的四角上出現(xiàn)一些黑色的小方塊，這些小方塊就是該模塊的關(guān)鍵點，拖動這些黑色小方塊可以改變模塊的大小。7.2 Simulink模塊操作 一、模塊的編輯處理 2模塊位置和外形的調(diào)整（1）調(diào)整模

11、塊位置（2）調(diào)整模塊大?。?）調(diào)整模塊方向（4）調(diào)整模塊顏色和效果3模塊名的處理（1）改變模塊名（2）隱藏模塊名（3）改變模塊名的位置4復(fù)制和刪除模塊（1）復(fù)制模塊（2）刪除模塊2模塊位置和外形的調(diào)整1模塊參數(shù)的設(shè)置 Simulink中幾乎所有模塊的參數(shù)（Parameters）都允許用戶進(jìn)行設(shè)置，雙擊要設(shè)置的模塊就可以打開模塊參數(shù)設(shè)置對話框，不同模塊參數(shù)設(shè)置對話框的項目會不同。 例如“Source”模塊庫的“Step”模塊參數(shù)對話框，如圖7-8所示。7.2.2 模塊屬性和參數(shù)的設(shè)置1模塊參數(shù)的設(shè)置7.2.2 模塊屬性和參數(shù)的設(shè)置圖7-8 模塊參數(shù)設(shè)置對話框圖7-8 模塊參數(shù)設(shè)置對話框2模塊屬性

12、的設(shè)置 選定要設(shè)置屬性的模塊，然后單擊“Diagram”“Properties”選項；或單擊右鍵在彈出的菜單中選擇“Properties”選項，將得到如圖7-9所示的屬性設(shè)置對話框。該對話框有三個標(biāo)簽Description（說明）、Priority（優(yōu)先級）和Tag（標(biāo)記）。根據(jù)需要設(shè)定的基本屬性如下：1）Description（說明）：對該模塊在模型中的用法進(jìn)行說明。2）Priority（優(yōu)先級）：規(guī)定該模塊在模型中相對于其它模塊的優(yōu)先順序，優(yōu)先級的數(shù)值必須是整數(shù)（可以是負(fù)數(shù)），該數(shù)值越小，優(yōu)先級越高。3）Tag（標(biāo)記）：用戶為模塊添加的文本格式的標(biāo)記。2模塊屬性的設(shè)置 選定要設(shè)置屬性的模塊

13、，然后單圖7-9 模塊屬性設(shè)置對話框圖7-9 模塊屬性設(shè)置對話框三、模塊間的連線 1連接兩個模塊 從一個模塊的輸出端連到另一個模塊的輸入端是Simulink仿真最基本的操作。方法是先移動鼠標(biāo)指向模塊的輸出端，鼠標(biāo)的箭頭會變成十字形光標(biāo)，這時按住鼠標(biāo)左鍵，拖動鼠標(biāo)到另一個模塊的輸入端，當(dāng)十字形光標(biāo)出現(xiàn)“重影”時，釋放鼠標(biāo)即完成了連接。2模塊間連線的調(diào)整 用鼠標(biāo)單擊連線，可以選中該連線。這時會看到線上的一些黑色小方塊，這些是連線的關(guān)鍵點。用鼠標(biāo)按住關(guān)鍵點，拖動即可以改變連線的方向。 三、模塊間的連線 1連接兩個模塊3連線的分支 仿真時經(jīng)常會碰到需要把信號輸送到不同的接收端的情況，這時就需要分支結(jié)構(gòu)

14、的連線?？梢韵冗B好一條線，然后把鼠標(biāo)移到支線的起點位置，先按下Ctrl鍵，然后按住鼠標(biāo)，將連線拖到目標(biāo)模塊，松開鼠標(biāo)和Ctrl鍵即可。4刪除連線 如果想要刪除某條連線，可單擊要刪除的連線，此時連線上出現(xiàn)標(biāo)記點，表示該連線已經(jīng)被選中，然后單擊模型窗口工具欄中的【剪切】按鈕或者直接單擊鍵盤上的Delete鍵，即可刪除該連線。3連線的分支5標(biāo)注連線 雙擊需要標(biāo)注的連線，可以看到一個文本框，在里面輸入標(biāo)注文字單擊回車鍵確定即可。用鼠標(biāo)拖動文本框還可以改變標(biāo)識的位置。另外，在FormatPort/Signal Displays菜單下還有與連線標(biāo)注有關(guān)的命令：1）Sample Time Color：將采樣

15、時間不同的模塊和連線用不同的顏色顯示。2）Port Data Types：在連線上顯示傳輸數(shù)據(jù)的類型，如double、int32等。5標(biāo)注連線7.3 仿真模型的參數(shù)設(shè)置 在仿真系統(tǒng)設(shè)計過程中，事先還必須對仿真算法、輸出模式等各種模型參數(shù)進(jìn)行設(shè)置。 單擊模型窗口菜單欄“Simulation”“Model Configuration Parameters”選項，將出現(xiàn)仿真參數(shù)配置窗口，如圖7-11所示。7.3 仿真模型的參數(shù)設(shè)置 在仿真系統(tǒng)設(shè)計過程中圖7-11 仿真模型參數(shù)配置窗口圖7-11 仿真模型參數(shù)配置窗口 仿真參數(shù)配置窗口主要分為7個選項卡： Solver（解題器）、Data Import

16、/Export（數(shù)據(jù)輸入/輸出）、Math and Data Types（數(shù)學(xué)和數(shù)據(jù)類型）、Diagnostics（診斷）、Hardware Implemention（硬件工具）、Model Referencing（模型引用）和Simulation Target（仿真輸出文件的格式）， 其中Solver、Data Import/Export和Diagnostics三項經(jīng)常用到。 仿真參數(shù)配置窗口主要分為7個選項卡：7.3.1 Solver選項卡（1）Simulation time（仿真時間） 設(shè)置仿真起始時間和停止時間。在Star time 和Stop time兩個輸入框內(nèi)，直接輸入數(shù)值，時間

17、單位是秒。但要注意，這里的時間只是計算機(jī)對時間的一種表示，執(zhí)行一次仿真所需的時間依賴于很多因素，包括計算機(jī)的時鐘頻率、模型的復(fù)雜程度、解題器及其步長等。（2）Solver selection（算法選擇） 仿真解題器的選擇操作。Type（算法類型）分為：Variable-step（變步長算法）和Fixed-step（固定步長算法）兩種。7.3.1 Solver選項卡（1）Simulation 7.3.2 Data Import/Export選項卡7.3.2 Data Import/Export選項卡（1）Load from workspace 從工作空間中載入數(shù)據(jù)。在仿真過程中，如果模型中有輸入

18、端口（In模塊），可從工作空間直接把數(shù)據(jù)載入到輸入端口。（2）Save to workspace or file 將輸出保存到工作空間或文件，一般情況選擇的輸出選項有：Time（時間），States（狀態(tài)），Output（輸出端口）和Final state（最終狀態(tài)）。（1）Load from workspace7.3.3 Diagnostics選項卡7.3.3 Diagnostics選項卡 在選項卡的中間列出了仿真過程中可能出現(xiàn)的一些事件，用戶可以在相應(yīng)事件右邊的下拉框中根據(jù)需要選擇系統(tǒng)的反應(yīng)（即采取的操作）。反應(yīng)的類型有以下幾種：1）None：不做任何反應(yīng)，不影響程序運行。2）Warnin

19、g：顯示警告信息，不影響程序的運行。3）Error：顯示錯誤信息，中止運行的程序。 在選項卡的中間列出了仿真過程中可能出現(xiàn)的一些事件MATLAB基礎(chǔ)及應(yīng)用（第5版）MATLAB基礎(chǔ)及應(yīng)用（第5版）第8章 MATLAB綜合實訓(xùn)本章要點 數(shù)字圖像的幾何操作 圖像增強(qiáng) 圖像置亂與恢復(fù) 數(shù)字水印算法的實現(xiàn) 第8章 MATLAB綜合實訓(xùn)本章要點8.1 數(shù)字圖像的幾何操作 一、項目說明 1項目要求（1）圖像的文件操作（2）圖像的剪裁操作（3）圖像的大小調(diào)整（4）圖像的旋轉(zhuǎn)（5）圖像的插值操作（6）拓展要求 8.1 數(shù)字圖像的幾何操作 一、項目說明 2實施步驟1）討論、研究項目要求，明確項目內(nèi)容；2）學(xué)習(xí)項

20、目設(shè)計提示，練習(xí)相關(guān)函數(shù)的用法；3）編寫程序，實現(xiàn)項目功能；4）圖形用戶界面的設(shè)計與實現(xiàn)；5）項目演示。 2實施步驟二、項目設(shè)計提示 函數(shù)功能函數(shù)格式說明讀入圖像A=imread（文件名，圖像格式)從圖像文件中讀入圖像數(shù)據(jù)到變量A 中。圖像格式包括bmp、tif、jpg、png等A=imread（文件名.擴(kuò)展名)圖像顯示imshow（圖像變量名）顯示數(shù)據(jù)類型為uint8的灰度圖像imshow（圖像變量名，map）顯示數(shù)據(jù)類型為double的索引圖像。map為顏色圖圖像剪裁變量名=imcrop（圖像變量名）交互式剪裁。拖動鼠標(biāo)選定剪裁區(qū)域變量名=imcrop（圖像變量名，x y w h）非交互式

21、剪裁。x和y為剪裁起點；w為寬、h度為高度二、項目設(shè)計提示 函數(shù)功能函數(shù)格式說明讀入圖像A=imrea二、項目設(shè)計提示 函數(shù)功能函數(shù)格式說明圖像大小調(diào)整變量名=imresize（圖像變量名，調(diào)整系數(shù)，參數(shù)）比例調(diào)整。調(diào)整系數(shù)為“放大”或“縮小”的倍數(shù)；參數(shù)表示插值方法，包括nearest（默認(rèn)）、bilinear、bicubic等變量名=imresize（圖像變量名，m n，參數(shù)）拉伸調(diào)整。得到尺寸為mn的圖像，參數(shù)表示插值方法圖像旋轉(zhuǎn)變量名=imrotate（圖像變量名，角度，參數(shù)）角度為正，逆時針旋轉(zhuǎn)；角度為負(fù)，順時針旋轉(zhuǎn)。參數(shù)表示插值方法變量名=imrotate（圖像變量名，角度，參數(shù)，

22、crop）crop可以將旋轉(zhuǎn)后超出原圖像的部分剪裁掉，得到與原圖像大小相同的部分圖像插值變量名=interp2（圖像變量名，n，參數(shù)）n為插值次數(shù)。每插值一次，在兩個元素間插入一個點。參數(shù)表示插值方法二、項目設(shè)計提示 函數(shù)功能函數(shù)格式說明圖像大小調(diào)整變量名=i【例8-1】對一幅圖像進(jìn)行剪裁、放大2倍和旋轉(zhuǎn)15的處理。clearclcX=imread(cameraman.tif); %調(diào)入原圖像subplot(221);imshow(X); title(原圖像); %顯示原圖像X1=imcrop(X,60 40 100 90); %60 40 100 90為剪裁區(qū)域subplot(222);im

23、show(X1); title(原圖像剪裁塊);【例8-1】對一幅圖像進(jìn)行剪裁、放大2倍和旋轉(zhuǎn)15的處理。X2=imresize(X,2,bilinear); %參數(shù) bilinear為雙線性插值subplot(223);imshow(X2); title(拉伸調(diào)整的原圖像);X3=imrotate(X,15,bilinear,crop); %參數(shù) crop為返回同樣大小的圖像subplot(224);imshow(X3); title(旋轉(zhuǎn)的原圖像);X2=imresize(X,2,bilinear); MATLAB基礎(chǔ)及應(yīng)用課件第7-8章三、項目評價 項目評價是在教師的主持下，通過項目負(fù)責(zé)

24、人的講解演示，評估項目的完成情況，評價內(nèi)容如下：1）圖像的視覺效果；2）處理圖像與原始圖像之間的偏離程度；3）圖形用戶界面是否友好。三、項目評價 項目評價是在教師的主持下，通過項目負(fù)責(zé)人的講解8.2 數(shù)字圖像增強(qiáng) 圖像增強(qiáng)就是對圖像進(jìn)行加工，以得到對具體應(yīng)用來說視覺效果更“好”，更“有用”的圖像。 一、 項目說明 1項目要求 （1）圖像測試 （2）圖像添加噪聲8.2 數(shù)字圖像增強(qiáng) 圖像增強(qiáng)就是對圖（3）圖像的對比度調(diào)整（4）圖像求反（5）圖像平滑（6）拓展要求實現(xiàn)圖像增強(qiáng)操作的圖形用戶界面設(shè)計，包括“圖像測試”、“添加噪聲”、“對比度調(diào)整”、“圖像求反”和“圖像平滑”等基本功能。（3）圖像的對

25、比度調(diào)整2實施步驟1）討論、研究項目要求，明確項目內(nèi)容；2）學(xué)習(xí)項目設(shè)計提示，分析算法；3）仿真算法，完成項目；4）項目演示、講解設(shè)計方案，完成項目評價 2實施步驟二、項目設(shè)計提示 （1）圖像測試 先將教師提供的圖像讀入工作空間，記下圖像尺寸和數(shù)據(jù)類型，利用hist函數(shù)繪制概率分布圖，注意橫坐標(biāo)的數(shù)值范圍和步長；使用統(tǒng)計函數(shù)mean計算圖像均值。（2）圖像添加噪聲 噪聲信號可以由隨機(jī)矩陣函數(shù)實現(xiàn)，直接與圖像相加即可，注意相加后的值不要超過圖像的數(shù)據(jù)范圍，也不要小于0；如果加噪的效果不明顯，可增大隨機(jī)矩陣的元素數(shù)值，例如乘以一個大于1的系數(shù)。二、項目設(shè)計提示 （1）圖像測試（3）圖像的對比度調(diào)整

26、（4）圖像求反 （5）圖像平滑 將圖像中除四個邊（第1行、末行、第1列、末列）的元素外，每個元素都用其8鄰域與其自身相加之和的平均值代替，這種平滑方法稱為均值濾波；如果是用這9個元素的中間值代替就稱為中值濾波。元素的8鄰域指在該元素上、下、左、右、對角和反對角位置上的八個元素。 （3）圖像的對比度調(diào)整【例8-2】使用均值濾波方法的圖像平滑clearclcI=imread(cameraman.tif); %讀入圖像，I為uin8類型subplot(131);imshow(I); title(原圖像); A=round(randn(256,256)*15); %隨機(jī)矩陣放大15倍，生 成噪聲I1=

27、I+uint8(A); %加入噪聲，A由double轉(zhuǎn)為uin8subplot(132);imshow(I1); title(加噪聲后的圖像);【例8-2】使用均值濾波方法的圖像平滑clearI2=double(I1);for i=2:255 for j=2:255 temp=0; for m=1:3 %均值濾波 for n=1:3 temp=temp+I2(i+m-2,j+n-2); end end I2(i,j)=round(temp/9); endendsubplot(133);imshow(uint8(I2); title(平滑操作后的圖像);I2=double(I1);MATLAB基

28、礎(chǔ)及應(yīng)用課件第7-8章 三、項目評價 圖像處理的效果與圖像本身有關(guān)，不同的圖像對同一種算法會有不同的表現(xiàn)?？蓮囊韵聨讉€方面評價：（1）主觀評價 觀察圖像處理的視覺效果，綜合多個觀察者的評價意見。（2）算法評價 主要評價算法是否正確、實現(xiàn)功能是否符合項目要求、有無功能擴(kuò)展、程序可讀性如何、算法是否簡練、編寫是否規(guī)范、程序運行效率如何等方面。（3）演示過程評價 主要評價演示效果如何、對算法理解程度如何、回答問題是否準(zhǔn)確、語言是否流暢等方面。 三、項目評價 圖像處理的效果與圖像本身有關(guān)，不同的8.3 數(shù)字圖像置亂 圖像置亂是對圖像元素的位置按照某種規(guī)律（密鑰）重新排列，達(dá)到隱藏圖像真實內(nèi)容的目的，置

29、亂后的圖像可以根據(jù)同樣的規(guī)律（密鑰）恢復(fù)。一、項目說明1項目要求（1）圖像的客觀評價（2）圖像置亂（3）置亂圖像恢復(fù)8.3 數(shù)字圖像置亂 圖像置亂是對圖像元素的位置按照某種2實施步驟1）討論、研究項目要求，明確項目內(nèi)容；2）學(xué)習(xí)項目設(shè)計提示，分析算法；3）仿真算法，討論、評價、修改算法，并確定一個項目實現(xiàn)方案；4）實現(xiàn)方案，完成項目；5）撰寫項目報告；6）項目演示、講解設(shè)計方案，完成項目評價。2實施步驟二、項目設(shè)計提示 （1）圖像的客觀評價 相關(guān)系數(shù)可以使用二維相關(guān)函數(shù)corr2實現(xiàn)；峰值信噪比的計算公式。式中： 大小為的原圖像 處理后的圖像二、項目設(shè)計提示 （1）圖像的客觀評價（2）圖像置亂

30、方案一：類似于撲克洗牌。將圖像的全部偶數(shù)行提出，構(gòu)成圖像1；再提出全部奇數(shù)行構(gòu)成圖像2，連接圖像1和圖像2，構(gòu)成與原圖像同樣尺寸的新圖像；對新圖像的列進(jìn)行同樣的操作。重復(fù)n次，滿足要求即可，n為密鑰。方案二：先將圖像變成單列矩陣A，使每一個圖像元素都有一個唯一的地址（單列矩陣的行號），再利用randperm（隨機(jī)排列整數(shù)矩陣）函數(shù)生成一個同樣長度的單列矩陣B，該矩陣與原圖像的尺寸信息作為密鑰保存；建立一個與A、B等長的單列矩陣C，C中存放以矩陣B的元素值為地址的矩陣A的元素。（3）置亂圖像恢復(fù) 需要知道置亂密鑰和置亂方案才能恢復(fù)出原圖像，用不同方案置亂的圖像需要用相應(yīng)的算法恢復(fù)。 （2）圖像置

31、亂【例8-3】圖像置亂算法（方案2）clearclcX=imread(cameraman.tif);subplot(121)imshow(X)title(原圖像) %顯示原圖像【例8-3】圖像置亂算法（方案2）clearA=X(:);B=randperm(65536); %生成隨機(jī)排列整數(shù) 矩陣，作為置亂密鑰C=zeros(65536,1);for i=1:65536 C(i)=A(B(i); %置亂，用密鑰重新排列 圖像矩陣endA=X(:); D=zeros(256);k=1;for i=1:256 %將單列矩陣轉(zhuǎn)換成與原圖像同樣尺寸的矩陣 for j=1:256 D(j,i)=C(k);

32、 k=k+1; endend D=zeros(256);subplot(122)imshow(uint8(D) %uint8(D)強(qiáng)制轉(zhuǎn) 換D的數(shù)據(jù)格式title(置亂后的圖像)imwrite(uint8(D),d:image.bmp,bmp) % 保存置亂后的圖像fid=fopen(d:test.bin,wb) %保存密鑰fwrite(fid,B,double)fclose(fid)subplot(122)MATLAB基礎(chǔ)及應(yīng)用課件第7-8章【例8-4】方案2的圖像恢復(fù)程序clearclcRD=imread(d:image.bmp,bmp); %讀入 置亂后的圖像subplot(121)im

33、show(RD)title(置亂后的圖像)fid=fopen(d:test.bin,rb) %讀入密鑰RB=fread(fid,65536,double);fclose(fid)【例8-4】方案2的圖像恢復(fù)程序clearRC=RD(:); %生成單列矩陣RA=zeros(65536,1);for i=1:65536 RA(RB(i)=RC(i); %恢復(fù)，用密鑰重 新排列置亂的矩陣EndX=zeros(256);k=1;RC=RD(:); for i=1:256 %將單列矩陣轉(zhuǎn)換成與原圖像 同樣尺寸的矩陣 for j=1:256 X(j,i)=RA(k); k=k+1; endendX=uin

34、t8(X);subplot(122)imshow(X)title(恢復(fù)的圖像)for i=1:256 MATLAB基礎(chǔ)及應(yīng)用課件第7-8章三、項目評價（1）客觀評價 利用峰值信噪比和相關(guān)系數(shù)等指標(biāo)評價圖像置亂效果和恢復(fù)效果。 （2）算法評價 主要評價算法實現(xiàn)方案的優(yōu)勢和缺陷，安全性如何？運算時間是多少？有無實用價值等方面。（3）項目開發(fā)文檔評價 主要評價文檔內(nèi)容是否完整、分析是否全面、結(jié)構(gòu)是否合理、語句是否通順、編輯排版是否規(guī)范等方面。（4）論述答辯過程評價 主要評價答辯態(tài)度如何、思路是否清晰、回答是否準(zhǔn)確、語言是否流暢、對算法不足方面有無認(rèn)識等方面。 三、項目評價8.4 數(shù)字水印技術(shù) 數(shù)字水

35、印就是將版權(quán)（或認(rèn)證）信息嵌入到多媒體數(shù)據(jù)中，但不影響原始數(shù)據(jù)的正常使用，目的是鑒別非法復(fù)制或盜用的數(shù)字圖像產(chǎn)品，主要用于數(shù)字產(chǎn)品的知識產(chǎn)權(quán)保護(hù)、產(chǎn)品防偽等方面。一、項目說明1項目背景介紹（1）數(shù)字水印分類1）按水印特性劃分2）按水印所附載的媒體劃分3）按水印的檢測過程劃分4）按水印的內(nèi)容劃分 8.4 數(shù)字水印技術(shù) 數(shù)字水印就是將版權(quán)（或認(rèn)證）（2）數(shù)字水印的評價1）隱蔽性是指原始圖像在嵌入水印后的改變程度?？梢杂迷紙D像與嵌入水印后的圖像進(jìn)行比較來確定其隱蔽性。2）魯棒性是指含有水印的圖像在經(jīng)過各種線性和非線性處理后，提取的水印與原始水印之間的差別。差別的大小常用來判斷待測圖像中是否存在水印

36、。3）水印容量是指水印的信息量。水印容量和魯棒性之間是相互矛盾的，水印容量的增加會帶來魯棒性的下降，對隱蔽性也有同樣的影響。 （2）數(shù)字水印的評價1）隱蔽性是指原始圖像在嵌入水印后的改變2項目要求 以一幅256256的256級灰度圖像woman作為原始圖像，以一幅256256的黑白圖像text.png的局部作為水印圖像。水印檢測不需要原始圖像，只需要密鑰；水印具有一定的隱藏性和魯棒性。（1）隱蔽性要求含有水印的圖像與原始圖像在視覺上無明顯改變、峰值信噪比大于30、相關(guān)系數(shù)大于0.95。提取的水印圖像與嵌入的水印圖像相關(guān)系數(shù)大于0.9，視覺上無明顯差別。（2）魯棒性要求含有水印的圖像經(jīng)過剪切、加噪聲、濾波、壓縮等處理后，提取的水印圖像在視覺上變化不大，與原水印的峰值信噪比大于30、相關(guān)系數(shù)大于0.85。（3）水印容量在滿足隱蔽性和魯棒性指標(biāo)的前提下，水印容量越多越好。2項目要求 以一幅256256的256級灰度圖像w3實施步驟1）討論、研究項目要求，明確項目內(nèi)容；2）檢索、閱讀參考資料，學(xué)習(xí)項目設(shè)計提示；3）仿真算法，討論、評價、修改算法，并確定一個項目實現(xiàn)方案；4）實現(xiàn)方案，完成項目；5）撰寫項目報告；6）項目演示、講解設(shè)計方案，完成項目評價。3實施步驟二、項目設(shè)計提示 （1）數(shù)字水印嵌入模型 （2）數(shù)字水印提取模型 二、項目設(shè)計提示 （1）數(shù)字水印嵌入模