03-醫(yī)學圖形圖像處理識別技術與醫(yī)學信息分析

1、醫(yī)學圖形圖像處理識別技術與醫(yī)學信息分析劉燕E-Mail: Phone No. 87331856-811目錄內容容序醫(yī)醫(yī)學信息息學簡述述第一講醫(yī)醫(yī)學信信號處理理概述第二講醫(yī)醫(yī)學圖圖形信號號的采集集、存儲儲與處理理第三講Mathlab 與圖圖形信號號的處理理識別第四講醫(yī)醫(yī)學圖圖像信號號的采集集、存儲儲格式、處理方方法第五講Mathlab、Vtk、Mtk與醫(yī)醫(yī)學圖像像的處理理第六講醫(yī)醫(yī)學信信息挖掘掘第七講用用于教教學的醫(yī)醫(yī)學圖片片處理技技術作業(yè)第三講Mathlab與與圖形信信號的處處理識別別MATLAB軟軟件的使使用介紹紹圖形信號號的處理理與識別別程序介介紹臨床電生生理數(shù)據(jù)據(jù)分析實實驗模擬擬軟件介介

2、紹實驗 一、MATLAB軟件件的使用用介紹1.MATLAB的的簡介(1)MATLAB高級級交互式式軟件包包(2)MATLAB典型應應用數(shù)學計算算算法推導導建模和仿仿真模擬擬數(shù)據(jù)分析析和結果果的可視視化工程圖形形繪制應應用程序序開發(fā)（包括用用戶圖形形界面的的建立）(3)MATLAB產(chǎn)生的的歷史背背景調用EISPACK和和LINPACK的FORTRAN子程序序庫EISPACK和LINPACK的的接口程程序MATLAB發(fā)展展成為適適合多學學科，多多種工作作平臺的的功能強強大的大大型軟件件。美國MathWorks公司司設計的的“科學學和工程程計算的的高級交交互式軟軟件包”CleveMoler、NewM

3、exico、JohnLittleMATLAB軟軟件_續(xù)1(3)MATLAB的的語言特特點語言簡潔潔緊湊，使用方方便靈活活，庫函函數(shù)極其其豐富。運算符豐豐富。MATLAB既既具有結結構化的的控制語語句，又又有面向向對象編編程的特特性。程序限制制不嚴格格，程序序設計自自由度大大。程序的可可移植性性很好，基本上上不做修修改就可可以在各各種型號號的計算算機和操操作系統(tǒng)統(tǒng)上運行行。MATLAB的的圖形功功能強大大。MATLAB的的缺點是是程序的的執(zhí)行速速度較慢慢。功能強大大的工具具箱。源程序的的開放性性。如：解方方程組Ax=bx=AbMATLAB軟軟件_續(xù)22.MATLA的安安裝進入MATLAB軟件件安

4、裝程程序目錄錄，然后后單擊“Setup”程序，在安裝裝程序的的引導下下，逐步步完成以以下操作作:(1)選選擇“Install”安裝裝或”Updatelicense”升級級許可證證，單擊擊Next；然后后輸入姓姓名、公公司與序序列號，單擊Next。(2)選選擇接受受許可證證協(xié)議，單擊“Yes”單單擊Next，然然后選擇擇“Typical”典型安安裝，單單擊Next。(3)選選擇或輸輸入安裝裝目的路路徑，單單擊Next，后后單擊“Yes”(4)單單擊Install，系系統(tǒng)開始始安裝。MATLAB軟軟件_續(xù)3(5)在在提示插插入第2張CD時，插插入CD2，然然后單擊擊OK，接接著在文文件擴展展名警示

5、示對話框框內單擊擊Yestoall。(6)單單擊Next，接接著單擊擊Finish。(7)如果在單單擊Finish按鈕之之前，選選擇了“Start”選項，則此時時進入MATLAB工工作環(huán)境境。MATLAB軟軟件_續(xù)43.MATLAB啟啟動與退退出啟動MATLAB，進進入MATLAB工作作環(huán)境啟動MATLABMATLAB軟軟件_續(xù)5結束MATLAB，退退出MATLAB工作作環(huán)境鍵入exit鍵入quit直接關閉閉MATLAB的命令令視窗退出MATLABMATLAB軟軟件_續(xù)6-14.MATLAB變變量與基基本輸入入輸出變量命名名的規(guī)則則：第一一個字母母必須是是英文字字母；字母間間不可留留空格；最多

