03-醫(yī)學(xué)圖形圖像處理識(shí)別技術(shù)與醫(yī)學(xué)信息分析-閱讀頁
2025-02-25 21:45本頁面
【正文】 展名）。 ? 舉例： clear all。 save x ascii % 將 x以 ASCII格式存至 load % 載入 who % 列出工作空間中的變量 Your variables are: testfile x ? 注意：上述過程由于是以 ASCII格式儲(chǔ)存與載入，所以產(chǎn)生了一個(gè)與文件名稱相同的變量 testfile，此變量的值和原變量 x完全相同。 n = length(x) ? 2) linspace 產(chǎn)生線性間隔向量 ? y = linspace(a,b) 在 [a,b]之間產(chǎn)生一組 100點(diǎn)的線性間隔的行向量。 generates a row vector y of n points linearly spaced between and including a and b. ? 舉例： t=linspace(,*length(y1), length(y1))。 ? 舉例 :畫出一條正弦曲線 close all。 % 1產(chǎn)生 00個(gè)點(diǎn)的 x座標(biāo) y=sin(x)。 關(guān)于 MATLAB基本繪圖函數(shù)的應(yīng)用說明 : ? 繪圖刻度 ? plot: x軸和 y軸均為線性刻度（ Linear scale） ? loglog: x軸和 y軸均為對(duì)數(shù)刻度（ Logarithmic scale） ? 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))。c39。g39。 ? 若要同時(shí)改變顏色及圖線型態(tài)（ Line style），也是在座標(biāo)對(duì)後面加上相關(guān)字串 ? plot(x, sin(x), 39。, x, cos(x), 39。)。 MATLAB軟件 _續(xù) 30 ? 4) 將一圖形窗口劃分為多個(gè)顯示區(qū)，使之同時(shí)畫出數(shù)個(gè)小圖形於同一個(gè)視窗之中。 ? 舉例： ? subplot(2,2,1)。 ? subplot(2,2,2)。 ? subplot(2,2,3)。 ? subplot(2,2,4)。 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. 峰值程序 ? for i=1:length(u) ? if u(i)==max(u) ? m=i。 %宮縮單波數(shù)據(jù) ? h=y1(k:j)。 處理程序 _續(xù) 3 ? 3. 面積程序 ? t=linspace(,*length(h),length(h))。 %心率 時(shí)間圖 。 % 宮縮 時(shí)間圖 。目前具有如下功能： ? （ 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ù)。 ? ③計(jì)算峰值或瞬時(shí)心率時(shí)，按著 [Alt]鍵，然后單擊曲線的選取點(diǎn)，接著單擊相應(yīng)的操作按鈕 ? （ 3）結(jié)束軟件使用 ? 單擊窗口的關(guān)閉按鈕。D:\matlab_dat39。 %清除變量 ? 心率： ? y1=load(39。)。39。 ? ④曲線作圖 ? 建立 XY坐標(biāo)軸與繪圖時(shí)間數(shù)組點(diǎn) ? t=linspace(,*length(y1), length(y1))。 %心率 時(shí)間圖 。 % 宮縮 時(shí)間圖； 實(shí)驗(yàn) _續(xù) 2 ? 宮縮局部加權(quán)平滑濾波 (七點(diǎn)平滑 ) ? n=length(y2) ? for i=4:n3 ? y2_7s(i)= (2*y2(i3)+3*y2(i2)+6*y2(i1)+7*y2(i)2*y2(i+3)+3*y2(i+2)+6*y2(i+1))/21。 ? y2_7s([1 2 3 n2 n1 n])=[y2(1) y2(2) y2(3) y2(n2) y2(n1) y2(n) ] ? 再次繪制心率與宮縮圖，觀察曲線變化的情況 ? subplot(2,1,1),plot(t,y1)。 ? subplot(2,1,2),plot(t,y2_7s)。 ? h=y1(k:j)。 ? subplot(2,1,1),plot(t,h)。 ? subplot(2,1,2),plot(t,u)。 ? ⑧計(jì)算宮縮面積 ? area=*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)行文件 ? 計(jì)算持續(xù)時(shí)間 ? t0=(length(u)1)*。觀察其與前面單步執(zhí)行的結(jié)果是否一致，給出你的結(jié)論。D:\matlab_dat39。39。 ? y2=load(39。)。 ? subplot(2,1,1),plot(t,y1)。 實(shí)驗(yàn) _續(xù) 5 ? 2. ? （ 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類型文件 。 ） ? （ 5） 再對(duì)上述選定區(qū)域的數(shù)據(jù)求取平均心率 、 時(shí)距 、面積 ， 然后人工選取峰值點(diǎn) ， 計(jì)算峰值和瞬時(shí)心率 。 比較其結(jié)
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1