医学图形图像处理识别技术与医学信息分析

1、 作业作业 第三讲第三讲 mathlab 与图形信号的处理与图形信号的处理 识别识别 nmatlab软件的使用介绍软件的使用介绍 n图形信号的处理与识别程序介绍图形信号的处理与识别程序介绍 n临床电生理数据分析实验模拟软件介绍临床电生理数据分析实验模拟软件介绍 n实验实验 一、一、matlab软件的使用介绍软件的使用介绍 n1. matlab的简介的简介 (1) matlab高级交互式软件包高级交互式软件包 (2) matlab典型应用典型应用 数学计算数学计算 算法推导算法推导 建模和仿真模拟建模和仿真模拟 数据分析和结果的可视化数据分析和结果的可视化 工程图形绘制应用程序开发（包括用户图形

2、界面的建立）工程图形绘制应用程序开发（包括用户图形界面的建立） (3) matlab产生的历史背景产生的历史背景 调用调用eispack和和linpack的的fortran子程序库子程序库 eispack和和linpack的接口程序的接口程序matlab 发展成发展成 为适合多学科，多种工作平台的功能强大的大型软件。为适合多学科，多种工作平台的功能强大的大型软件。 美国美国mathworks公司设计的公司设计的“科学科学 和工程计算的高级交互式软件包和工程计算的高级交互式软件包” cleve moler、new mexico、john little matlab软件软件_续续1 (3) mat

3、lab的语言特点的语言特点 n语言简洁紧凑，使用方便灵活，库函数极其丰富。语言简洁紧凑，使用方便灵活，库函数极其丰富。 n运算符丰富。运算符丰富。 nmatlab既具有结构化的控制语句，又有面向对象编程的既具有结构化的控制语句，又有面向对象编程的 特性。特性。 n程序限制不严格，程序设计自由度大。程序限制不严格，程序设计自由度大。 n程序的可移植性很好，基本上不做修改就可以在各种型号程序的可移植性很好，基本上不做修改就可以在各种型号 的计算机和操作系统上运行。的计算机和操作系统上运行。 nmatlab的图形功能强大。的图形功能强大。 nmatlab的缺点是程序的执行速度较慢。的缺点是程序的执行

4、速度较慢。 n功能强大的工具箱。功能强大的工具箱。 n源程序的开放性。源程序的开放性。 如：解方程组如：解方程组ax=b 4 3 2 1 35713498 65542343 12857923 67451332 x x=ab matlab软件软件_续续2 n2. matla的安装的安装 进入进入matlab软件安装程序目录，然后单击软件安装程序目录，然后单击“setup” 程序，在安装程序的引导下，逐步完成以下操作程序，在安装程序的引导下，逐步完成以下操作: (1)选择选择“install”安装或安装或”update license”升级许升级许 可证，单击可证，单击next；然后输入姓名、公司

5、与序列号，单；然后输入姓名、公司与序列号，单 击击next。 (2)选择接受许可证协议，单击选择接受许可证协议，单击“yes”单击单击next，然，然 后选择后选择“typical”典型安装，单击典型安装，单击next。 (3)选择或输入安装目的路径，单击选择或输入安装目的路径，单击next，后单击，后单击 “yes” (4)单击单击install，系统开始安装。，系统开始安装。 matlab软件软件_续续3 (5)在提示插入第在提示插入第2张张cd时，插入时，插入cd2，然后单，然后单 击击ok，接着在文件扩展名警示对话框内单击，接着在文件扩展名警示对话框内单击 yes to all。 (6

6、)单击单击next，接着单击，接着单击finish。 (7)如果在单击如果在单击finish按钮之前，选择了按钮之前，选择了 “start”选项，则此时进入选项，则此时进入matlab工作环境。工作环境。 matlab软件软件_续续4 n3. matlab启动与退出启动与退出 启动启动matlab，进入，进入matlab工作环境工作环境 启动启动matlab matlab软件软件_续续5 结束结束matlab，退出，退出matlab工作环境工作环境 n键入键入exit n键入键入quit n直接关闭直接关闭matlab的命令视窗的命令视窗 退出退出matlab matlab软件软件_续续6-1