6、多只能有有19個個字母，MATLAB會忽略略多余字字母。變量定義義與輸入入輸出向量與標標量矩陣行和列列組成的的數(shù)據(jù)表表行向量矩陣陣只有一一行時，稱為行行向量。列向量矩陣陣只有一一列時，稱為列列向量。標量矩陣只只有1個個元素時時，稱為為標量。MATLAB軟軟件_續(xù)6-2變量定義義var=expt;expt可以是是數(shù)字、變量、操作符符和函數(shù)數(shù)等組成成。例如：X=7A=1 23;45 6; 789B=11121314C=21;22;23;24D=X*B注意：用方括號號“”確確定的的范圍內內逐行給給定元素素，相同同行元素素由空格格“”或逗號號“,”分隔，不同行行數(shù)據(jù)由由分號“;”或或回車鍵鍵分隔。語句

7、結尾尾沒有分分號的命命令在屏屏幕顯示示結果；有分號號的命令令只執(zhí)行行而不顯顯示。MATLAB軟軟件_續(xù)75.MATLAB常常用命令令(1)MATLAB常用的的基本數(shù)數(shù)學函數(shù)數(shù)abs(x)：純量的的絕對值值或向量量的長度度sqrt(x)：開平平方round(x)：四四舍五入入至最近近整數(shù)fix(x)：無論正正負，舍舍去小數(shù)數(shù)至最近近整數(shù)sign(x)：符號號函數(shù)(Signumfunction)exp(x)：自然指指數(shù)pow2(x)：2的的指數(shù)log(x)：以e為為底的對對數(shù)log10(x)：以以10為為底的對對數(shù)舉例:4x2+58x-89=0（a=4,b=58,c=89）MATLAB軟軟件_續(xù)8

