醫(yī)學(xué)圖形圖像處理識(shí)別技術(shù)與醫(yī)學(xué)信息分析.ppt

《醫(yī)學(xué)圖形圖像處理識(shí)別技術(shù)與醫(yī)學(xué)信息分析.ppt》由會(huì)員分享，可在線閱讀，更多相關(guān)《醫(yī)學(xué)圖形圖像處理識(shí)別技術(shù)與醫(yī)學(xué)信息分析.ppt（48頁珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

醫(yī)學(xué)圖形圖像處理識(shí)別技術(shù)與醫(yī)學(xué)信息分析 目錄內(nèi)容 序醫(yī)學(xué)信息學(xué)簡述第一講醫(yī)學(xué)信號(hào)處理概述第二講醫(yī)學(xué)圖形信號(hào)的采集 存儲(chǔ)與處理第三講Mathlab與圖形信號(hào)的處理識(shí)別第四講醫(yī)學(xué)圖像信號(hào)的采集 存儲(chǔ)格式 處理方法第五講Mathlab Vtk Mtk與醫(yī)學(xué)圖像的處理第六講醫(yī)學(xué)信息挖掘第七講用于教學(xué)的醫(yī)學(xué)圖片處理技術(shù) 作業(yè) 第三講Mathlab與圖形信號(hào)的處理識(shí)別 MATLAB軟件的使用介紹圖形信號(hào)的處理與識(shí)別程序介紹臨床電生理數(shù)據(jù)分析實(shí)驗(yàn)?zāi)M軟件介紹實(shí)驗(yàn) 一 MATLAB軟件的使用介紹 1 MATLAB的簡介 1 MATLAB 高級(jí)交互式軟件包 2 MATLAB典型應(yīng)用數(shù)學(xué)計(jì)算算法推導(dǎo)建模和仿真模擬數(shù)據(jù)分析和結(jié)果的可視化工程圖形繪制應(yīng)用程序開發(fā) 包括用戶圖形界面的建立 3 MATLAB產(chǎn)生的歷史背景調(diào)用EISPACK和LINPACK的FORTRAN子程序庫 EISPACK和LINPACK的接口程序 MATLAB 發(fā)展成為適合多學(xué)科 多種工作平臺(tái)的功能強(qiáng)大的大型軟件 美國MathWorks公司設(shè)計(jì)的 科學(xué)和工程計(jì)算的高級(jí)交互式軟件包 CleveMoler NewMexico JohnLittle MATLAB軟件 續(xù)1 3 MATLAB的語言特點(diǎn)語言簡潔緊湊 使用方便靈活 庫函數(shù)極其豐富 運(yùn)算符豐富 MATLAB既具有結(jié)構(gòu)化的控制語句 又有面向?qū)ο缶幊痰奶匦?程序限制不嚴(yán)格 程序設(shè)計(jì)自由度大 程序的可移植性很好 基本上不做修改就可以在各種型號(hào)的計(jì)算機(jī)和操作系統(tǒng)上運(yùn)行 MATLAB的圖形功能強(qiáng)大 MATLAB的缺點(diǎn)是程序的執(zhí)行速度較慢 功能強(qiáng)大的工具箱 源程序的開放性 如 解方程組Ax b x A b MATLAB軟件 續(xù)2 2 MATLA的安裝進(jìn)入MATLAB軟件安裝程序目錄 然后單擊 Setup 程序 在安裝程序的引導(dǎo)下 逐步完成以下操作 1 選擇 Install 安裝或 Updatelicense 升級(jí)許可證 單擊 Next 然后輸入姓名 公司與序列號(hào) 單擊 Next 2 選擇接受許可證協(xié)議 單擊 Yes 單擊 Next 然后選擇 Typical 典型安裝 單擊 Next 3 選擇或輸入安裝目的路徑 單擊 Next 后單擊 Yes 4 單擊 Install 系統(tǒng)開始安裝 MATLAB軟件 續(xù)3 5 在提示插入第2張CD時(shí) 插入CD2 然后單擊 OK 接著在文件擴(kuò)展名警示對(duì)話框內(nèi)單擊 Yestoall 6 單擊 Next 接著單擊 Finish 7 如果在單擊 Finish 按鈕之前 選擇了 Start 選項(xiàng) 則此時(shí)進(jìn)入MATLAB工作環(huán)境 MATLAB軟件 續(xù)4 3 MATLAB啟動(dòng)與退出啟動(dòng)MATLAB 進(jìn)入MATLAB工作環(huán)境 啟動(dòng)MATLAB MATLAB軟件 續(xù)5 結(jié)束MATLAB 退出MATLAB工作環(huán)境鍵入exit鍵入quit直接關(guān)閉MATLAB的命令視窗 退出MATLAB MATLAB軟件 續(xù)6 1 4 MATLAB變量與基本輸入輸出變量命名的規(guī)則 第一個(gè)字母必須是英文字母 字母間不可留空格 最多只能有19個(gè)字母 MATLAB會(huì)忽略多余字母 變量定義與輸入輸出向量與標(biāo)量矩陣 行和列組成的數(shù)據(jù)表行向量 矩陣只有一行時(shí) 稱為行向量 列向量 矩陣只有一列時(shí) 稱為列向量 標(biāo)量 矩陣只有1個(gè)元素時(shí) 稱為標(biāo)量 MATLAB軟件 續(xù)6 2 變量定義var expt expt可以是數(shù)字 變量 操作符和函數(shù)等組成 例如 X 7A 123 456 789 B 11121314 C 21 22 23 24 D X B注意 用方括號(hào) 確定的范圍內(nèi)逐行給定元素 相同行元素由空格 或逗號(hào) 分隔 不同行數(shù)據(jù)由分號(hào) 或回車鍵分隔 語句結(jié)尾沒有分號(hào)的命令在屏幕顯示結(jié)果 有分號(hào)的命令只執(zhí)行而不顯示 MATLAB軟件 續(xù)7 5 MATLAB常用命令 1 MATLAB常用的基本數(shù)學(xué)函數(shù)abs x 純量的絕對(duì)值或向量的長度sqrt x 開平方round x 四舍五入至最近整數(shù)fix x 無論正負(fù) 舍去小數(shù)至最近整數(shù)sign x 符號(hào)函數(shù) Signumfunction exp x 自然指數(shù)pow2 x 2的指數(shù)log x 以e為底的對(duì)數(shù)log10 x 以10為底的對(duì)數(shù) 舉例 4x2 58x 89 0 a 4 b 58 c 89 MATLAB軟件 續(xù)8 2 適用于向量的常用函數(shù)min x 向量x的元素的最小值max x 向量x的元素的最大值mean x 向量x的元素的平均值median x 向量x的元素的中位數(shù)std x 向量x的元素的標(biāo)準(zhǔn)差diff x 向量x的相鄰元素的差sort x 對(duì)向量x的元素進(jìn)行排序 Sorting length x 向量x的元素個(gè)數(shù)sum x 向量x的元素總和 舉例 對(duì)130名正常成年男子紅細(xì)胞數(shù) RBC 萬 mm 的數(shù)據(jù)進(jìn)行描述性統(tǒng)計(jì) x zeros 130 1 平均 中位數(shù) 眾數(shù) 標(biāo)準(zhǔn)差 MATLAB軟件 續(xù)9 適用于向量的常用函數(shù) 續(xù)prod x 向量x的元素總乘積cumsum x 向量x的累計(jì)元素總和cumprod x 向量x的累計(jì)元素總乘積dot x y 向量x和y的內(nèi)積cross x y 向量x和y的外積 MATLAB軟件 續(xù)10 一些常用的操作函數(shù)format語句是控制輸出的顯示格式 句法如下 format設(shè)置默認(rèn)顯示格式5位的定點(diǎn)數(shù)formattype設(shè)置按指定的type格式輸出顯示 例如 formatratRatioofsmallintegersformat type 與上面語句一樣 但以函數(shù)形式出現(xiàn)disp x 顯示矩陣函數(shù) 注釋語句zeros ones預(yù)先定義矩陣的大小 例如 h zeros 6 生成6 6的0矩陣 MATLAB軟件 續(xù)11 3 重復(fù)命令1 for循環(huán) for loop 其基本形式為 for變量 矩陣 運(yùn)算式 end變量值會(huì)被依次設(shè)定為矩陣的每一行 來執(zhí)行for和end之間的運(yùn)算式 因此 若中斷情況 運(yùn)算式執(zhí)行的次數(shù)等于矩陣的行數(shù) 其中 是終結(jié)內(nèi)語句的可抑制的重復(fù)顯示 舉例 產(chǎn)生一個(gè)長度為6的調(diào)和數(shù)列x zeros 1 6 x是一個(gè)1 6的矩陣fori 1 6 x i 1 i end MATLAB軟件 續(xù)12 2 while循環(huán) 其基本形式為 while條件式 運(yùn)算式 End實(shí)現(xiàn)功能 只要條件成立 就會(huì)一再執(zhí)行運(yùn)算式 舉例 用while循環(huán)產(chǎn)生一個(gè)長度為6的調(diào)和數(shù)列x zeros 1 6 x是一個(gè)1 6的零矩陣i 1 whilei 6 x i 1 i i i 1 endformatshort MATLAB軟件 續(xù)13 4 邏輯命令if end 其基本形式為 if條件式 運(yùn)算式 end實(shí)現(xiàn)功能 條件成立 執(zhí)行運(yùn)算式 舉例 產(chǎn)生一個(gè)隨機(jī)數(shù)ifrand 1 1 0 5disp Givenrandomnumberisgreaterthan0 5 elsedisp Givenrandomnumberislessthan0 5 end MATLAB軟件 續(xù)14 5 將多個(gè)命令寫入一個(gè)M文件若要一次執(zhí)行大量的MATLAB命令 可將這些命令存于一個(gè)文件類型為m的文件中 并在MATLAB提示號(hào)下鍵入此文件名即可執(zhí)行 包含MATLAB命令的文件以m為擴(kuò)展名 通稱M文件 M files 舉例 名為test m的M文件有以下內(nèi)容 ThisismyfirsttestM file RogerJang March3 1997fprintf Startoftest m n fori 1 3 fprintf i d i 3 d n i i 3 endfprintf Endoftest m n MATLAB軟件 續(xù)15 建立test m的M文件進(jìn)入MATLAB編輯器File New M File在編輯器內(nèi)輸入以下內(nèi)容 ThisismyfirsttestM file RogerJang March3 1997fprintf Startoftest m n fori 1 3 fprintf i d i 3 d n i i 3 endfprintf Endoftest m n 保存文件到指定或默認(rèn)的目錄下File Save 輸入文件名 test m MATLAB軟件 續(xù)16 執(zhí)行文件 test m 進(jìn)入文件所在目錄或查看文件 CD dir執(zhí)行test m testStartoftest m i 1 i 3 1i 2 i 3 8i 3 i 3 27Endoftest m MATLAB軟件 續(xù)17 注意 M檔案可分為命令集 Scripts 及函數(shù) Functions 前述的test m即為命令集 其效用和將命令逐一輸入完全一樣 因此在命令集直接使用工作空間的變量 則在命令集中設(shè)定的變量 可在工作空間中看得到 函數(shù)則需要用到輸入?yún)?shù) Inputarguments 和輸出參數(shù) Outputarguments 來傳遞資訊 例如 計(jì)算一個(gè)正整數(shù)的階乘 Factorial 可以寫一個(gè)如下的MATLAB函數(shù)并將之存檔于fact m functionoutput fact n FACTCalculatefactorialofagivenpositiveinteger output 1 fori 1 n output output i end其中fact是函數(shù)名 n是輸入?yún)?shù) output是輸出參數(shù) 而i則是此函數(shù)用到的暫時(shí)變量 MATLAB軟件 續(xù)18 要使用此函數(shù) 直接鍵入函數(shù)名及適當(dāng)輸入?yún)?shù)值即可 y fact 5 y 120注意 在執(zhí)行fact之前 你必須先進(jìn)入fact m所在的目錄 在執(zhí)行fact 5 時(shí) MATLAB會(huì)進(jìn)入一個(gè)下層的暫時(shí)工作空間 Temperaryworkspace 將變量n的值設(shè)定為5 然后進(jìn)行各項(xiàng)函數(shù)的內(nèi)部運(yùn)算 所有內(nèi)部運(yùn)算所產(chǎn)生的變量都存在此暫時(shí)工作空間中 運(yùn)算完畢后 MATLAB會(huì)將最后輸出參數(shù)output的值設(shè)定給上層的變量y 并將清除此暫時(shí)工作空間及其所含的所有變量 MATLAB軟件 續(xù)19 6 MATLAB的文件 目錄等操作指令顯示當(dāng)前目錄pwd顯示當(dāng)前工作目錄s pwd將當(dāng)前工作目錄賦予變量s改變工作目錄cd顯示當(dāng)前工作目錄w cd將當(dāng)前工作目錄賦予wcd directory 設(shè)置當(dāng)前工作目錄為directory指定的目錄cd 將當(dāng)前工作目錄改為它的上級(jí)目錄cddirectoryorcd 與上面兩條指令等同 MATLAB軟件 續(xù)20 MATLAB的文件 目錄等操作指令 續(xù)顯示目錄列表dir列示當(dāng)前工作目錄的文件dirname列示指定的文件files dir directory 將指定目錄的文件列表信息返回到files中 例如 afiles dir d work afiles 3 name afiles 3 date顯示文件內(nèi)容type filename typefilenamefilename是要被顯示的文件的路徑和名字 MATLAB軟件 續(xù)21 MATLAB的文件 目錄等操作指令 續(xù)瀏覽或改變MATLAB目錄的搜尋路徑path顯示當(dāng)前MATLAB搜索路徑 初始搜索路徑由toolbox local pathdef m文件定義的 path newpath 將搜索路徑改變?yōu)橐粋€(gè)新路徑 這里新路徑是一個(gè)目錄的字符串?dāng)?shù)組 path path newpath 將一個(gè)新路徑增加到當(dāng)前的搜索路徑后 path newpath path 將一個(gè)新路徑加到當(dāng)前的搜索路徑前 p path 將指定的路徑返回到一個(gè)字符串變量 MATLAB軟件 續(xù)22 瀏覽或改變MATLAB目錄的搜尋路徑 續(xù)舉例 pathMATLABPATHd matlab5 toolbox matlab generald matlab5 toolbox matlab ops d matlab5 toolbox simulink deed matlab5 toolbox local查詢某一命令是在搜尋路徑的何處 用which命令 whichtestc data mlbook test m將c data mlbook加入MATLAB的搜尋路徑 用path path path c data mlbook MATLAB軟件 續(xù)23 7 資料的儲(chǔ)存與載入 續(xù)1 將計(jì)算所得的儲(chǔ)存倒文件的命令savesave將工作空間的所有變量儲(chǔ)存到名為matlab mat的二進(jìn)制文件 savefilename將工作空間的所有變量儲(chǔ)存到名為filename mat的二進(jìn)制文件 savefilenamexyz將變量x y z儲(chǔ)存到名為filename mat的二進(jìn)制文件 savefilenamex ascii將變量x以八位數(shù)存到名為filename的ASCII文件 savefilenamex ascii double將變量x以十六位數(shù)存到名為filename的ASCII文件 注意 命令后不加任何選項(xiàng)時(shí) save會(huì)將變量以二進(jìn)制 Binary 的方式儲(chǔ)存至擴(kuò)展名為mat的文件 另一個(gè)選項(xiàng)是 tab 可將同一列相鄰的數(shù)目以定位鍵 tab 隔開 MATLAB軟件 續(xù)24 將計(jì)算所得的儲(chǔ)存倒文件的命令save 續(xù)舉例 who列出工作空間的變量Yourvariablesare Bhjyansixz savetestBy將變數(shù)B與y儲(chǔ)存至test mat dir列出現(xiàn)在目錄中的文件 2plotxy docfact msimulink doctest m 1basic doc 3plotxyz docfirst doctemp doctest mat1basic docbook dotgo mtemplate doctestfile dat deletetest mat刪除test mat MATLAB軟件 續(xù)25 2 將文件載入儲(chǔ)存得變量命令loadloadfilename尋找名稱為filename mat的文件 并以二進(jìn)制格式載入 若找不到filename mat 則尋找名稱為filename的文件 并以ASCII格式載入 loadfilename ascii尋找名稱為filename的文件 并以ASCII格式載入 若以ASCII格式載入 則變量名稱即為文件名稱 但不包含擴(kuò)展名 若以二進(jìn)制載入 則可保留原有的變量名稱 舉例 clearall 清除工作空間中的變數(shù) x 1 10 savetestfile datx ascii 將x以ASCII格式存至testfile dat文件 loadtestfile dat 載入testfile dat who 列出工作空間中的變量Yourvariablesare testfilex注意 上述過程由于是以ASCII格式儲(chǔ)存與載入 所以產(chǎn)生了一個(gè)與文件名稱相同的變量testfile 此變量的值和原變量x完全相同 MATLAB軟件 續(xù)26 8 基本xy平面繪圖命令1 length獲得一組數(shù)的個(gè)數(shù)或一個(gè)向量的長度 x ones 1 8 n length x 2 linspace產(chǎn)生線性間隔向量y linspace a b 在 a b 之間產(chǎn)生一組100點(diǎn)的線性間隔的行向量 generatesarowvectoryof100pointslinearlyspacedbetweenandincludingaandb y linspace a b n 在 a b 之間產(chǎn)生一組n點(diǎn)的線性間隔的行向量 generatesarowvectoryofnpointslinearlyspacedbetweenandincludingaandb 舉例 t linspace 0 005 0 005 length y1 length y1 MATLAB軟件 續(xù)27 基本xy平面繪圖命令 續(xù)3 plot繪制一維曲線的基本函數(shù) 使用此函數(shù)之前 需先定義曲線上每一點(diǎn)的x及y座標(biāo) 舉例 畫出一條正弦曲線 closeall x linspace 0 2 pi 100 1產(chǎn)生00個(gè)點(diǎn)的x座標(biāo) y sin x 生成對(duì)應(yīng)的y座標(biāo) plot x y 關(guān)于MATLAB基本繪圖函數(shù)的應(yīng)用說明 繪圖刻度plot x軸和y軸均為線性刻度 Linearscale loglog x軸和y軸均為對(duì)數(shù)刻度 Logarithmicscale semilogx x軸為對(duì)數(shù)刻度 y軸為線性刻度semilogy x軸為線性刻度 y軸為對(duì)數(shù)刻度 MATLAB軟件 續(xù)28 關(guān)于MATLAB基本繪圖函數(shù)的應(yīng)用說明 續(xù)畫出多條曲線 只需將座標(biāo)對(duì)依次放入plot函數(shù)plot x sin x x cos x 若要改變顏色 在座標(biāo)對(duì)後面加上相關(guān)字串plot x sin x c x cos x g 若要同時(shí)改變顏色及圖線型態(tài) Linestyle 也是在座標(biāo)對(duì)後面加上相關(guān)字串plot x sin x co x cos x g MATLAB軟件 續(xù)29 關(guān)于MATLAB基本繪圖函數(shù)的應(yīng)用說明 續(xù)plot繪圖函數(shù)的叁數(shù)說明在用plot a b s 繪圖時(shí) 可以使用各類線型 圖符 顏色 這里s是由下表中的任一元素或3列中所有元素組合的字符串 MATLAB軟件 續(xù)30 4 將一圖形窗口劃分為多個(gè)顯示區(qū) 使之同時(shí)畫出數(shù)個(gè)小圖形於同一個(gè)視窗之中 subplot m n p 實(shí)現(xiàn)功能 將圖形窗口劃分為m n矩陣的顯示區(qū) 同時(shí)激活第p個(gè)顯示區(qū) 舉例 subplot 2 2 1 plot x sin x subplot 2 2 2 plot x cos x subplot 2 2 3 plot x sinh x subplot 2 2 4 plot x cosh x MATLAB軟件 續(xù)31 5 其他各種二維繪圖函數(shù)bar長條圖errorbar圖形加上誤差范圍fplot較精確的函數(shù)圖形polar極座標(biāo)圖hist累計(jì)圖rose極座標(biāo)累計(jì)圖stairs階梯圖stem針狀圖fill實(shí)心圖 二 圖形信號(hào)的處理程序介紹 1 峰值程序fori 1 length u ifu i max u m i breakendendp h h m 處理程序 續(xù)1 求起點(diǎn)和終點(diǎn)的程序 judge startfori 5 length y2 7s 5 ify2 7s i 4 0 y2 7s i 1 0 y2 7s i 2 0 y2 7s i 3 0 y2 7s i 4 0y2 7s i 5 0k ibreakendend judge endfori 1 length y2 7s 6 ify2 7s i 0 y2 7s i 1 0 y2 7s i 2 0 y2 7s i 3 0 y2 7s i 4 0 y2 7s i 5 0 y2 7s i 6 0j ibreakendend 處理程序 續(xù)2 2 時(shí)距程序 judge start judge end u y2 7s k j 宮縮單波數(shù)據(jù) h y1 k j 對(duì)應(yīng)單波宮縮的心率數(shù)據(jù) t0 length u 1 0 005 處理程序 續(xù)3 3 面積程序 t linspace 0 005 0 005 length h length h subplot 2 1 1 plot t h 心率 時(shí)間圖 subplot 2 1 2 plot t u 宮縮 時(shí)間圖 area 0 005 sum u 15 1024 計(jì)算面積 三 臨床電生理數(shù)據(jù)分析實(shí)驗(yàn)?zāi)M軟件介紹 1 實(shí)驗(yàn)?zāi)M軟件簡介臨床電生理數(shù)據(jù)分析實(shí)驗(yàn)?zāi)M軟件Med exe是用于對(duì)臨床功能類的信號(hào)進(jìn)行初步分析的軟件 該軟件是用VisualBasic程序語言來編寫的 目前具有如下功能 1 可以瀏覽醫(yī)學(xué)儀器采集的電生理信號(hào)數(shù)據(jù) 人工干預(yù)信號(hào)截取 并以二進(jìn)制格式和文本格式來保存截取的信號(hào)數(shù)據(jù) 2 可以利用面板按扭工具 對(duì)胎心率信號(hào) 孕婦的宮縮信號(hào)進(jìn)行人工干預(yù)的信號(hào)處理和分析 包括 奇異點(diǎn)處理 平均心率 瞬時(shí)心率 宮縮峰值 宮縮時(shí)距和宮縮面積的計(jì)算 3 可以進(jìn)行信號(hào)數(shù)字平滑濾波 基線調(diào)整 自動(dòng)計(jì)算部分時(shí)域參數(shù) 2 軟件的使用說明 實(shí)驗(yàn)?zāi)M軟件介紹 續(xù)1 1 原始采樣數(shù)據(jù)回顯 臨床信號(hào)截取 選取數(shù)據(jù)源位置 確定位置后 單擊 OK 在輸入病歷號(hào)的輸入框內(nèi)輸入 A0012 在選擇文件框內(nèi)選擇項(xiàng)目 R0012A 5 單擊 開始 單擊 回放 2 進(jìn)行人工干預(yù)操作 在需要截波或人工分析時(shí) 單擊 凍結(jié) 截取波形 計(jì)算時(shí)距 面積和平均心率 或去除奇異項(xiàng)時(shí) 要用鼠標(biāo)在視圖瀏覽區(qū)內(nèi)拖出選區(qū) 然后分別單擊相應(yīng)的操作按鈕注意 如要取消選區(qū) 可以單擊瀏覽視圖的任意位置即可 計(jì)算峰值或瞬時(shí)心率時(shí) 按著 Alt 鍵 然后單擊曲線的選取點(diǎn) 接著單擊相應(yīng)的操作按鈕 3 結(jié)束軟件使用單擊窗口的關(guān)閉按鈕 四 實(shí)驗(yàn) 1 均值計(jì)算2 時(shí)距計(jì)算3 峰值計(jì)算4 面積計(jì)算5 作圖 平均值 平均幅值 實(shí)驗(yàn) 續(xù)1 1 MatLab實(shí)驗(yàn)內(nèi)容的介紹 1 利用matlab計(jì)算宮縮的峰值 宮縮時(shí)的平均心率 以及宮縮峰值對(duì)應(yīng)的順時(shí)心率 計(jì)算宮縮波的面積 宮縮持續(xù)時(shí)間實(shí)驗(yàn)步驟 啟動(dòng)matlab 增加路徑 path path D matlab dat 導(dǎo)入數(shù)據(jù) clearall 清除變量心率 y1 load r0012A0h txt 宮縮 y2 load r0012a0u txt 曲線作圖建立XY坐標(biāo)軸與繪圖時(shí)間數(shù)組點(diǎn) t linspace 0 005 0 005 length y1 length y1 繪制心率與宮縮圖 subplot 2 1 1 plot t y1 心率 時(shí)間圖 subplot 2 1 2 plot t y2 宮縮 時(shí)間圖 實(shí)驗(yàn) 續(xù)2 宮縮局部加權(quán)平滑濾波 七點(diǎn)平滑 n length y2 fori 4 n 3y2 7s i 2 y2 i 3 3 y2 i 2 6 y2 i 1 7 y2 i 2 y2 i 3 3 y2 i 2 6 y2 i 1 21 end y2 7s 123n 2n 1n y2 1 y2 2 y2 3 y2 n 2 y2 n 1 y2 n 再次繪制心率與宮縮圖 觀察曲線變化的情況 subplot 2 1 1 plot t y1 心率 時(shí)間圖 subplot 2 1 2 plot t y2 7s 宮縮 時(shí)間圖 判斷宮縮波的起始 終止點(diǎn) judge start judge end 運(yùn)行文件judge start mj和judge end mk 取宮縮單波數(shù)據(jù)以及對(duì)應(yīng)單波宮縮的心率數(shù)據(jù) u y2 7s k j h y1 k j 實(shí)驗(yàn) 續(xù)3 計(jì)算宮縮峰值 Pu max u 建立時(shí)間坐標(biāo) t linspace 0 005 0 005 length h length h subplot 2 1 1 plot t h 心率 時(shí)間圖 subplot 2 1 2 plot t u 宮縮 時(shí)間圖 計(jì)算宮縮面積 area 0 005 sum u 15 1024 計(jì)算對(duì)應(yīng)的平均心率 mean p sum h length h 判斷宮縮最強(qiáng)時(shí)對(duì)應(yīng)的瞬時(shí)心率 p h 宮縮最強(qiáng)時(shí)的瞬時(shí)心率用p h表示 運(yùn)行文件p h m計(jì)算持續(xù)時(shí)間 t0 length u 1 0 005 實(shí)驗(yàn) 續(xù)4 2 單擊MATLAB程序窗口左邊窗格的Workspace 工作空間 卡片 瀏覽該卡片窗格內(nèi)的所有變量 并執(zhí)行以下命令 clearall然后再瀏覽該卡片窗格內(nèi)的所有變量 你有何結(jié)論 3 建立M 文件 輸入以下內(nèi)容 保存為 original dat m 文件 并執(zhí)行該M 文件 觀察其與前面單步執(zhí)行的結(jié)果是否一致 給出你的結(jié)論 path path D matlab dat y1 load r0012A0h txt y2 load r0012a0u txt t linspace 0 005 0 005 length y1 length y1 subplot 2 1 1 plot t y1 subplot 2 1 2 plot t y2 實(shí)驗(yàn) 續(xù)5 2 Med exe實(shí)驗(yàn)內(nèi)容的簡要介紹 1 啟動(dòng)臨床電生理數(shù)據(jù)分析實(shí)驗(yàn)?zāi)M軟件 2 在輸入病歷號(hào)的輸入框內(nèi)輸入 A0012 3 在選擇文件框內(nèi)選擇項(xiàng)目 R0012A 5 4 對(duì)回放的第一個(gè)宮縮波曲線對(duì)應(yīng)的心率曲線去除奇異點(diǎn) 然后通過菜單命令 保存截取波形數(shù)據(jù) 將數(shù)據(jù)保存為txt類型文件 記下你在保存文件提示輸入的序號(hào)n 保存后你的文件的名字應(yīng)為R0012Anh txt和R0012Anh txt 5 再對(duì)上述選定區(qū)域的數(shù)據(jù)求取平均心率 時(shí)距 面積 然后人工選取峰值點(diǎn) 計(jì)算峰值和瞬時(shí)心率 6 嘗試進(jìn)入MATLAB 對(duì)上面截取的數(shù)據(jù)文件重復(fù)上面實(shí)驗(yàn)一的操作 比較其結(jié)果