CN111281396B - 胸腹表面呼吸运动信号超分辨方法 - Google Patents

胸腹表面呼吸运动信号超分辨方法 Download PDF

Info

Publication number
CN111281396B
CN111281396B CN202010075888.0A CN202010075888A CN111281396B CN 111281396 B CN111281396 B CN 111281396B CN 202010075888 A CN202010075888 A CN 202010075888A CN 111281396 B CN111281396 B CN 111281396B
Authority
CN
China
Prior art keywords
resistor
phase
output
square wave
super
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010075888.0A
Other languages
English (en)
Other versions
CN111281396A (zh
Inventor
赵烟桥
陈睿
胡亚欣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN202010075888.0A priority Critical patent/CN111281396B/zh
Publication of CN111281396A publication Critical patent/CN111281396A/zh
Application granted granted Critical
Publication of CN111281396B publication Critical patent/CN111281396B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/113—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing
    • A61B5/1135—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing by monitoring thoracic expansion
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Physiology (AREA)
  • Surgery (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Analogue/Digital Conversion (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明胸腹表面呼吸运动信号超分辨方法属于精密仪器和胸腹部放射治疗技术领域；该方法首先通过信号转换模块，将呼吸运动信号转换成周期信号，再通过相位差多输出模块，形成相位差为π/2的等相位差三输出，再通过电阻链多相位生成模块形成相位差为π/10的等相位差十输出，再通过多相位正弦方波转换模块形成方波，并通过多相位融合逻辑门模块形成准相位超分辨方波输出和错相位超分辨方波输出，最后二次超分辨模块利用准相位超分辨方波输出和错相位超分辨方波输出的相位差进行进一步超分辨，将一个周期的呼吸运动信号变成多个周期的方波信号。

Description

胸腹表面呼吸运动信号超分辨方法
技术领域
本发明胸腹表面呼吸运动信号超分辨方法属于精密仪器和胸腹部放射治疗技术领域。
背景技术
在胸腹部肿瘤放射治疗过程中，呼吸使得肿瘤区域逸出靶区或正常组织进入靶区，不仅降低了放疗效果，而且容易产生并发症。
为了解决呼吸运动给放射治疗带来的负面影响，临床上最早采用了呼吸保持和呼吸门控等方法，通过对呼吸的干预实现肿瘤位置的控制，虽然取得了一定的效果，但是患者耐受性差；为了提高患者耐受性，不干预患者呼吸，又出现了跟踪系统，通过监测肿瘤区域实现对肿瘤位置的跟踪，然而，由于这种方法属于滞后补偿方法，因此总会出现“慢半拍”的问题；为了提高跟踪精度，学者们采用了预测手段，期初，将呼吸运动视为一个呼吸周期的简单重复，然而，由于呼吸运动本身具有无明显规律的准周期特性，因此预测精度并不高，随着时间的推移，误差会越来越大；后来，学者们根据呼吸运动的历史规律，通过拟合，形成以时间为变量的呼吸运动模型，进而对未来呼吸运动进行预测，这种方法的效果较好；本课题组的研究生樊琪和史领采用高斯过程回归方法对呼吸运动进行预测，并以均值和方差的形式给出预测结果，为呼吸运动提供了一种全新的预测手段。
对于呼吸运动的宏观特性由幅值和频率来决定，由于幅值和频率之间又存在着一定的对应关系，因此只需要准确获得呼吸频率，即可实现呼吸运动预测，然而，对于呼吸频率的预测，至少要采集一个呼吸周期，这样就给呼吸频率快速分析带来了一定的困难。如果能够对呼吸运动信号进行超分辨，就可以实现在少于一个呼吸运动周期的时间获得呼吸频率，进而快速实现对呼吸运动进行预测与分析。然而，在本领域，还没有发现对呼吸运动信号进行超分辨的技术手段。
发明内容
为了实现对呼吸运动信号进行超分辨，本发明公开了一种胸腹表面呼吸运动信号超分辨电路与方法，能够将一个周期的呼吸运动信号变成多个周期的方波信号，这个结果所能够带来的进一步技术优势在于，在多个呼吸周期过程中，通过判断方波信号的频率变化即可判断呼吸频率的变化，更重要的是，由于超分辨方波的周期远小于呼吸运动周期，因此能够在更短时间内，即小于一个呼吸周期的时间范围内，判断出呼吸频率的变化。
本发明的目的是这样实现的：
胸腹表面呼吸运动信号超分辨电路，包括信号转换模块，相位差多输出模块，电阻链多相位生成模块，多相位正弦方波转换模块，多相位融合逻辑门模块和二次超分辨模块；
所述信号转换模块的输入为一个周期的呼吸运动信号f(t)，输出为：
Figure BDA0002378482180000021
其中：t0为呼吸运动信号f(t)开始采集的时间；ω为呼吸运动信号f(t)的角频率；T为呼吸运动信号f(t)的周期；
所述相位差多输出模块包括运算放大器U1-1和运算放大器U1-2，所述运算放大器U1-1的反相输入端通过电容C1-1连接信号转换模块的输出端，通过电阻R1-1连接运算放大器U1-1的输出端，运算放大器U1-1的同相输入端连接地；所述运算放大器U1-2的反相输入端通过电容C1-2连接运算放大器U1-1的输出端，通过电阻R1-2连接运算放大器U1-2的输出端，运算放大器U1-2的同相输入端连接地；信号转换模块的输出端作为所述相位差多输出模块的第一输出，运算放大器U1-1的输出端作为所述相位差多输出模块的第二输出，运算放大器U1-2的输出端作为所述相位差多输出模块的第三输出；
所述电阻链多相位生成模块包括18个电阻；相位差多输出模块的第一输出和相位差多输出模块的第二输出之间分别通过电阻R2-2和电阻R2-6的串联结构连接，电阻R2-2和电阻R2-6的阻值比为9/28，通过电阻R2-3和电阻R2-7的串联结构连接，电阻R2-3和电阻R2-7的阻值比为8/11，通过电阻R2-4和电阻R2-8的串联结构连接，电阻R2-4和电阻R2-8的阻值比为11/8，通过电阻R2-5和电阻R2-9的串联结构连接，电阻R2-5和电阻R2-9的阻值比为28/9；相位差多输出模块的第二输出和相位差多输出模块的第三输出之间分别通过电阻R2-11和电阻R2-15的串联结构连接，电阻R2-11和电阻R2-15的阻值比为9/28，通过电阻R2-12和电阻R2-16的串联结构连接，电阻R2-12和电阻R2-16的阻值比为8/11，通过电阻R2-13和电阻R2-17的串联结构连接，电阻R2-13和电阻R2-17的阻值比为11/8，通过电阻R2-14和电阻R2-18的串联结构连接，电阻R2-14和电阻R2-18的阻值比为28/9；所述相位差多输出模块的第一输出经过电阻R2-1后作为第一相位输出，电阻R2-2和电阻R2-6之间的抽头作为第二相位输出，电阻R2-3和电阻R2-7之间的抽头作为第三相位输出，电阻R2-4和电阻R2-8之间的抽头作为第四相位输出，电阻R2-5和电阻R2-9之间的抽头作为第五相位输出，所述相位差多输出模块的第二输出经过电阻R2-10后作为第六相位输出，电阻R2-11和电阻R2-15之间的抽头作为第七相位输出，电阻R2-12和电阻R2-16之间的抽头作为第八相位输出，电阻R2-13和电阻R2-17之间的抽头作为第九相位输出，电阻R2-14和电阻R2-18之间的抽头作为第十相位输出；
所述多相位正弦方波转换模块包括10个运算放大器，运算放大器U2-1的反相输入端连接电阻链多相位生成模块的第一相位输出；运算放大器U2-1的同相输入端连接地，运算放大器U2-1的输出端为第一方波输出；运算放大器U2-2的反相输入端连接电阻链多相位生成模块的第二相位输出；运算放大器U2-2的同相输入端连接地，运算放大器U2-2的输出端为第二方波输出；运算放大器U2-3的反相输入端连接电阻链多相位生成模块的第三相位输出；运算放大器U2-3的同相输入端连接地，运算放大器U2-3的输出端为第三方波输出；运算放大器U2-4的反相输入端连接电阻链多相位生成模块的第四相位输出；运算放大器U2-4的同相输入端连接地，运算放大器U2-4的输出端为第四方波输出；运算放大器U2-5的反相输入端连接电阻链多相位生成模块的第五相位输出；运算放大器U2-5的同相输入端连接地，运算放大器U2-5的输出端为第五方波输出；运算放大器U2-6的反相输入端连接电阻链多相位生成模块的第六相位输出；运算放大器U2-6的同相输入端连接地，运算放大器U2-6的输出端为第六方波输出；运算放大器U2-7的反相输入端连接电阻链多相位生成模块的第七相位输出；运算放大器U2-7的同相输入端连接地，运算放大器U2-7的输出端为第七方波输出；运算放大器U2-8的反相输入端连接电阻链多相位生成模块的第八相位输出；运算放大器U2-8的同相输入端连接地，运算放大器U2-8的输出端为第八方波输出；运算放大器U2-9的反相输入端连接电阻链多相位生成模块的第九相位输出；运算放大器U2-9的同相输入端连接地，运算放大器U2-9的输出端为第九方波输出；运算放大器U2-10的反相输入端连接电阻链多相位生成模块的第十相位输出；运算放大器U2-10的同相输入端连接地，运算放大器U2-10的输出端为第十方波输出；
所述多相位融合逻辑门模块包括八个异或门，异或门U3-1的输入为第三方波输出和第七方波输出，异或门U3-1的输出为第一异或输出；异或门U3-2的输入为第二方波输出和第十方波输出，异或门U3-2的输出为第二异或输出；异或门U3-3的输入为第一方波输出和第一异或输出，异或门U3-3的输出为第三异或输出；异或门U3-4的输入为第二异或输出和第六方波输出，异或门U3-4的输出为第四异或输出；异或门U3-5的输入为第四异或输出和第八异或输出，异或门U3-5的输出为第五异或输出；异或门U3-6的输入为第三异或输出和第七异或输出，异或门U3-6的输出为第六异或输出；异或门U3-7的输入为第五方波输出和第九方波输出，异或门U3-7的输出为第七异或输出；异或门U3-8的输入为第四方波输出和第八方波输出，异或门U3-8的输出为第八异或输出；其中，第五异或输出作为错相位超分辨方波输出，第六异或输出作为准相位超分辨方波输出；
所述二次超分辨模块包括准相位超分辨方波变形电路，错相位超分辨方波变形电路，第一八输入与或非门，第二八输入与或非门和最终与门；
所述准相位超分辨方波变形电路的输入为准相位超分辨方波输出，准相位超分辨方波变形电路的第一输出为准相位超分辨方波输出，准相位超分辨方波变形电路的第二输出为准相位超分辨方波输出非，准相位超分辨方波变形电路的第三输出为准相位超分辨方波输出与准相位超分辨方波输出非且延时的与运算，准相位超分辨方波变形电路的第四输出为准相位超分辨方波输出非与准相位超分辨方波输出延时的与运算；
所述错相位超分辨方波变形电路的输入为错相位超分辨方波输出，错相位超分辨方波变形电路的第一输出为错相位超分辨方波输出，错相位超分辨方波变形电路的第二输出为错相位超分辨方波输出非，错相位超分辨方波变形电路的第三输出为错相位超分辨方波输出与错相位超分辨方波输出非且延时的与运算，错相位超分辨方波变形电路的第四输出为错相位超分辨方波输出非与错相位超分辨方波输出延时的与运算；
所述第一八输入与或非门和第二八输入与或非门的输入均为准相位超分辨方波输出、准相位超分辨方波输出非、准相位超分辨方波输出与准相位超分辨方波输出非且延时的与运算、准相位超分辨方波输出非与准相位超分辨方波输出延时的与运算、错相位超分辨方波输出、错相位超分辨方波输出非、错相位超分辨方波输出与错相位超分辨方波输出非且延时的与运算和错相位超分辨方波输出非与错相位超分辨方波输出延时的与运算，其中，第一八输入与或非门按照如下逻辑进行运算：
准相位超分辨方波变形电路的第一输出&错相位超分辨方波变形电路的第三输出
或非
准相位超分辨方波变形电路的第二输出&错相位超分辨方波变形电路的第四输出
或非
准相位超分辨方波变形电路的第三输出&错相位超分辨方波变形电路的第二输出
或非
准相位超分辨方波变形电路的第四输出&错相位超分辨方波变形电路的第一输出
第二八输入与或非门按照如下逻辑进行运算：
准相位超分辨方波变形电路的第一输出&错相位超分辨方波变形电路的第四输出
或非
准相位超分辨方波变形电路的第二输出&错相位超分辨方波变形电路的第三输出
或非
准相位超分辨方波变形电路的第三输出&错相位超分辨方波变形电路的第一输出
或非
准相位超分辨方波变形电路的第四输出&错相位超分辨方波变形电路的第二输出
其中，&表示与运算；
所述最终与门的输入为第一八输入与或非门的输出和第二八输入与或非门的输出，最终与门的输出为胸腹表面呼吸运动超分辨信号。
上述胸腹表面呼吸运动信号超分辨电路，所述信号转换模块，由通带频率包含ω的滤波器和放大器来实现。
上述胸腹表面呼吸运动信号超分辨电路，所述二次超分辨模块中，所述延时通过电容和电阻组成的延时电路实现。
胸腹表面呼吸运动信号超分辨方法，包括以下步骤：
步骤a、在信号转换模块中，利用通带频率包含ω的滤波器和放大器，实现输入为周期的呼吸运动信号f(t)，输出为：
Figure BDA0002378482180000051
其中，
滤波器选择频率为ω的分量通过；
放大器将幅值调整为
Figure BDA0002378482180000052
步骤b、在相位差多输出模块中，利用由电阻、电容和运算放大器组成的两个微分电路，分别实现对信号转换模块的输出余弦信号进行一次微分和两次微分，进而实现相位差为公差π/2的等差数列三输出；
步骤c、在电阻链多相位生成模块中，为了方便分析，将相位差多输出模块的三输出分别设定为sinα、sin(α-π/2)和sin(α-π)；
由于电阻R2-2和电阻R2-6的阻值比为9/28，因此电阻R2-2和电阻R2-6之间的抽头电压为：
Figure BDA0002378482180000061
由于电阻R2-3和电阻R2-7的阻值比为8/11，因此电阻R2-3和电阻R2-7之间的抽头电压为：
Figure BDA0002378482180000062
由于电阻R2-4和电阻R2-8的阻值比为11/8，因此电阻R2-4和电阻R2-8之间的抽头电压为：
Figure BDA0002378482180000063
由于电阻R2-5和电阻R2-9的阻值比为28/9，因此电阻R2-5和电阻R2-9之间的抽头电压为：
Figure BDA0002378482180000064
由于电阻R2-11和电阻R2-15的阻值比为9/28，因此电阻R2-11和电阻R2-15之间的抽头电压为：
Figure BDA0002378482180000065
由于电阻R2-12和电阻R2-16的阻值比为8/11，因此电阻R2-12和电阻R2-16之间的抽头电压为：
Figure BDA0002378482180000066
由于电阻R2-13和电阻R2-17的阻值比为11/8，因此电阻R2-13和电阻R2-17之间的抽头电压为：
Figure BDA0002378482180000067
由于电阻R2-14和电阻R2-18的阻值比为28/9，因此电阻R2-14和电阻R2-18之间的抽头电压为：
Figure BDA0002378482180000071
能够看出，电阻链多相位生成模块的输出为相位公差为π/10的等差数列十输出；
步骤d、在多相位正弦方波转换模块中，电阻链多相位生成模块的每一路输出都经过过零比较器，将幅值调整为运算放大器的饱和电压值，相位保留；
步骤e、在多相位融合逻辑门模块中，利用多相位正弦方波转换模块的十个方波输出和多相位融合逻辑门模块的六个异或输出进行逻辑运算，实现错相位超分辨方波输出和准相位超分辨方波输出；错相位超分辨方波输出将一个周期的呼吸运动信号变成五个周期的方波信号，相位相差错相位超分辨方波输出信号的四分之一周期，准相位超分辨方波输出将一个周期的呼吸运动信号变成五个周期的方波信号，相位相差准相位超分辨方波输出信号的零周期；
步骤f、在二次超分辨模块中，利用错相位超分辨方波输出信号与准相位超分辨方波输出信号之间相差四分之一周期的特性，采用非门、延时电路和与门，实现对错相位超分辨方波输出信号和准相位超分辨方波输出信号每一个上升沿和下降沿的采集，再通过两个与或非门和最终与门，实现对错相位超分辨方波输出信号和准相位超分辨方波输出信号的二次超分辨，最终实现将一个周期的呼吸运动信号变成二十个周期的方波信号。
有益效果：
第一、在本发明中，信号转换模块，相位差多输出模块，电阻链多相位生成模块，多相位正弦方波转换模块，多相位融合逻辑门模块和二次超分辨模块作为一个整体，缺一不可，共同实现了将一个周期的呼吸运动信号变成多个周期的方波信号，即实现了对胸腹表面呼吸运动信号的超分辨，这个结果所能够带来的进一步技术优势在于，在多个呼吸周期过程中，通过判断方波信号的频率变化即可判断呼吸频率的变化，更重要的是，由于超分辨方波的周期远小于呼吸运动周期，因此能够在更短时间内，即小于一个呼吸周期的时间范围内，判断出呼吸频率的变化。
第二、在二次超分辨模块中，由于将准相位超分辨方波变形电路和错相位超分辨方波变形电路设计成电路结构完全一致的形式，同时配合后续的两个与或非门和最终与门，因此即使多相位融合逻辑门模块的两个输出与二次超分辨模块的两个输入反接，也对胸腹表面呼吸运动超分辨信号的输出结果没有任何影响，即可以实现多相位融合逻辑门模块与二次超分辨模块无方向插接的效果，类似Type-C接口优于mini USB或micro USB的效果。
第三、在本发明中，同本项目团队同日申请的另外四项专利《一种胸腹表面呼吸运动信号错相位超分辨电路》、《一种胸腹表面呼吸运动信号错相位超分辨方法》、《一种胸腹表面呼吸运动信号准相位超分辨电路》、《一种胸腹表面呼吸运动信号准相位超分辨方法》相比，由于设计有二次超分辨模块，并同时利用了错相位超分辨方波输出和准相位超分辨方波输出，因此能够进行二次超分辨，细节程度翻两番，判断呼吸频率用时能够缩短为另外四项专利的四分之一。
附图说明
图1是本发明胸腹表面呼吸运动信号超分辨电路逻辑框图。
图2是本发明胸腹表面呼吸运动信号超分辨电路中相位差多输出模块电路图。
图3是本发明胸腹表面呼吸运动信号超分辨电路中电阻链多相位生成模块电路图。
图4是本发明胸腹表面呼吸运动信号超分辨电路中多相位正弦方波转换模块电路图。
图5是本发明胸腹表面呼吸运动信号超分辨电路中多相位融合逻辑门模块电路图。
图6是本发明胸腹表面呼吸运动信号超分辨电路中二次超分辨模块电路图。
图7是本发明胸腹表面呼吸运动信号超分辨电路各模块输出波形图。
具体实施方式
下面结合附图对本发明具体实施方式作进一步详细描述。
具体实施方式一
本实施方式为胸腹表面呼吸运动信号超分辨电路的实施方式。
胸腹表面呼吸运动信号超分辨电路，逻辑框图如图1所示，该胸腹表面呼吸运动信号超分辨电路包括信号转换模块，相位差多输出模块，电阻链多相位生成模块，多相位正弦方波转换模块，多相位融合逻辑门模块和二次超分辨模块；其中，相位差多输出模块电路图如图2所示，电阻链多相位生成模块电路图如图3所示，多相位正弦方波转换模块电路图如图4所示，多相位融合逻辑门模块电路图如图5所示，二次超分辨模块电路图如图6所示，本发明胸腹表面呼吸运动信号超分辨电路各模块输出波形图如图7所示；
所述信号转换模块的输入为一个周期的呼吸运动信号f(t)，输出为：
Figure BDA0002378482180000081
其中：t0为呼吸运动信号f(t)开始采集的时间；ω为呼吸运动信号f(t)的角频率；T为呼吸运动信号f(t)的周期；
所述信号转换模块，由通带频率包含ω的滤波器和放大器来实现；
所述相位差多输出模块包括运算放大器U1-1和运算放大器U1-2，所述运算放大器U1-1的反相输入端通过电容C1-1连接信号转换模块的输出端，通过电阻R1-1连接运算放大器U1-1的输出端，运算放大器U1-1的同相输入端连接地；所述运算放大器U1-2的反相输入端通过电容C1-2连接运算放大器U1-1的输出端，通过电阻R1-2连接运算放大器U1-2的输出端，运算放大器U1-2的同相输入端连接地；信号转换模块的输出端作为所述相位差多输出模块的第一输出，运算放大器U1-1的输出端作为所述相位差多输出模块的第二输出，运算放大器U1-2的输出端作为所述相位差多输出模块的第三输出；
所述电阻链多相位生成模块包括18个电阻；相位差多输出模块的第一输出和相位差多输出模块的第二输出之间分别通过电阻R2-2和电阻R2-6的串联结构连接，电阻R2-2和电阻R2-6的阻值比为9/28，通过电阻R2-3和电阻R2-7的串联结构连接，电阻R2-3和电阻R2-7的阻值比为8/11，通过电阻R2-4和电阻R2-8的串联结构连接，电阻R2-4和电阻R2-8的阻值比为11/8，通过电阻R2-5和电阻R2-9的串联结构连接，电阻R2-5和电阻R2-9的阻值比为28/9；相位差多输出模块的第二输出和相位差多输出模块的第三输出之间分别通过电阻R2-11和电阻R2-15的串联结构连接，电阻R2-11和电阻R2-15的阻值比为9/28，通过电阻R2-12和电阻R2-16的串联结构连接，电阻R2-12和电阻R2-16的阻值比为8/11，通过电阻R2-13和电阻R2-17的串联结构连接，电阻R2-13和电阻R2-17的阻值比为11/8，通过电阻R2-14和电阻R2-18的串联结构连接，电阻R2-14和电阻R2-18的阻值比为28/9；所述相位差多输出模块的第一输出经过电阻R2-1后作为第一相位输出，电阻R2-2和电阻R2-6之间的抽头作为第二相位输出，电阻R2-3和电阻R2-7之间的抽头作为第三相位输出，电阻R2-4和电阻R2-8之间的抽头作为第四相位输出，电阻R2-5和电阻R2-9之间的抽头作为第五相位输出，所述相位差多输出模块的第二输出经过电阻R2-10后作为第六相位输出，电阻R2-11和电阻R2-15之间的抽头作为第七相位输出，电阻R2-12和电阻R2-16之间的抽头作为第八相位输出，电阻R2-13和电阻R2-17之间的抽头作为第九相位输出，电阻R2-14和电阻R2-18之间的抽头作为第十相位输出；
所述多相位正弦方波转换模块包括10个运算放大器，运算放大器U2-1的反相输入端连接电阻链多相位生成模块的第一相位输出；运算放大器U2-1的同相输入端连接地，运算放大器U2-1的输出端为第一方波输出；运算放大器U2-2的反相输入端连接电阻链多相位生成模块的第二相位输出；运算放大器U2-2的同相输入端连接地，运算放大器U2-2的输出端为第二方波输出；运算放大器U2-3的反相输入端连接电阻链多相位生成模块的第三相位输出；运算放大器U2-3的同相输入端连接地，运算放大器U2-3的输出端为第三方波输出；运算放大器U2-4的反相输入端连接电阻链多相位生成模块的第四相位输出；运算放大器U2-4的同相输入端连接地，运算放大器U2-4的输出端为第四方波输出；运算放大器U2-5的反相输入端连接电阻链多相位生成模块的第五相位输出；运算放大器U2-5的同相输入端连接地，运算放大器U2-5的输出端为第五方波输出；运算放大器U2-6的反相输入端连接电阻链多相位生成模块的第六相位输出；运算放大器U2-6的同相输入端连接地，运算放大器U2-6的输出端为第六方波输出；运算放大器U2-7的反相输入端连接电阻链多相位生成模块的第七相位输出；运算放大器U2-7的同相输入端连接地，运算放大器U2-7的输出端为第七方波输出；运算放大器U2-8的反相输入端连接电阻链多相位生成模块的第八相位输出；运算放大器U2-8的同相输入端连接地，运算放大器U2-8的输出端为第八方波输出；运算放大器U2-9的反相输入端连接电阻链多相位生成模块的第九相位输出；运算放大器U2-9的同相输入端连接地，运算放大器U2-9的输出端为第九方波输出；运算放大器U2-10的反相输入端连接电阻链多相位生成模块的第十相位输出；运算放大器U2-10的同相输入端连接地，运算放大器U2-10的输出端为第十方波输出；
所述多相位融合逻辑门模块包括八个异或门，异或门U3-1的输入为第三方波输出和第七方波输出，异或门U3-1的输出为第一异或输出；异或门U3-2的输入为第二方波输出和第十方波输出，异或门U3-2的输出为第二异或输出；异或门U3-3的输入为第一方波输出和第一异或输出，异或门U3-3的输出为第三异或输出；异或门U3-4的输入为第二异或输出和第六方波输出，异或门U3-4的输出为第四异或输出；异或门U3-5的输入为第四异或输出和第八异或输出，异或门U3-5的输出为第五异或输出；异或门U3-6的输入为第三异或输出和第七异或输出，异或门U3-6的输出为第六异或输出；异或门U3-7的输入为第五方波输出和第九方波输出，异或门U3-7的输出为第七异或输出；异或门U3-8的输入为第四方波输出和第八方波输出，异或门U3-8的输出为第八异或输出；其中，第五异或输出作为错相位超分辨方波输出，第六异或输出作为准相位超分辨方波输出；
所述二次超分辨模块包括准相位超分辨方波变形电路，错相位超分辨方波变形电路，第一八输入与或非门，第二八输入与或非门和最终与门；
所述准相位超分辨方波变形电路的输入为准相位超分辨方波输出，准相位超分辨方波变形电路的第一输出为准相位超分辨方波输出，准相位超分辨方波变形电路的第二输出为准相位超分辨方波输出非，准相位超分辨方波变形电路的第三输出为准相位超分辨方波输出与准相位超分辨方波输出非且延时的与运算，准相位超分辨方波变形电路的第四输出为准相位超分辨方波输出非与准相位超分辨方波输出延时的与运算；
所述错相位超分辨方波变形电路的输入为错相位超分辨方波输出，错相位超分辨方波变形电路的第一输出为错相位超分辨方波输出，错相位超分辨方波变形电路的第二输出为错相位超分辨方波输出非，错相位超分辨方波变形电路的第三输出为错相位超分辨方波输出与错相位超分辨方波输出非且延时的与运算，错相位超分辨方波变形电路的第四输出为错相位超分辨方波输出非与错相位超分辨方波输出延时的与运算；
所述延时通过电容和电阻组成的延时电路实现；
所述第一八输入与或非门和第二八输入与或非门的输入均为准相位超分辨方波输出、准相位超分辨方波输出非、准相位超分辨方波输出与准相位超分辨方波输出非且延时的与运算、准相位超分辨方波输出非与准相位超分辨方波输出延时的与运算、错相位超分辨方波输出、错相位超分辨方波输出非、错相位超分辨方波输出与错相位超分辨方波输出非且延时的与运算和错相位超分辨方波输出非与错相位超分辨方波输出延时的与运算，其中，第一八输入与或非门按照如下逻辑进行运算：
准相位超分辨方波变形电路的第一输出&错相位超分辨方波变形电路的第三输出
或非
准相位超分辨方波变形电路的第二输出&错相位超分辨方波变形电路的第四输出
或非
准相位超分辨方波变形电路的第三输出&错相位超分辨方波变形电路的第二输出
或非
准相位超分辨方波变形电路的第四输出&错相位超分辨方波变形电路的第一输出
第二八输入与或非门按照如下逻辑进行运算：
准相位超分辨方波变形电路的第一输出&错相位超分辨方波变形电路的第四输出
或非
准相位超分辨方波变形电路的第二输出&错相位超分辨方波变形电路的第三输出
或非
准相位超分辨方波变形电路的第三输出&错相位超分辨方波变形电路的第一输出
或非
准相位超分辨方波变形电路的第四输出&错相位超分辨方波变形电路的第二输出
其中，&表示与运算；
所述最终与门的输入为第一八输入与或非门的输出和第二八输入与或非门的输出，最终与门的输出为胸腹表面呼吸运动超分辨信号。
具体实施方式二
本实施方式为胸腹表面呼吸运动信号超分辨方法的实施方式。
胸腹表面呼吸运动信号超分辨方法，其特征在于，包括以下步骤：
步骤a、在信号转换模块中，利用通带频率包含ω的滤波器和放大器，实现输入为周期的呼吸运动信号f(t)，输出为：
Figure BDA0002378482180000121
其中，
滤波器选择频率为ω的分量通过；
放大器将幅值调整为
Figure BDA0002378482180000122
步骤b、在相位差多输出模块中，利用由电阻、电容和运算放大器组成的两个微分电路，分别实现对信号转换模块的输出余弦信号进行一次微分和两次微分，进而实现相位差为公差π/2的等差数列三输出；
步骤c、在电阻链多相位生成模块中，为了方便分析，将相位差多输出模块的三输出分别设定为sinα、sin(α-π/2)和sin(α-π)；
由于电阻R2-2和电阻R2-6的阻值比为9/28，因此电阻R2-2和电阻R2-6之间的抽头电压为：
Figure BDA0002378482180000123
由于电阻R2-3和电阻R2-7的阻值比为8/11，因此电阻R2-3和电阻R2-7之间的抽头电压为：
Figure BDA0002378482180000131
由于电阻R2-4和电阻R2-8的阻值比为11/8，因此电阻R2-4和电阻R2-8之间的抽头电压为：
Figure BDA0002378482180000132
由于电阻R2-5和电阻R2-9的阻值比为28/9，因此电阻R2-5和电阻R2-9之间的抽头电压为：
Figure BDA0002378482180000133
由于电阻R2-11和电阻R2-15的阻值比为9/28，因此电阻R2-11和电阻R2-15之间的抽头电压为：
Figure BDA0002378482180000134
由于电阻R2-12和电阻R2-16的阻值比为8/11，因此电阻R2-12和电阻R2-16之间的抽头电压为：
Figure BDA0002378482180000135
由于电阻R2-13和电阻R2-17的阻值比为11/8，因此电阻R2-13和电阻R2-17之间的抽头电压为：
Figure BDA0002378482180000136
由于电阻R2-14和电阻R2-18的阻值比为28/9，因此电阻R2-14和电阻R2-18之间的抽头电压为：
Figure BDA0002378482180000137
能够看出，电阻链多相位生成模块的输出为相位公差为π/10的等差数列十输出；
步骤d、在多相位正弦方波转换模块中，电阻链多相位生成模块的每一路输出都经过过零比较器，将幅值调整为运算放大器的饱和电压值，相位保留；
步骤e、在多相位融合逻辑门模块中，利用多相位正弦方波转换模块的十个方波输出和多相位融合逻辑门模块的六个异或输出进行逻辑运算，实现错相位超分辨方波输出和准相位超分辨方波输出；错相位超分辨方波输出将一个周期的呼吸运动信号变成五个周期的方波信号，相位相差错相位超分辨方波输出信号的四分之一周期，准相位超分辨方波输出将一个周期的呼吸运动信号变成五个周期的方波信号，相位相差准相位超分辨方波输出信号的零周期；
步骤f、在二次超分辨模块中，利用错相位超分辨方波输出信号与准相位超分辨方波输出信号之间相差四分之一周期的特性，采用非门、延时电路和与门，实现对错相位超分辨方波输出信号和准相位超分辨方波输出信号每一个上升沿和下降沿的采集，再通过两个与或非门和最终与门，实现对错相位超分辨方波输出信号和准相位超分辨方波输出信号的二次超分辨，最终实现将一个周期的呼吸运动信号变成二十个周期的方波信号。

Claims (1)

1.胸腹表面呼吸运动信号超分辨方法，其特征在于，包括以下步骤：
步骤a、在信号转换模块中，利用通带频率包含ω的滤波器和放大器，实现输入为周期的呼吸运动信号f(t)，输出为：
Figure FDA0003761846250000011
其中，t0为呼吸运动信号f(t)开始采集的时间；ω为呼吸运动信号f(t)的角频率；T为呼吸运动信号f(t)的周期；
滤波器选择频率为ω的分量通过；
放大器将幅值调整为
Figure FDA0003761846250000012
步骤b、在相位差多输出模块中，利用由电阻、电容和运算放大器组成的两个微分电路，分别实现对信号转换模块的输出余弦信号进行一次微分和两次微分，进而实现相位差为公差π/2的等差数列三输出；
步骤c、在电阻链多相位生成模块中，为了方便分析，将相位差多输出模块的三输出分别设定为sinα、sin(α-π/2)和sin(α-π)；
由于电阻R2-2和电阻R2-6的阻值比为9/28，因此电阻R2-2和电阻R2-6之间的抽头电压为：
Figure FDA0003761846250000013
由于电阻R2-3和电阻R2-7的阻值比为8/11，因此电阻R2-3和电阻R2-7之间的抽头电压为：
Figure FDA0003761846250000014
由于电阻R2-4和电阻R2-8的阻值比为11/8，因此电阻R2-4和电阻R2-8之间的抽头电压为：
Figure FDA0003761846250000021
由于电阻R2-5和电阻R2-9的阻值比为28/9，因此电阻R2-5和电阻R2-9之间的抽头电压为：
Figure FDA0003761846250000022
由于电阻R2-11和电阻R2-15的阻值比为9/28，因此电阻R2-11和电阻R2-15之间的抽头电压为：
Figure FDA0003761846250000023
由于电阻R2-12和电阻R2-16的阻值比为8/11，因此电阻R2-12和电阻R2-16之间的抽头电压为：
Figure FDA0003761846250000024
由于电阻R2-13和电阻R2-17的阻值比为11/8，因此电阻R2-13和电阻R2-17之间的抽头电压为：
Figure FDA0003761846250000025
由于电阻R2-14和电阻R2-18的阻值比为28/9，因此电阻R2-14和电阻R2-18之间的抽头电压为：
Figure FDA0003761846250000026
能够看出，电阻链多相位生成模块的输出为相位公差为π/10的等差数列十输出；
步骤d、在多相位正弦方波转换模块中，电阻链多相位生成模块的每一路输出都经过过零比较器，将幅值调整为运算放大器的饱和电压值，相位保留；
步骤e、在多相位融合逻辑门模块中，利用多相位正弦方波转换模块的十个方波输出和多相位融合逻辑门模块的六个异或输出进行逻辑运算，实现错相位超分辨方波输出和准相位超分辨方波输出；错相位超分辨方波输出将一个周期的呼吸运动信号变成五个周期的方波信号，相位相差错相位超分辨方波输出信号的四分之一周期，准相位超分辨方波输出将一个周期的呼吸运动信号变成五个周期的方波信号，相位相差准相位超分辨方波输出信号的零周期；
步骤f、在二次超分辨模块中，利用错相位超分辨方波输出信号与准相位超分辨方波输出信号之间相差四分之一周期的特性，采用非门、延时电路和与门，实现对错相位超分辨方波输出信号和准相位超分辨方波输出信号每一个上升沿和下降沿的采集，再通过两个与或非门和最终与门，实现对错相位超分辨方波输出信号和准相位超分辨方波输出信号的二次超分辨，最终实现将一个周期的呼吸运动信号变成二十个周期的方波信号。
CN202010075888.0A 2020-01-22 2020-01-22 胸腹表面呼吸运动信号超分辨方法 Active CN111281396B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010075888.0A CN111281396B (zh) 2020-01-22 2020-01-22 胸腹表面呼吸运动信号超分辨方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010075888.0A CN111281396B (zh) 2020-01-22 2020-01-22 胸腹表面呼吸运动信号超分辨方法

Publications (2)

Publication Number Publication Date
CN111281396A CN111281396A (zh) 2020-06-16
CN111281396B true CN111281396B (zh) 2022-09-06

Family

ID=71023293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010075888.0A Active CN111281396B (zh) 2020-01-22 2020-01-22 胸腹表面呼吸运动信号超分辨方法

Country Status (1)

Country Link
CN (1) CN111281396B (zh)

Family Cites Families (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1384628A (en) * 1972-03-17 1975-02-19 United Aircraft Corp Aircraft navigation computer
US3821652A (en) * 1972-08-14 1974-06-28 Cincinnati Milacron Inc Square to sine converter
DE2741344A1 (de) * 1977-09-14 1979-03-15 Helge Dr Med Ruettgers Anordnung zur automatischen steuerung einer infusions- oder transfusionspumpe bei medizinischen anwendungen
DD229266A3 (de) * 1983-04-08 1985-10-30 Medizin Labortechnik Veb K Schaltungsanordnung zur messung von atemfrequenz und atemzeitverhaeltnis
JP3363171B2 (ja) * 1992-02-27 2003-01-08 本田技研工業株式会社 並列運転機能を備えた発電機
US5615099A (en) * 1994-08-24 1997-03-25 Kabushiki Kaisha Toshiba Control system for single-phase PWM converter
TW366796U (en) * 1998-07-21 1999-08-11 Hong-Yu Su Cordless transmission apparatus for infrared ray heart-beat sensor
US20020094796A1 (en) * 2001-01-12 2002-07-18 Woods Donnie W. Phase detector circuit and method therefor
US6771712B2 (en) * 2001-07-27 2004-08-03 The Pulsar Network, Inc. System for extracting a clock signal and a digital data signal from a modulated carrier signal in a receiver
AU2003217032A1 (en) * 2002-03-27 2003-10-13 C-Lect Medical Ltd A respiratory monitor
JP3985956B2 (ja) * 2002-09-11 2007-10-03 オータックス株式会社 脈拍計
JP3691479B2 (ja) * 2002-11-25 2005-09-07 三洋電機株式会社 心拍／呼吸計測装置
JP4174825B2 (ja) * 2003-05-23 2008-11-05 株式会社テクノリンク 生体刺激装置
TWI221765B (en) * 2003-07-03 2004-10-11 Tu Shiue He Two-point-type device and method for measuring heart rate
CN2675047Y (zh) * 2004-02-20 2005-02-02 肖青 非线性胎儿心率检测仪
JP4365795B2 (ja) * 2005-02-18 2009-11-18 株式会社ルネサステクノロジ パルス発生器及びそれを用いた送信機
TWI293682B (en) * 2006-02-24 2008-02-21 Yuan Jay Wang Multipurpose micro cantilever sensing sensitivity control device
KR20070122324A (ko) * 2006-06-26 2007-12-31 주식회사 대우일렉트로닉스 맥박수 체크가 가능한 리모콘
JP5081001B2 (ja) * 2008-01-28 2012-11-21 パナソニック株式会社 高圧放電灯点灯装置、照明器具
US7791912B2 (en) * 2008-05-02 2010-09-07 Advanced Energy Industries, Inc. Protection method, system and apparatus for a power converter
JP4911733B2 (ja) * 2009-03-13 2012-04-04 オムロン株式会社 電力変換装置、パワーコンディショナ、および発電システム
EP2309632B1 (en) * 2009-10-12 2013-05-29 STMicroelectronics Srl Half bridge resonant DC-DC control device
CN202069586U (zh) * 2010-12-12 2011-12-14 邵明省 一种对运动员心率动态监测的装置
JP2013219569A (ja) * 2012-04-10 2013-10-24 Seiko Epson Corp トランスコンダクタンス調整回路、回路装置及び電子機器
DE102013201686A1 (de) * 2013-02-01 2014-08-21 Siemens Aktiengesellschaft Verfahren und Schaltungsanordnung zur Wandlung eines Sinussignals in ein Rechtecksignal, sowie Verwendung hierfür
CN105162356B (zh) * 2015-08-24 2018-08-07 上海沈德医疗器械科技有限公司 一种电压相位可调且易于并联的全桥谐振超声驱动电路
JP6598365B2 (ja) * 2015-10-20 2019-10-30 日本放送協会 帯域合成装置、帯域分割装置、解像度変換装置、超解像装置およびプログラム
CN106645952B (zh) * 2016-10-18 2019-06-25 上海华虹计通智能系统股份有限公司 一种信号相位差的检测方法及系统
CN106725408A (zh) * 2016-12-02 2017-05-31 重庆软汇科技股份有限公司 基于自适应数字滤波的心率估计方法及装置
CN207837537U (zh) * 2017-05-03 2018-09-11 苏州工业职业技术学院 一种智能心率计
CN107260174A (zh) * 2017-07-10 2017-10-20 中央军委后勤保障部军需装备研究所 一种基于呼吸相位差法的呼吸率检测系统及方法
CN208837927U (zh) * 2018-02-12 2019-05-10 西安科技大学 一种脉搏和体温检测仪
DE112019001050T5 (de) * 2018-02-28 2020-12-03 Fujifilm Corporation Vorrichtung, verfahren und programm zur superauflösungsverarbeitung

Also Published As

Publication number Publication date
CN111281396A (zh) 2020-06-16

Similar Documents

Publication Publication Date Title
Arbabi et al. Study of dynamics in post-transient flows using Koopman mode decomposition
US7248994B1 (en) Digital method and apparatus for sensing position with a linear variable differential transformer
CN107561367B (zh) 一种基于压缩感知理论的宽频谱阻抗测量方法
CN103983849B (zh) 一种实时高精度的电力谐波分析方法
CN105933005A (zh) 一种基于等效采样的时域交替模数转换器失配校准方法
CN101320060A (zh) 快速相位计
CN104655929A (zh) 一种时域信号的数字时频测量方法及相应的目标识别方法
CN103399204A (zh) 一种基于Rife-Vincent（II）窗插值FFT的谐波与间谐波检测方法
CN103969508B (zh) 一种实时高精密的电力谐波分析方法及装置
CN111281396B (zh) 胸腹表面呼吸运动信号超分辨方法
CN111135480A (zh) 胸腹表面呼吸运动信号超分辨电路
CN211513177U (zh) 胸腹表面呼吸运动信号超分辨电路
Dongfang Dongfang Com Quantum Equations of LIGO Signal
CN104901577B (zh) 一种三相逆变器死区时间在线调整及补偿方法
CN111135482B (zh) 一种胸腹表面呼吸运动信号准相位超分辨方法
CN211513176U (zh) 一种胸腹表面呼吸运动信号准相位超分辨电路
CN110208601A (zh) 基于fpga的瞬时测频方法及其数字接收机
CN102809688B (zh) 一种基于迭代傅里叶变换计算三相电压实时相位的方法
CN111135481B (zh) 一种胸腹表面呼吸运动信号错相位超分辨方法
CN211513175U (zh) 一种胸腹表面呼吸运动信号错相位超分辨电路
CN111135478A (zh) 一种胸腹表面呼吸运动信号准相位超分辨电路
CN111135479A (zh) 一种胸腹表面呼吸运动信号错相位超分辨电路
CN105337612A (zh) 滤除工频干扰的软件锁相环
WO2021087632A1 (zh) 一种有限信号的精确分解方程式构建及分解方法
CN111551785B (zh) 基于无迹卡尔曼滤波的频率与谐波检测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant