TWI834615B - 用於一個體之一身體上及/或身體內部運動之非接觸量測之方法及用於量測一個體之運動之系統 - Google Patents
用於一個體之一身體上及/或身體內部運動之非接觸量測之方法及用於量測一個體之運動之系統 Download PDFInfo
- Publication number
- TWI834615B TWI834615B TW107121019A TW107121019A TWI834615B TW I834615 B TWI834615 B TW I834615B TW 107121019 A TW107121019 A TW 107121019A TW 107121019 A TW107121019 A TW 107121019A TW I834615 B TWI834615 B TW I834615B
- Authority
- TW
- Taiwan
- Prior art keywords
- signal
- tag
- antenna
- motion
- ncs
- Prior art date
Links
- 230000033001 locomotion Effects 0.000 title claims abstract description 228
- 238000000034 method Methods 0.000 title abstract description 52
- 238000005259 measurement Methods 0.000 claims abstract description 78
- 230000008878 coupling Effects 0.000 claims abstract description 28
- 238000010168 coupling process Methods 0.000 claims abstract description 28
- 238000005859 coupling reaction Methods 0.000 claims abstract description 28
- 230000005540 biological transmission Effects 0.000 claims abstract description 14
- 230000036772 blood pressure Effects 0.000 claims description 29
- 230000000241 respiratory effect Effects 0.000 claims description 26
- 230000001427 coherent effect Effects 0.000 claims description 17
- 238000004891 communication Methods 0.000 claims description 9
- 239000004744 fabric Substances 0.000 claims description 9
- 238000003306 harvesting Methods 0.000 claims description 7
- 230000000968 intestinal effect Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000005855 radiation Effects 0.000 abstract description 8
- 230000029058 respiratory gaseous exchange Effects 0.000 description 63
- 238000002565 electrocardiography Methods 0.000 description 55
- 210000000707 wrist Anatomy 0.000 description 48
- 238000001514 detection method Methods 0.000 description 27
- 210000001519 tissue Anatomy 0.000 description 26
- 238000012544 monitoring process Methods 0.000 description 25
- 238000005070 sampling Methods 0.000 description 25
- 241001465754 Metazoa Species 0.000 description 21
- 238000010586 diagram Methods 0.000 description 20
- 238000001914 filtration Methods 0.000 description 18
- 238000004088 simulation Methods 0.000 description 18
- 230000035945 sensitivity Effects 0.000 description 15
- 238000001228 spectrum Methods 0.000 description 15
- 230000001976 improved effect Effects 0.000 description 12
- 230000007774 longterm Effects 0.000 description 12
- 238000006243 chemical reaction Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 11
- 230000008859 change Effects 0.000 description 10
- 230000003595 spectral effect Effects 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 9
- 238000002474 experimental method Methods 0.000 description 9
- 230000007958 sleep Effects 0.000 description 9
- 238000012360 testing method Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 8
- 238000013461 design Methods 0.000 description 8
- 210000004204 blood vessel Anatomy 0.000 description 7
- 230000005684 electric field Effects 0.000 description 7
- 238000002955 isolation Methods 0.000 description 7
- 239000000523 sample Substances 0.000 description 7
- 230000001360 synchronised effect Effects 0.000 description 7
- 238000009826 distribution Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000012549 training Methods 0.000 description 6
- 241000251468 Actinopterygii Species 0.000 description 5
- 241000726096 Aratinga Species 0.000 description 5
- 241000699800 Cricetinae Species 0.000 description 5
- 230000008901 benefit Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000027288 circadian rhythm Effects 0.000 description 5
- 238000000605 extraction Methods 0.000 description 5
- 230000036541 health Effects 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 238000000926 separation method Methods 0.000 description 5
- 238000012706 support-vector machine Methods 0.000 description 5
- 241000271566 Aves Species 0.000 description 4
- 206010028347 Muscle twitching Diseases 0.000 description 4
- 241000700159 Rattus Species 0.000 description 4
- 241000270708 Testudinidae Species 0.000 description 4
- 230000000747 cardiac effect Effects 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 4
- 210000003205 muscle Anatomy 0.000 description 4
- 230000000737 periodic effect Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 238000002604 ultrasonography Methods 0.000 description 4
- 241001461451 Testudo horsfieldii Species 0.000 description 3
- 230000003542 behavioural effect Effects 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 3
- 238000013480 data collection Methods 0.000 description 3
- 230000035487 diastolic blood pressure Effects 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 3
- 230000004424 eye movement Effects 0.000 description 3
- 210000004072 lung Anatomy 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000009611 phonocardiography Methods 0.000 description 3
- 238000002360 preparation method Methods 0.000 description 3
- 230000036387 respiratory rate Effects 0.000 description 3
- 230000002269 spontaneous effect Effects 0.000 description 3
- 101710148586 ADP,ATP carrier protein 1 Proteins 0.000 description 2
- 101710111394 ADP,ATP carrier protein 1, mitochondrial Proteins 0.000 description 2
- 101710148588 ADP,ATP carrier protein 2 Proteins 0.000 description 2
- 101710165307 ADP,ATP carrier protein 2, mitochondrial Proteins 0.000 description 2
- 101710102716 ADP/ATP translocase 1 Proteins 0.000 description 2
- 102100032533 ADP/ATP translocase 1 Human genes 0.000 description 2
- 101710102718 ADP/ATP translocase 2 Proteins 0.000 description 2
- 102100026396 ADP/ATP translocase 2 Human genes 0.000 description 2
- 101710102715 ADP/ATP translocase 3 Proteins 0.000 description 2
- 241000269799 Perca fluviatilis Species 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 230000002612 cardiopulmonary effect Effects 0.000 description 2
- 230000009365 direct transmission Effects 0.000 description 2
- 230000005672 electromagnetic field Effects 0.000 description 2
- 230000008451 emotion Effects 0.000 description 2
- 239000003292 glue Substances 0.000 description 2
- 210000002837 heart atrium Anatomy 0.000 description 2
- 238000009532 heart rate measurement Methods 0.000 description 2
- 230000036039 immunity Effects 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 208000015181 infectious disease Diseases 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 230000000116 mitigating effect Effects 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000004806 packaging method and process Methods 0.000 description 2
- 210000000006 pectoral fin Anatomy 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 230000004202 respiratory function Effects 0.000 description 2
- 230000000284 resting effect Effects 0.000 description 2
- 230000033764 rhythmic process Effects 0.000 description 2
- 238000011808 rodent model Methods 0.000 description 2
- 230000003860 sleep quality Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000010146 3D printing Methods 0.000 description 1
- 206010002091 Anaesthesia Diseases 0.000 description 1
- 208000027896 Aortic valve disease Diseases 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 241000046053 Betta Species 0.000 description 1
- 241000218645 Cedrus Species 0.000 description 1
- 241000270617 Cheloniidae Species 0.000 description 1
- 241000938605 Crocodylia Species 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 240000007472 Leucaena leucocephala Species 0.000 description 1
- 235000010643 Leucaena leucocephala Nutrition 0.000 description 1
- 241000124008 Mammalia Species 0.000 description 1
- 241000699678 Mesocricetus Species 0.000 description 1
- 241000699673 Mesocricetus auratus Species 0.000 description 1
- 240000008790 Musa x paradisiaca Species 0.000 description 1
- 235000018290 Musa x paradisiaca Nutrition 0.000 description 1
- 241000700157 Rattus norvegicus Species 0.000 description 1
- 208000037656 Respiratory Sounds Diseases 0.000 description 1
- 241000283984 Rodentia Species 0.000 description 1
- 229910021607 Silver chloride Inorganic materials 0.000 description 1
- 241000270666 Testudines Species 0.000 description 1
- 241000251539 Vertebrata <Metazoa> Species 0.000 description 1
- 206010047924 Wheezing Diseases 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 210000000577 adipose tissue Anatomy 0.000 description 1
- 230000037005 anaesthesia Effects 0.000 description 1
- 210000001765 aortic valve Anatomy 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 238000002555 auscultation Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 210000000078 claw Anatomy 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000036461 convulsion Effects 0.000 description 1
- 239000006071 cream Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013481 data capture Methods 0.000 description 1
- 238000003066 decision tree Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 210000003746 feather Anatomy 0.000 description 1
- 210000000245 forearm Anatomy 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 150000003278 haem Chemical class 0.000 description 1
- 210000004247 hand Anatomy 0.000 description 1
- 210000005003 heart tissue Anatomy 0.000 description 1
- 238000011905 homologation Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005404 monopole Effects 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 230000000474 nursing effect Effects 0.000 description 1
- 230000000414 obstructive effect Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 238000010827 pathological analysis Methods 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000000541 pulsatile effect Effects 0.000 description 1
- 238000000718 qrs complex Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- HKZLPVFGJNLROG-UHFFFAOYSA-M silver monochloride Chemical compound [Cl-].[Ag+] HKZLPVFGJNLROG-UHFFFAOYSA-M 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000007920 subcutaneous administration Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000035488 systolic blood pressure Effects 0.000 description 1
- 238000013334 tissue model Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 230000000283 vasomotion Effects 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
- 210000001835 viscera Anatomy 0.000 description 1
- 238000009528 vital sign measurement Methods 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
- 230000003245 working effect Effects 0.000 description 1
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本發明可體現為用於一個體之一身體上及/或身體內部運動之非接觸量測之方法及/或系統。在待量測之一運動之一近場耦合範圍內提供一感測信號。以此方式,可產生一量測信號作為藉由該運動調變之該感測信號。該感測信號可係一ID調變信號。在一些實施例中,該感測信號係提供一無線標籤之一反向散射RFID鏈路。可提供一下行鏈路信號以為該無線標籤充電。該感測信號可係該下行鏈路信號之一諧波。偵測該量測信號。基於該量測信號量測該運動。該量測信號可經偵測為透射穿過該運動之一源之後的遠場輻射。該量測信號可經偵測為從該運動之一源反射為天線反射。
Description
本發明係關於運動之偵測,且特定言之,在無需實體接觸之情況下偵測生命徵象。
生命徵象之監測(諸如例如心率、血壓、呼吸率及呼吸運作)係用於患者管理及病理記錄之一關鍵程序。基於身體電極、光學吸收、壓力或應變計及超音波或射頻(RF)反向散射之當前實踐具有對感測能力及取樣率之各自限制。量測程序亦可歸因於直接皮膚接觸或運動約束而使受測試個體不舒適或破壞個體之晝夜節律。此不便顯著增加護理者之工作負載且阻礙連續長期監測。
歸因於身體表面之強反射及所反射信號之幾何平均值,用於基於射頻(「RF」)之遠場生命徵象偵測之習知技術可容易地揀取呼吸運動,但此等技術難以區分小機械振動細節,諸如具有低RF頻率之心跳及腕部脈搏波形。儘管可在仔細濾波之後擷取心率,但血壓之估計及多個自由移動之人的同時監測仍無法實現。
仍長久企盼用於偵測生命徵象之低侵入式技術。
本發明提供一種RF近場同調感測(NCS)而無需直接皮膚接觸之方法。本發明方法可用於將身體上及/或其內部之即時機械運動直接調變至一多工無線電上,該多工無線電可含有一唯一識別(ID)。本文中描述兩項例示性實施例以提供部署及操作之靈活性:心臟及腕部附近之被動標籤及主動標籤。為減少部署及維護成本,被動RFID標籤可整合至胸部及腕部區域處之服裝中,其中可在讀取器處收集兩個多工遠場反向散射波形以擷取心率、呼吸率、呼吸運作(breath effort)及血壓。為改良讀取範圍及對室內運動之免疫力,主動標籤可放置於例如一前部口袋及一腕部袖口中以量測歸因於NCS之天線反射,且接著完全以數位格式取樣及傳輸生命信號以消除室內多路徑干擾。本發明揭示之生命徵象監測系統可同時用於多個個體,且可帶來護理設施中之成本效益自動化。此外,直接皮膚接觸及運動約束之消除將增強患者舒適性,此可實現改良病理分析之長期監測。
透過使用用於生命徵象監測之一例示性諧波RFID系統之NCS之一實施方案,可監測呼吸、心跳及腕部脈搏以導出呼吸率、心率及血壓。在一些實施例中,可將頻譜等化應用至NCS信號以復原時域中之精確心跳時間間隔。
針對NCS展示感測天線阻抗匹配效應,其中一組織匹配天線經展示以提供改良NCS信號品質及波形細節。分析生命徵象波形擷取之效能且論述設計策略。
在另一實施例中,一心跳信號之高頻分量可用於緩解身體移動對心率估計之干擾。藉由以ECG為基準,NCS經展示為對於正常身體運動之即時心率及心率可變性係足夠精確的。
在另一態樣中,本發明可體現為一種用於一個體之一身體上及/或身體內部運動之非接觸量測之方法。在待量測之一第一運動之一近場耦合範圍內提供一第一射頻(「RF」)感測信號。以此方式,可產生一第一量測信號作為藉由該第一運動調變之該第一感測信號。該第一感測信號可係一ID調變信號。該第一感測信號可係一主動無線電鏈路。該第一感測信號可係一反向散射RFID鏈路。該第一感測信號可由一無線標籤提供。可提供一下行鏈路信號以為該無線標籤充電。該第一感測信號可係該下行鏈路信號之一諧波。
偵測該第一量測信號。基於該第一量測信號量測該第一運動。該第一量測信號可經偵測為透射穿過該第一運動之一源之後的遠場輻射。該第一量測信號可經偵測為從該第一運動之一源反射為天線反射。量測該第一運動可進一步包括濾波該第一量測信號以透過時序及波形獲得一第一運動信號。
在一些實施例中,該方法進一步包含在待量測之一第二運動之一近場耦合範圍內提供一第二RF感測信號。以此方式,可產生一第二量測信號作為藉由該第二運動調變之該第二感測信號。偵測該第二量測信號。基於該第二量測信號量測該第二運動。可基於同步量測之該第一運動及該第二運動量測一導數值。
在另一態樣中,提供一種用於量測一個體之運動之系統。該系統包含用於產生一第一感測信號之一第一信號源。一第一天線與該第一信號源電通信。該第一感測信號可係一ID調變波。該第一感測信號可係一主動無線電鏈路。該第一感測信號可係一反向散射RFID鏈路。該第一天線經組態以安置於待量測之一第一運動之一近場耦合範圍內。該第一天
線可經組態以安置於一心臟運動、一脈搏、一呼吸運動、一腸道運動、一眼睛運動等之一耦合範圍內。以此方式,產生一第一量測信號作為藉由該第一運動調變之該第一感測信號。該系統包含用於偵測一第一量測信號之一第一接收器。該第一接收器可經組態以將該第一量測信號偵測為一透射信號。該第一接收器可經組態以將該第一量測信號偵測為一反射信號。
該第一信號源及該第一天線可經組態為一無線標籤。一標籤讀取器經組態以將一下行鏈路信號傳輸至該無線標籤。該接收器可係該標籤讀取器之一部分。該無線標籤可經組態以由該下行鏈路信號充電。該第一感測信號可具有一頻率,該頻率係該下行鏈路信號之一頻率之一諧波。例如,該第一感測信號可係該下行鏈路信號之一第二諧波。該無線標籤可使用一正交ID調變該下行鏈路信號,使得該第一感測信號係一CDMA信號。
在一些實施例中,該系統可進一步包含用於產生一第二感測信號之一第二信號源。一第二天線與該第二信號源電通信。該第二天線經組態以安置於待量測之一第二運動之一近場耦合範圍內以將一第二量測信號產生為藉由該第二運動調變之該第二感測信號。該接收器可經進一步組態以偵測該第二量測信號。
在一些實施例中,一濾波器與該接收器通信。該濾波器可經組態以解調變及濾波該第一及/或第二量測信號以獲得一對應第一及/或第二運動信號。該濾波器可係例如經程式化以取樣、解調變及濾波該第一及/或第二量測信號以導出運動之一處理器。在一些實施例中,一處理器經程式化以基於所偵測耦合信號及第二耦合信號量測一導數值。
100:方法
103:提供
106:偵測
109:量測
112:濾波
115:提供
118:偵測
121:量測
124:濾波
127:判定
為更充分理解本發明之性質及目的,應參考結合附圖進行之以下實施方式,其中:圖1A描繪藉由近場同調感測(NCS)透過無線電傳輸對生命徵象進行之一CST Microwave Studio模擬模型,其中場型展示用於心跳感測之一身軀模擬模型中之共極化電場之實際部分。
圖1B描繪藉由NCS透過無線電傳輸對生命徵象進行之一CST Microwave Studio模擬模型,其中場型展示用於脈搏感測之一腕部模擬模型中之共極化電場之實際部分。
圖2A描繪具有一鄰近血管之一概念皮膚及組織結構之一NCS模擬,其中發射1.85GHz、0dBm信號之一偶極天線放置於皮膚上方且無需直接皮膚接觸。血管上方之組織包含皮膚、脂肪及肌肉。脈搏血管之幾何改變將調變近場且改變遠場反向散射場型。
圖2B展示用於圖2A之準靜態血管之天線反射參數S11,其中橫剖面標記為t1、t2及t3。
圖3A展示其中比較圖1A及圖1B之模擬振動振幅與遠場點處之取樣的一圖表。
圖3B展示其中比較圖1A及圖1B之模擬振動振幅與由散射參數S11表示之天線反射處之取樣的一圖表。
圖4A展示用於心跳之模擬RF輻射場型,其中遠場樣本點在胸部前方1m。
圖4B展示用於一腕部脈搏之模擬RF輻射場型,其中遠場樣本點在腕部上方1m。
圖5A係根據本發明之一實施例之用於近場同調感測之一被動諧波
RFID標籤之一圖式。
圖5B係根據本發明之一實施例之一諧波RFID讀取器之一圖式。
圖6A展示從一例示性諧波RFID系統解調變之一原始呼吸信號及低通濾波之後的波形。
圖6B展示從一例示性諧波RFID系統解調變之一原始心跳信號及十秒移動窗內之一平均心跳。標記展示來自安裝於左前臂上之OMRON血壓監測器之量測。
圖6C展示在三分鐘之資料收集期間之一心跳之波形。
圖6D展示在三分鐘之資料收集期間之一腕部脈搏之波形。
圖6E展示動態時間扭曲(「DTW」)波形分析,其展示心跳之各波形之距離(插圖展示盒鬚分佈)。
圖6F展示DTW波形分析,其展示腕部脈搏之各波形之距離(插圖展示盒鬚分佈)。
圖6G展示與用於心跳之各自DTW模板相比較之中值距離波形及最大距離波形。
圖6H展示與用於腕部脈搏之各自DTW模板相比較之中值距離波形及最大距離波形。
圖7A展示從同步心跳及腕部脈搏波形估計之脈波傳遞時間(「PTT」)。插圖展示信號及所提取PTT之一個週期。
圖7B展示PTT在3分鐘內之概率密度分佈。
圖7C展示在受測試人員坐下時從PTT提取之血壓。星形標記展示從商用血壓監測器量測之血壓。
圖7D展示在受測試人員經歷一適度活動且站立時從PTT提取之血
壓。星形標記展示從商用血壓監測器量測之血壓。
圖8A係展示藉由反向散射與標籤互動之RFID讀取器之RF前端之一習知RFID系統之一圖式。
圖8B係展示諧波RFID反向反射之本發明揭示原理之一實施例之一圖式,其中Tx表示傳輸器,Rx表示接收器,LPF表示低通濾波器,且HPF表示高通濾波器。
圖9A係與織物上之刺繡天線整合之本發明之一實施例之一RFID感測器標籤晶片之一照片。
圖9B係根據本發明之一實施例之一諧波RFID標籤之一PCB原型之一照片。
圖10展示具有各種系統SNR下之CDMA(200個標籤)之一諧波RFID標籤之定位誤差。
圖11展示具有各種SNR下之CDMA之諧波反向散射之取樣率,其中標籤數目在讀取範圍內。當晶片碼長度與標籤數目成比例時,實線係來自CDMA。虛線表示具有更有效半正交碼之實施方案。
圖12A係展示根據本發明之一實施例之一天線反射系統之一圖式。分離器係雙向的且可由一寬頻循環器替代。諧波收發器及諧波產生器類似於圖5A及圖5B之被動標籤系統中使用之諧波收發器及諧波產生器。
圖12B係展示使用圖12A之天線反射系統之經量測心跳波形之一圖式。
圖13展示當根據本發明之一實施例之一近場同調感測天線放置於一非對稱軸直流電減速馬達附近時之校準結果。解調變正弦波形係2D週期性非對稱旋轉之1D相位投影。
圖14係用於耦合至心臟運動之天線近場之處於1.8GHz之共極化電場之電磁波模擬之一橫剖面。
圖15係展示藉由標籤1之所收集馬達旋轉原始資料及對資料進行帶通濾波後的曲線之一圖式。
圖16係展示藉由標籤2針對具有500Hz之取樣頻率之一心跳波形收集之原始資料之一圖式。
圖17係展示從圖16提取之1.4Hz至15Hz帶通濾波之後的心跳波形之一圖式。
圖18係從圖16提取之心跳波形之一個循環,其中繪示特性點。
圖19展示從一胸部標籤(實線)及腕部標籤(虛線)記錄之一信號。時序差異給出血壓之一估計。
圖20係從圖17提取之心跳波形及腕部脈搏波形之一個循環,其中繪示兩個C峰值點。峰值1與峰值2之間的延遲係約0.074s,其轉變成約88mmHg之舒張壓。
圖21係展示根據本發明之另一實施例之一方法之一流程圖。
圖22A係一諧波RFID反向散射標籤之一示意圖。
圖22B係在圖22A中描繪之標籤之一PCB原型之一照片。
圖23係展示一實驗NCS設定之一照片。一第一標籤在胸部區域中且第二標籤在左腕上。感測波形展示在螢幕上。胸帶及腕帶在此處僅為部署便利。無需皮膚接觸或帶張力。
圖24A至圖24B係從胸部標籤NCS信號提取之心跳及呼吸波形。(A)主要藉由心跳調變NCS信號之振幅。頂部曲線係來自正規化振幅,其穿過數位BPF(0.8Hz至15Hz)以給出底部曲線。(B)主要藉由呼吸運動調變
NCS信號之相位。頂部曲線係正規化原始相位,其穿過數位BPF(0.1Hz至1.2Hz)以給出底部曲線。正交接收器中之振幅與相位調變之間的隔離致使呼吸波形與心跳波形之間的明確分離。
圖25A至圖25C係(A)同步時域NCS(實線)及ECG(點劃線)心跳信號。NCS之取樣率係5,000Sps,且ECG最初具有512Sps但增加取樣至5,000Sps。(B)NCS信號之頻譜,其中強度經正規化至約1Hz之峰值且放大以清晰展示較低強度部分。(C)ECG信號之頻譜,其中強度亦經正規化至約1Hz之峰值且放大以展示較低強度部分。2Hz與8Hz之間的高頻分量比NCS中之高頻分量更突出。
圖26A至圖26B係(A)高頻等化之後的時域NCS信號。三角形標記展示用於心跳時間間隔提取之尖銳峰值特徵點。(B)NCS及ECG信號之心跳時間間隔。虛線、點劃線及實線分別來自原始NCS(NCS1)、ECG及等化NCS(NCS2)。
圖27A至圖27B係NCS天線效應:(A)CST Microwave Studio中之一胸部RFID標籤之EM模擬模型。(B)用於標籤天線之模擬S11參數。
圖28A至圖28B係耦合至人體身軀中之NCS信號之處於2GHz之模擬EM功率流,其中(A)具有較差阻抗匹配且(B)具有適當阻抗匹配。
圖29A至圖29B係各種頻率中之實驗結果。(A)NCS天線在空氣中操作(實曲線)及放置於胸部上(虛曲線)時之S11。(B)從圖4中之胸部標籤解調變之心跳信號波形,其等具有不同感測頻率。來自實曲線(阻抗匹配條件)、虛曲線及點劃曲線之NCS信號分別對應於圖5A中之頻率1、2及3。
圖30係藉由導電油墨使用織物天線在電子產品碼(EPC)協定上運行之一例示性標籤之一照片。
圖31展示使用NCS從一胸部標籤提取之心跳及呼吸波形之實驗結果。
圖32A至圖32B係在大身體移動從40s開始時之NCS信號。(A)時域中之NCS信號之振幅。插圖係在受測試人員靜坐時的前40s中之解調變心跳波形。(B)在無移動之10s至30s期間的NCS頻譜(深灰)及具有大身體移動之85s至105s期間的NCS頻譜(淺灰)。
圖33A至圖33C係經提取心率。(A)用於心跳計數之三次諧波NCS信號(實線)及用於約97s之心跳參考之ECG信號(虛線)。(B)藉由NCS(實線)及ECG(虛線)計數心跳。底部線指示誤差(標記為ERROR)。插圖係在大身體移動在40s之後發生時的25s至55s期間的結果。(C)藉由NCS(實線)及ECG(虛線)監測之即時心率曲線。
圖34係用於呼吸、心跳及運動偵測之NCS設定。在胸部區域中部署一被動諧波RFID標籤且將諧波頻率反向散射至讀取器天線。對藉由SDR實施之諧波讀取器執行即時解調變及資料分析。
圖35A係DC濾波振幅及相位資料之一圖表。
圖35B係展示呼吸之一圖表。
圖35C係展示心跳之一圖表。
圖35D係展示峰值偵測之心跳之一圖表。
圖36A至圖36C係展示使用NCS獲得之心跳及對應運動偵測之圖表。
圖37A至圖37B係用於有意識小動物之生命徵象之近場同調感測(NCS)之兩個設定。藉由諧波讀取器收集信號且使用場可程式化閘極陣列(FPGA)及微控制器(MCU)進行數位基頻帶處理。(A)藉由具有被動諧波感測標籤之諧波RFID系統實現無線NCS系統。(B)有線NCS系統將讀取器與
諧波標籤之間的無線鏈路替換為RF纜線以減小干擾且便於室內實驗室部署。
圖38A係使用同步NCS及ECG量測之一麻醉大鼠之一實驗設定之一照片。
圖38B展示來自圖38A之實驗設定之NCS及ECG之5分鐘資料記錄。插圖展示一選定半秒持續時間內之波形細節。
圖38C展示從NCS提取之一心跳時間間隔,其展示與ECG信號之一緊密匹配。
圖38D展示在約第三秒至第四秒從NCS信號之相位解調變之代表性呼吸信號。
圖39A至圖39L展示使用非侵入性NCS設定對有意識小動物進行之生命徵象監測。(A)用於倉鼠之實驗設定。(B)從NCS信號解調變之心跳及呼吸信號。插圖展示約第八秒之心跳波形細節。(C)約30秒之心跳時間間隔。(D)用於長尾鸚鵡之有線NCS設定。(E)從NCS信號解調變之心跳及呼吸信號。(F)約1.5分鐘之經提取心跳時間間隔。(G)用於類似於圖39D之四爪陸龜之NCS設定,其中天線在木屑地板下方。(H)三分鐘之NCS信號之正規化原始振幅,其歸因於殼體結構而表示呼吸及心跳兩者。插圖展示重疊信號之波形細節。淺陰影區段指示心跳,而深陰影區段指示呼吸。強呼吸信號將在重疊期間覆蓋心跳信號。(I)在與(H)相同之週期期間藉由連續小波波形(CWT)線上處理之信號展示用於精確率估計之心跳及呼吸之明確分離。插圖展示經提取波形細節,其清楚指示心跳及呼吸中之每一峰值。(J)用於斗魚之類似於圖39D之NCS設定。在靠近魚之水中部署Tx及Rx天線。(K)由胸鰭移動引起之解調變NCS相位信號,其中插圖展示波形細
節。(L)可能由心跳引起之解調變NCS量值信號。
圖40係描繪根據本發明之另一實施例之一系統之一圖式。
本發明係在能源部授予的第DE-AR0000528號合約下由政府支援而實現。政府具有對本發明之某些權力。
本申請案主張2017年6月16日申請之現審查中之美國臨時申請案第62/521,163號及2018年1月17日申請之現審查中之美國臨時申請案第62/618,352號之優先權,該等案之揭示內容以引用的方式併入本文中。
生命徵象不僅對於病理研究係重要的,而且亦可由穿戴式裝置廣泛用於推斷行為、情緒及健康。儘管良好建立且廣泛應用許多此等裝置,但當前裝置具有限制其等感測精確性或長期便利性之缺點。本發明可體現為一種用於近場同調感測(「NCS」)之方法,其將一個體之移動調變至射頻(「RF」)信號上,射頻信號可係多工RF信號。一個體之移動可包含例如關於生命徵象之移動,例如心跳、脈搏、呼吸等。本發明方法之實施例可將個體之身體表面上或身體內部之機械運動直接調變至近場範圍內之RF信號上。可使用唯一數位識別(「ID」)將運動調變至多工諧波RF識別(「RFID」)反向散射信號上。
一天線之「近場」係其中感應特性主導優勢高於輻射特性且未良好界定電場(E場)與磁場(H場)之間的關係之一區域。在本發明之實施例中,「近場」可指一天線之接近區域,其中角視場分佈取決於距天線
之距離。在實施例中,近場延伸至天線之一個波長(λ)內之區域。在其他實施例中,近場延伸至天線之λ/2、λ/3、λ/4或λ/2π內之區域,其中λ係天線之操作波長。憑藉本發明之益處,熟習此項技術者將明白其他實施例。
參考圖21,本發明可體現為一種用於一個體之一身體運動(例如,身體上或身體內部運動)之非接觸量測之方法100。個體可係例如一人類或一非人類動物。所偵測運動可係例如一心臟運動、一脈搏、一呼吸運動、一腸道運動、一眼睛運動或如鑑於本發明將明白之其他身體運動。本發明方法100之實施例可將個體之身體表面上或身體內部之機械運動直接調變至與一唯一數位ID整合之多工無線電信號上。在待量測之一第一運動之一近場耦合範圍內提供103一第一射頻(「RF」)感測信號以產生一量測信號。所提供103第一感測信號可係一ID調變信號。在一些實施例中,第一感測信號係一主動無線電鏈路。在一些實施例中,第一感測信號係一反向散射RFID鏈路。例如,一天線可在一主動無線電鏈路或一反向散射RFID(無線電識別)鏈路中發射一信標或ID調變感測信號。將藉由第一運動調變第一感測信號,藉此產生一第一量測信號。方法100包含偵測106第一量測信號。在一些實施例中,偵測106可在遠場處完成,例如偵測傳輸通過個體之身體之第一量測信號。在一些實施例中,偵測106具有一反射信號,例如使用近場天線。
基於第一量測信號量測109第一運動。如上文提及,在NCS中,將比先前技術更多之能量引導至身體組織中,因此隱含地放大來自內臟之反向散射信號。而且,身體組織內之較短波長將一小機械運動呈現為一相對大相位變化。個體之身體內之較短波長自然地增大信雜比(「SNR」)。身體內信號之差分性質可隔離大表面移動。此亦可增大敏感
度,從而實現一弱運動信號(諸如例如一腕部脈搏)之量測。由於內部機械運動調變給出類似於一干擾計之一差分信號,故可藉由濾波容易地低降由外部移動導致之共同信號(例如,見圖2A)。憑藉在身體內部之機械運動之近場耦合範圍內之一天線,傳播或反射波可以一同調方式容易地偵測且將含有機械運動之即時幾何平均資訊。可藉由濾波112第一量測信號以獲得一運動信號而量測運動。在ID調變波之情況中,可以一同步方式同時讀取多個機械運動。可在被動反向散射或主動無線電傳輸中使用多工技術以促進在多個點及/或對多個人員之同時感測。NCS提供具有舒適性、便利性及低成本之生命徵象監測之新機會。
應注意,本發明揭示之NCS技術直接量測一身體內部或身體上之機械運動,而非間接藉由例如感測感應機械運動之電神經信號或由機械運動感應之電信號。因此,本發明NCS技術可提供比一習知心電圖更豐富之資訊。例如,藉由在其中可感受一脈搏之腕部或頸部區域上使用一第二標籤,來自一心臟標籤之波形差異可用於獲得血壓之一精確估計,且此可對房間內之多人進行而無歧義。使用本發明技術,可以一全新方式管理臨床領域:穿戴根據本發明之(若干)標籤之所有人可監測其等之ID、位置、心率、呼吸率、血壓等。另外,由於無需皮膚接觸(例如,相較於(若干)ECG襯墊),故可完成一個體之長期監測。
在一些實施例中,方法100可進一步包含在待量測之一第二運動之一近場耦合範圍內提供115一第二RF感測信號。以此方式,第二運動耦合至第二RF感測信號以產生一第二量測信號。偵測118第二量測信號,且基於第二量測信號量測121第二運動。可藉由濾波124第二量測信號以獲得一第二運動信號而量測第二運動。可基於同步量測之運動及第二
運動判定127一導數值。例如,在第一運動係一心跳(在胸部附近量測)且第二運動係一脈搏(在腕部附近量測)之情況下,導數值可係基於心跳及脈搏判定127之一血壓。
在另一態樣中,本發明可體現為一種用於量測一個體之運動之系統10(例如,見圖40)。系統10包含用於產生一第一感測信號之一第一信號源12。一第一天線14與第一信號源12電通信。第一天線14經組態以安置於待量測之一第一運動之一近場耦合範圍內。例如,第一天線14可經組態以安置於一心臟運動、一脈搏、一呼吸運動、一腸道運動、一眼睛運動等之一耦合範圍內。以此方式,由藉由第一運動調變之第一感測信號產生一第一量測信號。第一感測信號可係一ID調變波。例如,EM波可係一主動無線電鏈路或一反向散射RFID鏈路。
系統包含用於偵測第一量測信號(與第一運動耦合(藉由第一運動調變)之第一感測信號)之一接收器16。接收器16可經組態以將第一量測信號偵測為一透射信號,即,遠場輻射。接收器可經組態以將第一量測信號偵測為一反射信號,即,天線反射。系統可包含與接收器通信之一濾波器,其中濾波器經組態以解調變及濾波第一量測信號以獲得一運動信號。濾波器可係例如經程式化以取樣、解調變及/或濾波第一量測信號以導出運動信號之一處理器(諸如一數位信號處理器(「DSP」))。
在一些實施例中,一系統10可包含用於產生一第二感測信號之一第二信號源22。在此等實施例中,一第二天線24與第二信號源22電通信。第二天線24經組態以安置於待量測之一第二運動之一近場耦合範圍內。以此方式,可產生一第二量測信號作為藉由第二運動調變之第二感測信號。在一特定實例中,第一運動係一心跳且第二運動係一脈搏。在此
一實例中,第一天線可經組態以接近於一個體之胸部安置,且第二天線經組態以接近於個體之腕部安置。接收器16經進一步組態以偵測第二量測信號。系統10可包含用於基於所偵測第一量測信號及第二量測信號量測一導數值之一處理器30。在一心跳及脈搏之特定實例中,導數值可係例如個體之一血壓。
在一些實施例中,無線標籤(諸如被動(即,不具有諸如一電池之本端電源)RFID標籤)可整合至其中待量測生命徵象之區域附近之服裝中。此等RFID標籤可提供具有低部署及維護成本之一NCS實施方案。此等RFID標籤可提供ID調變信號,其中各標籤之一唯一ID有助於區分其信號與來自其他標籤及周圍信號之干擾。接著使用頻譜等化處理標籤反向散射信號以放大高頻分量,以不僅復原最初浸沒於低頻分量中之波形細節,而且亦復原具有改良峰值偵測確定性之精確心跳時間間隔之尖銳峰值。與同步ECG相比較,所導出心跳時間間隔展現改良穩定性。
NCS之實施例利用一EM場至一身體內部或身體表面上之機械運動的近場耦合。使用CST Microwave Studio進行電磁模擬以繪示NCS之操作原理。如圖1A及圖1B中展示,建構一男性身軀(1A)及一左下臂(1B)之EM模擬。基於Zubal Phantom建構EM模擬模型,Zubal Phantom具有3.6×3.6×3.6mm3之立體像素解析度且由磁共振成像(「MRI」)及電腦斷層掃描(「CT」)建立。各立體像素由人體組織之3D座標連同指數一起指示。在CST匯入身軀模型之後,Visual Basic for Applications(「VBA」)巨集語言用於使用CST生物庫映射EM性質。接著採用有限積分技術(「FIT」)以包含胸部區域附近之RFID標籤天線。如
從模擬可見,大量RF能量歸因於近場耦合而耦合於身軀內部。由於人體組織之高介電常數,波長相應地變短,此進一步增大NCS敏感度。
本發明NCS方法使用電磁場之振幅及相位兩者。由於相位對RF源與接收器之間的距離非常敏感,故因此可藉由相位評估一人員呼吸時之外部胸部移動。可容易地擷取呼吸率且可使用相位變化進一步解釋呼吸運作。與相位資訊相比較,電磁場之振幅對小距離變化不那麼敏感,此意謂呼吸或其他外部身體移動將改變相位而非振幅,從而提供待適當感測之身體內部之其他信號之良好隔離。在NCS中,干擾量測類結構將內臟/組織移動轉換為RF信號之振幅調變。
對於吾人之模擬(圖2A),當人體模型(human phantom)面向接收器時(圖4A),胸部上天線發射RF載波,其中藉由局部近場區域界定天線特性。根據天線指向性,將朝向接收器直接發射RF能量之部分,而其他部分將歸因於近場效應而耦合於身體內部。直觀地,可考量藉由心臟組織之機械移動調變來自心臟之反向散射RF信號且接著使用直接發射進行干擾,從而導致振幅改變。從干擾計類比,身體內部之移動係一「差分模式」調變,而身體表面移動係一「共同模式」調變。
不僅可從遠場而且亦可從如圖3B中之散射參數S11展示之天線反射記錄運動。使用天線反射,一NCS信號可使用一行動裝置直接記錄且因此對身體移動及一擁擠房間內之室內多路徑問題更具免疫力。由於NCS運用天線之近場區域內之組織運動操作,故幾何變化將影響天線反射S11,其中天線可被視為感測器之一部分。生命徵象將調變在天線之S11參數上且因此藉由反射信號擷取。
在傳輸天線接近於皮膚之情況下,一NCS裝置可透過一無
線電信號調變運動信號(生命徵象)。然而,習知微波傳輸器在局部振盪器及功率放大器中消耗顯著功率,且此等傳輸器可需要用於行動裝置之一電池。另外,身體上傳輸器與遠場接收器之間的同步亦將使系統設計更複雜。在一些實施例中,可使用被動諧波RF識別(RFID)標籤實施NCS,其中在諧波反向散射連同標籤ID上調變生命信號。除超低成本以外,被動標籤之簡單及穩健封裝亦實現與準備衣物之直接織物整合。在圖9A中展示與織物上之一刺繡天線整合之RFID感測器標籤晶片之一實例。在圖8A及圖8B中概述諧波反向散射優於習知RFID之益處。由於習知RFID讀取器及相位雜訊邊緣之高傳輸功率,來自非預期周圍物件之自洩漏、天線反射及反向散射皆促成雜訊且使反向散射標籤信號之SNR嚴重降級。然而,諧波反向散射可使用一大頻率分離來隔離下行鏈路(讀取器至標籤)與上行鏈路(標籤至讀取器),此增大SNR及敏感度兩者。標籤仍為一被動反向散射器,其可容易地遵循當前RF協定。在圖5A中展示諧波標籤之一示意圖(印刷電路板(PCB)原型之一照片見圖9B)。諧波標籤從讀取器接收處於f之下行鏈路RF信號,該信號通過標籤天線1(ANT 1)且分成兩個部分。一個部分藉由能量採集而為標籤電路提供DC功率,且另一部分在2f下饋送至被動諧波產生中以從天線2(ANT 2)重新發射,天線2用作為NCS傳輸器。諧波產生器前方之RF開關可藉由開關鍵控(OOK)而調變數位資訊,類似於習知RFID操作。數位資訊可包含標籤ID以及來自標籤上感測器之額外資訊。
在圖5B中展示作為一同調收發器之一例示性諧波RFID讀取器之一示意圖。相同數位時脈(虛線)經饋送至處於f及2f之兩個頻率合成器中以進行處於2f之同調解調變。數位模組執行CDMA協定。從讀取器至
標籤之下行鏈路命令藉由數位轉類比轉換器(「DAC」)調變且接著藉由混頻器升頻轉換為處於f之載波。諧波標籤反向散射至處於2f之讀取器,f藉由處於2f之同調局部振盪器降頻轉換為基頻帶且藉由正交類比轉數位轉換器(「ADC」)取樣。使用一軟體定義無線電(SDR)進行一諧波讀取器之一測試實施例之硬體。
相位對相對於讀取器之標籤實體位置更敏感。因此,當一標籤之天線2放置於一個體之胸部上時,可從正交方案中之相位導出呼吸資訊,如在圖6A中使用原始及低通濾波波形展示。基於反向散射相位資訊,可憑藉毫米解析度計算多個標籤之位置,其等可進一步導出呼吸運作。儘管由胸部移動導致之相位變化遠強於心跳及腕部脈搏之內部移動,但其係身體內部之組織運動之NCS之一「共同分量」(如上文進一步描述)。在使用本系統之實施例執行之實驗期間,NCS心跳信號對由個體之呼吸導致之移動具免疫力。多個頻率、改良信號處理及反射結構(圖12A及圖12B)之使用可進一步減輕嚴重多路徑干擾。首先,從瞬時週期(實線)及從10秒內之計數(虛線)擷取圖6B中之心率。從一商用血壓監測器(OMRON BP760N)量測星形標記。應注意,從正交解調變獨立地導出呼吸及心跳資訊,且無需如習知微波反向散射中之特殊濾波或圖案辨識。
由於從NCS擷取內部生命徵象,似干擾計結構(interferometer-like structure)顯著增大敏感度以實現運動波形之收集,類似於一心衝擊圖(「BCG」)。同時記錄來自胸部及腕部標籤之資料達3分鐘。使用PCB標籤進行試驗且讀取器天線距離受測試人員~1.5m至2m。藉由標籤轉換之諧波信號在1.9GHz(2f)下為約-20dBm。為分析波形
變化,疊置各週期以獲得如圖6C及圖6D中展示之心跳及腕部脈搏波形之平均值及盒鬚偏差。波形經正規化至記錄資料之第90個百分位數。應用動態時間扭曲(「DTW」)以分類波形以導出詳細特徵。圖6E及圖6F展示心臟及脈搏波形之DTW距離,且插圖展示變化。圖6G及圖6H展示所提取模板波形與最大距離波形及中值距離波形之比較,中值距離波形仍非常類似於模板且保持大部分主要特徵,諸如腕部脈搏中之反衝峰值。詳細運動波形分析可用作例如用於心律不整及主動脈瓣疾病之一心電圖候選。
CDMA協定不僅能夠同時監測多個人員,而且亦能夠同時監測相同人員上之多個點。CDMA標籤之可允許數目受限於基頻帶資料率且在圖11中展示。來自不同身體位置之波形時序之比較透過脈波傳遞時間(「PTT」)提供血壓(「BP」)之估計,其可從近端動脈波形及遠端動脈波形之特徵點提取。提出之非接觸式血壓感測呈現優於基於直接壓力之方法之顯著優勢,基於直接壓力之方法尤其針對年老患者之長期監測導致不適性且破壞晝夜節律。記錄胸部標籤信號(近端波形)及左腕標籤信號(遠端波形)之各者達三分鐘,如圖7A中展示。可從兩個波形之主峰容易地提取PTT。插圖展示一特定週期之詳細波形。圖7B展示PTT在3分鐘記錄期間的概率密度。PTT之分佈會受到取樣抖動及波形失真影響。可針對各心跳獲得一個PTT取樣,且可容易地應用移動平均值或其他信號處理方法以使PTT變化最小化。圖7C及圖7D展示從PTT計算之血壓以及來自商用血壓監測器(OMRON BP760N)之比較點(星形標記)。實線係每一心跳之收縮壓及舒張壓。虛線係來自約10秒之取樣點之移動平均值。圖7C中之資料係在受測試人員坐在一椅子上約30分鐘時收集,而圖7D中之資料係在一適度活動之後收集。
使用CST Microwave Studio進行電磁模擬。Zubal Phantom用於建構介電模型。使用來自電腦斷層掃描(CT)及磁共振成像(MRI)之資料校準組織幾何資訊。立體像素之解析度係3.6mm×3.6mm×3.6mm。使用CST生物庫映射各種組織之微波性質。首先逐層將Zubal Phantom資料預處理為組織幾何座標連同組織指數之檔案結構。CST接著匯入檔案且使用三維座標及組織性質自動建立每一立體像素,以建立藉由CST內建Visual Basic for Applications(VBA)巨集語言之描述性語言控制之介電模型。程序類似於三維印刷,但實際上僅在CST軟體中進行。心跳及腕部脈搏之動態模擬藉由幾何變化實現,其中心臟及腕部血管之幾何性質根據用作為實況之預設尺寸而變化。
藉由一自訂PCB原型化被動諧波反向散射標籤,從無線識別及感測平台(WISP)修改該自訂PCB。標籤上之諧波產生器經設計有一非線性傳輸線(NLTL),其包含電感器及變容器之一梯形結構。NLTL可提供具有低輸入功率之高轉換效率,此對於被動反向散射標籤設計係必要的。諧波RFID讀取器及天線反射系統經建立於國家儀器Ettus軟體定義無線電(SDR)B210之平台上。為實現同調諧波解調變,接收器之局部振盪器(LO)需要直接從傳輸器LO之第二諧波頻率予以導出。用LabVIEW撰寫即時控制及解調變軟體。運用零差調變方案,操作頻率係f=950MHz(處於2f之二次諧波=1.9GHz)。下行鏈路類比基頻帶係10kHz且諧波轉換之後的上行鏈路類比基頻帶係20kHz。數位轉類比轉換及類比轉數位轉換兩者皆按每秒106個樣本(Sps)而操作。原始數位信號接著經濾波、數位地降頻轉換為D.C.頻帶、且使用CDMA演算法解碼以區分來自各標籤之資
訊。接著藉由500Sps之取樣率減少取樣來自各標籤之信號。藉由具有0.8Hz之一截止頻率之一低通濾波器處理呼吸信號。藉由介於0.9Hz與15Hz之間的帶通濾波器處理心跳及脈搏信號。被動標籤之當前操作範圍係~1.5m,受限於WISP平台。範圍可根據習知RFID系統在相同頻帶中之操作而延伸朝向10m。
在另一模擬中,在接近於其中可感受脈搏之心臟及左腕處部署天線。信號源應在天線之近場區內,但無需藉由天線進行直接皮膚接觸。圖2A展示簡化下臂結構及CST Microwave Studio中之近場區域中之電場。天線經組態以將更多能量耦合至組織中以獲得較大信雜比(「SNR」)。在圖2B中展示天線反射參數S11。中心頻率係約1.85GHz。由於高組織電容率,天線頻寬更寬。圖2A中之模擬結果展示耦合至皮膚層、脂肪層及肌肉層以及附近血管中之電場。
為模擬耦合至EM場上之機械運動,將一小振動引入至心臟及腕部血管之幾何標度中,且組織模型中之血管橫剖面在t1、t2及t3之時間戳記中準靜態地變化以表示脈搏變化。在圖3A及圖3B中使用粗實線展示正規化振動振幅。圖4A及圖4B展示在CST Microwave Studio中模擬之遠場場型。陰影斜坡指示相位,而形狀表示共極化電場之振幅輪廓。記錄且在圖3A中展示遠場信號(其中取樣點在胸部(圖4A)前方及腕部(圖4B)上方1m)。可見,解調變心臟信號(在圖3A中使用一虛線描繪)及腕部脈搏信號(使用一點劃線描繪)與已知振動良好地匹配。
遠場可被視為兩個近場分量之干擾結果:來自天線之直接傳播波(圖2A之黑色箭頭)及來自內部組織之散射信號(白色箭頭)。當一心跳導致血管振動時,將歸因於干擾直接傳播而調變散射信號之相位。由於
從差分干擾導出解調變血管運動信號,故由呼吸或其他身體移動導致之天線運動可視為可被拒絕之共同模式。替代地,散射信號可耦合回至相同天線以進行同調解調變,此經指示為圖2A中之灰色箭頭及圖2B中之所得S11。
圖2A繪示如何使來自天線之EM場能量之一部分直接輻射至遠場(如由黑色箭頭指示)而EM能量之其他部分耦合至多層組織中直至感測目標(此處係動脈血管)。由於介電常數差異,來自脈搏之動脈血管之機械運動將調變反向散射信號,由白色箭頭指示。此信號亦連同直接輻射信號一起傳播至遠場。從EM場之觀點而言,此兩個信號源自相同源但通過不同路徑,即,此兩個信號係同調的且組合信號之振幅將歸因於相位差異而改變。可因此明白血管之機械運動將如何導致振幅調變。此操作原理類似於干擾計,且近場調變可被視為干擾計類比中之「差分分量」。另一方面,表面或整個身體之移動(諸如呼吸及身體運動)將同時改變兩個信號路徑之相位,此可被視為干擾計結構中之「共同模式」。因此,NCS不僅利用干擾計結構之敏感度以增強效能,而且亦隔離兩個不同調變:來自具有近場耦合之身體內部之一個調變及來自具有直接發射之表面運動之另一調變。
在一例示性實施例中,可藉由被動諧波RF識別(RFID)標籤實施NCS,其中在諧波反向散射連同標籤ID上調變生命徵象信號。除超低成本以外,此一被動標籤實施例之簡單及穩健封裝亦實現與準備衣物之直接織物整合。使用Ettus ResearchTM軟體定義無線電(「SDR」)B210平台建立一例示性諧波RFID讀取器及天線反射系統。
在圖5A中展示一例示性諧波標籤之一示意圖,且圖9B係此一標籤之一PCB原型之一照片。為易於協定存取,在圖8B中基於無線識別及感測平台(「WISP」)對諧波標籤之PCB原型執行所有例示性NCS操作。原型RFID感測器標籤晶片與織物上之一刺繡天線整合(圖9A中展示)。織物RFID標籤用於示範本發明揭示方法之服裝整合之可行性。根據電子生產碼(「EPC」)協定驗證標籤ID及溫度感測之正常RFID資料異動。
用LabVIEW撰寫即時解調變軟體。運用零差調變方案,操作頻率係f=950MHz(處於2f之二次諧波=1.9GHz)。下行鏈路類比基頻帶係10kHz且諧波轉換之後的上行鏈路類比基頻帶係20kHz。數位轉類比轉換及類比轉數位轉換兩者皆按每秒106個樣本(Sps)而操作。所提出之被動標籤操作範圍係約1.5m,受限於WISP平台。範圍可根據習知RFID系統在相同頻帶中之操作而延伸朝向10m。
在圖8B中展示諧波RFID反向散射之原理。在一實施例中,一讀取器傳輸處於基本頻率f之一下行鏈路信號,其將供電給讀取器之一範圍內之任何諧波標籤。諧波標籤從讀取器接收處於f之下行鏈路RF信號,該下行鏈路信號通過標籤天線1且分成兩個部分:一個部分用於能量採集以提供DC功率給標籤電路,且另一部分經饋送至處於2f之被動諧波產生中以從天線2作為上行鏈路載波重新發射以消除讀取器自干擾及來自附近物件之反射。以此方式,天線2用作為NCS傳輸器。傳輸器(Tx)處之一低通濾波器(「LPF」)及接收器(Rx)處之一高通濾波器(「HPF」)進一步隔離處於f及2f之兩個載波。運用同調諧波反向散射,Rx具有一非常低雜訊底限。因此,可增大Rx敏感度,此使系統能夠區分在上行鏈路上
調變之弱生命信號。諧波產生器(圖5A)前方之一RF開關可藉由開關鍵控(「OOK」)而調變數位資訊,類似於習知RFID操作。數位資訊可包含標籤ID以及來自標籤上感測器之額外資訊。
可透過使用諸如例如分碼多重存取(「CDMA」)協定之一協定而實現多標籤存取。此一協定之使用可提供諧波反向散射中之較佳標籤間同步、較高頻道效率及較高功率效率。CDMA標籤之可允許數目受限於基頻帶資料率且在圖11中展示。實線係CDMA碼,其中晶片碼長度與標籤數目成線性比例。然而,當數目變大時,半正交碼長度可近似對數,展示為虛線。為復原波形細節,取樣率應高於500Hz,此可在具有CDMA之諧波反向散射系統之例示性實施例中容易地實現。
藉由額外路徑且因此相關機械移動(此處係動脈及靜脈脈搏)調變反向散射信號之相位。在遠場處,兩個部分如一干擾計般彼此干擾,即,腕部脈搏將調變接收器處之EM波之振幅。同時,外部移動(諸如例如擺手)調變兩個部分之相位且將導致兩個信號之共同調變。簡言之,手部運動將為共同模式且腕部脈搏將為差分模式。圖1A展示部署於胸部區域附近以監測心跳之感測天線。歸因於近場效應,EM場耦合至身軀中,且藉由心跳調變反向散射信號。藉由遠場處之接收器接收直接傳輸及反向散射信號。兩個部分係來自相同源且具有不同路徑,且其等將彼此干擾,其中心跳可從振幅解調變。由呼吸導致之胸部移動改變兩個信號之相位,此係NCS之共同模式。類似於腕部標籤中之手部運動情況,呼吸及心跳信號歸因於與後者濾波無關之共同及差分模式而良好地隔離。
為將NCS實施為室內環境中之一穿戴式裝置,使用具有分碼多重存取(CDMA)協定之一例示性諧波RFID系統。在圖8B中展示RF示
意圖。讀取器傳輸處於頻率f之下行鏈路(讀取器至標籤)信號。諧波標籤接收下行鏈路信號且將處於2f之二次諧波反向散射為上行鏈路(標籤至讀取器)信號。傳輸器(Tx)及接收器(Rx)頻譜經分離以增大信雜比(SNR)及接收器敏感度。CDMA協定提供更佳標籤同步、較高取樣率、較低抖動及較低功率消耗。在圖22A中展示諧波標籤之示意圖。天線1(ANT 1)接收下行鏈路RF信號,其中標籤採集RF能量、供電給邏輯電路且解調變下行鏈路資訊。下行鏈路RF信號之部分耦合至非線性傳輸線(NLTL)以產生二次諧波信號,藉由RF開關調變該二次諧波且透過天線2(ANT 2)傳輸回至讀取器。反向散射二次諧波信號執行NCS功能且接著被讀取器之Rx天線接收。在圖22B中展示用於以下實驗之PCB原型,其具有約57×20mm之尺寸。標籤可經進一步整合為一被動晶片且使用服裝直接封裝,其中天線可由圖30中之導電油墨及由圖9A中之刺繡實施。
圖10展示在結合CDMA協定應用基於相位之諧波反向散射定位方法時之標籤定位誤差。模擬結果展示各種SNR下之定位誤差之累積概率函數。由於在諧波反向散射系統中已消除自干擾及直接反射,故雜訊底限可非常低以容易地達成20dB之SNR。
在圖12A中展示實驗天線反射系統且在圖12B中展示經量測心跳波形。圖12A中之示意圖可由安裝於心臟附近之服裝上之一主動標籤實現,其中標籤天線執行NCS功能。對於直接比較,一諧波產生器及一諧波收發器經選擇為類似於圖5A及圖5B之被動標籤系統以建立天線反射單元。可使用其他設計,包含例如在一單一頻率中使用高隔離度循環器。實驗系統中使用之信號分離器及諧波產生器兩者皆為雙向的。來自處於f之Tx之信號耦合至諧波產生器且接著耦合至處於f及2f兩者之天線。天線
反射再次饋送通過諧波產生器及分離器,但僅約2f之信號被帶通濾波器(「BPF」)選擇用以饋送回至Rx以進行同調解調變。對超過2f之生命信號取樣且以完全數位格式發送至遠端裝置,其中可藉由具有低資料率之標準技術容易地消除由多路徑及占有人運動導致之室內符號間干擾。圖12B展示解調變心跳信號(應用0.9Hz至15Hz濾波器),正規化至資料之90百分位數。運用天線反射方案,藉由NCS偵測到之生命信號對由周圍移動導致之嚴重多路徑干擾(諸如一擁擠房間之情況)更具免疫力。
作為無需直接機械接觸之NCS操作之一獨立驗證程序,使用具有已知旋轉速率之一非對稱軸馬達(ASLONG JGB37-520)(圖13)。被動標籤之天線經放置接近於非對稱旋轉軸,且透過相位調變及解調變將機械旋轉轉換為一正弦波。
在另一實驗中,將根據本發明之一實施例之一第一標籤(標籤1)放置於具有一已知旋轉之一已知馬達外部,且將一第二標籤(標籤2)放置於受測試之一個體之襯衫外部之心臟區域內。在圖15中展示來自馬達(標籤1)之所收集波形,其用作為驗證或校準。在圖16中展示來自胸部(標籤2)之所收集波形。在應用介於1.4Hz至15Hz之間的一簡單帶通濾波器之後,獲得圖17中之心臟運動波形。在圖18中放大約三秒之一特定波形,其中藉由包含特性點A至F而擷取主要運動特徵。例如,點C指示主要收縮運動,且點E指示閉合主動脈瓣之後的衝擊波反衝。
在另一實驗中,標籤1放置於胸部區域附近且標籤2在相同人員之腕部附近,且在圖19中展示所得經濾波波形。在圖20中放大約三秒之片段,其中兩個波形之主峰(點C)之時間延遲可用於使用已知技術給出血壓之一估計。峰1與峰2之間存在35個取樣點,且時間間隔係約0.074
秒且具有470Hz之一取樣率,從而導致約88mmHg之一舒張壓。對所判定血壓估計(例如在每一心跳處獲得)取平均可改良估計之可靠性。
生命徵象(包含心率、血壓、呼吸率及呼吸運作)之連續監測對於eHealth係關鍵的。先前方法受限於精確性、便利性及感測能力。一些方法(諸如心電圖(ECG))歸因於直接皮膚接觸(其使穿戴者不舒適、限制身體運動且破壞晝夜節律)之要求而難以用於長期應用。同時,心跳間時序特徵對於健康監測及人類情感研究係重要的。
對於現有系統,心電圖(ECG)係用於心跳監測之最流行方法,其使用身體電極以透過由心跳感應之小皮膚電流而聚集身體電位。為達成良好信號品質,電極需要藉由不適導電膠及脫毛而進行直接皮膚接觸。為進一步減少雜訊,需要具有不可極化Ag/AgCl之大電極。再者,當頻繁脫毛在動物測試中係不切實際時,插入皮下電極,此不僅造成額外量測困難而且亦造成感染顧慮。感染可改變體內環境恆定條件,此可使生命徵象量測產生嚴重偏差。光學體積掃描(PPG)廣泛用於醫烷及穿戴式裝置兩者中之脈搏率。PPG利用氧合血紅素位準之光學吸收之週期性改變以調變半導體雷射之強度,且需要繁重處理其信號以在詳細資訊經常丟失時擷取清晰血管脈搏波形。此外,如受限於雷射穿透深度,心率及心率可變性(HRV)之精確量測仍具挑戰性,且甚至相對於雷射束之小相對身體移動導致嚴重偏差。聲學方法(諸如基於聽診器之心音圖(PCG)及基於超音波之心回波圖)具有類似問題。轉換器大小破壞舒適性且限制穿戴式方法之連續監測。組織中之聲波之失真儘管可由體脂肪相依後處理補償,但亦使信號品質降級。通常藉由壓力或應變計量測諸如呼吸率/運作及血壓之其
他生命徵象,但來自腕帶或袖口之不適性阻礙長期使用。
一習知RF方法將一RF射束照射至胸部區域,且藉由身體表面反向散射遠場電磁(EM)波。在藉由接收器偵測之EM波上調變呼吸及心跳信號。遠場處之感測導致一些缺點:(1)空氣與人體組織之間的介電常數差異導致皮膚表面上之強反射,此意謂呼吸信號遠強於來自身體內部之信號。強呼吸信號可超過心跳信號,心跳信號之信雜比(SNR)及波形細節可歸因於有限量之能量及小幾何平均值而嚴重降級。(2)習知方法通常缺乏多頻道或多點感測,此限制其實際適用性。(3)從RF轉換器觀點而言,傳輸及接收頻帶彼此重疊,其中自干擾可使系統效能劣化。歸因於上述問題,習知RF生命徵象感測系統難以精確量測心跳波形或時間間隔。
相比之下,本發明揭示之NCS方法將RF近場範圍內之身體表面上及皮膚組織下方之機械運動調變至具有唯一數位識別(ID)之多工諧波RFID(RF識別)反向散射信號上。在NCS中,呼吸信號與具有與跳動心臟耦合之較高RF能量之心跳信號良好地隔離。Tx及Rx信號藉由諧波反向散射而廣泛分離,此改良SNR且提供標籤ID之多工。憑藉心跳之此改良信號品質以及高頻分量之進一步頻譜等化以減少取樣抖動,可以可靠地量測精確心跳時間間隔。
為針對心跳實驗性地實施NCS,感測天線經放置接近於胸部區域。具有全功能收發器之習知RF詢答器將需要來自電池或電力線之有功功率,此歸因於大小及便利性之考量而限制長期監測能力。一被動且無需維護之穿戴式裝置係較佳的。UHF RFID系統係一良好候選。在圖8A中展示RFID方案。讀取器傳輸處於頻率f之下行鏈路信號,其中傳輸器
(Tx)信號亦通過循環器至讀取器天線。藉由標籤接收之下行鏈路信號遵循電子產品碼(EPC)協定。標籤無需任何有功電源,因為其等使用下行鏈路RF信號中之能量來供電。標籤上之電荷泵採集藉由標籤天線接收之小RF能量,接著邏輯電路調變反向散射上行鏈路信號。上行鏈路信號接著被讀取器接收器(Rx)天線接收且通過循環器。然而,幾個問題將限制用於NCS目的之系統效能。(1)存在直接從下行鏈路至上行鏈路之強自干擾。如圖8A中展示,由於循環器、直接Tx天線反射及來自附近物件之周圍反射之不充分隔離,強Tx至Rx洩漏仍然存在。歸因於被動標籤之低調變率,Rx信號經受來自Tx相位雜訊邊緣之大雜訊,此限制SNR及讀取器敏感度。(2)習知RFID系統通常採用循環器作為Tx/Rx雙工器,其具有有限頻寬且因此易受室內多路徑干擾。(3)習知EPC協定利用分時多重存取(TDMA)以處置標籤衝突,其中標籤在Aloha方案中選擇一隨機延遲時間。此隨機延遲將額外孔徑抖動引入至感測信號,此進一步使SNR降級且導致信號失真。
為解決所有此等問題,可使用一諧波RFID系統。如圖8B中展示,Tx信號係處於f,從而傳輸通過低通濾波器(LPF)至Tx天線。諧波標籤由下行鏈路RF信號供電給諧波標籤,且透過Rx天線及高通濾波器(HPF)在上行鏈路中反向散射處於2f之諧波信號。Tx及Rx頻譜良好地分離,且Rx信號幾乎不受Tx相位雜訊邊緣影響。讀取器RF前端使用LPF及HPF作為雙工器,且頻寬可較寬以用於室內多路徑免疫。同時,諧波標籤亦可在CDMA協定上運行,藉由讀取器同步諧波標籤以提供具有最小孔徑抖動之多標籤存取。
在圖22A中展示諧波標籤之示意圖。天線1從讀取器接收下
行鏈路信號。信號之一部分經遞送至能量採集模組以供電給標籤控制電路,且標籤接收器解調變下行鏈路命令以進行適當空中通訊協定(air protocol)及對應邏輯運算。RF能量之其他部分耦合至非線性傳輸線(NLTL)以在上行鏈路中產生二次諧波信號,其藉由RF開關調變以輻射通過天線2。圖22B展示一例示性實施例之一PCB原型,從WISP平台修改其上之數位邏輯部分。
Ettus B200軟體定義無線電(SDR)經程式化為諧波RFID讀取器。Tx局部振盪器(LO)經設定在1GHz且Rx LO設定在2GHz,其等具有相同時脈源以維持同調諧波收發器方案。
圖23中展示一例示性生命徵象監測設定。一第一標籤在胸部區域中以感測心跳及呼吸,且第二標籤在左腕上以感測脈搏。胸帶及腕帶僅用於服裝上之便利標籤部署且無需皮膚接觸或彈性帶張力。標籤從SDR讀取器接收下行鏈路信號且將上行鏈路信號反向散射至讀取器天線。解調變NCS信號由LabVIEW處理且展示於螢幕上,其中白色曲線追蹤心跳且紅色曲線追蹤腕部脈搏。藉由具有0.8Hz至15Hz之帶通濾波器(BPF)處理波形。從心臟信號至脈搏信號之延遲係脈波傳遞時間(PPT),其可用於估計血壓。圖24A及圖24B展示從胸部RFID標籤收集之心跳及呼吸之解調變NCS信號。圖24A中之上曲線係NCS信號之正規化振幅部分。可見,主要調變係心跳,然而稍微耦合呼吸。正交接收器中之振幅與相位調變之間的隔離給出呼吸波形與心跳波形之間的明確分離,如干擾計類比中繪示。在數位BPF(0.8Hz至15Hz)之後,心跳信號經清晰擷取為下曲線。圖24B係NCS信號之正規化相位部分。原始相位波形經展示為上曲線,其中主要調變係呼吸,然而亦稍微耦合心跳。藉由數位BPF在一不同範圍
(0.1Hz至1.2Hz)中處理原始相位資料以擷取深紅色曲線中之經濾波呼吸信號。替代僅依靠先前工作中之頻帶濾波以分離呼吸信號及心跳信號,NCS可藉由正交方案區分外部及內部機械運動。
為量測心跳時間間隔,可精確偵測波形及相關聯峰間時間。時序精確性因此取決於特徵點之時間解析度及清晰度。藉由頻譜中之高頻分量更佳地反映時域中之尖銳邊緣或峰。然而,來自呼吸及主心室運動之低頻分量自然地具有非常大振幅且因此可超過用於精確時間間隔偵測之尖峰。如圖25A中展示,實曲線係NCS心跳信號且點劃曲線係ECG信號,ECG信號在SDR讀取器中同步為主鏈路。NCS之取樣率係每秒5000個樣本(Sps),且ECG之取樣率係512Sps,最初受限於電極雜訊。ECG信號經升頻轉換至5000Sps以用於同步及圖形顯示。為藉由ECG量測心跳時間間隔,可使用各週期中之QRS複合波之峰。由於峰已為尖銳的,故可容易地定位其時序。然而,NCS信號直接量測機械運動且因此較平滑而不具有明顯峰值,此減小峰值偵測精確性。為擷取精確時序,需要尖峰特徵點。因此,必須進一步處理原始NCS信號以改良峰值偵測確定性。圖25B及圖25C展示同步NCS及ECG信號之頻譜。兩個頻譜經正規化至約1Hz之各自峰值強度值。與NCS信號相比較,ECG具有更強高頻(2Hz至8Hz)分量且因此具有更尖峰。然而,由於NCS中之直接運動調變、諧波反向散射設計及較高取樣率,雜訊底限遠低於ECG之雜訊底限。然而,歸因於相同心跳源之同步量測,兩個頻譜之頻率分量經合理地對準且共用類似分佈。
NCS信號具有約4.5Hz至6Hz之豐富高頻分量,但強度振幅比約1Hz之主要心跳信號低至少10dB。相比之下,ECG亦具有該等分
量但具有更高強度及雜訊。因此,可應用一直接頻譜等化:將一4Hz至7Hz帶通有限脈衝回應(FIR)濾波器應用至原始信號,且接著將濾波器之輸出放大13dB以加回至原始信號。為使此兩個信號時域對準,濾波器可經設計為零相位濾波器以消除相位延遲。FIR結構有利地用於保持線性相頻回應。在圖26A中展示經處理信號,其與各心跳循環之原始信號相比具有更多特徵點及尖峰。各心跳循環之最高峰特徵點經選擇為時序指示符,展示為圖26A中之三角形標記。
心跳時間間隔經界定為一個循環與下一循環中之特徵點之間的時間。在圖26B中展示心跳時間間隔。針對原始NCS(虛線)、ECG(點劃線)及等化NCS(實線)信號計算心跳時間間隔。對於原始NCS及ECG信號,特徵點經選擇為最高正峰。由於基於心跳間時間計算心跳時間間隔,故水平軸係心跳指數而非時間。垂直軸係心跳時間間隔,且其倒數可被視為心率。所有三個曲線展示相同整體下降趨勢或等效地增大心率,此係因為量測係在身體鍛煉之後進行的。儘管此處未提供實況,但仍可進行以下觀察。與等化NCS心跳時間間隔相比,原始NCS歸因於不精確峰值偵測而具有較大變化;ECG信號歸因於512Sps之低取樣率而具有較大變化。在幾個心跳內之等化NCS中之更穩定心跳時間間隔變化具有生理意義。
使用相關聯天線設計策略分析天線對NCS效能之效應。當達成較佳天線阻抗匹配時,可改良能量耦合效率及SNR以偵測心電圖波形細節。
為示範NCS之操作原理,在CST Microwave Studio中建立一人體身軀電磁(EM)模擬模型,如圖27A中展示。從Zubal Phantom提取身軀之器官幾何性質及組織性質。感測天線經附接接近於胸部區域但無需接觸皮膚。所發射EM波部分傳輸至遠場,展示為灰色箭頭。RF能量之另一部分歸因於近場效應而耦合至身體中,展示為虛線箭頭。由於辯證邊界之移動,藉由心跳調變反向散射信號之相位。展示為黑色箭頭之經調變反向散射RF信號干擾直接傳輸(灰色箭頭)且被RF接收器接收。由於此兩個信號係來自相同源但具有不同路徑,故操作類似於一干擾計結構。內部心跳給出兩個信號之差分調變,且可從遠場RF量值解調變。同時,由呼吸或身體運動導致之胸部外部移動將改變兩個信號(灰色及黑色箭頭)之相位延遲(此可被稱為共同模式調變)且可從遠場RF相位解調變。可藉由正交混頻器容易且精確地分離接收器處之量值及相位資訊。
NCS之一部分係接近於外部及內部身體運動之感測天線。由於天線之近場區域中之人體組織之高介電常數,其S參數將顯著改變。從設計觀點而言,天線可被視為RF電路阻抗(通常約50Ω)與由自由空間及身軀構成之連結區域之阻抗之間的匹配組件,後者與天線相比將僅在自由空間中改變阻抗、頻率回應及輻射場型。對於圖27B中之CST模擬,當2GHz偶極天線在自由空間中操作時,天線之S11經展示為實線。然而,當該天線附接在胸部區域附近時,S11將偏移至點劃線,此意謂2GHz附近之頻率歸因於高反射而不再具有良好發射效率。存在兩個可能簡單解決策略:將操作頻率偏移至具有低S11之頻帶,或重新設計天線幾何性質以擬合原始2GHz頻帶。虛線係重新設計之天線,其與所存在身軀匹配。處於2GHz之反射從-3dB大幅改良至-18dB。
為進一步調查天線匹配效應之效能,圖28展示在一天線放置於圖28B中之胸部上時處於2GHz之功率流。在圖28A中,原始天線1具有點劃線,且圖28B之經修改天線2具有虛線。人體身軀之橫剖面經繪示有左(L)肺、右(R)肺及心臟。在相同驅動信號強度及顏色輪廓標度下,可觀察到從圖28A中之非匹配天線耦合至身軀之能量遠小於圖28B中之匹配狀況之能量。較強能量耦合亦增大反向散射信號之總強度。因此,在相同雜訊底限之情況下,將改良SNR及感測敏感度兩者。
在一例示性實施例中,為在一便利且高效能感測平台上實施NCS,利用一諧波反向散射RFID系統。NCS天線係諧波標籤之部分,諧波標籤之示意圖在圖22A中展示。標籤經設計為從讀取器採集下行鏈路RF能量以供電給標籤電路之被動裝置。藉由電荷泵採集由天線A(Ant.A)接收之RF能量之部分以操作接收器及微控制器單元(MCU)。RF能量之其他部分耦合至用於二次諧波產生之非線性傳輸線(NLTL)中以透過天線B(Ant.B)反向散射至讀取器。RF開關使用開關鍵控(OOK)將諧波信號調變為上行鏈路基頻帶。對標籤MCU執行分碼多重存取(CDMA)協定以達成多標籤案例中之較佳同步及效能。圖22B展示用於例示性諧波標籤之印刷電路板(PCB)原型,其中NCS感測天線安裝至天線B連接器。來自天線B之上行鏈路RF信號亦藉由呼吸及心跳予以調變且接著由讀取器予以接收及解調變。藉由軟體定義無線電(SDR,Ettus B200)實施諧波讀取器。讀取器傳輸器(Tx)之局部振盪器(LO)經設定在基本頻率f,且接收器(Rx)LO經設定在二次諧波頻率2f。藉由相同頻率參考驅動Tx及Rx LO之兩個合成
器,因此讀取器經組態為諧波同調收發器。使用諧波反向散射(而非習知RFID系統)之一主要益處係下行鏈路與上行鏈路之寬頻率分離(broad frequency separation),因此來自標籤之弱反向散射信號並不經受來自Tx洩漏之高相位雜訊邊緣。因此,敏感度及SNR可遠高於習知RFID方案。
在圖23中展示NCS量測設定。一個標籤在胸部區域中以獲得心跳及呼吸,且另一標籤在左腕上以獲得脈搏。胸帶及腕帶用於服裝上之便利標籤部署且無需皮膚接觸或彈性帶張力。標籤從SDR讀取器接收下行鏈路信號且將上行鏈路信號反向散射至讀取器。解調變NCS信號由LabVIEW中之0.8Hz至15Hz之帶通濾波器(BPF)處理且展示於螢幕上,其中白色曲線追蹤心跳且灰色曲線追蹤腕部脈搏。從心臟信號至脈搏信號之延遲係脈波傳遞時間(PPT),其可用於估計血壓。讀取器天線當前可距人1.5米至3米,此受限於RF功率採集及被動諧波轉換損耗。可藉由讀取器及標籤之自訂設計容易地改良該範圍。
一單極天線用作用於NCS之例示性實施例之感測天線,天線之S11在圖29A中展示為在天線在自由空間中操作時具有2.1GHz之中心頻率之實線。當NCS天線放置於胸部區域附近時,S11回應偏移至圖32A中之具有1.9GHz之中心頻率之虛線。約2.1GHz之反射現在非常高。為達成NCS中之大能量耦合及高SNR,下行鏈路信號可改變為950MHz,因此感測天線在1.9GHz下良好匹配。如圖29B中展示,實線係藉由處於1.9GHz(展示為圖29A中之頻率1之標記)之上行鏈路獲取之心跳信號。具有處於2GHz(圖29A中之頻率2)之上行鏈路之虛線及處於2.1GHz(圖29A中之頻率3)之點劃線歸因於耦合至身軀中之較少能量而具有非常弱NCS信
號。應注意,處於頻率2及頻率3之RF輻射效率仍相當高,分別在-10dB及-7dB之S11下估計為90%及80%,且因此NCS信號中之大部分降級可歸因於組織耦合之減小。在相同系統設定及提取程序(惟不同頻率除外)下進行所有量測。在相同系統雜訊底限下,較高NCS信號強度可增大SNR以復原更清晰波形細節。可沿著實線之良好匹配條件中之減小斜率觀察到可重複細節,此在具有處於2.0GHz及2.1GHz之上行鏈路之非匹配條件中幾乎不可見。
在此段落中,一例示性NCS生命徵象監測系統具備基於高頻心跳分量之一減輕方法以對抗來自身體移動之干擾。藉由使用同步ECG校準,在即時心率量測中實驗性地示範一低誤差概率。
如上文提及,心電圖(ECG)、光學體積掃描(PPG)及聲學方法(諸如聽診器心音圖及超音波心回波圖)係用於量測心率及其可變性之當前技術,其中ECG經最佳建立為臨床標準。然而,此等方法皆在穿戴式感測系統中受到大量關注。為達成ECG中之良好信號品質,電極需要藉由不適導電膠及脫毛而進行直接皮膚接觸。儘管已嘗試穿戴式ECG服裝,但仍難以用於日常服飾。當前PPG裝置廣泛用於智慧型手錶、腕帶、貼片及腕部或手指上之夾具上,但其需要一緊密接觸以避免來自周圍光之失真及相對運動。PPG信號依靠繁重處理以獲得相當清晰脈搏波形,且低頻率範圍及高頻率範圍兩者中之詳細資訊經常丟失。氧含量之感測深度受到限制,使得PPG最常應用至具有反射模式中之淺血流或透射模式中之高血濃度之身體區域,而直接心臟運動量測係難以達成的。聽診器心音圖易受內
部雜訊(呼吸及語音)及外部干擾(聲音及振動)且僅在一控制實驗室環境中係實用的。超音波心回波圖難以實施為一穿戴式裝置,此係因為其塊體轉換器需要使用阻抗匹配凝膠直接接觸皮膚表面。另外,當亦可期望同時量測呼吸率/運作及血壓時,通常需要應變計,但來自腕帶及袖口張力之不適性阻礙長期使用,且亦可破壞睡眠或晝夜節律。亦已提出若干基於RF之方法,其中作為遠場之一RF射束輻射至胸部區域以由人體反向散射。呼吸及心跳信號皆在RF載波上調變且接著被接收器天線接收。若藉由一RFID進一步調變反向散射信號,則認為即將發生之皮膚接觸對於可靠地擷取呼吸及心跳之微小皮膚移動係重要的。若未藉由無標籤情況中之一數位ID(識別)調變反向散射信號,則來自任何其他非特定反向散射之干擾可係不利的。在具有或不具有個人標籤之情況下,身體運動將導致基於RF之心率估計之嚴重不精確性。
在此段落中,一例示性NCS生命徵象監測系統具備基於高頻心跳分量之一減輕方法以對抗來自身體移動之干擾。藉由使用同步ECG校準,在即時心率量測中實驗性地示範一低誤差概率。
使用圖23中展示之實驗設定測試NCS技術之一實施例。一軟體定義無線電(SDR,Ettus B200)經程式化以用作為一諧波讀取器。在胸部及左腕區域處部署被動諧波標籤。胸帶及腕帶用於服裝上之便利標籤部署且無需皮膚接觸或彈性帶張力。NCS(量測)信號由讀取器Rx天線接收且由SDR諧波讀取器解調變。所得解調變心跳及腕部脈搏在螢幕上分別展示為淺色曲線及深色曲線。腕部脈搏相對於心跳之延遲係脈波傳遞時間
(「PTT」),其可用於估計血壓。圖31展示來自胸部標籤之解調變信號。心跳信號係來自帶通濾波(BPF:0.8Hz至15Hz)之後的量測信號之振幅。從0.1Hz至1.2Hz之BPF之後的量測信號之相位導出呼吸信號。
儘管NCS可隔離外部機械移動與內部移動(諸如來自心跳之呼吸及來自腕部脈搏之手部運動),但在身體運動之特性頻譜分量接近於生命徵象之特性頻譜分量時,量測信號仍可受到大身體移動干擾。此額外耦合類似於藉由振幅調變(AM)與頻率調變(FM)之混合之信號污染,尤其在AM及FM邊帶之一者非常大時。當包含大運動時,在圖32A中展示原始NCS量測信號。在前40秒中,受測試人員靜坐以建立清晰心跳信號之基線以供參考,如插圖中展示。在介於40s與70s之間的週期期間,人員猛烈地揮動其等之手部。在介於70s與100s之間的週期期間,人員移動其身體且在第97秒站起,接著繼續左右移動其身體。在第120秒,人員復原坐姿15s。在約第135秒,人員再次站起直至在第145秒再次靜坐。在圖32B中展示具有及不具有移動之20秒窗之頻譜。頻譜經正規化至約1.5Hz之主峰(不具有移動)之強度且放大至較低強度部分。濾出低於0.8Hz及高於15Hz之分量。淺灰色曲線係在從不具有身體移動之10s至30s之時間週期期間的頻譜,且深灰色曲線係在心跳信號受到大身體移動干擾時的85s至105s期間。不具有身體移動之頻譜不僅展示約1.5Hz之主峰,而且亦清晰展示較高頻率諧波。身體移動通常在較低頻率下具有強特性頻率。心跳之主譜峰被具有僅高於0dB之信號干擾比(「SIR」)之身體移動嚴重扭曲。然而,處於較高頻率之分量較少被身體移動污染。峰A、B及C之SIR分別為3.4dB、5.9dB及7.6dB。
由於心跳信號之高頻分量具有較佳SIR,故其等可用於心跳計數而具有來自身體移動之較少影響。NCS信號通過4Hz至5.5Hz之一BPF以擷取三次諧波之B峰,且當人員在第97秒站起時在圖33A(實線)中使用95s至99s期間的同步心電圖(ECG)(虛線)展示。三次諧波之選擇係歸因於其合理SIR及大量值。高頻分量亦與心率相乘以減小計數誤差。如在圖38中之前兩個心跳中展示,NCS中存在六個最大峰及最小峰,因此十二個總峰用於心率估計。然而,在96s至98s期間,人員站起且原始NCS信號大幅改變,同時在使用直接凝膠塗漿護皮墊之情況下之ECG信號亦失真,然而其QRS特徵歸因於其高頻特性而保持相當清晰。存在根據ECG信號之三個心跳及來自NCS之三次諧波之十七個峰。NCS遺漏一個峰,但心跳計數之誤差在此三個心跳內僅為5.6%且針對圖33A中之範圍係2.4%。
為分析即時心率,可在時域處均勻地重新取樣心跳計數。心跳計數係基於最大峰及最小峰之數目,但各峰之時間資訊並非均勻分佈為一離散信號。因此,心跳資料與三次樣條擬合且均勻分佈有0.05s之時間解析度。在圖33B中展示心跳計數對時間。實曲線及虛曲線分別為NCS及ECG計數結果,且底部曲線(標記為ERROR)係誤差。最大誤差發生在人員揮動手部時的第40秒之後。在圖33B之插圖中展示用於此週期之心跳計數曲線。在整個160s計數之後,ECG之結果係243(91.1beats/min之平均心率),且NCS係241(90.4beats/min之平均心率)。誤差係-0.8%。在圖33C中展示藉由NCS(實線)及ECG(虛線)之三次諧波監測之心率曲線。藉由0.5Hz之一低通濾波器處理曲線。在下文表I中展示基於二次諧波、三次諧波及四次諧波之峰計數(圖32B中之A、B、C)之平均心率及誤差。使用藉由ECG獲得之心率曲線及NCS之各諧波計算同調係數。NCS之二次諧
波給出歸因於低SIR之最大誤差。
與其他心跳去雜訊或計數方法相比較,當前使用之提取方法可提供即時心率,而無需計算需要資料之一長持續時間之整個頻譜以擷取合理頻率解析度。運算負載亦為小的,此可容易地在具有一基本微控制器之穿戴式裝置上執行。
長期睡眠計分在臨床設定中係非常重要的以監測患者之復原且在家中監測兒童及成人。以一具成本效益之方式,通常可藉由上身移動連同心跳及呼吸監測評估睡眠品質。取代歸因於感測器及電極之皮膚接觸而引起不適之習知多導睡眠圖(PSG),此段落呈現使用胸部區域中之一單一被動射頻(RF)識別(RFID)標籤之RF近場同調感測(NCS)而無需皮膚接觸之一例示性睡眠監測系統及方法,其中可同步提取心率、呼吸節律及運動偵測。運動分類係基於具有半監督式學習之支援向量機(SVM)。可在91.06%之測試情況中正確辨識突然身體抽搐、搖動及轉身。亦在運動人為誤差校正之後改良心率偵測精確性。
本實施例試圖具有同時監測生命徵象及身體運動之能力以按一低成本且非侵入性方式進行睡眠計分,此段落呈現可基於具有近場同調感測(NCS)之被動UHF RFID系統之一睡眠計分系統,其增強耦合至身
體中之RF能量且因此增加來自身體上及身體內部之心跳、呼吸及運動之反向散射信號。可在具有唯一數位識別(ID)之無線電信號上調變導致動態介電邊界之機械移動,此可容易地延伸至藉由一單一RFID讀取器監測多個標籤及人員。在一些實施例中,待研究人員可僅穿戴可整合於織物中之一單一胸部標籤且無需皮膚接觸或運動約束,如圖34中展示。
已針對PPG及ECG信號中之運動及其他人為誤差偵測開發基於生命徵象波形之變形之信號異常偵測演算法,其中採用用於SVM之分類之時間及頻譜特性、多層(MLP)感知器、決策樹或其他分類器。運動偵測演算法係基於偵測NCS波形在搖動、轉身及突然身體抽搐之上身移動期間的改變。一SVM分類器僅使用休息時之信號特徵進行訓練且可偵測運動期間的異常特徵。
可在運動源之近場區內部署NCS之天線時執行NCS。為評估睡眠品質,將諧波RFID標籤放置於胸部區域中,此在一無線電識別及感測平台(WISP)上原型化。使用國家儀器Ettus軟體定義無線電(SDR)B200實施諧波讀取器。按每秒106個樣本實行類比轉數位及數位轉類比轉換以擷取精確基頻帶波形。讀取器傳輸器頻率f係1GHz且具有10kHz類比基頻帶;對應讀取器接收器頻率係2GHz且具有20kHz基頻帶。
執行一受控資料收集以模擬睡眠期間的以下案例:
1)靜止狀態:對象有意識地保持靜止。
2)身軀運動:5秒之輕微身軀運動。
3)身體抽搐:0.5秒至2秒之身軀及手臂之一快速、高能運動。
4)轉身:對象在1秒至2秒內左轉或右轉。
收集資料1小時,其中NCS信號取樣率係每秒500個樣本。按1分鐘之時間間隔執行模擬移動。室內多路徑效應可被忽略,因為延遲擴展遠低於NCS取樣時間。
身體之表面及內部移動(主要由測試期間的呼吸及心跳信號組成)在標籤反向散射上調變且可被解調變數位資料之信號處理擷取。NCS資料之振幅及相位兩者擷取具有不同權重之呼吸、心跳及運動資訊。相位對整個標籤運動最敏感,而振幅對小標籤運動較不敏感但係由耦合至內部運動之近場導致之天線特性之一強函數。圖35A展示休息時之原始DC濾波振幅及相位資料。
藉由簡單濾波從相位提取呼吸節律。在圖35B中,一20階巴特威士(Butterworth)濾波器用於使介於0.01Hz與1Hz之間的頻率通過。在圖35C中,藉由移除含有自主呼吸資訊之低於0.6Hz之頻率而從NCS振幅提取心跳。一線性相位濾波器用於移除高於10Hz之所有頻譜含量。
心跳信號可在輕微喘息之事件期間受到呼吸信號嚴重干擾。簡單濾波亦可在阻塞性呼吸之情況中使信號失真。
使用精確心跳及呼吸資訊進行睡眠計分之運動偵測未直接藉由簡單濾波達成。運動影響心跳及呼吸兩者之波形特徵,然而可藉由校正受運動影響之信號而擷取精確生命信號。執行一逐心跳分段以達成運動之更精細時間解析度。在具有多個峰值之非靜止心跳波形中難以進行精確
峰值偵測。本方法開始於多層級一維小波分解,使用多西貝10級小波(Daubechies 10-tap wavelet)識別在重構之後導致與心跳之最大相關性之係數。在圖35D中展示含有心跳波形之主要分量之8層級(NCS-d8)細節係數之重構波形。可在此階段執行精確峰值偵測,其中藉由每分鐘200個心跳之一最大心率判定最小峰值距離之額外約束。
基於受運動影響之波形與在休息時獲得之波形之間的差異識別運動特徵。而且,特徵有利地係穩健的以考量呼吸、心跳及信號振幅隨時間之變化。相對心跳時間間隔及相對心跳均方根(RMS)係基於不期望心跳時間間隔及RMS在心跳間顯著變化之假設之特徵。計算統計平均值、方差、偏度及峰度以擷取運動與休息時之波形之間的主要差異,其中應用五個心跳之一啟發窗。在相同窗內之0.6Hz至10Hz之範圍內之正規化頻譜功率用作除心跳諧波以外的另一特徵運動。
運動分類係基於上述七個特徵。已採用具有徑向基函數核之SVM來偵測運動。半監督式學習用於訓練模型,即,訓練僅使用休息時之自主呼吸收集之資料執行。排除用於訓練之運動資料減少來自冗餘學習之運動之過度擬合及一般化問題以及減少在訓練期間執行移動常式之不便性。圖36A及圖36B展示NCS心跳及對應運動偵測。藉由人工註釋結果評估運動偵測演算法之效能。
表II展示各運動類別之正確肯定(TP)及錯誤否定(FN)情況之數目。錯誤肯定(FP)表示在未實際發生之情況下用信號傳輸運動,此主要歸因於不規則呼吸模式及心跳偵測誤差。一逐心跳運動分類給出97.58%之精確性、88.28%之敏感度及98.10%之特異性。表III展示具有各
類別中之心跳數目之混淆矩陣。圖36C展示移除運動人為誤差之後的改良心率估計之一實例。使用30拍窗長度之一移動平均值估計心率。
可憑藉良好精確性偵測身體抽搐運動,而緩慢轉身可被錯誤分類為休息。精確性可取決於訓練資料及相關演算法。使用休息時之自主呼吸(包含具有偶爾深呼吸之規則呼吸)訓練實驗實施例。使用一規則呼吸模式訓練增大對運動之敏感度且導致不規則呼吸之情況中之增大FP。
量測動物生命徵象之當前方法通常涉及複雜且侵入性準備程序且導致受測試動物之主要不適以致其等通常需要被麻醉。例如,用於心跳波形之心電圖(ECG)需要具有良好電接觸之身體電極且因此難以應用至具有厚毛皮之哺乳動物、具有鱗甲或殼之爬行動物、具有羽毛之鳥類及
具有鱗之魚類。諸如腳底及嘴唇之裸露皮膚區域通常具有不足電信號且可對觸摸敏感。通常可使用麻醉執行肌肉電極。類似地,身體表面條件亦導致光學體積掃描(PPG)設定中之複雜性,此限制其應用至動物。聽診及超音波需要緊密皮膚接觸或阻抗匹配凝膠以獲得清晰信號,此需要高程度之動物處置。基於微小皮膚運動之多普勒(Doppler)遠場反向散射之射頻(RF)方法具有非特定無線頻道且易受任何周圍運動干擾。呼吸通常係主導信號且變為對精確心跳偵測之一主要干擾。基於傳輸線模型之RF方法再次需要皮膚電極之良好阻抗匹配。小動物歸因於有限信號敏感度而造成大部分此等RF方法之進一步挑戰。
藉由多工無線電信號上之天線特性之近場調變之本發明NCS方法無需皮膚接觸且提供有意識小動物之生命徵象之長期監測之一有效解決方案。在圖37A及圖37B中展示適合於搭配小動物使用之兩個例示性NCS設定之示意圖。圖37A中之無線感測利用諧波RFID(無線電頻率識別)架構,其使被動感測標籤變得廉價且無需維護,但需要一特定多工讀取器。此版本可適於在具有耐風雨被動標籤之自然生境中部署,且從操作者或固定裝置上之一附近讀取器收集生命徵象。諧波系統之使用藉由分離傳輸器(Tx)及接收器(Rx)之頻帶而減少自擾,改良信號敏感度以及信雜比(SNR),且減少測試RF Tx信號之所需強度,所有此等對感測小動物之生命徵象係有利的。照射於動物身體上之NCS信號可遠低於先前技術(例如0.1mW/cm2及0.15W/kg)以遵循囓齒動物模型中之健康及安全標準。在圖37A中,Ettus X310軟體定義無線電(SDR)之諧波讀取器透過讀取器Tx天線傳輸處於約950MHz之f之下行鏈路信號。下行鏈路信號供電給被動諧波RFID標籤且接著轉換為處於2f之二次諧波頻率以作為目標動物身體
之近場範圍內之NCS感測信號。只要生命信號在感測標籤天線之此近場範圍內(通常約所採用波長的三分之一)而無需任何皮膚接觸,動物身體上及動物身體內部之運動可耦合至反向散射信號以由讀取器Rx天線接收。NCS中之高SNR允許微小內部運動(諸如人類腕部脈搏波形)之精確量測,此對小動物之生命徵象係有利的。相對於讀取器之主要整體標籤天線運動將主要反映為相位調變,其可與由於相對於近場中之標籤天線之介電邊界運動之量值調變自然地分離。源自感測標籤之信號可含有一唯一識別(ID)碼以達成分碼多重存取,此改良頻道隔離以防非特定干擾且實現多個感測標籤之同時讀取。被動感測標籤之製造類似於習知RFID標籤,其除無需維護之便利部署以外亦給出低成本生產及靈活基板選擇。
替代地,圖37B中之設定用RF纜線取代讀取器至標籤頻道,此減少干擾且可適於便利地部署於具有繁重操作者訊務或具有其他干擾源之一室內實驗室中。讀取器Tx天線直接傳輸處於2f之NCS感測信號且在受測試動物之近場範圍內。接著由讀取器Rx天線接收藉由生命徵象調變之NCS信號,此可根據所考量應用進行部署。
為比較小動物心臟波形上之NCS生命徵象,首先對一囓齒動物模型執行同步NCS及ECG量測,此不僅對於臨床試驗係重要的,而且亦具有非常微弱且快速心跳,從而提供一更具挑戰測試。麻醉一朗伊萬斯(Long-Evans)實驗室大鼠(溝鼠)(編碼為#110),以藉由剃刀及Veet凝膠霜完成腹部脫毛以進行ECG電極部署,如圖38A中展示。嘗試在兩個爪子上使用彈簧夾及導電貼紙墊,但ECG信號非常弱且雜訊多。藉由一橡膠綜將具有圖38B中之架構之感測標籤天線放置於頸部區域之後方附近而無需脫毛。亦使用合理NCS信號嘗試包含沿著尾部、在胸部前方及沿著後腿之其
他天線放置。在圖38B中展示一代表性5分鐘記錄,且插圖展示第三分鐘內之半秒內之波形細節。圖38B中之ECG及NCS之心電圖波形呈現如圖38C中展示之非常類似心跳間時間間隔,但詳細特徵時序(例如NCS特徵相對於ECG之S及T特徵點之間的時序之位置)可需要進一步特性化。根據時序比較,NCS可取代ECG以進行基於心率變化(HRV)之行為研究。在圖38D中藉由具有額外0.5Hz至2.5Hz低通濾波之NCS同步收集呼吸波形,從而允許無法單獨藉由ECG達成之進一步心肺分析。
NCS量測藉由近場幾何改變調變之標籤天線特性,而ECG量測藉由微小皮膚電流感應之身體電位差,微小皮膚電流進一步由心臟電刺激及血流感應。從此觀點,NCS具有類似於心衝擊圖(BCG)之一波形且係比ECG更直接之一心跳運動量測。NCS亦較少經受來自皮膚條件及準備步驟之變化。
在確認麻醉大鼠上之心肺信號之後,接著示範圖39中之可能非侵入性NCS設定,其等對於ECG及其他習知技術對有意識小動物係非常困難的(若有可能)。名為「Timo」之一寵物金倉鼠(金倉鼠(Mesocricetus))在其籠睡眠區中受到監測,如圖39A中展示。在籠外應用圖37A及圖37B中之無線及有線NCS架構兩者。對於無線版本,藉由下行鏈路信號供電給被動諧波標籤,下行鏈路信號之部分經轉換為二次諧波以作為耦合至倉鼠身體中之NCS感測信號。讀取器Rx天線距標籤約1.5m。對於有線版本,感測天線經安裝於睡眠區之右側上,其直接傳輸NCS感測信號。在圖39B中,在倉鼠不注意之情況下獲取其呼吸及心跳之生命徵象。生命徵象之振幅經正規化至整個記錄週期之最大值。插圖展示約第八秒之心跳波形細節。波形特徵不僅類似於記錄期間的各心跳,而且亦類似
於圖38B中之麻醉大鼠中之波形特徵。在應用具有20之窗大小之移動平均值之後,針對心跳時間間隔提取波形之最小負值,如圖39C中展示。倉鼠心跳時間間隔比大鼠心跳時間間隔長約20%。為示範對小鳥之NCS適用性,量測名為「Banana」之一寵物觀賞長尾鸚鵡(虎皮鸚鵡,亦稱為相思鸚鵡),如圖39D至圖39D中展示。有意識鳥類之生命徵象監測可實現行為研究以及鳥類健康檢查中之新能力。圖39D示範有線實驗設定。諧波Tx天線傳輸NCS感測信號,NCS感測信號經耦合至長尾鸚鵡身體中。諧波Rx天線亦整合至棲木且能夠獲取具有圖3F中之經提取心跳時間間隔之如圖3E中展示之心跳及呼吸兩者之詳細特徵。天線部署對於有線NCS設定係非常便利的。在長尾鸚鵡行為之簡單觀察之後,識別鳥類通常站立之棲木上之若干位置。選擇一個最頻繁位置且安裝天線。在NCS量測期間,不干擾寵物長尾鸚鵡之晝夜節律。
在名為「Blimp」之一俄羅斯寵物龜(四爪陸龜,亦稱為霍斯菲爾德陸龜(Horsfield’s tortoise)或中亞陸龜)上進行下一示範,如圖39G至圖39I中展示。俄羅斯陸龜係達到月球軌道之第一個脊椎動物物種且歸因於其等之亞冬眠能力而對長太空旅行之研究係重要的,此使生命徵象之連續長期記錄變得甚至更具科學意義。海龜具有由小氣隙分離之一硬殼及軟身體組織之一身體結構。藉由ECG及超音波之生命徵象獲取僅在使用頸部或肌肉探針時係可行的,此兩者導致動物之主要不適性。殼之生理特徵及外溫性連同適當診斷方法之缺乏因此使龜類急救護理變得非常困難。鳥類及囓齒動物之心臟具有含有兩個心房及兩個心室之四個腔(類似於人類),且因此圖38B、圖39B及圖39E中之NCS波形具有類似特徵。相比之下,陸龜具有含有兩個心房及一個心室之三腔心臟,且NCS心電圖具
有不同特徵,然而仍可精確擷取心率。在圖39G中展示NCS設定。圖39D中之一類似天線對經放置於玻璃籠外部之雪松木片下方。歸因於身體結構,呼吸及心跳信號皆嵌入圖39H中之NCS振幅之原始資料中,此係因為呼吸期間不存在殼表面移動。心跳信號可清晰(如插圖中之淺陰影區域中展示),但在呼吸週期期間,其將被歸因於肺之相對較大體積之強呼吸信號(深陰影區域)覆蓋。為分離兩個重疊信號以進行精確率估計,採用連續小波變換(CWT)以提取圖39I中之兩個波形之峰值特徵。亦已嘗試將天線直接放置於椎殼及腹甲上,其中亦可清晰記錄NCS心跳及呼吸信號。
最後但非最不重要的,RF信號可在水中之一短範圍內工作,尤其近場耦合。儘管ECG遙測術係可行的,但儀器植入程序僅適合於在可執行無偏量測之前具有長達一周之復原週期之較大魚類。魚類學者因此長期探索其他非侵入性解決方案以擷取小魚中之行為研究及進化生物學之生理資訊。在圖39J至圖39L中示範名為「Glee」之一寵物五彩博魚(鬪魚,亦稱為暹羅鬪魚)之生命徵象量測。解調變NCS相位信號經解釋為圖39K中之胸鰭移動,其具有含有藉由非同步視訊記錄確認之頻率之一週期性波形。據信,圖39L中之解調變NCS量值信號表示源自心跳之週期性波形。
如圖37A及圖37B中展示,藉由一軟體定義無線電(SDR)執行諧波讀取器。使用Ettus USRP X310及UBX 160MHz RF子板之組合。在吾人之實驗期間,下端Ettus USRP B200/210亦可經組態為用於此應用之諧波讀取器,但X310具有用於波形細節之高解析度之較高資料取樣率。為運作為圖37A及圖37B中之一同調諧波收發器,Tx鏈及Rx鏈共用相
同RF時脈源。Rx鏈中之合成器將局部振盪器(LO)頻率組態為Tx鏈中之OL頻率之兩倍。藉由場可程式化閘極陣列(FPGA)產生Tx基頻帶信號,且將Rx基頻帶饋送至FPGA中進行解調變以進行資料記錄及顯示。
圖37A及圖37B中展示之無線及有線版本兩者可應用於NCS。在圖37A中之無線版本中,在修改無線識別及感測平台(WISP)及非線性傳輸線(NLTL)之情況下設計被動諧波標籤。NLTL係電感器及變容器之梯形結構,因此使用變容器符號來表示諧波標籤及諧波產生器。在所採集RF功率下操作之諧波標籤可容易地大量部署於許多感測目標上或相同目標之多個點上。同時,對於單點監測,圖37B中之有線版本較易於應用至室內動物實驗室。此處之諧波產生器仍可係如以前的NLTL以提供高轉換效率。NLTL之輸入埠處之低通濾波器隔離從NLTL至讀取器之直接諧波反射,且NLTL之輸出埠處之高通濾波器將基本頻率信號衰減至Tx天線。然而,在不具有如被動標籤中之嚴功率約束之情況下,具有適當頻率回應之任何主動或被動倍頻器可取代有線版本中之NLTL。在實驗期間,嘗試將自訂二極體及商用被動倍頻器(CRYSTEK CPPD-0.85-2)作為諧波產生器,此兩者提供滿意效能。亦可使用讀取器CDMA技術將有線NCS系統延伸至多個點,但系統成本將按比例增大。在吾人之先前工作中示範諧波系統之益處。作為一簡要概述,諧波系統提供Tx鏈與Rx鏈之間的非常好隔離,因此Rx雜訊底限可非常低以改良SNR及讀取器Rx敏感度兩者。繼而,Tx功率亦可非常低以仍維持生命徵象感測所需之SNR,此消除關於照射至活體組織之RF功率之任何進一步健康顧慮。
藉由具有LabVIEW介面之電腦控制SDR。DAC(數位轉類
比轉換器)及ADC(類比轉數位轉換器)之取樣率經組態為皆處於10MSps。來自DAC之基頻帶輸出之頻率係1MHz。當Tx LO頻率係950MHz時,從Tx輸出之信號係951MHz。在諧波轉換之後,Rx信號之中心頻率將為1902MHz。處於1900MHz之Rx LO經設定為Tx LO之兩倍。因此,Rx基頻帶頻率係2MHz以由ADC進行取樣。數位化基頻帶信號經降頻轉換且將LabVIEW中之取樣減少至5kSps之NCS取樣率。
儘管已關於一或多個特定實施例描述本發明,但將理解,可在不脫離本發明之精神及範疇的情況下製作本發明之其他實施例。因此,本發明被視為僅受隨附發明申請專利範圍及其合理解釋限制。
Claims (20)
- 一種用於量測一個體之運動之近場同調感測系統(near-field coherent sensing system),其包括:一第一標籤,其包括用於使用由該第一標籤接收之一下行鏈路信號產生一第一感測信號之一第一信號源及與該第一信號源電通信之一第一近場天線,其中該第一近場天線經組態以經安置於待量測之一第一運動之一第一近場耦合範圍內以產生一第一量測信號作為藉由該第一運動調變之該第一感測信號;及一接收器,其用於偵測該第一量測信號,其中該第一近場耦合範圍約小於該第一近場天線之一個波長,及其中由該第一信號源所產生之該第一感測信號之一頻率係由該第一標籤所接收之該下行鏈路信號之一頻率之一諧波(harmonic)。
- 如請求項1之系統,其中該第一運動係一心臟運動(heart motion)、一脈搏、一呼吸運動、一腸道運動、或一眼睛運動。
- 如請求項1之系統,其中該第一感測信號係一ID調變波(ID-modulated wave)、一主動無線電鏈路(active radio link)或一反向散射RFID鏈路(backscattering RFID link)。
- 如請求項1之系統,其中該接收器係該第一標籤之一部分。
- 如請求項1之系統,其中該第一感測信號之該頻率係該下行鏈路信號之該頻率之第一諧波。
- 如請求項1之系統,其中該第一標籤使用一正交ID調變該下行鏈路信號,使得該第一感測信號係一CDMA信號。
- 如請求項1之系統,其進一步包括與該接收器通信之一濾波器,其中該濾波器經組態以解調變(demodulate)及濾波該第一量測信號以獲得一第一運動信號。
- 如請求項7之系統,其中該濾波器係經程式化以取樣、解調變、及濾波該第一量測信號以導出運動之一處理器。
- 如請求項1之系統,其進一步包括:一第二標籤,其包括用於產生一第二感測信號之一第二信號源及與該第二信號源電通信之一第二近場天線,其中該第二近場天線經組態以經安置於待量測之一第二運動之一第二近場耦合範圍內以產生一第二量測信號作為藉由該第二運動調變之該第二感測信號;及其中該接收器進一步經組態以偵測該第二量測信號,其中該第二近場耦合範圍約小於該第二近場天線之一個波長。
- 如請求項9之系統,其進一步包括用於基於所偵測之該第一量測信號及該第二量測信號量測一導數值(derivative value)之一處理器。
- 如請求項10之系統,其中該導數值係一血壓。
- 如請求項1之系統,其中該第一感測信號之該頻率係該下行鏈路信號之該頻率之第二諧波或該下行鏈路信號之該頻率之第三諧波。
- 如請求項1之系統,其中該第一標籤經組態以整合至一織物(fabric)。
- 如請求項1之系統,其中該接收器經組態以自該第一量測信號及該第一感測信號隔離(isolate)一差分模式(differential mode)及一共同模式(common mode)。
- 如請求項1之系統,其中該第一標籤經組態以安置於水中。
- 如請求項1之系統,其中該第一近場天線包括(i)經組態以接收來自一讀取器之該下行鏈路信號之該第一標籤之一第一天線,及(ii)經組態以將一上行鏈路信號反向散射至該第一標籤之該讀取器之一第二天線。
- 如請求項16之系統,其中該第一標籤之該第一天線經組態以自該下行鏈路信號採集(harvest)能量而為該第一標籤提供功率。
- 如請求項17之系統,其中該第一標籤進一步包括: 一微控制器;一電荷泵(charge pump),其經組態以自該下行鏈路信號採集能量之一第一部分以操作一RF接受器及該微控制器;及一非線性傳輸線(nonlinear transmission line),其經組態以耦合來自該下行鏈路信號之能量之一第二部分以用於諧波產生並經由該第一標籤之該第二天線反向散射。
- 如請求項16之系統,其中該第一標籤之該第一天線及該第一標籤之該第二天線經安裝至該第一標籤之一印刷電路板。
- 如請求項1之系統,其進一步包括一讀取器,其包括一遠場天線(far-field antenna),經安置於自該第一標籤之一遠場距離處,經組態以接收來自該第一標籤之該第一近場天線之一諧波信號。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201762521163P | 2017-06-16 | 2017-06-16 | |
| US62/521,163 | 2017-06-16 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| TW201909038A TW201909038A (zh) | 2019-03-01 |
| TWI834615B true TWI834615B (zh) | 2024-03-11 |
Family
ID=66590281
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW107121019A TWI834615B (zh) | 2017-06-16 | 2018-06-19 | 用於一個體之一身體上及/或身體內部運動之非接觸量測之方法及用於量測一個體之運動之系統 |
Country Status (1)
| Country | Link |
|---|---|
| TW (1) | TWI834615B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR102770857B1 (ko) * | 2019-05-31 | 2025-02-24 | 아사히 가세이 가부시키가이샤 | 계측 장치, 계측 방법 및 프로그램 |
| TWI750504B (zh) * | 2019-08-29 | 2021-12-21 | 鉅怡智慧股份有限公司 | 活體偵測的方法及相關裝置 |
| TWI796828B (zh) | 2021-11-10 | 2023-03-21 | 財團法人工業技術研究院 | 天線陣列校正裝置與其方法 |
| CN116327161A (zh) * | 2023-03-29 | 2023-06-27 | 清华大学 | 一种脑灌注成像方法与装置 |
| TWI885929B (zh) * | 2024-06-03 | 2025-06-01 | 廣達電腦股份有限公司 | 血壓量測裝置 |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100060424A1 (en) * | 2008-03-19 | 2010-03-11 | Checkpoint Systems, Inc. | Range Extension and Multiple Access in Modulated Backscatter Systems |
| US20100152600A1 (en) * | 2008-04-03 | 2010-06-17 | Kai Sensors, Inc. | Non-contact physiologic motion sensors and methods for use |
| CN102307518A (zh) * | 2009-02-03 | 2012-01-04 | 雅培糖尿病护理公司 | 分析物传感器和用于传感器插入的设备 |
| US8286887B2 (en) * | 2009-03-10 | 2012-10-16 | Wal-Mart Stores, Inc. | RFID tag sensors and methods |
| US20130030257A1 (en) * | 2010-05-14 | 2013-01-31 | Kai Medical, Inc. | Systems and methods for non-contact multiparameter vital signs monitoring, apnea therapy, apnea diagnosis, and snore therapy |
| CN103198283A (zh) * | 2013-04-23 | 2013-07-10 | 复旦大学 | 一种谐波射频识别系统 |
| US20130237798A1 (en) * | 2006-09-21 | 2013-09-12 | Noninvasive Medical Technologies, Inc. | Method of processing thoracic reflected radio interrogation signals |
| US8665098B2 (en) * | 2010-09-20 | 2014-03-04 | Industrial Technology Research Institute | Non-contact motion detection apparatus |
| US9314648B2 (en) * | 2011-12-23 | 2016-04-19 | Texas Tech University System | System, method and apparatus for tracking targets during treatment using a radar motion sensor |
| US9477812B2 (en) * | 2007-07-12 | 2016-10-25 | University Of Florida Research Foundation, Inc. | Random body movement cancellation for non-contact vital sign detection |
| CN106102566A (zh) * | 2014-01-10 | 2016-11-09 | 马尔西奥·马克·阿布雷乌 | 用于测量abreu脑热通道的红外输出的装置 |
| US20160338798A1 (en) * | 2015-05-21 | 2016-11-24 | Drexel University | Passive rfid based health data monitor |
-
2018
- 2018-06-19 TW TW107121019A patent/TWI834615B/zh active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130237798A1 (en) * | 2006-09-21 | 2013-09-12 | Noninvasive Medical Technologies, Inc. | Method of processing thoracic reflected radio interrogation signals |
| US9477812B2 (en) * | 2007-07-12 | 2016-10-25 | University Of Florida Research Foundation, Inc. | Random body movement cancellation for non-contact vital sign detection |
| US20100060424A1 (en) * | 2008-03-19 | 2010-03-11 | Checkpoint Systems, Inc. | Range Extension and Multiple Access in Modulated Backscatter Systems |
| US20100152600A1 (en) * | 2008-04-03 | 2010-06-17 | Kai Sensors, Inc. | Non-contact physiologic motion sensors and methods for use |
| CN102307518A (zh) * | 2009-02-03 | 2012-01-04 | 雅培糖尿病护理公司 | 分析物传感器和用于传感器插入的设备 |
| US8286887B2 (en) * | 2009-03-10 | 2012-10-16 | Wal-Mart Stores, Inc. | RFID tag sensors and methods |
| US20130030257A1 (en) * | 2010-05-14 | 2013-01-31 | Kai Medical, Inc. | Systems and methods for non-contact multiparameter vital signs monitoring, apnea therapy, apnea diagnosis, and snore therapy |
| US8665098B2 (en) * | 2010-09-20 | 2014-03-04 | Industrial Technology Research Institute | Non-contact motion detection apparatus |
| US9314648B2 (en) * | 2011-12-23 | 2016-04-19 | Texas Tech University System | System, method and apparatus for tracking targets during treatment using a radar motion sensor |
| CN103198283A (zh) * | 2013-04-23 | 2013-07-10 | 复旦大学 | 一种谐波射频识别系统 |
| CN106102566A (zh) * | 2014-01-10 | 2016-11-09 | 马尔西奥·马克·阿布雷乌 | 用于测量abreu脑热通道的红外输出的装置 |
| US20160338798A1 (en) * | 2015-05-21 | 2016-11-24 | Drexel University | Passive rfid based health data monitor |
Also Published As
| Publication number | Publication date |
|---|---|
| TW201909038A (zh) | 2019-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US12290338B2 (en) | Methods and systems for electromagnetic near-field coherent sensing | |
| US12226190B2 (en) | System and method for monitoring cardiorespiratory parameters | |
| TWI834615B (zh) | 用於一個體之一身體上及/或身體內部運動之非接觸量測之方法及用於量測一個體之運動之系統 | |
| CN101489478B (zh) | 用于监视生理症状的装置、系统和方法 | |
| Güler et al. | Theory and applications of biotelemetry | |
| US20170258366A1 (en) | Apparatus and method for continuous noninvasive measurement of respiratory function and events | |
| Hafner et al. | Fish heart rate monitoring by body-contact Doppler radar | |
| Zhou et al. | A wireless wearable RF sensor for brumation study of chelonians | |
| Wang et al. | Human physiology and contactless vital signs monitoring using camera and wireless signals | |
| Conroy et al. | Heart id: Biometric identification using wearable mimo rf heart sensors | |
| Hui et al. | Mitigation of body movement interference in near-field coherent sensing for heartrate monitoring | |
| Hui et al. | Accurate extraction of heartbeat intervals with near-field coherent sensing | |
| US20240306994A1 (en) | Wearable health monitoring device | |
| Hui | Harmonic RF Sensing: From Indoor Localization to Vital-Sign Monitoring | |
| Zhou | Radio Near-field Motion Sensing: Theory and Applications | |
| Nie et al. | Is Continuous Heart Rate Monitoring of Livestock a Dream or Is It Realistic? |