8、(2)適適用于于向量的的常用函函數(shù)min(x): 向量量x的元元素的最最小值max(x): 向量量x的元元素的最最大值mean(x):向向量x的的元素的的平均值值median(x): 向量量x的元元素的中中位數(shù)std(x): 向量量x的元元素的標標準差diff(x):向向量x的的相鄰元元素的差差sort(x):對對向量x的元素素進行排排序（Sorting）length(x): 向量量x的元元素個數(shù)數(shù)sum(x): 向量量x的元元素總和和舉例:對對130名正常常成年男男子紅細細胞數(shù)(RBC ,萬萬/mm)的數(shù)數(shù)據(jù)進行行描述性性統(tǒng)計。(x=zeros(130,1)平均、中中位數(shù)、眾數(shù)、標準差差MA

9、TLAB軟軟件_續(xù)9適用于向向量的常常用函數(shù)數(shù)_續(xù)prod(x):向向量x的的元素總總乘積cumsum(x): 向量量x的累累計元素素總和cumprod(x):向向量x的的累計元元素總乘乘積dot(x,y): 向量量x和y的內積積cross(x,y):向向量x和y的的外積MATLAB軟軟件_續(xù)10一些常用用的操作作函數(shù)format語語句是控控制輸出出的顯示示格式，句法如如下：format設設置默默認顯示示格式5位的定定點數(shù)formattype設設置置按指定定的type格格式輸出出顯示。例如：formatratRatio of small integersformat(type)與與上上面語句句

10、一樣，但以函函數(shù)形式式出現(xiàn)disp(x)顯顯示矩陣陣函數(shù)%注注釋語語句zeros、ones 預先先定義矩矩陣的大大小。例如：h =zeros(6);%生成66的0矩陣MATLAB軟軟件_續(xù)11(3)重重復命命令1）for循環(huán)環(huán)（for-loop），其其基本形形式為：for變變量=矩矩陣；運算式；end變量值會會被依次次設定為為矩陣的的每一行行，來執(zhí)執(zhí)行for和end之之間的運運算式。因此，若中斷斷情況，運算式式執(zhí)行的的次數(shù)等等于矩陣陣的行數(shù)數(shù)。其中中“;”是終結結內語句句的可抑抑制的重重復顯示示。舉例:產(chǎn)產(chǎn)生一個個長度為為6的調調和數(shù)列列x =zeros(1,6);% x是是一個16的的矩陣f

11、ori= 1:6,x(i) =1/i;endMATLAB軟軟件_續(xù)122）while循環(huán)，其基本本形式為為：while條條件式；運算式；End實現(xiàn)功能能：只要要條件成成立，就就會一再再執(zhí)行運運算式。舉例:用用while循循環(huán)產(chǎn)生生一個長長度為6的調和和數(shù)列x =zeros(1,6);%x是一個個16的零矩矩陣i =1;whilei 0.5disp(Givenrandomnumber is greaterthan0.5.);elsedisp(Givenrandomnumber is lessthan0.5.);endMATLAB軟軟件_續(xù)14(5)將多多個命令令寫入一一個M文文件若要一次次執(zhí)行大

12、大量的MATLAB命命令，可可將這些些命令存存于一個個文件類類型為m的文件件中，并并在MATLAB提示示號下鍵鍵入此文文件名即即可執(zhí)行行。包含含MATLAB命令的的文件以以m為擴擴展名，通稱M文件（M-files）。舉例:名名為test.m的M文件有有以下內內容：% ThisismyfirsttestM-file.% Roger Jang, March 3, 1997fprintf(Startoftest.m!n);fori= 1:3,fprintf(i =%d-i3=%dn,i,i3);endfprintf(Endoftest.m!n);MATLAB軟軟件_續(xù)15建立test.m的的M文件件

13、進入MATLAB編輯輯器FileNewM-File在編輯器器內輸入入以下內內容% ThisismyfirsttestM-file.% Roger Jang, March 3, 1997fprintf(Startoftest.m!n);fori= 1:3,fprintf(i =%d-i3=%dn,i,i3);endfprintf(Endoftest.m!n);保存文件件到指定定或默認認的目錄錄下FileSave輸入文文件名“test.m”MATLAB軟軟件_續(xù)16執(zhí)行文文件“test.m”進入文件件所在目目錄或查查看文件件CDdir執(zhí)行test.mtestStartoftest.m!i =1-

14、i3= 1i =2- i3= 8i =3- i3= 27Endoftest.m! MATLAB軟軟件_續(xù)17注意：M檔案案可分為為命令集(Scripts)及及函數(shù)(Functions)。前前述的test.m即即為命令令集，其其效用和和將命令令逐一輸輸入完全全一樣，因此在在命令集集直接使使用工作作空間的的變量，則在命命令集中中設定的的變量，可在工工作空間間中看得得到。函數(shù)則需需要用到到輸入?yún)?shù)(Inputarguments)和輸出參數(shù)數(shù)(Output arguments)來來傳遞資資訊。例如：計計算一個個正整數(shù)數(shù)的階乘乘（Factorial），可以寫寫一個如如下的MATLAB函函數(shù)并將將之存檔

15、檔于fact.m：functionoutput=fact(n)% FACTCalculate factorialofa given positive integer.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ù)數(shù)名及適適當輸入入?yún)?shù)值值即可。y=fact(5)y =120注意：在執(zhí)行fact之前，你必須須先進入入fact.m所在的的目錄。在執(zhí)行fact(5)時，MATLAB會會進入一一個下層層的暫

16、時時工作空空間（Temperaryworkspace），將將變量n的值設設定為5，然后后進行各各項函數(shù)數(shù)的內部部運算，所有內內部運算算所產(chǎn)生生的變量量都存在在此暫時時工作空空間中。運算完完畢后，MATLAB會將最最后輸出出參數(shù)output的的值設定定給上層層的變量量y，并并將清除除此暫時時工作空空間及其其所含的的所有變變量。MATLAB軟軟件_續(xù)19(6)MATLAB的的文件、目錄等等操作指指令顯示當前前目錄pwd顯顯示當前前工作目目錄s =pwd將將當前前工作目目錄賦予予變量s改變工作作目錄cd顯顯示當當前工作作目錄w =cd將將當當前工作作目錄賦賦予wcd(directory)設設置置當前

