第一篇：硫化氫處理實(shí)驗(yàn)報告

`

硫化氫處理實(shí)驗(yàn)報告

目

錄

一、現(xiàn)場簡介 ...........................................................................................................................2

二、現(xiàn)場實(shí)驗(yàn) ...........................................................................................................................3

2.1 第一階段實(shí)驗(yàn) 數(shù)據(jù)........................................................................................................3

2.1.1 設(shè)備安裝..................................................................................................................4

2.1.2 檢測方法.................................................................................................................4

2.1.3 第一階段實(shí)驗(yàn)數(shù)據(jù).................................................................................................4

2.2 第二階段實(shí)驗(yàn)數(shù)據(jù)........................................................................................................7

2.2.1 各檢測點(diǎn)實(shí)驗(yàn)數(shù)據(jù).................................................................................................7

2.2.2 藥劑室內(nèi)實(shí)驗(yàn).........................................................................................................8

2.3 第三階段實(shí)驗(yàn)數(shù)據(jù)........................................................................................................9

三、總結(jié) .................................................................................................................................13

一、現(xiàn)場簡介 勝臨商 56 聯(lián)合站隸屬于東勝公司勝臨分公司，位于商河?xùn)|南，屬于油水一體化處理綜合聯(lián)合站。污水量約為 500m 3 /d，現(xiàn)場硫化氫超標(biāo)嚴(yán)重，對商 56 站各節(jié)點(diǎn)硫化氫進(jìn)行了測試，兩相分離器出口檢測點(diǎn) A 的硫化氫濃度為 1700ppm，商 56-116 井排油水來液檢測點(diǎn) B 的硫化氫濃度 1600ppm，一次沉降油罐檢測點(diǎn)C 的硫化氫濃度為 2000ppm，外輸油檢測點(diǎn) D 的硫化氫濃度 200ppm，一次除油

水罐污水進(jìn)罐處檢測點(diǎn) E 的硫化氫濃度 2000ppm（檢測方法：用 5L 取樣桶，取1L 采出液，放置 30 秒，檢測液面以上氣體中的硫化氫濃度）。各節(jié)點(diǎn)的硫化氫濃度均大于安全臨界濃度 20ppm，硫化氫含量超標(biāo)嚴(yán)重。

圖 1 現(xiàn)場工藝流程及檢測點(diǎn)

二、現(xiàn)場實(shí)驗(yàn)

自 10 月 22 號進(jìn)行了現(xiàn)場實(shí)驗(yàn)，現(xiàn)場實(shí)驗(yàn)分三個階段，第一階段 10 月22 日 9 時至 10 月 27 日 14 時，在井排匯管處投加硫化氫吸收劑 1 號，第二階段10 月 27 日 14 時至 10 月 30 日 15 時，在井排匯管處投加硫化氫吸收劑 2 號，第三階段 10 月 30 日 15 時至 11 月 13 日 10 時，在商 56-116 計量站投加硫化氫吸收劑 2 號。加藥后對各硫化氫檢測點(diǎn)的硫化氫濃度進(jìn)行檢測。

2.1 第一階段實(shí)驗(yàn)數(shù)據(jù) 計量站來液 加熱爐 兩相分離器 一 次 沉 降油罐 外輸油罐 外輸油罐

外輸油泵 一 次 除 油水罐 回注水 卸油 硫化氫檢測點(diǎn) C 硫化氫檢測 A硫化氫檢測點(diǎn) B 硫化氫檢測點(diǎn) E 硫化氫檢測點(diǎn) D 原油 污水 加藥點(diǎn)

2.1.1 設(shè)備安裝

現(xiàn)場加藥設(shè)備

兩相分離器 現(xiàn)場管道內(nèi)壓力為 0.4Mpa，使用加藥泵出口壓力為 5Mpa，能夠?qū)⑺巹┘尤雲(yún)R管內(nèi)。加藥口選擇在計量站來液的管匯處，藥劑通過加熱爐，藥劑可以和聯(lián)合站所有來液充分接觸，去除油水采出液中的高濃度硫化氫。

2.1.2 檢測方法 實(shí)驗(yàn)設(shè)備：便攜式 H 2 S 氣體測試儀（相關(guān)配件）、燒杯、塑料桶（5L）。

硫化氫吸收劑劑：LQ 硫化氫吸收劑、測試管、注射器。

結(jié)合氣體測試的特性，本實(shí)驗(yàn)采用密閉式氣體收集裝置，保證采集樣本中氣體不擴(kuò)散，不流失，科學(xué)地得出實(shí)驗(yàn)數(shù)據(jù)，具體實(shí)驗(yàn)方法如下：

使用 5L 塑料桶，在內(nèi)蓋先開好僅適合檢測管進(jìn)入的小孔，墊入一次薄膜后蓋好內(nèi)蓋，形成密封效果，測試時用檢測管插入刺穿薄膜。此設(shè)計相對氣密性較好，采樣后至檢測前無氣體泄漏。并導(dǎo)入 1L 液體后標(biāo)注液位，確保每次取樣均為 1L 液體。

使用硫化氫檢測設(shè)備為：便攜式硫化氫測試儀、測試管、注射器。

儀器測試方法：將硫化氫檢測軟管插入取樣桶中檢測液面以上,檢測氣相中硫化氫含量，通過硫化氫測試儀測出硫化氫濃度。

測試管測試方法：用注射器抽取 50ml 液面上層氣體，然后緩慢注入測試管中，50ml 氣體注射完畢后，進(jìn)行讀數(shù)。

2.1.3 第一階段實(shí)驗(yàn)數(shù)據(jù) 10 月 22 日至 10 月 27 日在聯(lián)合站加藥點(diǎn)投加硫化氫吸收劑 1 號，每天加藥

量為 72kg/d（加藥泵最大排量），然后檢測各硫化氫檢測點(diǎn)的硫化氫濃度，數(shù)據(jù)如下：

兩相分離器檢測點(diǎn) A 硫化氫濃度數(shù)據(jù)

檢測點(diǎn)名稱 加藥量 kg 硫化氫濃度 ppm 日期 兩相分離器 0 1600（儀器）1800（測試管）

兩相分離器 6.6 1600（儀器）1800（測試管）

兩相分離器 78 553（儀器）600（測試管）

兩相分離器 81.1 450（儀器）

兩相分離器 85.5 416（儀器）

兩相分離器 90 193（儀器）

兩相分離器 95 600（儀器）

兩相分離器 98 549（儀器）

兩相分離器 142 620（儀器）1000（測試管）

兩相分離器 148 450（儀器）500（測試管）

兩相分離器 164 400（儀器）450（測試管）

兩相分離器 170 500（儀器）800（測試管）

兩相分離器 175 470（儀器）750（測試管)

兩相分離器 218 700（儀器）750（測試管)

兩相分離器 240 650（儀器）700（測試管)

兩相分離器 258 500（儀器）650（測試管)

兩相分離器 287 600（儀器）900（測試管)

兩相分離器 303 350（儀器）400（測試管)

兩相分離器 315 560(儀器)

兩相分離器 320 312（儀器）

兩相分離器 368 760（儀器）

兩相分離器 370 560（儀器）

兩相分離器 385 600（儀器）

通過以上數(shù)據(jù)，加入硫化氫吸收劑 24h 后，硫化氫濃度從 1600ppm 降至553ppm。之后的繼續(xù)沖擊加藥過程中，檢測點(diǎn)的硫化濃度在 300ppm-700ppm 之間波動。由于加藥點(diǎn)至檢測點(diǎn)，管輸路程短，藥劑的的反應(yīng)時間在 5 分鐘左右，藥劑反應(yīng)時間短，硫化氫濃度下降不明顯。

一次沉降油罐檢測點(diǎn) C

檢測點(diǎn)名稱 硫化氫濃度 ppm 日期 一次沉降油罐 2000（儀器）

一次沉降油罐 2000（儀器）

一次沉降油罐 2000（儀器）

一次沉降油罐 2000（儀器）

一次沉降油罐 2000（儀器）

一次沉降油罐 2000（儀器）

在一次沉降油罐的罐頂部進(jìn)行了硫化氫濃度檢測，每次檢測數(shù)據(jù)均為2000ppm（2000ppm 為儀器的最大量程，實(shí)際硫化氫濃度要超出 2000ppm），加藥后，罐頂?shù)牧蚧瘹錆舛纫恢略?2000ppm，由于一次沉降罐為 2000m 3，罐體內(nèi)部積存大量的硫化氫和 SRB 細(xì)菌，完全去除罐體內(nèi)部的硫化氫，需要較長時間的藥劑沖擊。

檢測點(diǎn) E 一次除油水罐數(shù)據(jù)

檢測點(diǎn)名稱 硫化氫濃度 ppm 日期 一次除油水罐 2000（儀器）2000（測試管）

一次除油水罐 1900（儀器）2000（測試管）

一次除油水罐 2000（儀器）

一次除油水罐 2000（儀器）

一次除油水罐 1700（儀器）

一次除油水罐 1700（儀器）

通過一次除油水罐污水進(jìn)罐處檢測點(diǎn)數(shù)據(jù)，加入藥劑后第四天，污水硫化氫濃度從 2000ppm 降至 1700ppm，污水中硫化氫濃度有下降趨勢。

綜上數(shù)據(jù)，藥劑自 10 月 22 日加入后，兩相分離器處檢測點(diǎn)和一次除油水罐污水進(jìn)罐處檢測點(diǎn)硫化氫濃度均有下降趨勢。兩相分離器檢測點(diǎn) A 硫化氫濃度降至300ppm-700ppm之間。一次除油水罐進(jìn)罐處檢測點(diǎn)E硫化氫濃度由2000ppm降至 1700ppm。一次沉降油罐硫化氫濃度檢測均為 2000ppm，由于一次沉降油罐，罐體內(nèi)硫化氫濃度高，硫化氫濃度超出儀器檢測范圍，沒有檢測到硫化氫下降趨勢。

劑 加入硫化氫吸收劑 1 號后，各檢測點(diǎn)硫化氫濃度下降趨勢不明顯。加藥量

為 為 72kg/d，已是加藥泵的最大量程，為了提高藥劑每天沖擊加藥量，盡快降低各檢測點(diǎn)的硫化氫濃度，劑 準(zhǔn)備用高濃度的硫化氫吸收劑 2 號進(jìn)行沖擊加藥。

2.2 第二階段實(shí)驗(yàn)數(shù)據(jù) 2.2.1 各檢測點(diǎn)實(shí)驗(yàn)數(shù)據(jù) 自 10 月 27 日更換硫化氫吸收劑 2 號進(jìn)行沖擊加藥，每天加藥量為 72kg/d（加藥泵最大排量），然后檢測各硫化氫檢測點(diǎn)的硫化氫濃度，數(shù)據(jù)如下：

兩相分離器檢測點(diǎn) A 硫化氫濃度數(shù)據(jù)

檢測點(diǎn)名稱 加藥量 kg 硫化氫濃度 ppm 日期 兩相分離器 3 560（儀器）

兩相分離器 51 400（儀器）

兩相分離器 60 265（儀器）

兩相分離器 69 182（儀器）

兩相分離器 75 80（儀器）

兩相分離器 123 140（儀器）

兩相分離器 132 100（儀器）

兩相分離器 141 80（儀器）

兩相分離器 147 120（儀器）

兩相分離器 217 234（儀器）

通過實(shí)驗(yàn)數(shù)據(jù)看出，加入高濃度藥劑后，硫化氫濃度降至 80ppm-200ppm 之間，采出液油、水、氣三相中的硫化氫濃度仍然沒有完全消除，檢測點(diǎn)的硫化氫濃度沒有降至 0ppm。

一次沉降油罐檢測點(diǎn) C

檢測點(diǎn)名稱 硫化氫濃度 ppm 日期 一次沉降油罐 2000（儀器）

一次沉降油罐 2000（儀器）

一次沉降油罐 2000（儀器）

在一次沉降油罐的罐頂部進(jìn)行了硫化氫濃度檢測，每次檢測數(shù)據(jù)均為

2000ppm（2000ppm 為儀器的最大量程，實(shí)際硫化氫濃度要超出 2000ppm），藥劑沖擊后，一次沉降油罐硫化氫濃度沒有下降趨勢。

外輸油檢測點(diǎn) D

檢測點(diǎn)名稱 硫化氫濃度 ppm 日期 外輸油 200（儀器）

外輸油 150（儀器）

外輸油 145（儀器）

外輸油 145（儀器）

使用硫化氫吸收劑 2 號沖擊加藥后，外輸油硫化氫濃度不斷下降，停止加藥后，硫化氫濃度降至 145ppm 左右。

檢測點(diǎn) E 一次除油水罐數(shù)據(jù)

檢測點(diǎn)名稱 硫化氫濃度 ppm 日期 一次除油水罐 1700（儀器）

一次除油水罐 1300（儀器）

一次除油水罐 1450（儀器）

一次除油水罐 1450（儀器）

通過試驗(yàn)數(shù)據(jù)，高濃度藥劑沖擊加藥后，檢測點(diǎn) E 的硫化氫濃度從 1700ppm降至 1450ppm，硫化氫濃度有下降趨勢。

綜上數(shù)據(jù)，高濃度硫化氫吸收劑 2 號沖擊加藥后，兩相分離器檢測點(diǎn) A 的硫化氫濃度由 1600ppm 降至 80ppm-200ppm 之間，一次沉降油罐的罐頂部檢測點(diǎn) C 的硫化氫濃度檢測數(shù)據(jù)均為 2000ppm，藥劑沖擊后，一次沉降油罐硫化氫濃度沒有下降趨勢。外輸油檢測點(diǎn) D 的硫化氫濃度由 200ppm 降至 145ppm，一次除油水罐進(jìn)罐處檢測點(diǎn) E 的硫化氫濃度由 2000ppm 降至 1450ppm。

2.2.2 藥劑室內(nèi)實(shí)驗(yàn) 在兩相分離器檢測點(diǎn) A 處取得油水采出液仍有硫化氫殘存，硫化氫沒有完

全被藥劑吸收，我們?nèi)〔沙鲆核畼雍螅?40 攝氏度的條件下，進(jìn)行了放置檢測，在取樣桶內(nèi)延長藥劑的反應(yīng)時間，測試加入的高濃度藥劑能否完全吸收采出液中高濃度的硫化氫，實(shí)驗(yàn)數(shù)據(jù)如下：

樣品名稱 H 2 S 濃度 0min 30min 40min 60min 商 56-116 計量站來液 1600 1600 1600 1596 檢測點(diǎn) A 采出液水樣 1 212 4 0 0 檢測點(diǎn) A 采出液水樣 2 84 0 0 0 檢測點(diǎn) A 采出液水樣 3 147 0 0 0

通過以上實(shí)驗(yàn)數(shù)據(jù)，商 56-116 計量站來液沒有加入硫化氫吸收劑，放置 1h后，取樣桶內(nèi)的硫化氫濃度基本沒有變化，兩相分離器處檢測點(diǎn) A 的采出液水樣，放置 40min 后，均降為 0ppm，因此加入的硫化氫吸收劑能夠完全吸收采出液中的高濃度硫化氫，檢測點(diǎn) A 硫化氫濃度未降至 0ppm，原因是到達(dá)檢測點(diǎn) A藥劑和采出液中硫化氫反應(yīng)時間短，無法完全去采出液中的硫化氫。

點(diǎn) 由于現(xiàn)在的加藥點(diǎn)距離檢測點(diǎn) A 管程太短，藥劑反應(yīng)時間不夠，檢測點(diǎn) E點(diǎn) 管程距離太長，藥劑消耗嚴(yán)重，檢測點(diǎn) C 油罐內(nèi)含有大量的硫化氫和 SRB 細(xì)菌無法短時間內(nèi)消除。取 兩相分離點(diǎn) 器處檢測點(diǎn) A 的采出液水樣，放置 40min 后，降為 為 0ppm。現(xiàn)場工藝流程中無法選取藥劑和采出液反應(yīng)時間為 1h 左右的檢測點(diǎn)。

加藥點(diǎn)的選取不夠科學(xué)，無法短時間內(nèi)確定使硫化氫濃度降至安全范圍內(nèi)需要的 最佳商 加藥量。因此，更換加藥點(diǎn)，在商 56-116 計量站進(jìn)行藥劑的投加。

2.3 第三階段實(shí)驗(yàn)數(shù)據(jù) 自 10 月 30 號開始在商 56-116 計量站處進(jìn)行加藥，在聯(lián)合站內(nèi) S56-116 來液井排處進(jìn)行檢測，加藥點(diǎn)和檢測點(diǎn)兩點(diǎn)之間初步估算 2500m 左右，計量站每天的液量為 200m 3 /d 左右，理論上加藥點(diǎn)到檢測點(diǎn)藥劑的反應(yīng)時間為 6h 左右，10 月 30 日 16：00 進(jìn)行加藥，藥劑為硫化氫吸收劑 1 號，加藥量為 160kg/d，數(shù)據(jù)如下：

計量站加藥點(diǎn)

井排來液處檢測點(diǎn)

檢測點(diǎn) B 商 56-116 來液井排處

檢測點(diǎn)名稱 加藥量 kg 硫化氫濃度 ppm 日期 站內(nèi)商 56-116 井排 0 1600（儀器）1800（測試管）

站內(nèi)商 56-116 井排 116 1600（儀器）1700（測試管）

站內(nèi)商 56-116 井排 150 1600（儀器）

站內(nèi)商 56-116 井排 162 1560（儀器）

站內(nèi)商 56-116 井排 268 1600（儀器）

站內(nèi)商 56-116 井排 310 1540（儀器）

站內(nèi)商 56-116 井排 326 1660（儀器）

站內(nèi)商 56-116 井排 430 286（儀器）175（測試管）

站內(nèi)商 56-116 井排 448 246（儀器）

站內(nèi)商 56-116 井排 474 170（儀器）

站內(nèi)商 56-116 井排 490 220（儀器）120(測試管）

站內(nèi)商 56-116 井排 598 165（儀器）

站內(nèi)商 56-116 井排 1124 65（儀器）

站內(nèi)商 56-116 井排 1446 95（儀器）

站內(nèi)商 56-116 井排 1726 30（儀器）

站內(nèi)商 56-116 井排 1760 40（儀器）

站內(nèi)商 56-116 井排 1874 20（儀器）

站內(nèi)商 56-116 井排 1920 10（儀器）

站內(nèi)商 56-116 井排 2036 0（儀器）

站內(nèi)商 56-116 井排 2070 0（儀器）0（測試管）

站內(nèi)商 56-116 井排 2204 0（儀器）0（測試管）

站內(nèi)商 56-116 井排 2240 0（儀器）0（測試管）

站內(nèi)商 56-116 井排 2400 0（儀器）0（測試管）

通過以上數(shù)據(jù)，在商 56-116 計量站沖擊加藥后，井排處檢測點(diǎn)的硫化氫濃度不斷下降。沖擊第 13 天至 11 月 12 日檢測點(diǎn)硫化氫濃度降為 0ppm，由于管道內(nèi)有殘余的硫化物、SRB 細(xì)菌和其他活性物，它們會和硫化氫吸收劑反應(yīng)，因此需要沖擊一段時間，井排處檢測點(diǎn)的硫化氫濃度降至 0ppm。之后檢測點(diǎn)的硫化氫濃度仍為 0ppm，數(shù)據(jù)穩(wěn)定。

.兩相分離器檢測點(diǎn) A 硫化氫濃度數(shù)據(jù)

檢測點(diǎn)名稱

硫化氫濃度 ppm 兩相分離器 234（儀器）

兩相分離器 1700（儀器）

兩相分離器 1600（儀器）

兩相分離器 1700（儀器）

兩相分離器 1600（儀器）

兩相分離器 1750（儀器）

兩相分離器 224（儀器）

兩相分離器 186（儀器）

兩相分離器 120（儀器）

兩相分離器 86（儀器）

通過實(shí)驗(yàn)數(shù)據(jù)，在商 56-116 計量站加藥后，檢測點(diǎn) A 的硫化氫濃度先上升后下降的趨勢。由于聯(lián)合站內(nèi)停止加藥后，藥劑在計量站投加，短時間內(nèi)藥劑無法沖擊到兩相分離器檢測點(diǎn)，造成檢測點(diǎn)硫化氫濃度升高。藥劑沖擊到一定時間后，兩相分離器檢測點(diǎn) A 的硫化氫濃度不斷下降，停止加藥時檢測點(diǎn) A 的硫化氫濃度降至 100ppm 左右。

一次沉降油罐檢測點(diǎn) C

檢測點(diǎn)名稱 硫化氫濃度 ppm 一次沉降油罐 2000（儀器）

一次沉降油罐 2000（儀器）

在一次沉降油罐的罐頂部進(jìn)行了硫化氫濃度檢測，每次檢測數(shù)據(jù)均為2000ppm（2000ppm 為儀器的最大量程，實(shí)際硫化氫濃度要超出 2000ppm），藥劑沖擊后，一次沉降油罐硫化氫濃度沒有下降趨勢。

外輸油檢測點(diǎn) D

檢測點(diǎn)名稱 硫化氫濃度 ppm 日期 外輸油 125（儀器）

外輸油 98（儀器）

外輸油 62（儀器）

外輸油 0（儀器）0(測試管）

外輸油 0（儀器）

外輸油 0（儀器）0(測試管）

外輸油 0（儀器）

在計量站沖擊加藥后，外輸油檢測點(diǎn)硫化氫濃度不斷下降，11 月 8 號，外輸油檢測點(diǎn)硫化氫濃度降至 0ppm，之后數(shù)據(jù)穩(wěn)定。

檢測點(diǎn) E 一次除油水罐數(shù)據(jù)

檢測點(diǎn)名稱 硫化氫濃度 ppm 日期 一次除油水罐 1700（儀器）

一次除油水罐 1300（儀器）

一次除油水罐 1450（儀器）

一次除油水罐 1450（儀器）

一次除油水罐 1500（儀器）

一次除油水罐 324（儀器）

一次除油水罐 220（儀器）

一次除油水罐 220（儀器）

一次除油水罐 220（儀器）

在計量站沖擊加藥后，最后一次水罐的硫化氫濃度降至 220ppm 左右，可能因?yàn)榱鞒梯^長，水相中硫化氫含量較高，無法完全去除采出液中高濃度硫化氫。

綜上數(shù)據(jù)，在計量站沖擊加藥后，沖擊第 13 天，商 56-116 計量站井排來液處檢測點(diǎn) B 硫化氫濃度降至 0ppm，之后檢測，檢測點(diǎn)硫化氫濃度穩(wěn)定在 0ppm，沒有上升趨勢。兩相分離器處檢測點(diǎn) A 硫化氫濃度降至 100ppm 左右，一次沉降油罐罐頂檢測點(diǎn) C 硫化氫濃度仍為 2000ppm，短時間內(nèi)，無法將硫化氫濃度降至安全范圍內(nèi)。外輸油檢測點(diǎn) D 硫化氫濃度，11 月 8 日降至 0ppm，之后檢測穩(wěn)定，一次除油水罐進(jìn)罐處檢測點(diǎn) E，硫化氫濃度降至 220ppm 左右。商 56-116 計量站井排處來液處檢測點(diǎn) B 硫化氫濃度降至 0ppm，證明硫化氫吸收劑 1 號能夠完全處理采出液油、水、氣三相中的高濃度硫化氫。后期硫化氫治理，藥劑可以在各計量站投加，藥劑和采出液中硫化氫有充足的反應(yīng)時間，最后完全消除采出液中高濃度硫化氫，使聯(lián)合站內(nèi)各檢測點(diǎn)的硫化氫濃度降至安全范圍，保證站內(nèi)人員安全。

三、總結(jié) 1、對商 56 站各節(jié)點(diǎn)硫化氫進(jìn)行了測試，兩相分離器出口檢測點(diǎn) A 的硫化氫濃度為 1700ppm，商 56-116 井排油水來液檢測點(diǎn) B 的硫化氫濃度 1600ppm，一次沉降油罐檢測點(diǎn) C 的硫化氫濃度為 2000ppm，外輸油檢測點(diǎn) D 的硫化氫濃度 200ppm，一次除油水罐污水進(jìn)罐處檢測點(diǎn) E 的硫化氫濃度 2000ppm（檢測方法：用 5L 取樣桶，取 1L 采出液，放置 30 秒，檢測液面以上氣體中的硫化氫濃度）。各節(jié)點(diǎn)的硫化氫濃度均大于安全臨界濃度 20ppm，硫化氫含量超標(biāo)嚴(yán)重。

2、實(shí)驗(yàn)第一階段，藥劑自 10 月 22 日加入后，兩相分離器處檢測點(diǎn)和一次除油水罐污水進(jìn)罐處檢測點(diǎn)硫化氫濃度均有下降趨勢。兩相分離器檢測點(diǎn) A 硫化氫濃度降至 300ppm-700ppm 之間。一次除油水罐進(jìn)罐處檢測點(diǎn) E 硫化氫濃度由 2000ppm 降至 1700ppm。一次沉降油罐硫化氫濃度檢測均為 2000ppm，由于一次沉降油罐，罐體內(nèi)硫化氫濃度高，硫化氫濃度超出儀器檢測范圍，沒有檢測到硫化氫下降趨勢。加入硫化氫吸收劑 1 號后，各檢測點(diǎn)硫化氫濃度下降趨勢不明顯。加藥量為 72kg/d，已是加藥泵的最大量程，為了提高藥劑每天沖擊加藥量，盡快降低各檢測點(diǎn)的硫化氫濃度，準(zhǔn)備用高濃度的硫化氫吸收劑 2 號進(jìn)行沖擊加藥。

3、實(shí)驗(yàn)第二階段，高濃度硫化氫吸收劑 2 號沖擊加藥后，兩相分離器檢測點(diǎn) A 的硫化氫濃度由 1600ppm 降至 80ppm-200ppm 之間，一次沉降油罐的罐頂

部檢測點(diǎn) C 的硫化氫濃度檢測數(shù)據(jù)均為 2000ppm，藥劑沖擊后，一次沉降油罐硫化氫濃度沒有下降趨勢。外輸油檢測點(diǎn) D 的硫化氫濃度由 200ppm 降至145ppm，一次除油水罐進(jìn)罐處檢測點(diǎn) E 的硫化氫濃度由 2000ppm 降至 1450ppm由于現(xiàn)在的加藥點(diǎn)距離檢測點(diǎn) A 管程太短，藥劑反應(yīng)時間不夠，檢測點(diǎn) E 管程距離太長，藥劑消耗嚴(yán)重，檢測點(diǎn) C 油罐內(nèi)含有大量的硫化氫和 SRB 細(xì)菌無法短時間內(nèi)消除。取兩相分離器處檢測點(diǎn) A 的采出液水樣，放置 40min 后，將為0ppm。現(xiàn)場工藝流程中無法選取藥劑和采出液反應(yīng)時間為 1h 左右的檢測點(diǎn)。加藥點(diǎn)的選取不夠科學(xué)，無法短時間內(nèi)確定使硫化氫濃度降至安全范圍內(nèi)需要的加藥量。因此，更換加藥點(diǎn)，在商 56-116 計量站進(jìn)行藥劑的投加。

4、高濃度硫化氫吸收劑 2 號沖擊加藥后，取兩相分離器檢測點(diǎn) A 處的采出液樣品進(jìn)行了室內(nèi)實(shí)驗(yàn)，兩相分離器處檢測點(diǎn) A 的采出液水樣，放置 40min 后，均將為 0ppm，因此加入的硫化氫吸收劑能夠完全吸收采出液中的高濃度硫化氫，檢測點(diǎn) A 硫化氫濃度未降至 0ppm，原因是到達(dá)檢測點(diǎn) A 藥劑和采出液中硫化氫反應(yīng)時間短，無法完全去采出液中的硫化氫。由于現(xiàn)在的加藥點(diǎn)距離檢測點(diǎn) A管程太短，藥劑反應(yīng)時間不夠，檢測點(diǎn) E 管程距離太長，藥劑消耗嚴(yán)重，檢測點(diǎn) C 油罐內(nèi)含有大量的硫化氫和 SRB 細(xì)菌無法短時間內(nèi)消除，加藥點(diǎn)的選取不夠科學(xué)，無法短時間內(nèi)確定使硫化氫濃度降至安全范圍內(nèi)需要的加藥量。因此，更換加藥點(diǎn)，在商 56-116 計量站進(jìn)行藥劑的投加。

5、實(shí)驗(yàn)第三階段，改變加藥點(diǎn)，在商 56-116 計量站內(nèi)進(jìn)行藥劑投加。在計量站沖擊加藥后，沖擊第 13 天，商 56-116 計量站井排來液處檢測點(diǎn) B 硫化氫濃度降至 0ppm，之后之后檢測，檢測點(diǎn)硫化氫濃度穩(wěn)定在 0ppm，沒有上升趨勢。兩相分離器處檢測點(diǎn) A 硫化氫濃度降至 100ppm 左右，一次沉降油罐罐頂檢測點(diǎn)C 硫化氫濃度仍為 2000ppm，短時間內(nèi)，無法將硫化氫濃度降至安全范圍內(nèi)。外輸油檢測點(diǎn) D 硫化氫濃度，11 月 8 日降至 0ppm，之后檢測穩(wěn)定，一次除油水罐進(jìn)罐處檢測點(diǎn) E，硫化氫濃度降至 220ppm 左右。商 56-116 計量站井排處來液處檢測點(diǎn) B 硫化氫濃度降至 0ppm，證明硫化氫吸收劑 1 號能夠完全處理采出液油、水、氣三相中的高濃度硫化氫。后期硫化氫治理，藥劑可以在各計量站投加，藥劑和采出液中硫化氫有充足的反應(yīng)時間，最后完全消除采出液中高濃度硫化氫，使聯(lián)合站內(nèi)各檢測點(diǎn)的硫化氫濃度降至安全范圍，保證站內(nèi)人員安全。

6、實(shí)驗(yàn)分三個階段，投加硫化氫吸收劑 1 號 2785kg，投加硫化氫吸收劑 2號 217kg，藥劑投加后能夠使兩相分離器處檢測點(diǎn) A 的硫化氫濃度降至 100ppm，商 56-116 計量站井排來液處檢測點(diǎn) B 的硫化氫濃度降至 0ppm，一次沉降油罐罐頂檢測點(diǎn) C 的硫化氫濃度 2000ppm，外輸油檢測點(diǎn) D 的硫化氫濃度降至 0ppm，一次除油水罐進(jìn)罐處檢測點(diǎn) E 的硫化氫濃度降至 220ppm。檢測點(diǎn) B 和檢測點(diǎn) D的硫化氫濃度降至 0ppm，證明三嗪硫化氫吸收劑能夠處理現(xiàn)場采出液中高濃度硫化氫，使油、水、氣三相中的硫化氫濃度降至安全范圍內(nèi)。

第二篇：數(shù)字信號處理實(shí)驗(yàn)報告

南京郵電大學(xué)

實(shí) 驗(yàn) 報 告

實(shí)驗(yàn)名稱_____熟悉MATLAB環(huán)境 ___ 快速傅里葉變換及其應(yīng)用 ____IIR數(shù)字濾波器的設(shè)計_ FIR數(shù)字濾波器的設(shè)計

課程名稱 數(shù)字信號處理A

班級學(xué)號_______09002111___________ 姓 名 王都超

開課時間 2011/2012學(xué)年，第 二 學(xué)期

實(shí)驗(yàn)一

熟悉MATLAB環(huán)境

一、實(shí)驗(yàn)?zāi)康?/p>

（1）熟悉MATLAB的主要操作命令。（2）學(xué)會簡單的矩陣輸入和數(shù)據(jù)讀寫。（3）掌握簡單的繪圖命令。

（4）用MATLAB編程并學(xué)會創(chuàng)建函數(shù)。（5）觀察離散系統(tǒng)的頻率響應(yīng)。

二、實(shí)驗(yàn)內(nèi)容

(1)數(shù)組的加、減、乘、除和乘方運(yùn)算。輸入A=[1 2 3 4]，B=[3,4,5,6]，求

C=A+B，D=A-B，E=A.*B,F=A./B,G=A.^B。并用stem語句畫出A、B、C、D、E、F、G。

D =

-2 E =

F =

0.3333

0.5000

0.6000

0.6667 G =

243

4096(2)用MATLAB實(shí)現(xiàn)下列序列： a)x(n)?0.8n 0?n?1

5n=0:1:15;x1=0.8.^n;a=(0.2+3*i)*n;stem(x1)b)x(n)?e(0.2?3j)n 0?n?15

n=0:1:15;x2=exp(a);a=(0.2+3*i)*n;stem(x2)

c)x(n)?3cos(0.125?n?0.2?)?2sin(0.25?n?0.1?)

0?n?15

(4)繪出下列時間函數(shù)的圖形，對x軸、y軸以及圖形上方均須加上適當(dāng)?shù)臉?biāo)注： a)x(t)?sin(2?t)0?t?10s

b)x(t)?cos(100?t)sin(?t)0?t?4s t=0:0.01:4;x=cos(100*pi*t).*sin(pi*t);plot(t,x, 'r-');xlabel('t'),ylabel('x(t)'),title('cos')

(6)給定一因果系統(tǒng)H(z)?(1?頻響應(yīng)和相頻響應(yīng)。

2z?1?z?2)/(1?0.67z?1?0.9z?2)，求出并繪制H(z)的幅

(7)計算序列{8-2-1 2 3}和序列{2 3-1-3}的離散卷積，并作圖表示卷積結(jié)果。

(8)求以下差分方程所描述系統(tǒng)的單位脈沖響應(yīng)h(n), 0?n?50

y(n)?0.1y(n?1)?0.06y(n?2)?x(n)?2x(n?1)

實(shí)驗(yàn)過程與結(jié)果（含實(shí)驗(yàn)程序、運(yùn)行的數(shù)據(jù)結(jié)果和圖形）； clear all;N=50;a=[1-2];b=[1 0.1-0.06];x1=[1 zeros(1,N-1)];n=0:1:N-1;h=filter(a,b,x1);stem(n,h)axis([-1 53-2.5 1.2])

實(shí)驗(yàn)二

快速傅里葉變換及其應(yīng)用

一、實(shí)驗(yàn)?zāi)康?/p>

(1)在理論學(xué)習(xí)的基礎(chǔ)上，通過本實(shí)驗(yàn)，加深對FFT的理解，熟悉MATLAB中的有關(guān)函數(shù)。(2)應(yīng)用FFT對典型信號進(jìn)行頻譜分析。

(3)了解應(yīng)用FFT進(jìn)行信號頻譜分析過程中可能出現(xiàn)的問題，以便在實(shí)際中正確應(yīng)用FFT。(4)應(yīng)用FFT實(shí)現(xiàn)序列的線性卷積和相關(guān)。

二、實(shí)驗(yàn)內(nèi)容

實(shí)驗(yàn)中用到的信號序列 a)高斯序列

??(n?p)q?xa(n)??e?0?20?n?15 其他

b)衰減正弦序列

?e?ansin(2?fn)xb(n)??0?0?n?15其他

c)三角波序列 ?n?xc(n)??8?n?0?0?n?34?n?7 其他

d)反三角波序列

?4?n?xd(n)??n?4?0?0?n?34?n?7 其他

(1)觀察高斯序列的時域和幅頻特性，固定信號xa(n)中參數(shù)p=8，改變q的值，使q分別等于2，4，8，觀察它們的時域和幅頻特性，了解當(dāng)q取不同值時，對信號序列的時域幅頻特性的影響；固定q=8，改變p，使p分別等于8，13，14，觀察參數(shù)p變化對信號序列的時域及幅頻特性的影響，觀察p等于多少時，會發(fā)生明顯的泄漏現(xiàn)象，混疊是否也隨之出現(xiàn)？記錄實(shí)驗(yàn)中觀察到的現(xiàn)象，繪出相應(yīng)的時域序列和幅頻特性曲線。

(3)觀察三角波和反三角波序列的時域和幅頻特性，用N=8點(diǎn)FFT分析信號序列xc(n)和觀察兩者的序列形狀和頻譜曲線有什么異同？繪出兩序列及其幅頻特性xd(n)的幅頻特性，曲線。

在xc(n)和xd(n)末尾補(bǔ)零，用N=32點(diǎn)FFT分析這兩個信號的幅頻特性，觀察幅頻特性發(fā)生了什么變化？兩種情況的FFT頻譜還有相同之處嗎？這些變化說明了什么？

(5)用FFT分別實(shí)現(xiàn)xa(n)（p＝8，q＝2）和xb(n)（a＝0.1，f＝0.0625）的16點(diǎn)循環(huán)卷積和線性卷積。

n=0:15;p=8;q=2;

xa=exp(-(n-p).^2/q);subplot(2,3,1);stem(n,xa,'.');title('xa波形');

Xa=fft(xa,16);subplot(2,3,4);stem(abs(Xa),'.');

title('Xa(k)=FFT[xa(n)]的波形 ');A=1;f=0.0625;a=0.1;

xb=exp(-a*n).*sin(2*pi*f*n);subplot(2,3,2);stem(n,xb,'.');title('xb波形');Xb=fft(xb,16);subplot(2,3,5);stem(abs(Xb),'.');

title('Xb(k)=FFT[xb(n)]的波形 ');

實(shí)驗(yàn)過程與結(jié)果（含實(shí)驗(yàn)程序、運(yùn)行的數(shù)據(jù)結(jié)果和圖形）；

實(shí)驗(yàn)三 IIR數(shù)字濾波器的設(shè)計

一、實(shí)驗(yàn)?zāi)康?/p>

(1)掌握雙線性變換法及脈沖響應(yīng)不變法設(shè)計IIR數(shù)字濾波器的具體設(shè)計方法及其原理，熟悉用雙線性變換法及脈沖響應(yīng)不變法設(shè)計低通、高通和帶通IIR數(shù)字濾波器的計算機(jī)編程。

(2)觀察雙線性變換及脈沖響應(yīng)不變法設(shè)計的濾波器的頻域特性，了解雙線性變換法及脈沖響應(yīng)不變法的特點(diǎn)。

(3)熟悉巴特沃思濾波器、切比雪夫?yàn)V波器和橢圓濾波器的頻率特性。

二、實(shí)驗(yàn)內(nèi)容（1）P162 例4.4 設(shè)采樣周期T=250?s(采樣頻率fs=4kHz)，分別用脈沖響應(yīng)不變法和雙線性變換法設(shè)計一個三階巴特沃思低通濾波器，其3dB邊界頻率為fc=1kHz。

脈沖響應(yīng)不變法： fc=1000;fs=4000;OmegaC=2*pi*fc;[B,A]=butter(3, OmegaC,'s');[num1,den1]=impinvar(B,A,fs);[h1,w]=freqz(num1,den1);f = w/pi*fs/2;plot(f,abs(h1));

雙線性變換法： fc=1000;fs=4000;

OmegaC=2*fs*tan(pi*fc/fs);[B,A]=butter(3, OmegaC,'s');[num2,den2]=bilinear(B,A,fs);[h2,w]=freqz(num2,den2);f = w/pi*fs/2;plot(f,abs(h2));

同一圖中畫兩條曲線： fc=1000;fs=4000;OmegaC=2*pi*fc;[B,A]=butter(3, OmegaC,'s');[num1,den1]=impinvar(B,A,fs);[h1,w]=freqz(num1,den1);f = w/pi*fs/2;

OmegaC=2*fs*tan(pi*fc/fs);[B,A]=butter(3, OmegaC,'s');[num2,den2]=bilinear(B,A,fs);[h2,w]=freqz(num2,den2);f = w/pi*fs/2;plot(f,abs(h1),'r-.');hold on;plot(f,abs(h2),'g-');

(選做)(2)fc=0.2kHz,?=1dB,fr=0.3kHz,At=25dB,T=1ms；分別用脈沖響應(yīng)不變法及雙線性變換法設(shè)計一巴特沃思數(shù)字低通濾波器，觀察所設(shè)計數(shù)字濾波器的幅頻特性曲線，記錄帶寬和衰減量，檢查是否滿足要求。比較這兩種方法的優(yōu)缺點(diǎn)。

實(shí)驗(yàn)過程與結(jié)果（含實(shí)驗(yàn)程序、運(yùn)行的數(shù)據(jù)結(jié)果和圖形）；

實(shí)驗(yàn)四

FIR數(shù)字濾波器的設(shè)計

一、實(shí)驗(yàn)?zāi)康?/p>

(1)掌握用窗函數(shù)法，頻率采樣法及優(yōu)化設(shè)計法設(shè)計FIR濾波器的原理及方法，熟悉相應(yīng)的計算機(jī)編程；

(2)熟悉線性相位FIR濾波器的幅頻特性和相頻特性；

(3)了解各種不同窗函數(shù)對濾波器性能的影響。

二、實(shí)驗(yàn)內(nèi)容

(1)生成一個長度為20的矩形窗，畫出其時域和幅頻特性曲線。n=0:1:19;N=20;win(1:20)=1;[H,w]=freqz(win,1);subplot(2,1,1);stem(n,win)subplot(2,1,2);plot(w,abs(H));

(2)用矩形窗設(shè)計一個21階的線性相位低通FIR數(shù)字濾波器，截止頻率Wc＝0.25π，求出濾波器系數(shù)，并繪出濾波器的幅頻特性。修改程序，分別得到階次為N＝41，61的濾波器，并顯示其各自的幅頻曲線。

a)在上面所得的幾幅圖中，在截止頻率兩邊可以觀察到幅頻響應(yīng)的擺動行為。請問波紋的數(shù)量與濾波器脈沖響應(yīng)的長度之間有什么關(guān)系？

b)最大波紋的高度與濾波器脈沖響應(yīng)的長度之間有什么關(guān)系？

實(shí)驗(yàn)過程與結(jié)果（含實(shí)驗(yàn)程序、運(yùn)行的數(shù)據(jù)結(jié)果和圖形）； 21階的線性相位低通FIR數(shù)字濾波器： Wc=0.25*pi;N=21;M=(N-1)/2;

%位移量

for n=0:(N-1)

if(n== fix(M))

%中間的點(diǎn)單獨(dú)算

hd(n+1)=Wc/pi;

else

hd(n+1)=sin(Wc*(n-M))/(pi*(n-M));end;end;win=boxcar(N);%%%不同窗函數(shù)

h=hd.*win';[H,w]=freqz(h,1);n=0:1:N-1;subplot(3,1,1);stem(n,h)subplot(3,1,2);plot(w,abs(H));subplot(3,1,3);plot(w,angle(H));

41階的線性相位低通FIR數(shù)字濾波器： Wc=0.25*pi;N=41;M=(N-1)/2;

%位移量

for n=0:(N-1)

if(n== fix(M))

%中間的點(diǎn)單獨(dú)算

hd(n+1)=Wc/pi;

else

hd(n+1)=sin(Wc*(n-M))/(pi*(n-M));end;end;win=boxcar(N);%%%不同窗函數(shù)

h=hd.*win';[H,w]=freqz(h,1);n=0:1:N-1;subplot(3,1,1);stem(n,h)subplot(3,1,2);plot(w,abs(H));subplot(3,1,3);plot(w,angle(H));

61階的線性相位低通FIR數(shù)字濾波器： Wc=0.25*pi;N=61;M=(N-1)/2;

%位移量

for n=0:(N-1)

if(n== fix(M))

%中間的點(diǎn)單獨(dú)算

hd(n+1)=Wc/pi;

else

hd(n+1)=sin(Wc*(n-M))/(pi*(n-M));end;end;win=boxcar(N);%%%不同窗函數(shù)

h=hd.*win';[H,w]=freqz(h,1);n=0:1:N-1;subplot(3,1,1);stem(n,h)subplot(3,1,2);plot(w,abs(H));subplot(3,1,3);plot(w,angle(H));

數(shù)字信號處理實(shí)驗(yàn)小結(jié)及心得體會：

通過這次實(shí)驗(yàn)，我對MATLAB語言有了一定的認(rèn)識，雖然還不能完全用MATLAB獨(dú)立編寫程序，但對這種語言環(huán)境有了新的了解。我知道了一般的加減乘除在MATLAB中不同的意義。知道輸入、輸出語句怎么形成。通過快速傅里葉變換及其應(yīng)用的實(shí)驗(yàn)，加深了我對FFT的理解，還有對各典型信號的頻譜分析，改變參數(shù)后時域和幅頻特性的變化。IIR數(shù)字濾波器的設(shè)計讓我知道了巴特沃思濾波器和切比雪夫?yàn)V波器的頻率特性，還有雙線性變換及脈沖響應(yīng)不變法設(shè)計的濾波器的頻率特性。做這個實(shí)驗(yàn)的時候程序有點(diǎn)困難，很多細(xì)節(jié)問題不能考慮清楚，導(dǎo)致圖形出不來。FIR數(shù)字濾波器的設(shè)計出來的是三種窗的圖形，通過三種窗的比較，我了解了他們各自的特點(diǎn)，幅頻和相頻特性。我在這次實(shí)驗(yàn)中的收獲很大，接觸了很多新的知識，但在實(shí)驗(yàn)寫程序時，我發(fā)現(xiàn)自己還有很多不足。很多程序?qū)懖煌耆＿@是自己今后要加強(qiáng)的地方。

第三篇：數(shù)字信號處理實(shí)驗(yàn)報告

JIANGSU

UNIVERSITY OF TECHNOLOGY

數(shù)字信號處理實(shí)驗(yàn)報告

學(xué)院名稱： 電氣信息工程學(xué)院

專 業(yè)：

班 級： 姓 名： 學(xué) 號： 指導(dǎo)老師： 張維璽（教授）

2013年12月20日

實(shí)驗(yàn)一 離散時間信號的產(chǎn)生

一、實(shí)驗(yàn)?zāi)康?/p>

數(shù)字信號處理系統(tǒng)中的信號都是以離散時間形態(tài)存在的，所以對離散時間信號的研究是數(shù)字信號的基本所在。而要研究離散時間信號，首先需要產(chǎn)生出各種離散時間信號。使用MATLAB軟件可以很方便地產(chǎn)生各種常見的離散時間信號，而且它還具有強(qiáng)大繪圖功能，便于用戶直觀地處理輸出結(jié)果。

通過本實(shí)驗(yàn)，學(xué)生將學(xué)習(xí)如何用MATLAB產(chǎn)生一些常見的離散時間信號，實(shí)現(xiàn)信號的卷積運(yùn)算，并通過MATLAB中的繪圖工具對產(chǎn)生的信號進(jìn)行觀察，加深對常用離散信號和信號卷積和運(yùn)算的理解。

二、實(shí)驗(yàn)原理

離散時間信號是指在離散時刻才有定義的信號，簡稱離散信號，或者序列。離散序列通常用x（n）來表示，自變量必須是整數(shù)。常見的離散信號如下：（1）單位沖激序列δ（n）

如果δ（n）在時間軸上延遲了k個單位，得到δ（n-k），即長度為N的單位沖激序列δ（n）可以通過下面的MATLAB命令獲得。

n=-（N-1）：N-1 x=[zeros(1,N-1)1 zeros(1,N-1)]； stem(n,x)延遲K個采樣點(diǎn)的長度為N的單位沖激序列δ（n-k）（k

n=0：N-1 y=[zeros(1,M)1 zeros(1,N-M-1)]； stem（n，y）

（2）單位階躍序列u（n）

如果u（n）在時間軸上延遲了k個單位，得到u（n-k），即長度為N的單位階躍序列u（n）可以通過下面的MATLAB命令獲得。

n=-（N-1）：N-1 x=[zeros(1,N-1)ones(1,N)]； stem(n,x)延遲的單位階躍序列可以使用類似于單位沖激序列的方法獲得。（3）矩形序列

矩形序列有一個重要的參數(shù)，就是序列的寬度N。矩形序列與u（n）之間的關(guān)系為矩形序列等= u（n）— u（n-N）。

因此，用MATLAB表示矩形序列可利用上面的單位階躍序列組合而成。（4）正弦序列x（n）

這里，正弦序列的參數(shù)都是實(shí)數(shù)。與連續(xù)的正弦信號不同，正弦序列的自變量n必須為整數(shù)。可以證明，只有當(dāng)2π/w為有理數(shù)時，正弦序列具有周期性。

長度為N的正弦序列x（n）可以通過下面的MATLAB命令獲得。n=0：N-1 x=A*cos（2*pi*f*n/Fs+phase）（5）單邊實(shí)指數(shù)序列x（n）

長度為N的實(shí)指數(shù)序列x（n）可以通過下面的MATLAB命令實(shí)現(xiàn)。n=0：N-1 x=a.^n stem(n,x)單邊指數(shù)序列n的取值范圍為n>=0。當(dāng)|a|>1時，單邊指數(shù)序列發(fā)散；當(dāng)|a|<1時，單邊指數(shù)序列收斂。當(dāng)a>0時，該序列均取正值；當(dāng)a<0時，序列在正負(fù)擺動。

（6）負(fù)指數(shù)序列x（n）

當(dāng)a=0時，得到虛指數(shù)序列x（n）。

與連續(xù)負(fù)指數(shù)信號一樣，我們將負(fù)指數(shù)序列實(shí)部和虛部的波形分開討論，得到如下結(jié)論：

1）當(dāng)a>0時，負(fù)指數(shù)序列x（n）的實(shí)部和虛部分別是按指數(shù)規(guī)律增長的正弦振蕩序列；

2)當(dāng)a<0時，負(fù)指數(shù)序列x（n）的實(shí)部和虛部分別是按指數(shù)規(guī)律衰減的正弦振蕩序列；

3)當(dāng)a=0時，負(fù)指數(shù)序列x（n）即為虛指數(shù)序列，其實(shí)部和虛部分別是等幅的正弦振蕩序列；

長度為N的實(shí)指數(shù)序列x（n）可以通過下面的MATLAB命令實(shí)現(xiàn)。n=0：N-1 x=exp（（a.+j*w）*n）stem(n,real(x))或

stem(n,imag(x))

三、實(shí)驗(yàn)內(nèi)容及分析

?1n?01、編制程序產(chǎn)生單位沖激序列??n???“?并繪出其圖及??n?”學(xué)號后兩位0n?0?形。程序:(1)N=4;

n=-(N-1):N-1;

x=[zeros(1,N-1)1 zeros(1,N-1)];stem(n,x);

title('單位沖激序列');

grid on;

(2)N=6;

M=1;%學(xué)號01 n=-(N-1):N-1;

y=[zeros(1,N-M+1)1 zeros(1,N-M-1)];stem(n,y);

title('單位沖激序列');grid on;

分析:在上圖的基礎(chǔ)上向右平移了1個單位。

?1n?02、編制程序產(chǎn)生單位階躍序列u?n???、u?n?“學(xué)號后兩位”?及

0n?0?u?n??u?n?“學(xué)號后兩位”?，并繪出其圖形。程序: 4

(1)N=5;

n=-(N-1):N-1;

x=[zeros(1,N-1)ones(1,N)];stem(n,x);

title('單位階躍序列');grid on;

(2)N=6;

M=1;%學(xué)號01 n=-(N-1):N-1;

x=[zeros(1,N-M+1)ones(1,N-M)];stem(n,x);

title('單位階躍序列');grid on;

分析:在上圖的基礎(chǔ)上平移了1個單位.(3)N=6;

M=1;%學(xué)號01 n=-(N-1):N-1;

x=[zeros(1,N-1)ones(1,N)];y=[zeros(1,N-M+1)ones(1,N-M)];z=x-y;stem(n,z);

title('單位階躍序列');grid on;

2??

3、編制程序產(chǎn)生正弦序列x?n??cos?2?n?、x?n??cos??n?及

?學(xué)號后兩位?x?n??sin?2n?并繪出其圖形。

程序:(1)N=5;

A=1;

w=2*pi;phi=0;n=0:0.05:N-1;x=A*cos(w*n+phi);stem(n,x);title('余弦信號');grid on;

分析:該序列具有周期性,且輸出為余弦信號.(2)N=5;

A=1;

w=2*pi/1;%學(xué)號01 phi=0;n=0:0.05:N-1;x=A*cos(w*n+phi);stem(n,x);title('余弦信號');grid on;

;

分析:該序列具有周期性,且輸出為余弦信號.(3)N=5;

A=1;

w=2*pi;phi=0;

n=0:0.05:N-1;x=A*sin(w*n+phi);stem(n,x);title('正弦信號');grid on;

分析:該序列具有周期性,且輸出為正弦信號.4、編制程序產(chǎn)生復(fù)正弦序列x?n??e(2?j學(xué)號后兩位）n，并繪出其圖形。N=3;

n=0:0.2:N-1;

w=1;%學(xué)號01 x=exp((2+j*w)*n);subplot(2,1,1)

stem(n,real(x)),title('實(shí)部');grid on;subplot(2,1,2)

stem(n,imag(x)),title('虛部');grid on;

5、編制程序產(chǎn)生指數(shù)序列x?n??an，并繪出其圖形。其中a=學(xué)號后兩位、a=1/“學(xué)號后兩位”。

(1)N=10;

n=0:N-1;

a=1;%學(xué)號01 x=a.^n;stem(n,x);title('指數(shù)序列');grid on;

(2)N=10;

n=0:N-1;

a=1;%學(xué)號01 x=a.^(-n);stem(n,x);title('指數(shù)序列');grid on;

實(shí)驗(yàn)三 離散時間信號的頻域分析

一、實(shí)驗(yàn)?zāi)康?/p>

信號的頻域分析是信號處理中一種有效的工具。在離散信號的頻域分析中，通常將信號表示成單位采樣序列的線性組合，而在頻域中，將信號表示成復(fù)變量或的線性組合。通過這樣的表示，可以將時域的離散序列映射到頻域以便于進(jìn)一步的處理。

在本實(shí)驗(yàn)中，將學(xué)習(xí)利用MATLAB計算離散時間信號的DTFT和DFT，并加深對其相互關(guān)系的理解。

二、實(shí)驗(yàn)原理

（１）DTFT和DFT的定義及其相互關(guān)系。

（２）使用到的MATLAB命令有基于DTFT離散時間信號分析函數(shù)以及求解序列的DFT函數(shù)。

三、實(shí)驗(yàn)內(nèi)容及分析

（１）編程計算并畫出下面DTFT的實(shí)部、虛部、幅度和相位譜。

X(e)?jw0.0518?0.1553e1?1.2828ex(n)?cos?jw?jw?0.1553e?j2w?1.0388e?j2w?0.0518e?j3w?0.3418e?j3w

（２）計算32點(diǎn)序列

5?n16，0≦n≦31的32點(diǎn)和64點(diǎn)DFT，分別繪出幅度譜圖形，并繪出該序列的DTFT圖形。

３－１

clear;

x=[0.0518,-0.1553,0.1553,0.0518];y=[1,1.2828,1.0388,0.3418];w=[0:500]*pi/500 H=freqz(x,y,w);

magX=abs(H);angX=angle(H);realX=real(H);imagX=imag(H);subplot(221);plot(w/pi,magX);grid;

xlabel('frequency in pi unit');ylabel('magnitude');title('幅度 part');axis([0 0.9 0 1.1]);

subplot(223);plot(w/pi,angX);grid;

xlabel('frequency in pi unit');ylabel('radians');title('相位 part');axis([0 1-3.2 3.2]);

subplot(222);plot(w/pi,realX);grid;

xlabel('frequency in pi unit');ylabel('real part');title('實(shí)部 part');axis([0 1-1 1]);

subplot(224);plot(w/pi,imagX);grid;

xlabel('frequency in pi unit');ylabel('imaginary');title('虛部 part');axis([0 1-1 1.1]);

３－２

N=32;n=0:N-1;

xn=cos(5*pi*n/16);k=0:1:N-1;Xk=fft(xn,N);subplot(2,1,1);stem(n,xn);subplot(2,1,2);stem(k,abs(Xk));title('32點(diǎn)');figure N=64;n=0:N-1;

xn=cos(5*pi*n/16);k=0:1:N-1;Xk=fft(xn,N);subplot(2,1,1);stem(n,xn);subplot(2,1,2);stem(k,abs(Xk));title('64點(diǎn)');

（１）

（２）

實(shí)驗(yàn)四 離散時間LTI系統(tǒng)的Z域分析

一、實(shí)驗(yàn)?zāi)康?/p>

本實(shí)驗(yàn)通過使用MATLAB函數(shù)對離散時間系統(tǒng)的一些特性進(jìn)行仿真分析，以加深對離散時間系統(tǒng)的零極點(diǎn)、穩(wěn)定性，頻率響應(yīng)等概念的理解。學(xué)會運(yùn)用MATLAB分析離散時間系統(tǒng)的系統(tǒng)函數(shù)的零極點(diǎn)；學(xué)會運(yùn)用MATLAB分析系統(tǒng)函數(shù)的零極點(diǎn)分布與其時域特性的關(guān)系；學(xué)會運(yùn)用MATLAB進(jìn)行離散時間系統(tǒng)的頻率特性分析。

二、實(shí)驗(yàn)原理

離散時間系統(tǒng)的系統(tǒng)函數(shù)定義為系統(tǒng)零狀態(tài)響應(yīng)的Z變化與激勵的Z變化之比。

在MATLAB中系統(tǒng)函數(shù)的零極點(diǎn)可通過函數(shù)roots得到，也可借助函數(shù)tf2zp得到，tf2zp的語句格式為

[Z,P,K]=tf2zp(B,A)其中，B與A分別表示H（z）的分子與分母多項式的系數(shù)向量。它的作用是將H（z）的有理分式表示式轉(zhuǎn)換為零極點(diǎn)增益形式。

若要獲得系統(tǒng)函數(shù)H（z）的零極點(diǎn)分布圖，可直接應(yīng)用zplane函數(shù)，其語句格式為

Zplane（B，A）

其中，B與A分別表示H（z）的分子和分母多項式的系數(shù)向量。它的作用是在z平面上畫出單位圓、零點(diǎn)與極點(diǎn)。

離散系統(tǒng)中z變化建立了時域函數(shù)h（n）與z域函數(shù)H（z）之間的對應(yīng)關(guān)系。因此，z變化的函數(shù)H（z）從形式可以反映h（n）的部分內(nèi)在性質(zhì)。可根據(jù)系統(tǒng)的傳遞函數(shù)H（z）求單位沖激響應(yīng)h（n）的函數(shù)impz、filter等。

利用系統(tǒng)的頻率響應(yīng)，可以分析系統(tǒng)對各種頻率成分的響應(yīng)特性，并推出系統(tǒng)的特性（高通、低通、帶通、帶阻等）。

MATLAB提供了求離散時間系統(tǒng)頻響特性的函數(shù)freqz，調(diào)用freqz的格式主要有兩種。一種形式為

[H,w]= reqz(B,A,N)其中，B與A分別表示H（z）分子和分母多項式的系數(shù)向量；N為正整數(shù)，默認(rèn)值為512；返回值w包含[0，π]范圍內(nèi)的N個頻率等分點(diǎn)；返回值H則是離散時間系統(tǒng)頻率響應(yīng)在0~π范圍內(nèi)N個頻率處的值。另一種形式為

[H，w]= freqz（B，A，N，‘whole’）

與第一種方式不同之處在于角頻率的范圍由[0，π]擴(kuò)展到[0，2π]。

三、實(shí)驗(yàn)內(nèi)容與結(jié)果分析

已知LTI離散時間系統(tǒng)，要求由鍵盤實(shí)現(xiàn)系統(tǒng)參數(shù)輸入，并繪出幅頻和相頻響應(yīng)曲線和零極點(diǎn)分布圖，進(jìn)而分析系統(tǒng)的濾波特性和穩(wěn)定性。

（一）程序

b=[0.0528,0.797,0.1295,0.1295,0.797,0.0528];

a=[1,-1.8107,2.4947,-1.8801,0.9537,-0.2336];w=[0:20:500]*pi/500;

x1=0.0528+0.797*exp(-1*j*w)+0.1295*exp(-2*j*w)+0.1295*exp(-3*j*w)+0.797*exp(-4*j*w)+0.0528*exp(-5*j*w);

x2=1-1.8107*exp(-1*j*w)+2.4947*exp(-2*j*w)+1.8801*exp(-3*j*w)+0.9537*exp(-4*j*w)+0.2336*exp(-5*j*w);x22=x2+(x2==0)*eps;x=x1./x22;magx=abs(x);

angx=angle(x).*180/pi;

subplot(2,2,3);zplane(b,a);title('零極點(diǎn)圖');subplot(2,2,2);stem(w/pi,magx);title('幅度部分');ylabel('振幅');subplot(2,2,4);stem(w/pi,angx);

xlabel('以pi為單位的頻率');title('相位部分');ylabel('相位');

（二）波形圖

圖4-1 幅頻、相頻響應(yīng)曲線、零極點(diǎn)分布圖

實(shí)驗(yàn)六 IIR數(shù)字濾波器的設(shè)計

一、實(shí)驗(yàn)?zāi)康?/p>

從理論上講，任何的線性是不變（LTI）離散時間系統(tǒng)都可以看做一個數(shù)字濾波器，因此設(shè)計數(shù)字濾波器實(shí)際就是設(shè)計離散時間系統(tǒng)。數(shù)字濾波器你包括IIR（無限沖激響應(yīng)）和FIR（有限沖激響應(yīng)）型，在設(shè)計時通常采用不同的方法。

本實(shí)驗(yàn)通過使用MATLAB函數(shù)對數(shù)字濾波器進(jìn)行設(shè)計和和實(shí)現(xiàn)，要求掌握IIR數(shù)字巴特沃斯濾波器、數(shù)字切比雪夫?yàn)V波器的設(shè)計原理、設(shè)計方法和設(shè)計步驟；能根據(jù)給定的濾波器指標(biāo)進(jìn)行濾波器設(shè)計；同時也加深學(xué)生對數(shù)字濾波器的常用指標(biāo)和設(shè)計過程的理解。

二、實(shí)驗(yàn)原理

在IIR濾波器的設(shè)計中，常用的方法是：先根據(jù)設(shè)計要求尋找一個合適的模擬原型濾波器，然后根據(jù)一定的準(zhǔn)則將此模擬原型濾波器轉(zhuǎn)換為數(shù)字濾波器。

IIR濾波器的階數(shù)就等于所選的模擬原型濾波器的階數(shù)，所以其階數(shù)確定主要是在模擬原型濾波器中進(jìn)行的。

IIR數(shù)字濾波器的設(shè)計方法如下：（1）沖激響應(yīng)不變法。（2）雙線性變化法。

一般來說，在要求時域沖激響應(yīng)能模仿模擬濾波器的場合，一般使用沖激響應(yīng)不變法。沖激響應(yīng)不變法一個重要特點(diǎn)是頻率坐標(biāo)的變化是線性的，因此如果模擬濾波器的頻率響應(yīng)帶限于折疊頻率的話，則通過變換后濾波器的頻率響應(yīng)可不失真地反映原響應(yīng)與頻率的關(guān)系。

與沖激響應(yīng)不變法比較，雙線性變化的主要優(yōu)點(diǎn)是靠頻率的非線性關(guān)系得到s平面與z平面的單值一一對應(yīng)關(guān)系，整個值對應(yīng)于單位圓一周。所以從模擬傳遞函數(shù)可直接通過代數(shù)置換得到數(shù)字濾波器的傳遞函數(shù)。

MATLAB提供了一組標(biāo)準(zhǔn)的數(shù)字濾波器設(shè)計函數(shù)，大大簡化了濾波器的設(shè)計工程。

（1）butter。

（2）cheby1、cheby2。

三、實(shí)驗(yàn)內(nèi)容及分析

利用MATLAB編程方法或利用MATLAB中fdatool工具設(shè)計不同功能的IIR數(shù)字濾波器。

1、基于chebyshev I型模擬濾波器原型使用沖激不變轉(zhuǎn)換方法設(shè)計數(shù)字濾波器，要求參數(shù)為通帶截止頻率?p?0.4?；通帶最大衰減Ap?1dB；阻帶截止頻率?s?0.4?；阻帶最小衰減As?35dB。

程序:

wp=0.2*pi;

%通帶邊界頻率

ws=0.4*pi;

%阻帶截止頻率 rp=1;

%通帶最大衰減 rs=35;

%阻帶最小衰減

Fs=1000;

%?ùéè3é?ù??3?1000hz

[N,Wn]=cheb1ord(wp,ws,rp,rs,'s');

[Z,P,K]=cheby1(N,rp,Wn,'s');[H,W]=zp2tf(Z,P,K);

figure(1);freqs(H,W);[P,Q]=freqs(H,W);figure(2);plot(Q*Fs/(2*pi),abs(P));grid on;

xlabel('頻率/Hz');ylabel('幅度');

2、基于Butterworth型模擬濾波器原型使用雙線性變換方法設(shè)計數(shù)字濾波器的，要求參數(shù)為截止頻率?p?0.4?；通帶最大衰減Ap?1dB；阻帶截止頻率?s?0.25?；阻帶最小衰減AS?40dB。程序: wp=0.4*pi;ws=0.25*pi;rp=1;rs=40;fs=500;ts=1/fs;wp1=wp*ts;ws1=ws*ts;

wp2=2*fs*tan(wp1/2);ws2=2*fs*tan(ws1/2);

[N,Wn]=buttord(wp2,ws2,rp,rs,'s');[Z,P,K]=buttap(N);[Bap,Aap]=zp2tf(Z,P,K);[b,a]=lp2lp(Bap,Aap,Wn);[bz,az]=bilinear(b,a,fs);[H,W]=freqz(bz,az);subplot(2,1,1);plot(W/pi,abs(H));grid on;xlabel('頻率')ylabel('幅度')subplot(2,1,2);

plot(W/pi,20*log10(abs(H)));grid on;xlabel('頻率');ylabel('幅度(dB)');

實(shí)驗(yàn)七 FIR數(shù)字濾波器的設(shè)計

一、實(shí)驗(yàn)?zāi)康?/p>

掌握用窗函數(shù)設(shè)計FIR數(shù)字濾波的原理及其設(shè)計步驟；熟悉線性相位數(shù)字濾波器的特性。學(xué)習(xí)編寫數(shù)字濾波器的設(shè)計程序的方法，并能進(jìn)行正確編程；根據(jù)給定的濾波器指標(biāo)，給出設(shè)計步驟。

二、實(shí)驗(yàn)原理

如果系統(tǒng)的沖激響應(yīng)h（n）為已知，則系統(tǒng)的輸入輸出關(guān)系為

y（n）=x（n）*h（n）

對于低通濾波器，只要設(shè)計出低通濾波器的沖激響應(yīng)函數(shù)，就可以由式得到系統(tǒng)的輸出了。

但是將h（n）作為濾波器的脈沖響應(yīng)有兩個問題：一是它是無限長的；二是它是非因果的。對此，采取兩項措施：一是將h（n）截短；二是將其右移。

設(shè)計時，要根據(jù)阻帶的最小衰減和過渡帶寬度來選擇恰當(dāng)?shù)拇昂瘮?shù)類型和窗口長度N。常用的窗函數(shù)有矩形窗、海明窗和布萊克曼窗等。

窗函數(shù)設(shè)計FIR濾波器步驟如下：

（1）給定理想頻率響應(yīng)的幅頻特性和相頻特性；

（2）求理想單位脈沖響應(yīng)，在實(shí)際計算中，可對理想頻率響應(yīng)采樣。（3）根據(jù)過渡帶寬度和阻帶最小衰減，確定窗函數(shù)類型和窗口長度N；（4）求FIR濾波器單位脈沖響應(yīng)；

（5）分析幅頻特性，若不滿足要求，可適當(dāng)改變窗函數(shù)形式或長度N，重復(fù)上述設(shè)計過程，以得到滿意的結(jié)果。

三、實(shí)驗(yàn)內(nèi)容及分析

1、分別用海明窗和布萊克曼窗設(shè)計一個48階的FIR帶通濾波器，通帶為Wn??0.450.55?。程序1:海明窗設(shè)計

N=48;

Window=hamming(N+1);w1=0.45;w2=0.55;ws=[w1,w2];

b=fir1(N,ws/pi,Window);freqz(b,1,512);title('海明窗');grid on;

程序2:萊克曼窗設(shè)計

N=48;

Window=blackman(N+1);w1=0.45;w2=0.55;ws=[w1,w2];

b=fir1(N,ws/pi,Window);freqz(b,1,512);title('布萊克曼窗');grid on;

2、用矩形窗設(shè)計一個線性相位高通濾波器。其中He??jw?e?j????????00.3?????

0???0.3?程序: N=9;

alpha=(N-1)/2;Wc=0.7*pi;n=(0:8);i=n-alpha;i=i+(i==0)*eps;

h=(-1).^n.*sin((i).*Wc)./((i).*pi);%矩形窗函數(shù)設(shè)計的系統(tǒng)脈沖響應(yīng) w=(0:1:500)*2*pi/500;

H=h*exp(-j*n'*w);%矩形窗函數(shù)設(shè)計的頻響 magH=abs(H);% 矩形窗函數(shù)設(shè)計的振幅 subplot(211);stem(n,h);

axis([0,8,-0.4,0.4]);title('矩形窗設(shè)計h(n)');line([0,10],[0,0]);xlabel('n');ylabel('h');subplot(212);plot(w/pi,magH);

xlabel('以pi為單位的頻率');ylabel('H振幅');axis([0,2,0,1.7]);title('矩形窗設(shè)計振幅譜');

實(shí)驗(yàn)心得體會：

這次實(shí)驗(yàn)使我進(jìn)一步加深了對MATLAB軟件的使用。從上次的信號系統(tǒng)實(shí)驗(yàn)的初步使用到這一次的深入了解，有了更深刻的認(rèn)識。對這種語言環(huán)境也有了新的了解。

在實(shí)驗(yàn)的過程中，我對數(shù)字濾波器的整個過程有了很好的理解和掌握。IIR數(shù)字濾波器的設(shè)計讓我知道了巴特沃思濾波器和切比雪夫?yàn)V波器的頻率特性，還有雙線性變換及脈沖響應(yīng)不變法設(shè)計的濾波器的頻率特性。做這兩個實(shí)驗(yàn)的時候程序有點(diǎn)困難，但經(jīng)過細(xì)心的改寫圖形最終出來了。FIR數(shù)字濾波器的設(shè)計出來的是兩種窗的圖形，通過兩種窗的比較，我了解了他們各自的特點(diǎn)，幅頻和相頻特性。

最后，感謝張老師對我的諄諄教導(dǎo)！

第四篇：數(shù)字圖像處理實(shí)驗(yàn)報告

數(shù)字圖像處理

實(shí)驗(yàn)報告

目錄

1.數(shù)字圖像處理簡介

2.實(shí)驗(yàn)?zāi)康?.實(shí)驗(yàn)內(nèi)容

4.實(shí)驗(yàn)結(jié)果及代碼展示

5.算法綜述

6.Matlab優(yōu)勢

7.總結(jié)

8.存在問題

一、數(shù)字圖像處理簡介

圖像處理，是對圖像進(jìn)行分析、加工、和處理，使其滿足視覺、心理以及其他要求的技術(shù)。圖像處理是信號處理在圖像域上的一個應(yīng)用。目前大多數(shù)的圖像是以數(shù)字形式存儲，因而圖像處理很多情況下指數(shù)字圖像處理。此外，基于光學(xué)理論的處理方法依然占有重要的地位。

圖像處理是信號處理的子類，另外與計算機(jī)科學(xué)、人工智能等領(lǐng)域也有密切的關(guān)系。

傳統(tǒng)的一維信號處理的方法和概念很多仍然可以直接應(yīng)用在圖像處理上，比如降噪、量化等。然而，圖像屬于二維信號，和一維信號相比，它有自己特殊的一面，處理的方式和角度也有所不同。

二、實(shí)驗(yàn)?zāi)康?/p>

鞏固所學(xué)知識，提高所學(xué)能力

三、實(shí)驗(yàn)內(nèi)容

利用matlab的GUI程序設(shè)計一個簡單的圖像處理程序，并含有如下基本功能： 1.讀入一幅RGB圖像，變換為灰度圖像和二值圖像，并在同一個窗口內(nèi)分成三個子窗口來分別顯示RGB圖像和灰度圖像，注上文字標(biāo)題 2.對給定圖像進(jìn)行旋轉(zhuǎn)

3．對給定的圖像添加噪聲（椒鹽噪聲、高斯噪聲）

四、實(shí)驗(yàn)結(jié)果及代碼展示

1.軟件設(shè)計界面

2.各模塊功能展示以及程序代碼

（1）讀入一幅RGB圖像，變換為灰度圖像和二值圖像，并在同一個窗口內(nèi)分成三個子窗口來分別顯示RGB圖像和灰度圖像，注上文字標(biāo)題

效果展示：

代碼：

a = imread('C:Documents and SettingsAdministrator桌面數(shù)字圖像舞美.JPG');

i = rgb2gray(a);I = im2bw(a,0.5);

subplot(3,1,1);imshow(a);title('源圖像')subplot(3,1,2);imshow(i);title('灰度圖像')subplot(3,1,3);imshow(I);title('二值圖像')

（2）圖像旋轉(zhuǎn) 原圖

效果展示：

代碼：

clc;clear all;close all;

Img=imread('D:My DocumentsMy Pictures5.JPG');Img=double(Img);[h w]=size(Img);alpha=pi/4;

wnew=w*cos(alpha)+h*sin(alpha);hnew=w*sin(alpha)+h*cos(alpha);wnew=ceil(wnew);

hnew=ceil(hnew);u0=w*sin(alpha);

T=[cos(alpha),sin(alpha);-sin(alpha),cos(alpha)];Imgnew2=zeros(hnew,wnew);Imgnew1=zeros(hnew,wnew);for u=1:hnew

for v=1:wnew

tem=T*([u;v]-[u0;0]);x=tem(1);y=tem(2);if x>=1&&x<=h&&y>=1&&y<=w x_low=floor(x);x_up=ceil(x);y_low=floor(y);y_up=ceil(y);if(x-x_low)<=(x_up-x)x=x_low;

else

x=x_up;

end

if(y-y_low)<=(y_up-y)y=y_low;

else

y=y_up;

end

p1=Img(x_low,y_low);p2=Img(x_up,y_low);p3=Img(x_low,y_low);p4=Img(x_up,y_up);s=x-x_low;t=y-y_low;Imgnew1(u,v)=Img(x,y);

Imgnew2(u,v)=(1-s)*(1-t)*p1+(1-s)*t*p3+(1-t)*s*p2+s*t*p4;end

end end

figure;imshow(Imgnew2,[]);B=imrotate(Img,alpha/pi*180);figure;imshow(B,[]);

（3）對給定的圖像添加噪聲（斑點(diǎn)噪聲、高斯噪聲）效果展示：

代碼：

I= imread('D:My DocumentsMy Pictures5.JPG');figure,subplot(211);imshow(I);title('原圖');J1=imnoise(I,'gaussian',0,0.02);

subplot(223);imshow(J);title('添加高斯噪聲');J=imnoise(I,'speckle',0.04);

subplot(224);imshow(J);title('添加斑點(diǎn)噪聲');

五、算法綜述 灰度圖像：

一幅完整的圖像，是由紅色、綠色、藍(lán)色三個通道組成的。紅色、綠色、藍(lán)色三個通道的縮覽圖都是以灰度顯示的。用不同的灰度色階來表示“ 紅，綠，藍(lán)”在圖像中的比重。通道中的純白，代表了該色光在此處為最高亮度，亮度級別是255。

通道是整個Photoshop顯示圖像的基礎(chǔ)。色彩的變動，實(shí)際上就是間接在對通道灰度圖進(jìn)行調(diào)整。通道是Photoshop處理圖像的核心部分，所有的色彩調(diào)整工具都是圍繞在這個核心周圍使用的。

在計算機(jī)領(lǐng)域中，灰度數(shù)字圖像是每個像素只有一個采樣顏色的圖像。這類圖像通常顯示為從最暗黑色到最亮的白色的灰度，盡管理論上這個采樣可以任何顏色的不同深淺，甚至可以是不同亮度上的不同顏色。灰度圖像與黑白圖像不同，在計算機(jī)圖像領(lǐng)域中黑白圖像只有黑色與白色兩種顏色；灰度圖像在黑色與白色之間還有許多級的顏色深度。但是，在數(shù)字圖像領(lǐng)域之外，“黑白圖像”也表示“灰度圖像”，例如灰度的照片通常叫做“黑白照片”。在一些關(guān)于數(shù)字圖像的文章中單色圖像等同于灰度圖像，在另外一些文章中又等同于黑白圖像。灰度圖像經(jīng)常是在單個電磁波頻譜如可見光內(nèi)測量每個像素的亮度得到的。

用于顯示的灰度圖像通常用每個采樣像素 8 位的非線性尺度來保存，這樣可以有 256 級灰度。這種精度剛剛能夠避免可見的條帶失真，并且非常易于編程。在醫(yī)學(xué)圖像與遙感圖像這些技術(shù)應(yīng)用中經(jīng)常采用更多的級數(shù)以充分利用每個采樣 10 或 12 位的傳感器精度，并且避免計算時的近似誤差。在這樣的應(yīng)用領(lǐng)域每個采樣 16 位即 65536 級得到流行。

二值圖像：

是指每個像素不是黑就是白，其灰度值沒有中間過渡的圖像。二值圖像一般用來描述文字或者圖形，其優(yōu)點(diǎn)是占用空間少，缺點(diǎn)是，當(dāng)表示人物，風(fēng)景的圖像時，二值圖像只能描述其輪廓，不能描述細(xì)節(jié)。這時候要用更高的灰度級。

二值圖像是每個像素只有兩個可能值的數(shù)字圖像。人們經(jīng)常用黑白、B&W、單色圖像表示二值圖像，但是也可以用來表示每個像素只有一個采樣值的任何圖像，例如灰度圖像等。

二值圖像中所有的像素只能從0和1這兩個值中取，因此在MATLAB中，二值圖像用一個由0和1組成的二維矩陣表示。這兩個可取的值分別對應(yīng)于關(guān)閉和打開，關(guān)閉表征該像素處于背景，而打開表征該像素處于前景。以這種方式來操作圖像可以更容易識別出圖像的結(jié)構(gòu)特征。二值圖像操作只返回與二值圖像的形式或結(jié)構(gòu)有關(guān)的信息，如果希望對其他類型的圖像進(jìn)行同樣的操作，則首先要將其轉(zhuǎn)換為二進(jìn)制的圖像格式，可以通過調(diào)用MATLAB提供的im2bw（）來實(shí)現(xiàn)。

二值圖像經(jīng)常出現(xiàn)在數(shù)字圖像處理中作為圖像掩碼或者在圖像分割、二值化和dithering的結(jié)果中出現(xiàn)。一些輸入輸出設(shè)備，如激光打印機(jī)、傳真機(jī)、單色計算機(jī)顯示器等都可以處理二值圖像。

二值圖像經(jīng)常使用位圖格式存儲。

二值圖像可以解釋為二維整數(shù)格Z，圖像變形處理領(lǐng)域很大程度上就是受到這個觀點(diǎn)啟發(fā)。

圖像旋轉(zhuǎn)：

圖像旋轉(zhuǎn)是指圖像以某一點(diǎn)為中心旋轉(zhuǎn)一定的角度,形成一幅新的圖像的過程。當(dāng)然這個點(diǎn)通常就是圖像的中心。既然是按照中心旋轉(zhuǎn),自然會有這樣一個屬性：旋轉(zhuǎn)前和旋轉(zhuǎn)后的點(diǎn)離中心的位置不變.根據(jù)這個屬性,我們可以得到旋轉(zhuǎn)后的點(diǎn)的坐標(biāo)與原坐標(biāo)的對應(yīng)關(guān)系。由于原圖像的坐標(biāo)是以左上角為原點(diǎn)的,所以我們先把坐標(biāo)轉(zhuǎn)換為以圖像中心為原點(diǎn)。假設(shè)原圖像的寬為w,高為h,（x0,y0）為原坐標(biāo)內(nèi)的一點(diǎn),轉(zhuǎn)換坐標(biāo)后的點(diǎn)為（x1,y1）。那么不難得到： x1 = x0-w/2;y1 =-y0 + h/2;在新的坐標(biāo)系下,假設(shè)（x0,y0）距離原點(diǎn)的距離為r,點(diǎn)與原點(diǎn)之間的連線與x軸的夾角為b,旋轉(zhuǎn)的角度為a,旋轉(zhuǎn)后的點(diǎn)為（x1,y1）

噪聲：

是電路或系統(tǒng)中不含信息量的電壓或電流。在工業(yè)與自然界中，存在著各種干擾源（噪聲源），如大功率電力電子器件的接入、大功率用電設(shè)備的開啟與斷開、雷擊閃電等都會使空間電場和磁場產(chǎn)生有序或無序的變化，這些都是干擾源（或噪聲源）。這些源產(chǎn)生的電磁波或尖峰脈沖通過磁、電耦合或是通過電源線等路徑進(jìn)入放大電路，各種電氣設(shè)備，形成各種形式的干擾。

斑點(diǎn)噪聲：

斑點(diǎn)噪聲是SAR成像系統(tǒng)的一大特色，源自基本分辨單元內(nèi)地物的隨機(jī)散射，在圖像上表現(xiàn)為信號相關(guān)（如在空間上相關(guān)）的小斑點(diǎn)，它既降低了圖像的畫面質(zhì)量，又嚴(yán)重影響圖像的自動分割、分類、目標(biāo)檢測以及其它定量專題信息的提取。

SAR圖像斑點(diǎn)噪聲的去除一方面要抑制圖像均勻區(qū)域斑點(diǎn)噪聲，另一方面要保持圖像邊緣和紋理細(xì)節(jié)信息。SAR斑點(diǎn)噪聲的抑制可通過非相干多視處理，也可使用空間域?yàn)V波實(shí)現(xiàn)。非相干多視處理會降低圖像的地面分辨率。因此，涌現(xiàn)出了一系列空間域?yàn)V波方法，如均值濾波、中值濾波、Lee濾波、Kuan濾波、Frost濾波、Sigma濾波以及Gamma Map濾波等。但這類算法存在自身無法克服的矛盾：一方面為增強(qiáng)斑點(diǎn)去噪效果需選較大的濾波窗口，另一方面為保持圖像的實(shí)際分辨率要求所選的窗口較小。

高斯噪聲：

所謂高斯噪聲是指它的概率密度函數(shù)服從高斯分布（即正態(tài)分布）的一類噪聲。如果一個噪聲，它的幅度分布服從高斯分布，而它的功率譜密度又是均勻分布的，則稱它為高斯白噪聲。高斯白噪聲的二階矩不相關(guān)，一階矩為常數(shù)，是指先后信號在時間上的相關(guān)性。高斯白噪聲包括熱噪聲和散粒噪聲。

實(shí)驗(yàn)中是通過MATLAB自帶的函數(shù)產(chǎn)生噪聲，各函數(shù)如下： J1=imnoise(I,'salt & pepper',0.05);%添加椒鹽噪聲

J2=imnoise(I,'gaussian',0,0.03);

%添加均值為0，方差為0.03的高斯噪聲。

六、Matlab優(yōu)勢

MATLAB是一個包含大量算法的集合。其可以快捷的實(shí)現(xiàn)用戶所需的各種計算功能。函數(shù)中所使用的算法都是科研和工程計算中的最新研究成果，而前經(jīng)過了各種優(yōu)化和差錯處理。在通常情況下，可以用它來代替底層編程語言，如C和C++。在計算要求相同的情況下，使用MATLAB的編程工作量會大大減少。MATLAB的這些函數(shù)集包括從最簡單最基本的函數(shù)到諸如矩陣，特征向量、快速傅立葉變換的復(fù)雜函數(shù)。函數(shù)所能解決的問題其大致包括矩陣運(yùn)算和線性方程組的求解、微分方程及偏微分方程的組的求解、符號運(yùn)算、傅立葉變換和數(shù)據(jù)的統(tǒng)計分析、工程中的優(yōu)化問題、稀疏矩陣運(yùn)算、復(fù)數(shù)的各種運(yùn)算、三角函數(shù)和其他初等數(shù)學(xué)運(yùn)算、多維數(shù)組操作以及建模動態(tài)仿真等。圖形處理功能

圖形處理功能MATLAB自產(chǎn)生之日起就具有方便的數(shù)據(jù)可視化功能，以將向量和矩陣用圖形表現(xiàn)出來，并且可以對圖形進(jìn)行標(biāo)注和打印。高層次的作圖包括二維和三維的可視化、圖象處理、動畫和表達(dá)式作圖。可用于科學(xué)計算和工程繪圖。新版本的MATLAB對整個圖形處理功能作了很大的改進(jìn)和完善，使它不僅在一般數(shù)據(jù)可視化軟件都具有的功能（例如二維曲線和三維曲面的繪制和處理等）方面更加完善，而且對于一些其他軟件所沒有的功能（例如圖形的光照處理、色度處理以及四維數(shù)據(jù)的表現(xiàn)等），MATLAB同樣表現(xiàn)了出色的處理能力。同時對一些特殊的可視化要求，例如圖形對話等，MATLAB也有相應(yīng)的功能函數(shù)，保證了用戶不同層次的要求。另外新版本的MATLAB還著重在圖形用戶界面（GUI）的制作上作了很大的改善，對這方面有特殊要求的用戶也可以得到滿足

模塊集合工具箱

MATLAB對許多專門的領(lǐng)域都開發(fā)了功能強(qiáng)大的模塊集和工具箱。一般來說，它們都是由特定領(lǐng)域的專家開發(fā)的，用戶可以直接使用工具箱學(xué)習(xí)、應(yīng)用和評估不同的方法而不需要自己編寫代碼。目前，MATLAB已經(jīng)把工具箱延伸到了科學(xué)研究和工程應(yīng)用的諸多領(lǐng)域，諸如數(shù)據(jù)采集、數(shù)據(jù)庫接口、概率統(tǒng)計、樣條擬合、優(yōu)化算法、偏微分方程求解、神經(jīng)網(wǎng)絡(luò)、小波分析、信號處理、圖像處理、系統(tǒng)辨識、控制系統(tǒng)設(shè)計、LMI控制、魯棒控制、模型預(yù)測、模糊邏輯、金融分析、地圖工具、非線性控制設(shè)計、實(shí)時快速原型及半物理仿真、嵌入式系統(tǒng)開發(fā)、定點(diǎn)仿真、DSP與通訊、電力系統(tǒng)仿真等，都在工具箱（Toolbox）家族中有了自己的一席之地。

七、總結(jié)

運(yùn)用matlab軟件對圖像進(jìn)行處理，讓我鞏固了之前所學(xué)的知識，同時也在這次作業(yè)中更加了解到matlab語言在生活中的運(yùn)用環(huán)境和掌握這門語言的重要性

八、存在問題

1.在進(jìn)行圖像增強(qiáng)時要不要講圖像先進(jìn)行平滑處理？ 2.如何增加這個算法的準(zhǔn)確度

3.在此次作業(yè)中，為何添加椒鹽噪聲時無法顯示

第五篇：數(shù)字圖像處理實(shí)驗(yàn)報告

數(shù)字圖像處理

實(shí)驗(yàn)報告

班級:通信103 學(xué)號:201027201 姓名:計富威 指導(dǎo)教師:孫潔

實(shí)驗(yàn)一 MATLAB數(shù)字圖像處理初步

一、實(shí)驗(yàn)?zāi)康呐c要求

1．熟悉及掌握在MATLAB中能夠處理哪些格式圖像。2．熟練掌握在MATLAB中如何讀取圖像。

3．掌握如何利用MATLAB來獲取圖像的大小、顏色、高度、寬度等等相關(guān)信息。

4．掌握如何在MATLAB中按照指定要求存儲一幅圖像的方法。5．圖像間如何轉(zhuǎn)化。

二、實(shí)驗(yàn)內(nèi)容及步驟

1．利用imread()函數(shù)讀取一幅圖像，假設(shè)其名為”第一個.tif”，存入一個數(shù)組中；

>>I=imread('第一個.tif');2．利用whos命令提取該讀入圖像”第一個.tif”的基本信息； >>whos I 3．利用imshow()函數(shù)來顯示這幅圖像； >>imshow(I);

第一個.tif 4．利用imfinfo函數(shù)來獲取圖像文件的壓縮，顏色等等其他的詳細(xì)信息；

>>imfinfo('第一個.tif');5．利用imwrite()函數(shù)來壓縮這幅圖象，將其保存為一幅壓縮了像素的jpg文件,設(shè)為flower.jpg；語法：imwrite(原圖像，新圖像，‘quality’,q), q取0-100。

>>imwrite(I,'第一個.jpg','quality',50)6．同樣利用imwrite()函數(shù)將最初讀入的tif圖象另存為一幅bmp圖像，設(shè)為flower.bmp。>>imwrite(I,'第一個.bmp');7．用imread()讀入圖像：Lenna256.jpg 和camemaman.jpg； >>b=imread('lena256.bmp');>>c=imread('cameraman.tif');8．用imfinfo()獲取圖像Lenna256.jpg和camemaman.jpg 的大小； >>imfinfo('lena256.bmp');>>imfinfo('cameraman.tif');9．用figure,imshow()分別將Lenna256.jpg和camemaman.jpg顯示出來，觀察兩幅圖像的質(zhì)量。>>figure >>imshow(b);>>figure >>imshow(c);

（Lena256.jpg圖像截圖）

（cameraman.jpg圖像截圖）

10．用im2bw將一幅灰度圖像轉(zhuǎn)化為二值圖像，并且用imshow顯示出來觀察圖像的特征。>> d=im2bw(b);>>figure >>imshow(b);>>figure >>imshow(d);

（二值化截圖）

三、實(shí)驗(yàn)總結(jié)

通過實(shí)驗(yàn)MatLab軟件的基本使用有了基本的了解，學(xué)會了使用MatLab軟件來讀取一個特定格式的圖像，并通過相關(guān)的命令語句對圖像進(jìn)行格式轉(zhuǎn)換、圖像壓縮、二值化等的處理，掌握了利用MATLAB來獲取圖像的大小、顏色、高度、寬度等等相關(guān)信息，掌握在MATLAB中如何通過imshow()語句來讀取圖像等等。

第二 圖像基本運(yùn)算

一、實(shí)驗(yàn)?zāi)康?/p>

1．了解圖像的算術(shù)運(yùn)算在數(shù)字圖像處理中的初步應(yīng)用。2．體會圖像算術(shù)運(yùn)算處理的過程和處理前后圖像的變化。

二、實(shí)驗(yàn)原理

圖像的代數(shù)運(yùn)算是圖像的標(biāo)準(zhǔn)算術(shù)操作的實(shí)現(xiàn)方法，是兩幅輸入圖像之間進(jìn)行的點(diǎn)對點(diǎn)的加、減、乘、除運(yùn)算后得到輸出圖像的過程。如果輸入圖像為A(x,y)和B(x,y)，輸出圖像為C(x,y)，則圖像的代數(shù)運(yùn)算有如下四種形式：

C(x,y)= A(x,y)+ B(x,y)C(x,y)= A(x,y)-B(x,y)C(x,y)= A(x,y)* B(x,y)C(x,y)= A(x,y)/ B(x,y)

三、實(shí)驗(yàn)步驟 1．圖像的加法運(yùn)算

在MATLAB中，如果要進(jìn)行兩幅圖像的加法，或者給一幅圖像加上一個常數(shù)，可以調(diào)用imadd函數(shù)來實(shí)現(xiàn)。imadd函數(shù)將某一幅輸入圖像的每一個像素值與另一幅圖像相應(yīng)的像素值相加，返回相應(yīng)的像素值之和作為輸出圖像。imadd函數(shù)的調(diào)用格式如下：

Z = imadd（X，Y）首先讀入兩幅圖像

>>a=imread('第二個原圖1.jpg');>>b=imread('第二個原圖2.jpg')通過一個加法操作：>> c=imadd(a,b);

給圖像的每一個像素加上一個常數(shù)可以使圖像的亮度增加。如截圖

第一張為原圖，第二張為亮度加50，第三張為亮度減50 2．圖像的減法運(yùn)算

在MATLAB中，使用imsubtract函數(shù)可以將一幅圖像從另一幅圖像中減去，或者從一幅圖像中減去一個常數(shù)。imsubtract函數(shù)將一幅輸入圖像的像素值從另一幅輸入圖像相應(yīng)的像素值中減去，再將這個結(jié)果作為輸出圖像相應(yīng)的像素值。imsubtract函數(shù)的調(diào)用格式如下：

Z = imsubtract(X,Y)； 讀入一幅畫后通過減法 >>a3=imsubtract(a,50);

3.圖像的乘法運(yùn)算

在MATLAB中，使用immultiply函數(shù)實(shí)現(xiàn)兩幅圖像的乘法。immultiply函數(shù)將兩幅圖像相應(yīng)的像素值進(jìn)行元素對元素的乘法操作（MATLAB點(diǎn)乘），并將乘法的運(yùn)算結(jié)果作為輸出圖形相應(yīng)的像素值。immulitply函數(shù)的調(diào)用格式如下：

Z = immulitply(X,Y)讀入一幅圖后通過乘法操作 >> a=imread('cameraman.tif');>> b=immultiply(a,1.5);

4．圖像的除法運(yùn)算

在MATLAB中使用imdivide函數(shù)進(jìn)行兩幅圖像的除法。imdivide函數(shù)對兩幅輸入圖像的所有相應(yīng)像素執(zhí)行元素對元素的除法操作（點(diǎn)除），并將得到的結(jié)果作為輸出圖像的相應(yīng)像素值。imdivide函數(shù)的調(diào)用格式如下：

Z = imdivide(X,Y)讀入一幅圖后通過除法操作

四、實(shí)驗(yàn)總結(jié)

通過對圖像的四則運(yùn)算了結(jié)圖像的不同變化過程，對軟件的進(jìn)一步使用也有了更加深刻的認(rèn)識。

實(shí)驗(yàn)三 圖像增強(qiáng)—空域?yàn)V波

一、實(shí)驗(yàn)?zāi)康?/p>

進(jìn)一步了解MatLab軟件/語言，學(xué)會使用MatLab對圖像作濾波處理，使學(xué)生有機(jī)會掌握濾波算法，體會濾波效果。

了解幾種不同濾波方式的使用和使用的場合，培養(yǎng)處理實(shí)際圖像的能力，并為課堂教學(xué)提供配套的實(shí)踐機(jī)會。

二、實(shí)驗(yàn)設(shè)備與軟件

(1)IBM-PC計算機(jī)系統(tǒng)

(2)MatLab軟件/語言包括圖像處理工具箱(Image Processing Toolbox)

(3)實(shí)驗(yàn)所需要的圖片

三、實(shí)驗(yàn)內(nèi)容與步驟

a)調(diào)入并顯示原始圖像“原圖像.jpg”。>>I=imread('原圖像.jpg');b)利用imnoise命令在圖像“原圖像.jpg”上加入高斯(gaussian)噪聲

>>J = imnoise(I,'gauss',0.02);

%添加高斯噪聲 c)利用預(yù)定義函數(shù)fspecial命令產(chǎn)生平均(average)濾波器 d）分別采用3x3和5x5的模板，分別用平均濾波器以及中值濾波器，對加入噪聲的圖像進(jìn)行處理并觀察不同噪聲水平下，上述濾波器處理的結(jié)果；

>>ave1=fspecial('average',3);

%產(chǎn)生3×3的均值模版 >>ave2=fspecial('average',5);

%產(chǎn)生5×5的均值模版 >>K = filter2(ave1,J)/255;

%均值濾波3×3 >>L = filter2(ave2,J)/255;

%均值濾波5×5 e）選擇不同大小的模板，對加入某一固定噪聲水平噪聲的圖像進(jìn)行處理，觀察上述濾波器處理的結(jié)果。

>>M = medfilt2(J,[3 3]);

%中值濾波3×3模板 >>N = medfilt2(J,[4 4]);

%中值濾波4×4模板

f）利用imnoise命令在圖像Sample2-1.jpg 上加入椒鹽噪聲(salt & pepper)>>J = imnoise(I,'salt& pepper',0.02);

%添加椒鹽噪聲

四、實(shí)驗(yàn)總結(jié)

椒鹽噪聲是由圖像傳感器，傳輸信道，解碼處理等產(chǎn)生的黑白相間的亮暗點(diǎn)噪聲。椒鹽噪聲往往由圖像切割引起，去除脈沖干擾及椒鹽噪聲最常用的算法是中值濾波。椒鹽噪聲是指兩種噪聲，一種是鹽噪聲，另一種是胡椒噪聲。鹽=白色，椒=黑色。前者是高灰度噪聲，后者屬于低灰度噪聲。一般兩種噪聲同時出現(xiàn)，呈現(xiàn)在圖像上就是黑白雜點(diǎn)。這點(diǎn)我們通過實(shí)驗(yàn)結(jié)果可以明顯看到。中值濾波對于濾除圖像的椒鹽噪聲非常有效。

實(shí)驗(yàn)四圖像分割

一、實(shí)驗(yàn)?zāi)康?/p>

使用MatLab 軟件進(jìn)行圖像的分割。使學(xué)生通過實(shí)驗(yàn)體會一些主要的分割算子對圖像處理的效果，以及各種因素對分割效果的影響。

二、實(shí)驗(yàn)要求

要求學(xué)生能夠自行評價各主要算子在無噪聲條件下和噪聲條件下的分割性能。能夠掌握分割條件(閾值等)的選擇。完成規(guī)定圖像的處理并要求正確評價處理結(jié)果，能夠從理論上作出合理的解釋。

三、實(shí)驗(yàn)內(nèi)容與步驟

(1)使用Roberts 算子的圖像分割實(shí)驗(yàn)，使用的原圖是cameraman.jpg,截圖如下

(2)使用Prewitt 算子的圖像分割實(shí)驗(yàn) 截圖如下

(3)使用Sobel 算子的圖像分割實(shí)驗(yàn)

(4)使用LoG(拉普拉斯-高斯)算子的圖像分割實(shí)驗(yàn)

四、實(shí)驗(yàn)結(jié)果

對Roberts算子、Prewitt 算子、Sobel 算子、LoG(拉普拉斯-高斯)算子的運(yùn)算對圖像的結(jié)果有了基本的認(rèn)識，加深學(xué)習(xí)效果。

實(shí)驗(yàn)五 形態(tài)學(xué)運(yùn)算

1、實(shí)驗(yàn)?zāi)康?/p>

學(xué)習(xí)常見的數(shù)學(xué)形態(tài)學(xué)運(yùn)算基本方法，了解腐蝕、膨脹、開運(yùn)算、閉運(yùn)算取得的效果，培養(yǎng)處理實(shí)際圖像的能力，并為課堂教學(xué)提供配套的實(shí)踐機(jī)會。

2、實(shí)驗(yàn)要求

利用MatLab工具箱中關(guān)于數(shù)學(xué)形態(tài)學(xué)運(yùn)算的函數(shù)，計算本指導(dǎo)書中指定二值圖像進(jìn)行處理。

3、實(shí)驗(yàn)設(shè)備與軟件

1.LC-PC計算機(jī)系統(tǒng)

2.MatLab軟件/語言包括圖像處理工具箱(Image Processing Toolbox)3.實(shí)驗(yàn)所需要的圖片

4、實(shí)驗(yàn)內(nèi)容與步驟

1.調(diào)入并顯示圖像“原圖.jpg”； 2.調(diào)入并顯示圖像“原圖.jpg”；

3.選取合適的閾值，得到二值化圖像“原圖.jpg”； >>bw = im2bw(I,level);

%二值化 4.設(shè)置結(jié)構(gòu)元素；

5.對得到的二值圖像“原圖.jpg”進(jìn)行腐蝕運(yùn)算； >>BW2 = imerode(bw,SE1);

%腐蝕 6.對得到的二值圖像“原圖.jpg”進(jìn)行膨脹運(yùn)算； >>BW1 = imdilate(bw,SE);

%膨脹 7.對得到的二值圖像“原圖.jpg”進(jìn)行開運(yùn)算；

>>BW3 = bwmorph(bw, 'open');

%開運(yùn)算 8.對得到的二值圖像“原圖.jpg”進(jìn)行閉運(yùn)算； >>BW4 = bwmorph(bw, 'close');

%閉運(yùn)算 9.將兩種處理方法的結(jié)果作比較；

五、實(shí)驗(yàn)總結(jié)

通過本次實(shí)驗(yàn)，學(xué)習(xí)了常見的數(shù)學(xué)形態(tài)學(xué)運(yùn)算基本方法，了解腐蝕、膨脹、開運(yùn)算、閉運(yùn)算取得的效果，培養(yǎng)處理實(shí)際圖像的能力，通過自己動手的實(shí)驗(yàn)，對課本上的知識有了更加深刻的理解。