7、 n4. matlab变量与基本输入输出变量与基本输入输出 变量命名的规则变量命名的规则：第一个字母必须是英文字第一个字母必须是英文字 母；母；字母间不可留空格；字母间不可留空格；最多只能有最多只能有19个个 字母，字母，matlab会忽略多余字母。会忽略多余字母。 变量定义与输入输出变量定义与输入输出 n向量与标量向量与标量 矩阵矩阵行和列组成的数据表行和列组成的数据表 行向量行向量矩阵只有一行时，称为行向量。矩阵只有一行时，称为行向量。 列向量列向量矩阵只有一列时，称为列向量。矩阵只有一列时，称为列向量。 标量标量矩阵只有矩阵只有1个元素时，称为标量。个元素时，称为标量。 matlab软件

8、软件_续续6-2 n变量定义变量定义 var=expt; nexpt可以是数字、变量、操作符和函数等组成。可以是数字、变量、操作符和函数等组成。 例如：例如： x=7 a=1 2 3; 4 5 6; 7 8 9 b=11 12 13 14 c=21;22;23;24 d=x*b 注意：注意： n用方括号用方括号“ ” 确定的范围内逐行给定元素，相同行元素由空确定的范围内逐行给定元素，相同行元素由空 格格“ ”或逗号或逗号“,”分隔，不同行数据由分号分隔，不同行数据由分号“;”或回车键分隔。或回车键分隔。 n语句结尾没有分号的命令在屏幕显示结果；有分号的命令只执语句结尾没有分号的命令在屏幕显示结

9、果；有分号的命令只执 行而不显示。行而不显示。 matlab软件软件_续续7 n5. matlab常用命令常用命令 (1) matlab常用的基本数学函数常用的基本数学函数 nabs(x)：纯量的绝对值或向量的长度：纯量的绝对值或向量的长度 nsqrt(x)：开平方：开平方 nround(x)：四舍五入至最近整数：四舍五入至最近整数 nfix(x)：无论正负，舍去小数至最近整数：无论正负，舍去小数至最近整数 nsign(x)：符号函数：符号函数 (signum function) nexp(x)：自然指数：自然指数 npow2(x)：2的指数的指数 nlog(x)：以：以e为底的对数为底的对数

10、 nlog10(x)：以：以10为底的对数为底的对数 举例举例: 4x2+58x-89=0（a=4, b=58, c=89） a acbb x a acbb x 2 4 2 4 2 2 2 1 matlab软件软件_续续8 (2) 适用于向量的常用函数适用于向量的常用函数 nmin(x): 向量向量x的元素的最小值的元素的最小值 nmax(x): 向量向量x的元素的最大值的元素的最大值 nmean(x): 向量向量x的元素的平均值的元素的平均值 nmedian(x): 向量向量x的元素的中位数的元素的中位数 nstd(x): 向量向量x的元素的标准差的元素的标准差 ndiff(x): 向量向量