17、工工作目錄錄為directory指定定的目錄錄cd(.)將將當前前工作目目錄改為為它的上上級目錄錄cddirectoryorcd.與與上面兩兩條指令令等同MATLAB軟軟件_續(xù)20MATLAB的的文件、目錄等等操作指指令_續(xù)顯示目錄錄列表dir列列示示當前工工作目錄錄的文件件dirname列列示示指定的的文件files= dir(directory)將將指指定目錄錄的文件件列表信信息返回回到files中。例如：afiles=dir(d:/work/*.*)afiles(3).nameafiles(3).date顯示文件件內容type (filename)type filenamefilenam

18、e是要要被顯示示的文件件的路徑徑和名字字MATLAB軟軟件_續(xù)21MATLAB的的文件、目錄等等操作指指令_續(xù)瀏覽或改改變MATLAB目目錄的搜搜尋路徑徑path顯顯示當前前MATLAB 搜索索路徑。初始搜索索路徑由由toolbox/local/pathdef.m文文件定義義的。path(newpath)將將搜索路路徑改變變?yōu)橐粋€個新路徑徑，這里里新路徑徑是一個個目錄的的字符串串數(shù)組。path(path,newpath)將將一個個新路徑徑增加到到當前的的搜索路路徑后。path(newpath,path)將將一個個新路徑徑加到當當前的搜搜索路徑徑前。p =path(.)將將指定定的路徑徑返回到到

19、一個字字符串變變量。MATLAB軟軟件_續(xù)22瀏覽或改改變MATLAB目目錄的搜搜尋路徑徑_續(xù)舉例:pathMATLABPATHd:matlab5toolboxmatlabgenerald:matlab5toolboxmatlabopsd:matlab5toolboxsimulinkdeed:matlab5toolboxlocal查詢某一一命令是是在搜尋尋路徑的的何處，用which命令：which testc:datamlbooktest.m將c:datamlbook加加入MATLAB的搜搜尋路徑徑，用pathpath(path,c:datamlbook)MATLAB軟軟件_續(xù)23(7)資資料

20、的儲儲存與載載入_續(xù)1)將將計算所所得的儲儲存倒文文件的命命令savesave將將工作作空間的的所有變變量儲存存到名為為matlab.mat的二二進制文文件。save filename將將工工作空間間的所有有變量儲儲存到名名為filename.mat的二二進制文文件。save filename xyz將將變量量x、y、z儲儲存到名名為filename.mat的二二進制文文件。save filename xascii將將變變量x以以八位數(shù)數(shù)存到名名為filename的ASCII文件。save filename x-ascii double將將變量量x以十十六位數(shù)數(shù)存到名名為filename的AS

21、CII文件。注意:命令后后不加任任何選項項時，save會將變變量以二二進制（Binary）的方方式儲存存至擴展展名為mat的的文件；另另一個選選項是-tab，可將將同一列列相鄰的的數(shù)目以以定位鍵鍵（tab）隔隔開。MATLAB軟軟件_續(xù)24將計算所所得的儲儲存倒文文件的命命令save_續(xù)舉例：who列列出工工作空間間的變量量Your variablesare:B hjyansix zsavetestB y將將變變數(shù)B與與y儲存存至test.matdir列列出現(xiàn)在在目錄中中的文件件. 2plotxy.doc fact.msimulink.doctest.m $1basic.doc.3plotxy

22、z.docfirst.doc temp.doc test.mat1basic.doc book.dot go.mtemplate.doc testfile.datdeletetest.mat刪刪除除test.matMATLAB軟軟件_續(xù)252)將文文件載入入儲存得得變量命命令loadload filename尋尋找名稱稱為filename.mat的文文件，并并以二進進制格式式載入。若找不不到filename.mat，則則尋找名名稱為filename的文文件，并并以ASCII格式載載入。load filename ascii尋尋找名稱稱為filename的文件件，并以以ASCII格格式載入入。若

23、以以ASCII格格式載入入，則變變量名稱稱即為文文件名稱稱（但不不包含擴擴展名）。若以以二進制制載入，則可保保留原有有的變量量名稱。舉例：clearall;%清清除工作作空間中中的變數(shù)數(shù)x =1:10;savetestfile.dat x-ascii%將將x以ASCII格式式存至testfile.dat文文件loadtestfile.dat% 載入入testfile.datwho%列列出工作作空間中中的變量量Your variablesare:testfilex注意：上上述過程程由于是是以ASCII格式儲儲存與載載入，所所以產(chǎn)生生了一個個與文件件名稱相相同的變變量testfile，此變變量的值

24、值和原變變量x完完全相同同。MATLAB軟軟件_續(xù)26(8)基基本xy平面繪繪圖命令令1)length獲獲得得一組數(shù)數(shù)的個數(shù)數(shù)或一個個向量的的長度x=ones(1,8);n=length(x)2)linspace產(chǎn)產(chǎn)生線線性間隔隔向量y =linspace(a,b)在在a,b之之間產(chǎn)生生一組100點點的線性性間隔的的行向量量。generatesa rowvectoryof100points linearly spacedbetween andincludingaandb.y =linspace(a,b,n)在在a,b之之間產(chǎn)生生一組n點的線線性間隔隔的行向向量。generatesarowvec

25、tor yofnpointslinearlyspaced betweenand includinga andb.舉例：t=linspace(0.005,0.005*length(y1), length(y1);MATLAB軟軟件_續(xù)27基本xy平面繪繪圖命令令_續(xù)3）plot繪繪制制一維曲曲線的基基本函數(shù)數(shù)，使用用此函數(shù)數(shù)之前，需先定定義曲線線上每一一點的x及y座座標。舉例:畫畫出一條條正弦曲曲線close all;x=linspace(0,2*pi,100);% 1產(chǎn)產(chǎn)生00個點的的x座標標y=sin(x);%生成成對應的的y座標標plot(x,y);關于MATLAB基本本繪圖函函數(shù)的應應用

26、說明明 :繪圖刻度度plot:x軸和y軸均為為線性刻刻度（Linearscale）loglog: x軸軸和y軸軸均為對對數(shù)刻度度（Logarithmic scale）semilogx:x軸為為對數(shù)刻刻度，y軸為線線性刻度度semilogy:x軸為為線性刻刻度，y軸為對對數(shù)刻度度MATLAB軟軟件_續(xù)28關于MATLAB基本本繪圖函函數(shù)的應應用說明明_續(xù)畫出多條條曲線，只需將將座標對對依次放放入plot函函數(shù)plot(x, sin(x),x,cos(x);若要改變變顏色，在座標標對後面面加上相相關字串串plot(x, sin(x),c,x,cos(x), g);若要同時時改變顏顏色及圖圖線型態(tài)態(tài)

27、（Linestyle），也是是在座標標對後面面加上相相關字串串plot(x, sin(x),co,x,cos(x),g*);MATLAB軟軟件_續(xù)29關于MATLAB基本本繪圖函函數(shù)的應應用說明明_續(xù)plot繪圖函函數(shù)的叁叁數(shù)說明明在用plot(a,b,s)繪圖時時，可以以使用各各類線型型、圖符符、顏色色，這里里s是由由下表中中的任一一元素或或3列中中所有元元素組合合的字符符串。MATLAB軟軟件_續(xù)304)將將一圖形形窗口劃劃分為多多個顯示示區(qū)，使使之同時時畫出數(shù)數(shù)個小圖圖形於同同一個視視窗之中中。subplot(m,n,p)實現(xiàn)功能能：將圖圖形窗口口劃分為為mn矩陣的的顯示區(qū)區(qū)，同時時激活

28、第第p個顯顯示區(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ù)315)其其他各種種二維繪繪圖函數(shù)數(shù)bar長長條條圖errorbar圖圖形加加上誤差差范圍fplot較較精確確的函數(shù)數(shù)圖形polar極極座標標圖hist累累計圖rose極極座標累累計圖stairs階階梯梯圖stem針針狀圖fill實實心圖二、圖形形信號的的處理程程序介紹紹1.峰峰值程序序fori=1:length(u

29、)ifu(i)=max(u)m=i;breakendendp_h=h(m)處理程序序_續(xù)1求起點和和終點的的程序%judge_startfori=5:(length(y2_7s)-5)ify2_7s(i-4)=0&y2_7s(i-3)=0&y2_7s(i-2)=0&y2_7s(i-1)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)=

30、0&y2_7s(i+4)=0&y2_7s(i+5)=0&y2_7s(i+6)judge_startjudge_endu=y2_7s(k:j);%宮縮縮單波數(shù)數(shù)據(jù)h=y1(k:j);%對應應單波宮宮縮的心心率數(shù)據(jù)據(jù)t0=(length(u)-1)*0.005;處理程序序_續(xù)33.面面積程程序t=linspace(0.005,0.005*length(h),length(h);subplot(2,1,1),plot(t,h);%心率-時間圖圖;subplot(2,1,2),plot(t,u);%宮宮縮-時時間圖;area=0.005*sum(u)*15./1024%計計算面面積三、臨床床電生理理數(shù)

31、據(jù)分分析實驗驗模擬軟軟件介紹紹1.實驗驗模擬軟軟件簡介介臨床電生生理數(shù)據(jù)據(jù)分析實實驗模擬擬軟件Med.exe是用于于對臨床床功能類類的信號號進行初初步分析析的軟件件，該軟軟件是用用Visual Basic程序語語言來編編寫的。目前具具有如下下功能：（1）可可以瀏覽覽醫(yī)學儀儀器采集集的電生生理信號號數(shù)據(jù)，人工干干預信號號截取，并以二二進制格格式和文文本格式式來保存存截取的的信號數(shù)數(shù)據(jù)；（2）可可以利用用面板按按扭工具具，對胎胎心率信信號、孕孕婦的宮宮縮信號號進行人人工干預預的信號號處理和和分析，包括：奇異點點處理，平均心心率、瞬瞬時心率率、宮縮縮峰值、宮縮時時距和宮宮縮面積積的計算算；（3）可可

32、以進行行信號數(shù)數(shù)字平滑滑濾波、基線調調整、自自動計算算部分時時域參數(shù)數(shù)。2.軟軟件的使使用說明明實驗模擬擬軟件介介紹_續(xù)1（1）原原始采樣樣數(shù)據(jù)回回顯“臨床床信號截截取”“選取取數(shù)據(jù)源源位置”確定定位置后后，單擊擊OK在輸入入病歷號號的輸入入框內輸輸入：A0012在選擇擇文件框框內選擇擇項目：R0012A_5單擊開始單擊擊回放放（2）進進行人工工干預操操作在需要要截波或或人工分分析時，單擊凍結截取波波形，計計算時距距、面積積和平均均心率，或去除除奇異項項時，要要用鼠標標在視圖圖瀏覽區(qū)區(qū)內拖出出選區(qū)，然后分分別單擊擊相應的的操作按按鈕注意：如要取取消選區(qū)區(qū)，可以以單擊瀏瀏覽視圖圖的任意意位置即即

33、可。計算峰峰值或瞬瞬時心率率時，按按著Alt鍵，然然后單擊擊曲線的的選取點點，接著著單擊相相應的操操作按鈕鈕（3）結結束軟件件使用單擊窗口口的關閉閉按鈕。四、實驗驗1均值值計算2.時時距計算算3.峰峰值計算算4.面面積計算算5.作作圖平均值：平均幅值值：實驗_續(xù)11.MatLab實實驗內容容的介紹紹（1）利利用matlab計算算宮縮的的峰值、宮縮時時的平均均心率，以及宮宮縮峰值值對應的的順時心心率、計計算宮縮縮波的面面積、宮宮縮持續(xù)續(xù)時間實驗步驟驟啟動matlab；增加路路徑：path(path,D:matlab_dat)導入數(shù)數(shù)據(jù)clear all;%清除變變量心率：y1=load(r001

34、2A0h.txt);宮縮：y2=load(r0012a0u.txt);曲線作作圖建立XY坐標軸軸與繪圖圖時間數(shù)數(shù)組點t=linspace(0.005,0.005*length(y1),length(y1);繪制心率率與宮縮縮圖subplot(2,1,1),plot(t,y1);%心心率-時時間圖;subplot(2,1,2),plot(t,y2);%宮宮縮-時間圖圖；實驗_續(xù)2宮縮局部部加權平平滑濾波波(七點點平滑)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(12 3n-2 n-1n)=y2(1) y2(2) y2(3) y2(n-2)y2(n-1)y2(n)再次繪制制心率與與宮縮圖圖，觀察察曲線變變化的情情況subplot(2,1,1),plot(t,y1);%心心率-時時間圖;subplot(2,1,2),plot(t,y2_7s);%宮宮縮-時間