心率监测仪.doc

黄石理工学院课程设计（设计者：占东杰）摘 要七段数码管作为显示系统，以AT89S52单片机为核心，光电式脉搏波传感器采集信号，经信号处理电路后脉冲送入单片机，由数码管显示心率。心率监测仪通电后，数码管应显示0。把手指轻轻置于右下角的传感器中，以稍微有压迫感为宜。很快红色发光二极管会伴随你的脉搏而闪烁，并伴有蜂鸣器响，同时数码管显示出你的心率。如果偶尔出现显示不稳的情况，请按复位键对系统进行复位。本测试器适用于79V直流电压，正常工作电流约为100mA。关键词：AT89C2051、LM324、ULN2003 24 目 录1 前言12 系统总体设计2 2.1 医学常识22.2心率监测仪系统总框图 32.3心率监测仪工作原理 33 硬件设计43.1 主要芯片介绍4 3.1.1 AT89C2051内部结构与功能 4 3.1.2 ULN2003的功能5 3.1.3 LM324的功能53.2 传感器与信号处理电路63.2.1 光电式脉搏波传感器 63.2.2 前置放大与滤波电路 73.2.3 波形变换电路 83.3 显示电路93.4报警电路 93.5时钟电路 103.6复位电路 103.7单片机处理与显示部分电路图 114 软件设计 124.1程序流程图 124.1.1主程序流程图 124.1.2定时器中断程序流程图 124.1.3INT中断程序流程图 134.1.4显示程序流程图 134.2源程序 135调试与仿真 195.1 Proteus7.12与KeiluVision219结论 19致谢 19参考文献 20附录 21附录A 传感器与信号处理电路图 21附录 B 心率监测仪电气原理图21附录 C AT89C2051结构框图211 前言 随着人民生活水平的提高，心脏疾病的发病率呈上升趋势，已成为威胁人类身体健康的杀手之一。因为心脏病的发作具有突发性和随机性，为患者性实时监测发展的必然趋势。随着电子科技的不断发展，生命科学和信息科学的结合越来越紧密，许多研究人员都投身于人类的健康事业发展中。Holter的出现，使心电图机进入家庭变成了可能，但基于心电工作站的模式，使个别地区的患者因医院分析诊断系统的不健全，而变得不适用；基于嵌入式及DSP的心电监护仪功能强大，但又因芯片价格的高昂而有悖于我国基本国情，不利于家庭的普及。因此，一种性能优良，带有自动检测诊断功能，而又符合我国人均收入水平不高这一国情的心率监测仪的研制显得尤其重要。基于这一目的，这次课题是关于医疗器械的设计，我的课题是设计心率监测仪。本文所介绍的心率监测仪运用了单片机技术，与采用一般数字电路制作的心率计相比较，其体积更小巧，抗干扰能力更强，而且单片机的可编程性使得其适应性、灵活性大大增强。它使用方便，只需将手指端轻轻放在传感器上，即可实时显示出你的每分钟脉搏次数，特别适合体育训练和外出旅游等场合使用。采用红外光学检测法，摒弃了不便于运动状态下测量脉搏的听诊器和吸附在人体上的电极等老式测量方法。检测的基本原理是：随着心脏的搏动，人体组织半透明度随之改变：当血液送到人体组织时，组织的半透明度减小；当血液流回心脏，组织则半透明度增大。这种现象在人体组织较薄的手指尖、耳垂等部位最为明显。因此，本心率监测仪将红外发光二极管产生的红外线照射到人体的上述部位，并用装在该部位另一侧或旁边的红外光电管来检测机体组织的透明程度并把它转换成电信号。由于此信号的频率与人体每分钟的脉搏次数成正比，故只要把它转换成脉沖并进行整形、计数和显示，即可实时地测出脉搏的次数。2 系统总体设计2.1 医学常识心率(Heart Rate)：用来描述心动周期的专业术语，是指心脏每分钟跳动的次数,以第一声音为准。心电信号是一种非常弱且频率较低的信号，一般幅值在0.055mV，频率在0.05100Hz结构特征：心脏是一厚壁的肌性器官，由左右2心房和左右2心室 4个心腔组成。心脏的自动节律性收缩，推动血液在循环系统的各种血管中环流，使机体各组织、器官能不断地吐故纳新、新陈代谢。正常成年人安静时的心率有显著的个体差异，平均在75次/分左右(60100次/分之间)。心率可因年龄、性别及其它生理情况而不同。初生儿的心率很快，可达130次/分以上。在成年人中，女性的心率一般比男性稍快。同一个人，在安静或睡眠时心率减慢，运动时或情绪激动时心率加快，在某些药物或神经体液因素的影响下，会使心率发生加快或减慢。经常进行体力劳动和体育锻炼的人，平时心率较慢。健康成人的心率为60100次/分，大多数为6080次/分，女性稍快；3岁以下的小儿常在100次/分以上；老年人偏慢。成人每分钟心率超过100次（一般不超过 160次/分）或婴幼儿超过 150次/分者，称为窦性心动过速。常见于正常人运动、兴奋、激动、吸烟、饮酒和喝浓茶后。也可见于发热、休克、贫血、甲亢、心力衰竭及应用阿托品、肾上腺素、麻黄素等。如果心率在 160220次/分，常称为阵发性心动过速。心率低于60次/分者（一般在40次/分以上），称为窦性心动过缓。可见于长期从事重体力劳动和运动员；病理性的见于甲状腺机能低下、颅内压增高、阻塞性黄疸、以及洋地黄、奎尼丁或心得安类药物过量或中毒。如心率低于40次/分，应考虑有房室传导阻滞。心率过快超过160次/分，或低于40次/分，大多见于心脏病病人，病人常有心悸、胸闷、心前区不适，应及早进行详细检查，以便针对病因进行治疗。心电图如图1-1-1所示图1-1-1 心电图脉搏波：人体心室周期性的收缩和舒张导致主动脉的收缩和舒张,使血流压力以波的形式从主动脉根部开始沿着整个动脉系统传播,这种波称为脉搏波。脉搏波所呈现出的形态、强度、速率和节律等方面的综合信息,很大程度上反映出人体心血管系统中许多生理病理的血流特征。正常人脉搏数为6080次/分钟，婴儿为90100次/分钟，老人为100150次/分钟。脉搏波如图1-1-2所示 图1-1-2 脉搏波心率与脉搏的联系：心率与脉搏在身体正常的时候是相等的。在房颤等心脏疾病时候可出现不等。因此心率测量问题可以转化为脉搏的测量，而脉搏的测量有更容易实现特点，在实际应用中得到广泛运用。本监测仪的有效测量范围为50次199次/分钟。2.2心率监测仪系统总框图 心率监测仪由单片机AT89C2051、复位电路、时钟电路、传感器与信号处理电路、显示电路和报警电路等组成，如图2-1-1所示。AT89C2051单片机显 示 电 路传感器与信号处理电路复 位 电 路报 警 电 路时 钟 电 路 图2-1-1 心率监测仪系统总框图2.3心率监测仪工作原理光电式脉搏波传感器从手指获取脉搏信号，通过前置放大和滤波电路，再经过波形变换电路后，将脉搏信号转换成脉冲信号传给单片机。显示器将单片机处理后的信号显示出来，同时每次脉冲的到来均响铃，与脉搏同步。3 硬件设计3.1主要芯片介绍3.1.1 AT89C2051内部结构与功能一、AT89C2051主要性能 AT89C2051是ATMEL公司生产的带2K字节闪速可编程可擦除只读存储器(EEPROM)的8位单片机,它具有如下主要特性： 和MCS-51产品的兼容 15根可编程I/O引线2K字节可重编程闪速存储器 两个16位定时器/计数器耐久性：1,000写擦除周期 六个中断源2.7V6V的操作范围 可编程串行UART通道全静态操作：0Hz24MHz 直接LED驱动输出两级加密程序存储器 低功耗空载和掉电方式1288位内部RAM 二、AT89C2051的结构框图 AT89C2051是一带有2K字节闪速可编程可擦除只读存储体(EEPROM)的低电压,高性能8位CMOS微型计算机。它采用ATMEL的高密非易失存储技术制造并和工业标准MCS51指令集和引脚结构兼容。通过在单块芯片上组合通用的CPL1和闪速存储器,ATMEL的AT89C2051是一强劲的微型计算机,它对许多嵌入式控制应用提供一高度灵活和成本低的解决办法。此外,从AT89C2051内部结构图也可看出,其内部结构与8051内部结构基本一致（除模拟比较器外）,引脚RST、XTAL1、XTAL2的特性和外部连接电路也完全与51系列单片机相应引脚一致,但P1口、P3口有其独特之处。结构框图见附录C三、AT89C2051的引脚说明AT89C2051是一个有20个引脚的芯片,引脚如图3-1-1所示,与8051内部结构进行对比可发现,AT89C2051减少了两个对外端口（即P0、P2口）,使它最大可能地减少了对外引脚,因而芯片尺寸有所减少。 图3-1-1 AT89C2051引脚图AT89C2051芯片的20个引脚功能为： (1).Vcc：电源电压。 (2). GND：地。 (3). P1口：P1口是一8位双向I/O口。口引脚P1.2P1.7提供内部上拉电阻 P1.0和P1.1要求外部上拉电阻。P1.0和P1.1还分别作为片内精密模拟比较器的同相输入(AIN0)和反相输入（AIN1)。P1口输出缓冲器可吸收20mA电流并能直接驱动LED显示。当P1口引脚写入“1”时,其可用作输入端。当引脚P1.2P1.7用作输入并被外部拉低时,它们将因内部的上拉电阻而流出电流(IIL)。P1口还在闪速编程和程序校验期间接收代码数据。(4). P3口：P3口的P3.0P3.5、P3.7是带有内部上拉电阻的七个双向I/0引脚。P3.6用于固定输入片内比较器的输出信号并且它作为一通用I/O引脚而不可访问。P3口缓冲器可吸收20mA电流。当P3口引脚写入“1”时,它们被内部上拉电阻拉高并可用作输入端。用作输入时,被外部拉低的P3口引脚将用上拉电阻而流出电流(IIL)。P3口还用于实现AT89C2051的各种功能,如下表3.1.1所示。 P3口还接收一些用于闪速存储器编程和程序校验的控制信号。 (5). RST：复位输入。RST一旦变成高电平,所有的I/O引脚就复位到“1”。当振荡器正在运行时,持续给出RST引脚两个机器周期的高电平便可完成复位。每一个机器周期需12个振荡器或时钟周期。(6). XTAL1：作为振荡器反相放大器的输入和内部时钟发生器的输入。 (7). XTAL2：作为振荡器反相放大器的输出。 四、 表3-1-1 P3口的功能P3口引脚 功能P3.0RXD(串行输入端口) P3.1TXD(串行输出端口) P3.2INT0(外中断0) P3.3INT1(外中断1) P3.4 TO(定时器0外部输入) P3.5T1(定时器1外部输入) 3.1.2 ULN2003的功能ULN2003(国产型号为5G1413)是七路达林顿驱动器阵列，是个集电极开路(OC)输出的反向器.最大驱动电流可以达到500mA，并且能够在关态时承受 50V的电压，输出还可以在高负载电流并行运行。 采用DIP16或SOP16 塑料封装，如图3-1-3所示。图3-1-3 ULN2003引脚图、内部方框图3.1.3 LM324的功能LM324是四运放集成电路，它采用14脚双列直插塑料封装，外形如图所示。它的内部包含四组形式完全相同的运算放大器，除电源共用外，四组运放相互独立。具有电源电压范围宽，静态功耗小，可单电源使用，价格低廉等优点，因此被广泛应用在各种电路中，引脚如图3-1-4所示。图3-1-4 LM324引脚图3.2 传感器与信号处理电路3.2.1 光电式脉搏波传感器一、光电式脉搏波传感器的原理 其原理是利用光电信号来测量脉搏容量的变化。当血管内血容量变化时，组织对光的吸收程度相对发生变化，利用光电传感器可测出这种变化，该变化反映出血液动脉的基本参数情况（包括频率、幅度、脉搏波形状的改变）。根据郎伯特比尔（lambertbeer）定律,物质在一定波长处的吸光度和他的浓度成正比,当恒定波长的光照射到人体组织上时,通过人体组织吸收、反射衰减后测量到的光强在一定程度上反映了被照射部位组织的结构特征。光电式脉搏波传感器可分为反射型光电式脉搏波传感器和透射型光电式脉搏波传感器。可利用发光二极管做光源，光敏晶体管做光检测器件组成光电式脉搏波传感器。将发光二极管和光敏晶体管分别放在组织的两边（透射法）或同一侧（反射法），当被测处血管中的血液流动改变时，此处组织的透光率和反射率随之变化，光敏晶体管就可将此时引起的光线变化转换为相应的电信号。反射型光电式脉搏波传感器：光源和光敏元件处于被测部位的同一侧，光源所发处的光线经指尖组织的反射有一小部分可以透过指尖组织和血液照射到光敏元件所在的窗口中，再由光敏元件检测出其变化，因此血液的搏动情况可以被描记，通常称这种传感器为反射型光电式脉搏波传感器，其电路如图3-2-1所示。图3-2-1反射型光电式脉搏波传感器 光源和光敏元件分别处于被测部位的两侧，光源发出的光线可以经指尖部组织透射到光敏元件所在的窗口，从而由光敏元件检测出脉搏的波动信号，这样记录的波也有将其作为指尖容积波处理，通常称这种传感器为透射型光电式脉搏波传感器，其电路如图3-2-2所示。图3-2-2 透射型光电式脉搏波传感器本次设计选用透射型光电式脉搏波传感器，其电路如图3-2-3所示。图3-2-3透射型光电式脉搏波传感器电路图因为传感器输出信号的频率很低，如当脉搏为50次/分钟时，只有0.78Hz，200次/分钟时也只有3.33Hz，因此信号首先经R14、C8组成的低通滤波器滤除高频干扰，当传感器与手指断开或检测到较强的干扰光线时，输出端的直流电压会出现很大变化，用C6、C7背靠背串联组成的双极性耦合电容把它隔断，滤除直流成分。3.2.2 前置放大与滤波电路运放IC1a将信号放大200倍，并与R11、C5组成截止频率为10Hz左右的低通滤波器以进一步滤除残留的干扰。9v电源电压由R3、R4分压成4.5v，再经IC1b缓冲后用作IC1 a的参考电压，这样即使电池电电压降低到6v本电路也能正常工作。其电路如图3-2-4所示。 图3-2-4前置放大与滤波电路3.2.3 波形变换电路IC1a输出的信号是叠加有噪声的脉动脉搏波（图1-1-2），此信号由滞回比较器（施密特触发器）IC1d转换成方波。该比较器的阈值可用P2调定在脉搏波的幅值范围之内，但对P2的调定要求并不严格，因为IC1d的输出信号经C4、C3的微分后总是将正、负相间的尖脉冲加到单稳多谐振荡器IC1c的反相输入端，不会造成很大的触发误差。当有输入信号时，IC1c在滞回比较器IC1d输入信号的每个下降沿到来时输出高电平，使C3通过R6充电。大约持续20ms之后，因C3充电电流减小而使IC1c同相输入端的电位降低到低于反相输入端的电位（尖脉冲已过去很久），于是IC1c改变状态并再次输出低电平。这长的脉冲是与脉搏同步的，并由红色发光二极管D14的闪亮指示出来。同时，该脉冲电平通过R17送到单片机P3.3脚，进行对心率的计算和显示。9v电源电压由R3、R4分压成4.5v，再经IC1b缓冲后用作IC1 a、IC1d、 IC1c的参考电压，这样即使电池电电压降低到6v本电路也能正常工作。其电路如图3-2-5所示。图3-2-5 波形变换电路3.3 显示电路选用七段共阴数码管作为数据显示器，单片机AT89C2051的P1口接数码管的a、b、c、d、e、f、g和dp引脚。P3.0P3.2接ULN2003的1B、2B和3B，ULN2003的1C、2C和3C分别接三个数码管的公共端子COM2 、COM1 和COM0。当P3.0P3.2其中有一个变为高电平时，分别驱动7SEG2、7SEG1 和7SEG0，此时从P1口输出经过单片机处理后的数据数码管就可以显示心率的次数（每分钟送一次数据），其中P3.3是对脉冲计数的。8个220的排阻接电源后接P1口，作限流用，维持数码管正常显示，当单片机复位时，使P1口电平全为高电平。其电路如图3-3-1所示。图3-3-1 显示电路3.4报警电路 报警电路是对每次脉冲的到来均响铃，与脉搏同步。这样，就可以通过声光的形式形象地把脉搏的快慢显示出来，其电路如图3-4-1所示。图3-4-1 报警电路3.5时钟电路单片机内部时候方式如图3-5-1所示。在单片机内部有一振荡器电路，只要在单片机的XTAL1和XTAL2引脚外接石英晶体，就构成了自激振荡器并在单片机内部产生时钟脉冲信号。图中电容C9和C10的作用是稳定频率和快速起振，电容值在530pF。晶振X1的振荡频率范围为1.212MHz，典型值为12 MHz和6 MHz.图3-5-1 时钟电路3.6复位电路按键与上电复位电路如图3-6-1所示。开机瞬间RESET引脚获得高电平，随着电容 C11的充电，RESET引脚的高电平将逐渐下降。RESET引脚的高电平只要保持足够的时间（2个机器周期），单片机就可以进行复位操作；另外，在单片机运行期间，还可以利用按键来完成复位操作。图3-6-1 复位电路3.7单片机处理与显示部分电路图如图3-7-1，本部分运用89C2051单片机作核心元件，在这里运用了单片机能更快更准确地对数据进行运算，而且可根据实际情况进行编程，所用外围元件少，轻巧省电，故障率低。ULN2003为内置达林顿管集成电路，作动态扫描时的选通驱动用。220排阻作限流用，维持数码管正常显示。本部分电源电压为5V稳压。来自传感和整形输出电路的脉冲电平输入单片机89C2051的P3.3脚，单片机设为负跳变中断触发模式，故每次脉冲下降沿到达时触发单片机产生中断并进行计时；当下一次脉冲的下降沿到达时，单片机对两次脉冲间的时间进行运算得出心率，通过P1口把结果送到数码管显示出来。同时，对每次脉冲的到来均响铃，与脉搏同步。这样，就可以通过声光的形式形象地把脉搏的快慢显示出来。为避免干扰的影响，单片机对两个脉冲之间的时间间距进行检测，若发现有干扰则忽略该干扰而不显示。本测试器的有效测量范围为50次199次/分钟。图3-7-1 单片机处理与显示部分电路图4 软件设计4.1程序流程图4.1.1主程序流程图 4.1.2定时器中断程序流程图图4-4-1 主程序流程图 图4-4-2 定时中断程序流程图4.1.3 INT中断程序流程图 4.1.4显示程序流程图图4-4-3 INT中断程序流程图 图4-4-4 显示程序流程图4.2源程序 ORG 0000H AJMP START ORG 000BH AJMP TIME ORG 0013H AJMP INT_1TIMES EQU 30HMORE_200 EQU 31HMORE200 EQU 32HBELL EQU 33H ORG 00e0HSTART: MOV SP,#50H MOV TMOD,#01H;方式1，定时时间为6ms MOV TH0,#0E8H; MOV TL0,#90H SETB EA ;开总中断 SETB IT1;INT1 负跳变触发 SETB EX1;INT1 开中断 SETB ET0 ;SETB IE1 MOV R1,#00H MOV R2,#00H MOV R3,#02H SETB TR0 ORG 0100H MOV R3,#06H AJMP MAIN ORG 0110HMAIN: ACALL SHOW AJMP MAIN;* ORG 0200HTIME: CLR ex1 PUSH ACC PUSH PSW MOV TH0,#0E8H MOV TL0,#90H CLR C MOV A,R6 ;大于200既是小于50次每分;清r6，并置标志位MORE200为1 SUBB A,#0C8H JC T_TEMP1 MOV R6,#00H MOV MORE200,#01H AJMP T_TEMP2T_TEMP1: INC R6T_TEMP2: CLR C MOV A,BELL SUBB A,#01H;响铃18ms，与脉搏脉冲同步 JC STOPBELL ;SETB P3.7 AJMP SSSSSTOPBELL: ;CLR P3.7SSSS: SETB ex1 POP PSW POP ACC RETI NOP NOP LJMP 0100H;* ORG 0300HINT_1: PUSH PSW ; PUSH ACC MOV TIMES,R6;保存r6的时间 MOV R6,#00H;令r6为0，重新计时 MOV MORE_200,MORE200;保存MORE200标志位 MOV MORE200,#00H;令MORE200为0 MOV BELL,#03H;设响铃 MOV A,MORE_200 ;判断是否大于200，即是否小于50次 CJNE A,#00H,TOKEEP;大于200，即小于50次，;保持原来的显示 CLR C MOV A,TIMES;判断是否小于50，即是否大于200次 SUBB A,#32H; JC TOKEEP ;小于50，保持原来的显示 AJMP CALCULATETOKEEP: POP ACC POP PSW AJMP C_TORETICALCULATE:push psw PSW ACC MOV R1,times MOV R7,#0H MOV R3,#27H MOV R2,#10H CLR Csubs1: MOV A,R2 SUBB A,R1 MOV R2,A MOV A,R3 SUBB A,#0H INC R7 MOV R3,A CLR C CJNE R3,#0H,subs1subs2: MOV A,R2 SUBB A,R1 MOV R2,A JC chg2dec INC R7 AJMP subs2CHG2DEC: MOV A,R7; 由r7 的值转换成3位十进制，其中r1为百位,r2为十位，r3为个位 CLR C SUBB A,#0C7H;大于200次每分钟，令为199次每分钟 JC CHANGE MOV R7,#0C7HCHANGE: MOV A,R7 MOV R0,A MOV A,R0 CLR C ;减100 SUBB A,#64H JC NO_MINUS; MOV R0,A;大于等于100次每分钟，令百位为1 MOV R1,#01H AJMP divvNO_MINUS: MOV R1,#00Hdivv: MOV A,R0; MOV B,#0AH; DIV AB MOV R2,A MOV R3,B POP ACC POP PSW AJMP C_TORETIC_TORETI: RETI NOP NOP LJMP 0100H;*;SHOW的子程序 ;r1为百位数字，r2为十位，r3为个位 ;p3。0接百位，p3。1接十位，p3。2接个位 ORG 0500HSHOW: MOV A,R1; CJNE A,#00H,S_TEMP1; ACALL DELAY; AJMP S_TEMP2;S_TEMP1: MOV A,R1;显示百位 MOV DPTR,#TABLE MOVC A,A+DPTR MOV P1,A SETB P3.0 ACALL DELAY CLR P3.0S_TEMP2: MOV A,R1 CJNE A,#00H,S_TEMP2_2;百位为0，十位又为0时，;十位的0消隐 MOV A,R2 CJNE A,#00H,S_TEMP2_2; ACALL DELAY AJMP S_TEMP3S_TEMP2_2: MOV A,R2 ;显示十位 MOV DPTR,#TABLE MOVC A,A+DPTR MOV P1,A SETB P3.1 ACALL DELAY CLR P3.1S_TEMP3: MOV A,R3;显示个位 MOV DPTR,#TABLE MOVC A,A+DPTR MOV P1,A SETB P3.2 ACALL DELAY CLR P3.2 RET NOP NOP LJMP 0100H;* ORG 0600HDELAY: PUSH ACC MOV A,R5 PUSH ACC MOV R5,#0AH ;延时1msD2_1: MOV R4,#32HD2_2: DJNZ R4,D2_2 DJNZ R5,D2_1 POP ACC MOV R5,A POP ACC RET NOP NOP LJMP 0100H;* ORG 0700HTABLE: DB 03FH ;LED 显示码0，共阴 DB 006H ; 1 DB 05BH ; 2 DB 04FH ; 3 DB 066H ; 4 DB 06DH ; 5 DB 07DH ; 6 DB 007H ; 7 DB 07FH ; 8 DB 06FH ; 9 NOP NOP LJMP 0100H END5调试与仿真5.1 Proteus7.12与KeiluVision3（1）打开Keil uVision3,新建Keil项目，选择AT89C2051单片机作为CPU，新建汇编源程序文件，编写程序，并将其导入”Source Group1”中。“ptions for Target ”对话窗口中，选择“Output”选择卡中的“Create HEX File ”选项和“Debug”选项卡中的“Debug Use:Proteus VSW Simulator”选项。编译汇编源程序，改正程序中的错误。（2）在Proteus ISIS 中，选中AT89C2051并单击鼠标左键，打开“Edit Component对话窗口，设置单片机晶振频率为12MHZ，在此窗口中的“Program File”栏中，选择先前用Keil生成的. HEX文件。在Proteus ISIS的菜单栏中选择“File”/“Save Design”保存设计。打开“Debug”/Use Remote Debug Monitor以支持与Keil的联合调试。（3）在Keil的菜单栏中选择“Debug”/Start/Stop Debug Session选项，进入程序调试环境。按“F5”键，顺序运行