11、x的相邻元素的差的相邻元素的差 nsort(x): 对向量对向量x的元素进行排序（的元素进行排序（sorting） nlength(x): 向量向量x的元素个数的元素个数 nsum(x): 向量向量x的元素总和的元素总和 举例举例:对对130名正常成年男子红名正常成年男子红 细胞数细胞数(rbc ,万万/mm)的数据进的数据进 行描述性统计。行描述性统计。 (x=zeros(130,1) 平均、中位数、众数、标准差平均、中位数、众数、标准差 matlab软件软件_续续9 适用于向量的常用函数适用于向量的常用函数_续续 nprod(x): 向量向量x的元素总乘积的元素总乘积 ncumsum(x)

12、: 向量向量x的累计元素总和的累计元素总和 ncumprod(x): 向量向量x的累计元素总乘积的累计元素总乘积 ndot(x, y): 向量向量x和和y的内积的内积 ncross(x, y): 向量向量x和和y的外积的外积 matlab软件软件_续续10 一些常用的操作函数一些常用的操作函数 nformat语句是控制输出的显示格式，句法如下：语句是控制输出的显示格式，句法如下： format 设置默认显示格式设置默认显示格式5位的定点数位的定点数 format type 设置按指定的设置按指定的type格式输出显示。格式输出显示。 n例如：例如：format rat ratio of sma

13、ll integers format(type) 与上面语句一样，但以函数形式出现与上面语句一样，但以函数形式出现 ndisp(x) 显示矩阵函数显示矩阵函数 n% 注释语句注释语句 nzeros、ones 预先定义矩阵的大小。预先定义矩阵的大小。 例如：例如：h = zeros(6); %生成生成66的的0矩阵矩阵 matlab软件软件_续续11 (3) 重复命令重复命令 n1）for循环（循环（for-loop），其基本形式为：），其基本形式为： for 变量变量 = 矩阵；矩阵； 运算式；运算式； end 变量值会被依次设定为矩阵的每一行，来执行变量值会被依次设定为矩阵的每一行，来执行f

14、or和和end之之 间的运算式。因此，若中断情况，运算式执行的次数等于间的运算式。因此，若中断情况，运算式执行的次数等于 矩阵的行数。其中矩阵的行数。其中“;”是终结内语句的可抑制的重复显示。是终结内语句的可抑制的重复显示。 n举例举例:产生一个长度为产生一个长度为6的调和数列的调和数列 x = zeros(1,6); % x是一个是一个16的矩阵的矩阵 for i = 1:6, x(i) = 1/i; end matlab软件软件_续续12 n2）while循环，其基本形式为：循环，其基本形式为： while 条件式；条件式； 运算式；运算式； end 实现功能：只要条件成立，就会一再执行运

15、算式。实现功能：只要条件成立，就会一再执行运算式。 n举例举例:用用while循环产生一个长度为循环产生一个长度为6的调和数列的调和数列 x = zeros(1,6); % x是一个是一个16的零矩阵的零矩阵 i = 1; while i 0.5 disp(given random number is greater than 0.5.); else disp(given random number is less than 0.5.); end matlab软件软件_续续14 (5)将多个命令写入一个将多个命令写入一个m文件文件 n若要一次执行大量的若要一次执行大量的matlab命令，可将这

16、些命令命令，可将这些命令 存于一个文件类型为存于一个文件类型为m的文件中，并在的文件中，并在matlab提提 示号下键入此文件名即可执行。包含示号下键入此文件名即可执行。包含matlab命令命令 的文件以的文件以m为扩展名，通称为扩展名，通称m文件（文件（m-files）。）。 n举例举例:名为名为test.m的的m文件有以下内容：文件有以下内容： % this is my first test m-file. % roger jang, march 3, 1997 fprintf(start of test.m!n); for i = 1:3, fprintf(i = %d - i3 = %

17、dn, i, i3); end fprintf(end of test.m!n); matlab软件软件_续续15 建立建立test.m的的m文件文件 进入进入matlab编辑器编辑器 nfilenewm-file 在编辑器内输入以下内容在编辑器内输入以下内容 % this is my first test m-file. % roger jang, march 3, 1997 fprintf(start of test.m!n); for i = 1:3, fprintf(i = %d - i3 = %dn, i, i3); end fprintf(end of test.m!n); 保存文

18、件到指定或默认的目录下保存文件到指定或默认的目录下 nfilesave输入文件名输入文件名“test.m” matlab软件软件_续续16 执行文件执行文件“test.m” n进入文件所在目录或查看文件进入文件所在目录或查看文件 cd dir n执行执行test.m test start of test.m! i = 1 - i3 = 1 i = 2 - i3 = 8 i = 3 - i3 = 27 end of test.m! matlab软件软件_续续17 n注意注意：m档案可分为档案可分为命令集命令集(scripts)及及函数函数(functions)。前述。前述 的的test.m即为命

19、令集，其效用和将命令逐一输入完全一样，因即为命令集，其效用和将命令逐一输入完全一样，因 此在命令集直接使用工作空间的变量，则在命令集中设定的变此在命令集直接使用工作空间的变量，则在命令集中设定的变 量，可在工作空间中看得到。量，可在工作空间中看得到。函数则需要用到输入参数函数则需要用到输入参数(input arguments)和和输出参数输出参数(output arguments)来传递资讯。来传递资讯。 n例如：计算一个正整数的阶乘（例如：计算一个正整数的阶乘（factorial），可以写一个如下），可以写一个如下 的的matlab函数并将之存档于函数并将之存档于fact.m： functi

20、on output = fact(n) % fact calculate factorial of a given positive integer. output = 1; for i = 1:n, output = output*i; end 其中其中fact是函数名，是函数名，n是输入参数，是输入参数，output是输出参数，而是输出参数，而i则是此则是此 函数用到的暂时变量。函数用到的暂时变量。 matlab软件软件_续续18 要使用此函数，直接键入函数名及适当输入参数值即可。要使用此函数，直接键入函数名及适当输入参数值即可。 y = fact(5) y = 120 注意：注意： n在

21、执行在执行fact之前，你必须先进入之前，你必须先进入fact.m所在的目录。所在的目录。 n在执行在执行fact(5)时，时，matlab会进入一个下层的暂时工作会进入一个下层的暂时工作 空间（空间（temperary workspace），将变量），将变量n的值设定为的值设定为 5，然后进行各项函数的内部运算，所有内部运算所产，然后进行各项函数的内部运算，所有内部运算所产 生的变量都存在此暂时工作空间中。运算完毕后，生的变量都存在此暂时工作空间中。运算完毕后， matlab会将最后输出参数会将最后输出参数output的值设定给上层的变的值设定给上层的变 量量y，并将清除此暂时工作空间及其所

22、含的所有变量。，并将清除此暂时工作空间及其所含的所有变量。 matlab软件软件_续续19 (6)matlab的文件、目录等操作指令的文件、目录等操作指令 n显示当前目录显示当前目录 pwd 显示当前工作目录显示当前工作目录 s = pwd 将当前工作目录赋予变量将当前工作目录赋予变量s n改变工作目录改变工作目录 cd 显示当前工作目录显示当前工作目录 w = cd 将当前工作目录赋予将当前工作目录赋予w cd(directory) 设置当前工作目录为设置当前工作目录为directory指定指定 的目录的目录 cd(.) 将当前工作目录改为它的上级目录将当前工作目录改为它的上级目录 cd d

23、irectory or cd . 与上面两条指令等同与上面两条指令等同 matlab软件软件_续续20 matlab的文件、目录等操作指令的文件、目录等操作指令 _续续 n显示目录列表显示目录列表 dir 列示当前工作目录的文件列示当前工作目录的文件 dir name 列示指定的文件列示指定的文件 files = dir(directory) 将指定目录的文件列表信息返回将指定目录的文件列表信息返回 到到files中。中。 例如：例如： afiles=dir(d:/work/*.*) afiles(3).name afiles(3).date n显示文件内容显示文件内容 type (filen

24、ame) type filename nfilename是要被显示的文件的路径和名字是要被显示的文件的路径和名字 matlab软件软件_续续21 matlab的文件、目录等操作指令的文件、目录等操作指令 _续续 n浏览或改变浏览或改变matlab 目录的搜寻路径目录的搜寻路径 path 显示当前显示当前matlab 搜索路径。搜索路径。 n初始搜索路径由初始搜索路径由toolbox/local/pathdef.m文件定义的。文件定义的。 path(newpath) 将搜索路径改变为一个新路径，这里将搜索路径改变为一个新路径，这里 新路径是一个目录的字符串数组。新路径是一个目录的字符串数组。 p

25、ath(path,newpath) 将一个新路径增加到当前的搜索将一个新路径增加到当前的搜索 路径后。路径后。 path(newpath,path) 将一个新路径加到当前的搜索路将一个新路径加到当前的搜索路 径前。径前。 p = path(.) 将指定的路径返回到一个字符串变量。将指定的路径返回到一个字符串变量。 matlab软件软件_续续22 n浏览或改变浏览或改变matlab 目录的搜寻路径目录的搜寻路径 _续续 举例举例: path matlabpath d:matlab5toolboxmatlabgeneral d:matlab5toolboxmatlabops d:matlab5to

26、olboxsimulinkdee d:matlab5toolboxlocal n查询某一命令是在搜寻路径的何处，用查询某一命令是在搜寻路径的何处，用which命令：命令： which test c:datamlbooktest.m n将将c:datamlbook加入加入matlab的搜寻路径，用的搜寻路径，用path path(path, c:datamlbook) matlab软件软件_续续23 (7)资料的储存与载入资料的储存与载入 _续续 n1) 将计算所得的储存倒文件的命令将计算所得的储存倒文件的命令save save 将工作空间的所有变量储存到名为将工作空间的所有变量储存到名为mat

27、lab.mat的的 二进制文件。二进制文件。 save filename 将工作空间的所有变量储存到名为将工作空间的所有变量储存到名为 filename.mat的二进制文件。的二进制文件。 save filename x y z 将变量将变量x、y、z储存到名为储存到名为filename.mat 的二进制文件。的二进制文件。 save filename x ascii 将变量将变量x以八位数存到名为以八位数存到名为filename的的 ascii文件。文件。 save filename x -ascii double 将变量将变量x以十六位数存到名为以十六位数存到名为 filename的的as

28、cii文件。文件。 n注意注意:命令后不加任何选项时，命令后不加任何选项时，save会将变量以二进制会将变量以二进制 （binary）的方式储存至扩展名为）的方式储存至扩展名为mat的文件的文件 ；另一个选项另一个选项 是是-tab，可将同一列相邻的数目以定位键（，可将同一列相邻的数目以定位键（tab）隔开。）隔开。 matlab软件软件_续续24 n将计算所得的储存倒文件的命令将计算所得的储存倒文件的命令save _续续 举例：举例： who 列出工作空间的变量列出工作空间的变量 your variables are: b h j y ans i x z save test b y 将变数将

29、变数b与与y储存至储存至test.mat dir 列出现在目录中的文件列出现在目录中的文件 . 2plotxy.doc fact.m simulink.doc test.m $1basic.doc . 3plotxyz.doc first.doc temp.doc test.mat 1basic.doc book.dot go.m template.doc testfile.dat delete test.mat 删除删除test.mat matlab软件软件_续续25 n2)将文件载入储存得变量命令将文件载入储存得变量命令load load filename 寻找名称为寻找名称为filena

30、me.mat的文件，并以二进制格的文件，并以二进制格 式载入。若找不到式载入。若找不到filename.mat，则寻找名称为，则寻找名称为filename的文件，的文件， 并以并以ascii格式载入。格式载入。 load filename ascii 寻找名称为寻找名称为filename的文件，并以的文件，并以ascii 格式载入。若以格式载入。若以ascii格式载入，则变量名称即为文件名称（但不格式载入，则变量名称即为文件名称（但不 包含扩展名）。若以二进制载入，则可保留原有的变量名称。包含扩展名）。若以二进制载入，则可保留原有的变量名称。 举例：举例： clear all; % 清除工作空

31、间中的变数清除工作空间中的变数 x = 1:10; save testfile.dat x -ascii % 将将x以以ascii格式存至格式存至testfile.dat文件文件 load testfile.dat % 载入载入testfile.dat who % 列出工作空间中的变量列出工作空间中的变量 your variables are: testfile x n注意：上述过程由于是以注意：上述过程由于是以ascii格式储存与载入，所以产生了一个与格式储存与载入，所以产生了一个与 文件名称相同的变量文件名称相同的变量testfile，此变量的值和原变量，此变量的值和原变量x完全相同。完全

32、相同。 matlab软件软件_续续26 (8)基本基本xy平面绘图命令平面绘图命令 n1) length 获得一组数的个数或一个向量的长度获得一组数的个数或一个向量的长度 x = ones(1,8); n = length(x) n2) linspace 产生线性间隔向量产生线性间隔向量 y = linspace(a,b) 在在a,b之间产生一组之间产生一组100点的线性间隔点的线性间隔 的行向量。的行向量。generates a row vector y of 100 points linearly spaced between and including a and b. y = lins

33、pace(a,b,n) 在在a,b之间产生一组之间产生一组n点的线性间隔的点的线性间隔的 行向量。行向量。generates a row vector y of n points linearly spaced between and including a and b. 举例：举例： t=linspace(0.005,0.005*length(y1), length(y1); matlab软件软件_续续27 基本基本xy平面绘图命令平面绘图命令_续续 n3）plot 绘制一维曲线的基本函数，使用此函数之前，绘制一维曲线的基本函数，使用此函数之前， 需先定义曲线上每一点的需先定义曲线上每一点的

34、x及及y座标。座标。 举例举例:画出一条正弦曲线画出一条正弦曲线 close all; x=linspace(0, 2*pi, 100); % 1产生产生00个点的个点的x座标座标 y=sin(x); %生成对应的生成对应的y座标座标 plot(x,y); 关于关于matlab基本绘图函数的应用说明基本绘图函数的应用说明 : 绘图刻度绘图刻度 nplot: x轴和轴和y轴均为线性刻度（轴均为线性刻度（linear scale） nloglog: x轴和轴和y轴均为对数刻度（轴均为对数刻度（logarithmic scale） nsemilogx: x轴为对数刻度，轴为对数刻度，y轴为线性刻度轴

35、为线性刻度 nsemilogy: x轴为线性刻度，轴为线性刻度，y轴为对数刻度轴为对数刻度 matlab软件软件_续续28 关于关于matlab基本绘图函数的应用说明基本绘图函数的应用说明 _续续 画出多条曲线，只需将座标对依次放入画出多条曲线，只需将座标对依次放入plot函数函数 nplot(x, sin(x), x, cos(x); 若要改变颜色，在座标对後面加上相关字串若要改变颜色，在座标对後面加上相关字串 nplot(x, sin(x), c, x, cos(x), g); 若要同时改变颜色及图线型态（若要同时改变颜色及图线型态（line style），也是在座标），也是在座标 对後面

36、加上相关字串对後面加上相关字串 nplot(x, sin(x), co, x, cos(x), g*); matlab软件软件_续续29 关于关于matlab基本绘图函数的应用说明基本绘图函数的应用说明 _续续 plot绘图函数的叁数说明绘图函数的叁数说明 n在用在用plot(a,b,s)绘图时，可以使用各类线型、图符、颜绘图时，可以使用各类线型、图符、颜 色，这里色，这里s是由下表中的任一元素或是由下表中的任一元素或3列中所有元素组合列中所有元素组合 的字符串。的字符串。 matlab软件软件_续续30 n4) 将一图形窗口划分为多个显示区，使之同时画出将一图形窗口划分为多个显示区，使之同时

37、画出 数个小图形於同一个视窗之中数个小图形於同一个视窗之中。 subplot(m,n,p) 实现功能：将图形窗口划分为实现功能：将图形窗口划分为mn矩阵的显示区，同时激矩阵的显示区，同时激 活第活第p个显示区个显示区。 举例：举例： nsubplot(2,2,1); plot(x, sin(x); nsubplot(2,2,2); plot(x, cos(x); nsubplot(2,2,3); plot(x, sinh(x); nsubplot(2,2,4); plot(x, cosh(x); matlab软件软件_续续31 n5) 其他各种二维绘图函数其他各种二维绘图函数 bar 长条图长

38、条图 errorbar 图形加上误差范围图形加上误差范围 fplot 较精确的函数图形较精确的函数图形 polar 极座标图极座标图 hist 累计图累计图 rose 极座标累计图极座标累计图 stairs 阶梯图阶梯图 stem 针状图针状图 fill 实心图实心图 二、二、图形信号的处理程序介绍图形信号的处理程序介绍 n1. 峰值程序峰值程序 nfor i=1:length(u) n if u(i)=max(u) n m=i; n break n end nend np_h=h(m) 处理程序处理程序_续续1 n求起点和终点的程序求起点和终点的程序 %judge_start for i=5

39、:(length(y2_7s)-5) if y2_7s(i-4)=0 %宫缩单波数据宫缩单波数据 h=y1(k:j); %对应单波宫缩的心率数据对应单波宫缩的心率数据 t0=(length(u)-1)*0.005; 处理程序处理程序_续续3 n3. 面积程序面积程序 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 % 计算面积计算面积 三、临

40、床电生理数据分析实验三、临床电生理数据分析实验 模拟软件介绍模拟软件介绍 n1.实验模拟软件简介实验模拟软件简介 临床电生理数据分析实验模拟软件临床电生理数据分析实验模拟软件med.exe是用于对临是用于对临 床功能类的信号进行初步分析的软件，该软件是用床功能类的信号进行初步分析的软件，该软件是用 visual basic程序语言来编写的。目前具有如下功能：程序语言来编写的。目前具有如下功能： n（1）可以浏览医学仪器采集的电生理信号数据，人工干预信）可以浏览医学仪器采集的电生理信号数据，人工干预信 号截取，并以二进制格式和文本格式来保存截取的信号数据；号截取，并以二进制格式和文本格式来保存截

41、取的信号数据； n（2）可以利用面板按扭工具，对胎心率信号、孕妇的宫缩信）可以利用面板按扭工具，对胎心率信号、孕妇的宫缩信 号进行人工干预的信号处理和分析，包括：奇异点处理，平均号进行人工干预的信号处理和分析，包括：奇异点处理，平均 心率、瞬时心率、宫缩峰值、宫缩时距和宫缩面积的计算；心率、瞬时心率、宫缩峰值、宫缩时距和宫缩面积的计算； n（3）可以进行信号数字平滑滤波、基线调整、自动计算部分）可以进行信号数字平滑滤波、基线调整、自动计算部分 时域参数。时域参数。 n2. 软件的使用说明软件的使用说明 实验模拟软件介绍实验模拟软件介绍_续续1 （1）原始采样数据回显）原始采样数据回显 n“临床

42、信号截取临床信号截取”“选取数据源位置选取数据源位置”确定位置后，单击确定位置后，单击ok n在输入病历号的输入框内输入：在输入病历号的输入框内输入：a0012 n在选择文件框内选择项目：在选择文件框内选择项目：r0012a_5 n单击单击开始开始单击单击回放回放 （2）进行人工干预操作）进行人工干预操作 n在需要截波或人工分析时，单击在需要截波或人工分析时，单击冻结冻结 n截取波形，计算时距、面积和平均心率，或去除奇异项时，要用鼠截取波形，计算时距、面积和平均心率，或去除奇异项时，要用鼠 标在视图浏览区内拖出选区，然后分别单击相应的操作按钮标在视图浏览区内拖出选区，然后分别单击相应的操作按钮

43、 注意注意：如要取消选区，可以单击浏览视图的任意位置即可。：如要取消选区，可以单击浏览视图的任意位置即可。 n计算峰值或瞬时心率时，按着计算峰值或瞬时心率时，按着alt键，然后单击曲线的选取点，接键，然后单击曲线的选取点，接 着单击相应的操作按钮着单击相应的操作按钮 （3）结束软件使用）结束软件使用 n单击窗口的关闭按钮。单击窗口的关闭按钮。 四、实验四、实验 n1均值计算均值计算 n2. 时距计算时距计算 n3. 峰值计算峰值计算 n4. 面积计算面积计算 n5. 作图作图 平均值平均值 ： n i i x n x 1 1 平均幅值平均幅值 ： n i i x n x 1 1 tnnt st

44、artend )( )(max 1 i n i p xx end starti i txs 实验实验_续续1 n1. matlab实验内容的介绍实验内容的介绍 （1）利用）利用matlab计算宫缩的峰值、宫缩时的平均心率，以及宫缩峰值计算宫缩的峰值、宫缩时的平均心率，以及宫缩峰值 对应的顺时心率、计算宫缩波的面积、宫缩持续时间对应的顺时心率、计算宫缩波的面积、宫缩持续时间 实验步骤实验步骤 n启动启动matlab； n增加路径：增加路径： n path(path,d:matlab_dat) n导入数据导入数据 nclear all; %清除变量清除变量 n心率：心率： n y1=load(r0

45、012a0h.txt); n宫缩：宫缩： n y2=load(r0012a0u.txt); n曲线作图曲线作图 n建立建立xy坐标轴与绘图时间数组点坐标轴与绘图时间数组点 nt=linspace(0.005,0.005*length(y1), length(y1); n绘制心率与宫缩图绘制心率与宫缩图 nsubplot(2,1,1),plot(t,y1); %心率心率-时间图时间图; nsubplot(2,1,2),plot(t,y2); % 宫缩宫缩-时间图；时间图； 实验实验_续续2 n宫缩局部加权平滑滤波宫缩局部加权平滑滤波(七点平滑七点平滑) n n=length(y2) nfor i

46、=4:n-3 n y2_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; n end; ny2_7s(1 2 3 n-2 n-1 n)=y2(1) y2(2) y2(3) y2(n-2) y2(n-1) y2(n) n再次绘制心率与宫缩图，观察曲线变化的情况再次绘制心率与宫缩图，观察曲线变化的情况 nsubplot(2,1,1),plot(t,y1); %心率心率-时间图时间图; nsubplot(2,1,2),plot(t,y2_7s); % 宫缩宫缩-时间图；时间图； n判断宫缩波的起始、终止点判断宫缩波的起始、终止点 n ju