CH706854A2 - Verfahren und Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien. - Google Patents
Verfahren und Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien. Download PDFInfo
- Publication number
- CH706854A2 CH706854A2 CH14342012A CH14342012A CH706854A2 CH 706854 A2 CH706854 A2 CH 706854A2 CH 14342012 A CH14342012 A CH 14342012A CH 14342012 A CH14342012 A CH 14342012A CH 706854 A2 CH706854 A2 CH 706854A2
- Authority
- CH
- Switzerland
- Prior art keywords
- impedance
- current
- voltage
- materials
- monitoring
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 149
- 239000000463 material Substances 0.000 title claims abstract description 43
- 238000012360 testing method Methods 0.000 claims abstract description 55
- 238000002847 impedance measurement Methods 0.000 claims abstract description 30
- 238000004458 analytical method Methods 0.000 claims abstract description 20
- 239000006185 dispersion Substances 0.000 claims abstract description 11
- 239000007788 liquid Substances 0.000 claims abstract description 11
- 238000004611 spectroscopical analysis Methods 0.000 claims abstract description 8
- -1 tissues Substances 0.000 claims abstract description 7
- 239000003814 drug Substances 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 43
- 239000013598 vector Substances 0.000 claims description 41
- 230000000638 stimulation Effects 0.000 claims description 35
- 238000012544 monitoring process Methods 0.000 claims description 25
- 230000008569 process Effects 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000012512 characterization method Methods 0.000 claims description 16
- 238000005259 measurement Methods 0.000 claims description 15
- 238000011156 evaluation Methods 0.000 claims description 12
- 238000004891 communication Methods 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 11
- 230000010355 oscillation Effects 0.000 claims description 10
- 238000002593 electrical impedance tomography Methods 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 8
- 239000008280 blood Substances 0.000 claims description 7
- 210000004369 blood Anatomy 0.000 claims description 7
- 239000002775 capsule Substances 0.000 claims description 7
- 210000004027 cell Anatomy 0.000 claims description 7
- 238000005260 corrosion Methods 0.000 claims description 7
- 230000007797 corrosion Effects 0.000 claims description 7
- 239000000839 emulsion Substances 0.000 claims description 7
- 238000001453 impedance spectrum Methods 0.000 claims description 7
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 239000012528 membrane Substances 0.000 claims description 7
- 239000004753 textile Substances 0.000 claims description 7
- 239000012620 biological material Substances 0.000 claims description 6
- 230000033228 biological regulation Effects 0.000 claims description 6
- 239000003990 capacitor Substances 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 claims description 5
- 239000008103 glucose Substances 0.000 claims description 5
- 235000013372 meat Nutrition 0.000 claims description 5
- 239000000853 adhesive Substances 0.000 claims description 4
- 230000001070 adhesive effect Effects 0.000 claims description 4
- 230000010261 cell growth Effects 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 201000009030 Carcinoma Diseases 0.000 claims description 3
- 206010024612 Lipoma Diseases 0.000 claims description 3
- 206010039491 Sarcoma Diseases 0.000 claims description 3
- 230000023555 blood coagulation Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 208000031513 cyst Diseases 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000001727 in vivo Methods 0.000 claims description 3
- 210000004243 sweat Anatomy 0.000 claims description 3
- 230000002110 toxicologic effect Effects 0.000 claims description 3
- 231100000027 toxicology Toxicity 0.000 claims description 3
- 210000002700 urine Anatomy 0.000 claims description 3
- 238000011161 development Methods 0.000 claims description 2
- 230000018109 developmental process Effects 0.000 claims description 2
- 230000001360 synchronised effect Effects 0.000 claims description 2
- 238000003908 quality control method Methods 0.000 claims 8
- 239000000126 substance Substances 0.000 claims 8
- 239000000203 mixture Substances 0.000 claims 6
- 239000000047 product Substances 0.000 claims 6
- 241000251468 Actinopterygii Species 0.000 claims 4
- 230000000694 effects Effects 0.000 claims 4
- 235000013305 food Nutrition 0.000 claims 4
- 230000007246 mechanism Effects 0.000 claims 4
- 210000000056 organ Anatomy 0.000 claims 4
- 125000006850 spacer group Chemical group 0.000 claims 4
- 210000000130 stem cell Anatomy 0.000 claims 4
- 238000003860 storage Methods 0.000 claims 4
- 238000012384 transportation and delivery Methods 0.000 claims 4
- 244000005700 microbiome Species 0.000 claims 3
- 235000001674 Agaricus brunnescens Nutrition 0.000 claims 2
- 241000894006 Bacteria Species 0.000 claims 2
- 241000167854 Bourreria succulenta Species 0.000 claims 2
- 241000283153 Cetacea Species 0.000 claims 2
- 241001481833 Coryphaena hippurus Species 0.000 claims 2
- 241000195493 Cryptophyta Species 0.000 claims 2
- 241000192700 Cyanobacteria Species 0.000 claims 2
- 102000004190 Enzymes Human genes 0.000 claims 2
- 108090000790 Enzymes Proteins 0.000 claims 2
- 238000004566 IR spectroscopy Methods 0.000 claims 2
- 108090000862 Ion Channels Proteins 0.000 claims 2
- 102000004310 Ion Channels Human genes 0.000 claims 2
- 235000007688 Lycopersicon esculentum Nutrition 0.000 claims 2
- 241001465754 Metazoa Species 0.000 claims 2
- 241000234295 Musa Species 0.000 claims 2
- 206010030113 Oedema Diseases 0.000 claims 2
- 240000003768 Solanum lycopersicum Species 0.000 claims 2
- 210000004381 amniotic fluid Anatomy 0.000 claims 2
- 235000021015 bananas Nutrition 0.000 claims 2
- 230000004888 barrier function Effects 0.000 claims 2
- 235000013405 beer Nutrition 0.000 claims 2
- 239000002551 biofuel Substances 0.000 claims 2
- 238000011138 biotechnological process Methods 0.000 claims 2
- 210000004556 brain Anatomy 0.000 claims 2
- 230000024245 cell differentiation Effects 0.000 claims 2
- 230000019522 cellular metabolic process Effects 0.000 claims 2
- 238000006243 chemical reaction Methods 0.000 claims 2
- 235000019693 cherries Nutrition 0.000 claims 2
- 238000000502 dialysis Methods 0.000 claims 2
- 239000003651 drinking water Substances 0.000 claims 2
- 235000020188 drinking water Nutrition 0.000 claims 2
- 235000013601 eggs Nutrition 0.000 claims 2
- 238000001839 endoscopy Methods 0.000 claims 2
- 239000003344 environmental pollutant Substances 0.000 claims 2
- 238000006911 enzymatic reaction Methods 0.000 claims 2
- 239000004744 fabric Substances 0.000 claims 2
- 235000019688 fish Nutrition 0.000 claims 2
- 235000013312 flour Nutrition 0.000 claims 2
- 238000001506 fluorescence spectroscopy Methods 0.000 claims 2
- 235000011389 fruit/vegetable juice Nutrition 0.000 claims 2
- 239000004463 hay Substances 0.000 claims 2
- 235000012907 honey Nutrition 0.000 claims 2
- 229930195733 hydrocarbon Natural products 0.000 claims 2
- 150000002430 hydrocarbons Chemical class 0.000 claims 2
- 239000007943 implant Substances 0.000 claims 2
- 238000000338 in vitro Methods 0.000 claims 2
- 238000011065 in-situ storage Methods 0.000 claims 2
- 239000013067 intermediate product Substances 0.000 claims 2
- 150000002632 lipids Chemical class 0.000 claims 2
- 239000003550 marker Substances 0.000 claims 2
- 238000010327 methods by industry Methods 0.000 claims 2
- 238000000386 microscopy Methods 0.000 claims 2
- 239000008267 milk Substances 0.000 claims 2
- 235000013336 milk Nutrition 0.000 claims 2
- 210000004080 milk Anatomy 0.000 claims 2
- 239000002105 nanoparticle Substances 0.000 claims 2
- 230000003287 optical effect Effects 0.000 claims 2
- 210000000496 pancreas Anatomy 0.000 claims 2
- 230000008635 plant growth Effects 0.000 claims 2
- 231100000719 pollutant Toxicity 0.000 claims 2
- 108090000623 proteins and genes Proteins 0.000 claims 2
- 102000004169 proteins and genes Human genes 0.000 claims 2
- 239000002994 raw material Substances 0.000 claims 2
- 239000004460 silage Substances 0.000 claims 2
- 239000003381 stabilizer Substances 0.000 claims 2
- 235000013311 vegetables Nutrition 0.000 claims 2
- 238000012795 verification Methods 0.000 claims 2
- 239000002351 wastewater Substances 0.000 claims 2
- 235000014101 wine Nutrition 0.000 claims 2
- 241000251730 Chondrichthyes Species 0.000 claims 1
- 206010028980 Neoplasm Diseases 0.000 claims 1
- 241000282898 Sus scrofa Species 0.000 claims 1
- 241000700605 Viruses Species 0.000 claims 1
- 230000003044 adaptive effect Effects 0.000 claims 1
- 201000011510 cancer Diseases 0.000 claims 1
- 230000021164 cell adhesion Effects 0.000 claims 1
- 230000008878 coupling Effects 0.000 claims 1
- 238000010168 coupling process Methods 0.000 claims 1
- 238000005859 coupling reaction Methods 0.000 claims 1
- 238000005520 cutting process Methods 0.000 claims 1
- 230000002380 cytological effect Effects 0.000 claims 1
- 238000002347 injection Methods 0.000 claims 1
- 239000007924 injection Substances 0.000 claims 1
- 230000008611 intercellular interaction Effects 0.000 claims 1
- 201000004792 malaria Diseases 0.000 claims 1
- 235000015277 pork Nutrition 0.000 claims 1
- 238000001566 impedance spectroscopy Methods 0.000 abstract description 8
- 239000007787 solid Substances 0.000 abstract description 5
- 239000007789 gas Substances 0.000 abstract description 4
- 230000001419 dependent effect Effects 0.000 abstract description 3
- 239000010410 layer Substances 0.000 abstract description 2
- 239000011159 matrix material Substances 0.000 description 16
- 238000013461 design Methods 0.000 description 9
- 238000013459 approach Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000000157 electrochemical-induced impedance spectroscopy Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 210000001519 tissue Anatomy 0.000 description 5
- 238000003745 diagnosis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 238000010561 standard procedure Methods 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- 206010006187 Breast cancer Diseases 0.000 description 2
- 208000026310 Breast neoplasm Diseases 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 238000004146 energy storage Methods 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- 230000004199 lung function Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 240000008213 Brosimum alicastrum Species 0.000 description 1
- 206010009944 Colon cancer Diseases 0.000 description 1
- 206010027476 Metastases Diseases 0.000 description 1
- 206010060862 Prostate cancer Diseases 0.000 description 1
- 208000000236 Prostatic Neoplasms Diseases 0.000 description 1
- 206010037423 Pulmonary oedema Diseases 0.000 description 1
- 208000000453 Skin Neoplasms Diseases 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 239000000090 biomarker Substances 0.000 description 1
- 238000001574 biopsy Methods 0.000 description 1
- 210000001124 body fluid Anatomy 0.000 description 1
- 239000010839 body fluid Substances 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 208000029742 colonic neoplasm Diseases 0.000 description 1
- 238000011960 computer-aided design Methods 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 210000003722 extracellular fluid Anatomy 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000005534 hematocrit Methods 0.000 description 1
- 210000002977 intracellular fluid Anatomy 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 230000000144 pharmacologic effect Effects 0.000 description 1
- 239000006187 pill Substances 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 210000002307 prostate Anatomy 0.000 description 1
- 208000005333 pulmonary edema Diseases 0.000 description 1
- 235000005828 ramon Nutrition 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 210000003491 skin Anatomy 0.000 description 1
- 201000000849 skin cancer Diseases 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/028—Circuits therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/026—Dielectric impedance spectroscopy
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Die Erfindung betrifft ein spektroskopisches Verfahren und eine Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Testobjekten (TO) mittels Messung des frequenzabhängigen Wechselstromwiderstandes der zu charakterisierenden Materialien, wie Festkörper, Gewebe, Flüssigkeiten, Dispersionen, Schichten (auch Ionen-Schichten) oder Gase. Die Ersatzimpedanzen bzw. Ersatzwechselstromwiderstände der zu charakterisierenden Testobjekte bestehen aus bekannten, beliebig bzw. geeignet verschalteten, Grundelementen, wie Widerständen, Induktivitäten und Kapazitäten. Das Verfahren zur Impedanzmessung bzw. Impedanzspektrosokopie ist ein digitales modellgestütztes Messverfahren zur Charakterisierung der elektrischen Eigenschaften von Testobjekten, welches erfindungsgemäss eine wesentlich einfachere Hardware als die derzeit erforderliche ermöglicht sowie wesentlich kürzere Analysezeiten bei gleicher bzw. sogar höherer Genauigkeit als die bislang bekannten Verfahren. Die Erfindung betrifft auch Anwendungen des Verfahrens bzw. Verwendungen der Vorrichtung, insbesondere auf den Gebieten der Medizin und Biologie. Beim erfindungsgemässen Verfahren werden am Testobjekt (TO) über Elektroden Ströme angelegt. Der Spannungsabfall (u(t)) über dem Testobjekt (TO) und der das Testobjekt durchfliessende Strom (i(t)) werden vom Messgerät (3) bzw. den Messgeräten (2a, 2b) gemessen.
Description
[0001] Die Erfindung betrifft ein spektroskopisches Verfahren und eine Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Testobjekten mittels Messung des frequenzabhängigen Wechselstromwiderstandes der zu charakterisierenden Materialien, wie Festkörper, Gewebe, Flüssigkeiten, Dispersionen, Schichten (auch Ionen-Schichten) oder Gase. Die Ersatzimpedanzen bzw. Ersatzwechselstromwiderstände der zu charakterisierenden Testobjekte bestehen aus bekannten, beliebig bzw. geeignet verschalteten, Grundelementen, wie Widerständen, Induktivitäten und Kapazitäten. Die Ersatzimpedanzen werden im Folgenden vereinfachend als Impedanzen bezeichnet. Bekannt sind derartige Verfahren und Systeme zur Bestimmung der elektrischen Materialeigenschaften unter den Namen: Elektrische Impedanz-Spektroskopie (EIS), Elektrische Impedanz-Tomographie (EIT), sowie Dielektrische Spektroskopie. Diese Verfahren gehören zur Gruppe von Verfahren wie Röntgen oder Kernspintomographie und liefern vergleichbare Informationen zu Materialeigenschaften, ohne den menschlichen Körper mit Röntgenstrahlen zu belasten, wobei die Raum und Zeitauflösung noch wesentlich schlechter ist. Deshalb werden derartige Verfahren und Systeme verwendet, um biologische und nichtbiologische Materialien, die sog. Testobjekte, wie Festkörper, Gewebe, Flüssigkeiten, Dispersionen und Emulsionen in unterschiedlichsten Anwendungen zu charakterisieren [1–4]. Wie z.B. für die Analyse auf Korrosion, Charakterisierung von zellulären Wachstumsprozessen für die Zwecke der Biopsie sowie pharmakologischer, toxikologischer und allergologischer Tests, sowie die Differenzierung von geordneten und ungeordneten Zellwachstum, für nichtinvasive und minimalinvasive Diagnose von Pathologien wie Zysten, Lipomen, Karzinomen, Sarkomen, Metastasen, von Haut- Brust-, Prostata-, oder Colonkrebs sowie von Lungenwasser oder die Analyse von Körperflüssigkeiten wie z.B. Hämatokrit-Anteil im Blut oder die Blugerinnung, Glukosegehalt im Schweiss, Urin und Blut, Steuerungen von Reaktoren in der Prozesstechnik, wie beispielsweise für die Bioäthanol-Gewinnung. Ebenso gibt es zahlreiche kommerzielle Systeme, Geräte und Komponenten sowie patentierte Ideen, wie z.B. den Chip AD5933 von der Firma Analoge Devices, das Bioimpedanzspektrokop IMSPEC der Firma MEODAT, sowie das Konzept eines miniaturisierten Impedanzmessgerät in einer Pille nach [EP 115 997 A1]. Eine exemplarische Anwendung auf dem Gebiet der EIT ist beispielsweise die Analyse der Lungenfunktion von Neugeborenen [5, 6]. Dazu lässt man Ströme durch die Testobjekte fliessen, misst Strom und Spannung mit Analog/Digital-Konvertern und berechnet aus den Abtastwerten von Strom und Spannung mit einem Rechner den Wechselstromwiderstand, die sog. Impedanz.
[0002] Die Einspeisung der Ströme kann resistiv oder kapazitiv über Elektroden sowie induktiv über Spulen erfolgen. Sog. 2- bzw. 4 - Elektroden-Verfahren und -Systeme werden für EIS und Mehrelektroden Systeme und Verfahren für EIT angewendet. Der Strom wird in der Regel mit Stromquellen erzeugt, für die es ebenfalls unterschiedliche bekannte Schaltungen gibt. Das Verfahren und System zur Impedanzmessung bzw. Impedanzspektrosokopie gemäss Anspruch 1 ist ein digitales modellgestütztes Messverfahren zur Charakterisierung der elektrischen Eigenschaften von Testobjekten, welches erfindungsgemäss eine wesentlich einfachere Hardware als die derzeit erforderliche ermöglicht sowie wesentlich kürzere Analysezeiten bei gleicher bzw. sogar höherer Genauigkeit als die bislang bekannten Verfahren und Systeme. Für den Schutz von Elektroenergiesystemen werden bereits seit den 90er-Jahren erfolgreich digitale modellgestützte Verfahren zur schnellen, kostengünstigeren und sensitiven zur Bestimmung der Leitungsimpedanz eingesetzt. [EP 1 203 961 B1]. Die erfindungsgemässe Anwendung derartiger Verfahren für Schmal- und Breitband-Impedanzspektroskopie sowie -Tomographie zur Charakterisierung von biologischen und nicht-biologischen Materialien gemäss Anspruch 1–10 zur Verbesserung der Qualität der spektroskopischen und tomographischen Messungen bzw. Analyse hinsichtlich Sensitivität und Verarbeitungszeit in der Medizin, Pharmazie, Chemie, Bio-und Umwelttechnologie sowie der Materialwissenschaften nach Anspruch 16 und 17 ist somit Hauptgegenstand. Erfindungsgemäss wird durch das Verfahren und die Vorrichtung die Raum- und Zeitauflösung der Impedanz-spektroskopischen und Impedanz-tomographischen Verfahren und Vorrichtungen erhöht. Ausserdem ergeben sich erfindungsgemäss System-Vereinfachungen nach den Ansprüchen 11 und 12, welche insbesondere bei der Miniaturisierung und Low-Cost-Produkten wirksam werden. Exemplarisch seien einige Anwendungen gemäss den Ansprüchen 16 und 17 genannt, z.B.: zur kostengünstigen Bestimmung des Wassergehaltes, der Viskosität oder der Marker-freien Identifikation von Biomarkern. Die erfindungsgemässen Weiterentwicklungen nach den Ansprüchen 1–15 werden beispielsweise zur Impedanzmessung in einer kostengünstigen, autarken, verschluckbaren Kapsel als Mittel für die Zwecke der medizinischen Diagnose und Therapie benutzt.
Stand der Technik
[0003] Um die Impedanz und insbesondere die Parameter der in ihr verschalteten Widerstände, Kapazitäten und Induktivitäten zu bestimmen, werden in der Regel der durch das Testobjekt fliessende Strom und die über dem Testobjekt anliegende Spannung zeitaufgelöst gemessen und rechnerisch in die gesuchten Werte umgewandelt. Derartige Verfahren und Systeme zur Charakterisierung der elektrischen Materialeigenschaften sind grundsätzlich bekannt, kommerziell erhältlich und wurden bereits auch patentiert, sowie publiziert. Jedoch haben die zurzeit für die Impedanzspektroskopie verwendeten Verfahren und Systeme erhebliche Limitationen bezüglich Sensitivität, Zeitaufwand und Kosten, was im Folgenden detaillierter ausgeführt wird. Zur Kurzschlussdetektion auf Stromversorgungsleitungen werden effizientere Verfahren zur Impedanzbestimmung bezüglich Sensitivität, Zeitaufwand und Kosten entwickelt, wie z.B.: in [US 2002/158 641 A1]. Das zuletzt genannte Verfahren ermittelt direkt die Parameter der Impedanz, im speziellen Fall die Induktivität und den Widerstand. Ausserdem wird in [17] eine mathematische Methode zur Offsetkompensation für die Zwecke der Kurzschlussidentifikation in elektrischen Energiesystemen vorgestellt, welche ebenfalls für die elektrische Impedanzspektroskopie und Impedanztomographie zur Anwendung nach Anspruch 5 und 9 vorgeschlagen werden. Im Unterschied zur Impedanzspektroskopie werden im Prozess der Energieverteilung entstehenden Strom- und Spannungsverläufe zur Ermittlung der gesuchten Impedanzwerte verwendet und keine Ströme in das System eingespeist, um die gesuchten Werte zu bestimmen. Um die unterschiedlichen Wechselstromwiderstände bei unterschiedlichen Frequenzen zu ermitteln, werden die Ströme und Spannungen zeitaufgelöst gemessen und rechnerisch in die gesuchten Werte umgewandelt. Die digitale Standardberechnungsmethode zur Charakterisierung der elektrischen Materialeigenschaften mit Hilfe der Impedanz-Spektroskopie basiert grundsätzlich auf der Fouriertansformation. Im ersten Schritt werden die komplexen Zeiger von Strom, Ik, und Spannung, Uk, für die gewünschten Frequenzen k·Δf nach den wohl bekannten Formeln berechnet.
[0004] Wobei Δf=l/N·Ta ist. Im zweiten Schritt werden die Zeiger nach dem Ohmschen Gesetz komplex dividiert.
[0005] Diese Methode setzt stationäre sinusförmige Signale für Strom und Spannung voraus. Während in-stationärer Übergangsvorgänge ergeben sich falsche Ergebnisse. In der heutigen Praxis kann man in Geräten, wie z.B. von MEODAT, eine sog. Einschwingzeit einstellen, die so gross sein sollte, bis Strom und Spannung stationär sind. Diese Wartezeiten sind erfindungsgemäss nicht mehr erforderlich, was die Analysezeiten verkürzt, da das Verfahren auch Modelle für instationäre Prozesse und Signale impliziert. Die Genauigkeit des Verfahrens wird durch das Signal/Rausch-Verhältnis bestimmt. Insbesondere bei Multifrequenz-Signalen ist das nachteilig, da sich bei konstant angenommen Rauschen das Signal/Rausch-Verhältnis für jede einzelne Frequenz verschlechtert, je mehr Frequenzen das Multifrequenz-Signal enthält. Um ein optimales Signal/Rausch-Verhältnis zu erhalten, beaufschlagen deshalb Präzisions-Spektroskope die Testobjekte sequentiell mit einfrequenten Oszillationen mit der erforderlichen Frequenz. Die Messung kann erst erfolgen, wenn das System im stationären Zustand eingeschwungen ist. Hierfür sind in den qualitativ hochwertigen Impedanzspektroskopen Verzögerungs- bzw. Einschwingzeiten einzustellen. Da das Verfahren erfindungsgemäss ein Multifrequenzverfahren ist, welches nicht an stationäre sinusförmige Signalformen gebunden ist, kann die Analysezeit erfindungsgemäss durch «die Beaufschlagung des Testobjektes mit geeigneten Multifrequenz-Signalen weiter signifikant verkürzt werden.
[0006] Mit der Methode [2] kann man aus der komplexen Impedanz die einzelnen Komponenten, wie Widerstände, Kapazitäten, Induktivitäten etc., aus denen die Impedanz zusammengesetzt ist, ermitteln. Dafür wird das nichtlineare Verfahren nach Levenberg/Marquart verwendet, für welches die generelle Konvergenz nicht nachgewiesen werden kann, wodurch es sich für Echtzeitanwendungen nicht eignet.
[0007] Die Methode [1] basiert auf Signalmodellen, die über die komplexe Impedanz verknüpft sind. Die Bestimmung der Impedanz erfolgt über den Umweg der Korrelation mit Sinus- und Kosinus-Signalen. Bei in-stationären Stimulationsgrössen, versagt der Algorithmus prinzipbedingt, solange bis alle in-stationären Signalkomponenten ausser den Sinus- und Kosinusanteilen abgeklungen sind. Diese Wartezeiten bis zum Abklingen der in-stationären Signalkomponenten sind erfindungsgemäss nicht mehr erforderlich, da das Verfahren, wie bereits erwähnt, auch Modelle für in-stationäre Prozesse und Signale impliziert.
[0008] Die Methode [3] verwendet zur energieeffizienten Berechnung wenige Abtastwerte von speziellen stationären Strom- und Spannungsignalen. Die Methode [4] berechnet ebenfalls das Spektrum sehr effizient. Alle diese Methoden setzen stationäre Signalmodelle voraus mit den oben beschriebenen Nachteilen, die erfindungsgemäss durch das Verfahren behoben werden. Die kommerziell verfügbaren Impedanzspektroskope sowie das in [EP 115 997 A1] beschriebene miniaturisierte System verwenden die Standardmethode und Standardhardware, bei welcher der sog. Stimulationskreis aus Stromquellen besteht, welche mit unterschiedlichen Schaltungen aufgebaut sein können, wie z.B. Bipolar nach [8] oder nach Howland [7], oder auch [9–14]. Diesen Schaltungen ist gemeinsam, dass sie die multifrequenten Ströme möglichst breitbandig und genau erzeugen. Erfindungsgemäss können durch die Anwendung der Verfahren nach Anspruch 1–10 für die elektrische Impedanzspektroskopie einfachere Schaltungen nach Anspruch 11 und 12 zur breitbandigen Stimulation und Auswertung nach Anspruch 5 ohne Genauigkeitsverluste bei gleichzeitig kürzeren Analysezeiten benutzt werden.
[0009] In [21] werden Methoden vorgestellt, die auf Signalmodellen von u und i basieren, die über die Impuls- und Sprungantwort verknüpft sind. Die Bestimmung der Impedanz erfolgt über Regression, wobei unterschiedliche Zeitkonstanten ermittelt werden. Aufgrund numerischer Eigenschaften der eingesetzten Verfahren ist die Komplexität zu untersuchenden Objekte limitiert. Ein Prototyp eines vorgestellten Gerätes zur Ermittlung der Fleischqualität konnte sich bis heute nicht kommerziell durchsetzen. Die erfindungsgemässe Kombination von stationären und in-stationären Signalkomponenten sowie mit einem Term zur Offsetkompensation nach Anspruch 4 und 9 für die Zwecke der Ansprüche 16 und 17 ist bislang in den erfindungsgemäss angeführten Anwendungen nicht bekannt.
Kurze Beschreibung der Darstellungen
[0010] Fig. 1 . stellt die Vorrichtung mit der die erfindungsgemäss vorteilhaften Verfahren angewendet werden können, sowie die Verbindungen und Lage der einzelnen Komponenten
1 Stimulationseinheit, Auswerte- und Kommunikations- und Ausgabeeinheit, sowie die Energieversorgung
2a, 2b Strommessgeräte, welche Strom und Spannung zeitaufgelöst erfassen
3 Spannungsmessgerät, welches die Spannung zeitaufgelöst erfassen,
4 Testobjekt, welches ein Festkörper, eine Flüssigkeit, Dispersion, Gas biologisches oder nichtbiologisches Gewebe sein kann
5 Elektroden welche als Klebe- oder Nadelelektroden sowie mit Mikronadeln besetzte Flächen oder Flächen, Stab-, Zylinder-, Kugel, Halbkugel-oder Plattenelektroden sowie flexiblen Textilelektroden oder auch kapazitiv (mit isolierender Schicht) ausgeführt sein können
6 Stromsensor, der resistiv, z.B. als Shunt, aktiv mit einem Transimpedanzwandler oder induktiv mit einem Transformator, einer Rogowski- oder anders geformte Luftspule ausgeführt sein kann 7a,7b, 7c, 7d Schutz- und Filtereinheit, für die Einhaltung der elektromagnetischen Verträglichkeit sowie des Abtasttheoremes garantieren, sowie die Modifikation der Signalformen von Strom und Spannung
[0011] Fig. 2 zeigt die erfindunggemäss einfache Stimulations- und Messeinheit und ihre Verbindungen untereinander.
2 Strommessgerät, welches den Strom zeitaufgelöst erfasst,
3 Spannungsmessgerät, welches die Spannung zeitaufgelöst erfassen,
4 Testobjekt, welches ein Festkörper, eine Flüssigkeit, Dispersion, Gas biologisches oder nichtbiologisches Gewebe sein kann
-Strom- und Spannungsmessgeräte (Fig. 1 ( 2 , 3 )), welche Strom und Spannung zeitaufgelöst erfassen
5 Elektroden welche als Klebe- oder Nadelelektroden sowie mit Mikronadeln besetzte Flächen oder Flächen, Stab-, Zylinder-, Kugel, Halbkugel-oder Plattenelektroden sowie flexiblen Textilelektroden oder auch kapazitiv (mit isolierender Schicht) ausgeführt sein können
6 Stromsensor, der resistiv, z.B. als Shunt, aktiv mit einem Transimpedanzwandler oder induktiv mit einem Transformator, einer Rogoiwski- oder anders geformte Luftspule ausgeführt sein kann 8 Schaltelemente, welche mechanisch als einfacher Schalter, Mikroschalter, oder elektronisch mit Transistoren, FET’s oder auch als einfach FET-Stromquellen ausgeführt sein können
[0012] Fig. 3 zeigt die Integration der Vorrichtung in eine Tablettenkapsel.
1 Elektroden, die auch als Antennen genutzt werden können
2 Energiequelle z.B. Batterie oder Spule, über die Energie eingekoppelt werden kann
3 Kommunikationseinheit inkl. Antenne, die Informationen nach aussen kommuniziert
4 mC...MikroController
5 S...Stimulationseinheit, die als H-Brücke oder mit einfachen FET-Stromquellen ausgeführt sein können
[0013] Fig. 4 zeigt exemplarisch ausgewählte elektrische Ersatzschaltbilder, mit denen Testobjekte modelliert werden können.
a) RC-Parallelschaltung mit R...Parallelwiderstand, C-Parallelwiderstand
b) RC-Serienschaltung mit R...Serienwiderstand, C-Serienwiderstand
c) RCR-Schaltung mit Rs...Serienwiderstand, Rp...Polarisationswiderstand, CH-Helmholzkapazität
d) RCR-Schaltung mit Ri... Widerstand der intrazellulären Flüssigkeit, Ri... Wider stand der extrazellulären Flüssigkeit, CM...Membrankapazität
[0014] Fig. 5 zeigt eine einfache FET-Stromquelle.
[0015] Fig. 6 zeigt welche Komponenten benötigt werden, um mit einem kommerziellen Mikrokontroller ein digitales Breitbandspektroskop aufzubauen. Das Übersichtschaltbild nur die für die Implementation eines Impedanzspektroskopes notwendigen Komponenten und deren Verbindungen.
BO binäre Ausgänge
MUX...Multiplexer
Gain...Verstärker mit variabler Verstärkung
DAC Digital/ Analog-Konverter
µC...Mikrokontroller-Core
RAM, FLASH Speicher
VREf...interne Referenzspannung
fCPU…interner Oszillator
serial out... Serielle Interfaces z.B ZigBee oder USB
Mikrokontroller-externe Komponenten: Verstärker und Antialisfilter, Spannungsversorgung mit Transformator, Graetzbrücke und DC/DC-Wandler
[0016] Fig. 7 enthält das Übersichtschaltbild mit dem die Implementation eines digitalen Breitband-Impedanzspektroskopes auf einem einzigen Chip, der die Komponenten nach Fig. 6 enthält, für Lowcost-Anwendungen. Man benötigt
Rshunt... Shunt-Widerstand
Batterie
Kommunikationsleitungen oder Radio-Chip für die drahtlose Kommunikation
Ebenso ist das komplette, funktionelle elektrische Schema offengelegt, mit dem das digitale Breitbandspektroskop aufgebaut werden kann.
[0017] Fig. 8 zeigt unterschiedliche Elektrodenausführungen
a) 2-Elektroden-Konfiguration mit 2 Nadelelektroden
b) 2-Elektroden-Konfiguration mit einer bipolaren konzentrischen Nadelelektrode
c) 2-Elektroden-Konfiguration mit 2 Elektroden in sog. Interdigitalstruktur am Boden eines Gefässes zur Zellanalyse (Im Gefäss wachsen die zu analysierenden Zellen.)
[0018] Fig. 9 zeigt unterschiedliche Elektrodenausführungen
a) Mehr-Elektroden-Konfiguration mit 4 Elektroden am Rand des Gefässes zur räumlichen Auflösung
b) 2-Elektroden-Konfiguration mit einer Elektrode im Gefässzentrum und dem Gefäss als Gegenelektrode
[0019] Fig. 10 zeigt exemplarische unterschiedliche Elektrodenausführungen, bei denen die Elektroden an der Aussen- oder Innenwand des Gefässes angebracht sind
Darstellung der Erfindung:
[0020] Es stellt sich deshalb die Aufgabe, ein Verfahren und ein System der eingangs genannten Art bereitzustellen, welches eine einfache, schnelle und genaue Möglichkeit zur Bestimmung des Impedanz-Spektrums sowie der einzelnen, die Impedanz bildenden, Parameter, wie Widerstände, Kapazitäten und Induktivitäten, erlaubt. Diese Aufgabe wird in verschiedenen Verfahrens- und Vorrichtungsausprägungen gemäss den Ansprüchen 1 bis 9 gelöst. Dazu wird der durch das Testobjekt fliessende Strom und die über dem Testobjekt anliegende Spannung zeitaufgelöst gemessen und rechnerisch in die gesuchten Werte umgewandelt. Vorzugsweise wird eine grössere Anzahl von-Messwerten verwendet, so dass ein Gleichungssystem überbestimmt wird. Die Impedanz-Parameter können sodann mit Ausgleichsrechnung ermittelt werden. Vorzugsweise werden die Fehlerquadrate bzw. die Varianzen minimiert, was mittels einfacher Matrixinversion oder rekursiv bzw. iterativ möglich ist.
Anstelle der bislang verwendeten stationären Signalmodelle werden Systemmodelle und Signalmodelle für in-stationäre Signalkomponenten verwendet. Die Berechnung der Impedanz-Parameter, wie die Ortskurve bzw. Widerstände, Kapazitäten und Induktivitäten basiert auf einem angenommenen Modell-Zusammenhang zwischen Strom und Spannung im Zeitbereich, der durch eine lineare Differenzialgleichung beschrieben werden kann.
[0021] Aus dieser Gleichung lässt sich der Zusammenhang im z-Bereich zwischen Strom und Spannung nach unterschiedlichen Methoden ableiten [16]. Die unterschiedlichen Methoden haben unterschiedliche Eigenschaften bezüglich Genauigkeit, Stabilität versus einer bestimmten Abtastfrequenz. Favorisiert ist das Verfahren der bilinearen Substitution, wobei andere Entwurfsmethoden ebenfalls erfolgreich angewendet werden können.
[0022] Aus dieser Gleichung lässt sich der Zusammenhang zwischen den Abtastwerten von Strom und Spannung im Zeitbereich [16].
[0023] In Matrixnotation ergibt sich folgende Darstellung für den Modell-Zusammenhang zwischen Strom, Spannung und den Parametern der elektrischen Ersatzschaltung.
[0024] Mu& Milegen die Modellordnung, und die die Anzahl der gesuchten Parameter (M=MU+Mj) fest, welche durch die Anzahl der elektrischen Energiespeicher in der elektrischen Ersatzschaltung und deren Verschaltung bestimmt wird. Ausserdem ist mnder Messvektor, Θ der Parametervektor und ynder Modellausgang. In diesem Fall ist der Modellausgang z.B. der aktuelle Spannungsabtastwert. Das Gleichungssystem, welches den Modellausgang, den Parameter- und den Messvektor definiert, wird im Folgenden vereinfachend als Modell bezeichnet, dessen Eingang in diesem Fall der Strom und Ausgang die Spannung ist.
Die Impedanz ist durch den Parametervektor definiert. Alternativ kann aber auch die folgende Darstellung mit dem Strom als Ausgang bzw. dem aktuellen Stromabtastwert als Modellausgang und demzufolge die Spannung als Eingang.
[0025] Alternativ kann aber der Modellausgang auch aus einer Kombination von Abtastwerten gebildet werden. Um die Stimulations- und Messkreise zu vereinfachen, kann man eine Offsetkompensation rechnerisch in das Modell integrieren. Das Verfahren kombiniert das Systemmodell für den Zusammenhang zwischen Strom und Spannung mit dem Signalmodell für das Offset. yn=ξ·1
[0026] Das Offset wird explizit durch den eingefügten Parameter ξ berechnet, aber nicht in die Berechnung der Impedanz verwendet. Man verwendet nur die Parameter α und β nicht den Parameter ξ. Dadurch ist das Berechnungsergebnis der Parameter robust gegen das Auftreten von Offsets und dessen Drift. Aufwendige Kompensationshardware oder Methoden können entfallen. Die Berechnung des Parametervektors, der die Impedanz bestimmt, erfolgt durch ein Fehlerquadratminimierungsverfahren. Zum Beispiel durch die Methode der kleinsten Quadrate von Gauss, was exemplarisch für ein Modell mit einem Spannungsausgang, weiter ausgeführt ist.
[0027] Die Berechnung des Parametervektors erfolgt nach der Formel unten.
[0028] Es existieren mehrere Alternativen zur Berechnung des Parametervektors, die unterschiedliche Rand- bzw. Prozessbedingungen berücksichtigen. Alternativ lässt sich der Parametervektor sich auch mit der gewichteten Methode der kleinsten Quadrate ermitteln. Dabei werden die einzelnen Fehler zwischen Modellausgang, Mess- und Parametervektor mit dem Wichtungsfaktor wngewichtet, was vorzugsweise nach folgender Vorschrift erfolgen kann.
[0029] Anstelle der Wichtungsfaktoren wnkönnen aber auch Wichtungsmatrizen W, die von die Diagonalform abweichen, Anwendung finden. Diese Wichtungsmatrizen können vorzugsweise über Korrelationsfunktionen gebildet werden, oder wie bei der verallgemeinerten Methode der kleinsten Quadrate (GLS...Generalized Least Square), bei der aus dem korrelierten Fehler en(farbiges Rauschen) der nichtkorrelierte Fehler εndurch ein entsprechendes autoregressives Filter (AR-Filter) erzeugt wird.
[0030] Die Erzeugung des nichtkorrelierten Fehlers erfolgt sequentiell in mehreren Schritten, wobei zunächst die Systemparameter nach Gig. (10) bestimmt werden, wobei die Modelle nach Gig. (6), (7), oder (8) angewendet werden können. Danach werden die Parameter des Filter ebenfalls mit der Methode der kleinsten Quadrate bestimmt, wobei als Modell nach Gig. (12) zu verwenden ist. Danach werden die Ein- und Ausgangsgrössen mit dem Filter F verarbeitet, was neue Ein- und Ausgangsgrössen ergibt, auf die wiederum Glg. (10) angewendet werden kann. Diese Schritte werden solange angewendet, bis sich der Parametervektor nicht mehr signifikant ändert. Das Verfahren kann aber auch mit dem später angewendeten rekursiven Verfahren nach Gig. (17) implementiert werden.
Alternativ kann aber auch ein autoregressives Signalmodell mit gleitendem Mittel und endogenen Variablen (ARMAX) verwendet werden (ELS...Extended Least Square). Dieses Verfahren kann nichtrekursiv in Schritten oder rekursiv Anwendung finden. Dazu ist folgendes Modell zu wählen.
[0031] Aus dem Parametervektor ergibt sich die Impedanz für das Modell mit Spannung als Ausgang im z-Bereich wie folgt. Aus der Impedanzbeziehung im z-Bereich kann leicht die Impedanz des Testobjektes Zmeasberechnet werden, wenn man die Beziehung zwischen z-Bereich und Frequenzbereich, wie unten gezeigt, anwendet.
[0032] Aus dem Impedanz-Verlauf im Frequenzbereich und im z-Bereich können die einzelnen Parameter der Ersatzschaltung, wie Widerstand, Kapazität und Induktivität, bestimmt werden, da die Impedanzverläufe Funktionen von den Parametern R, L und C sind.
[0033] Durch Lösen der inversen Problemstellung erhält man eine Methode, mit aus dem Parametervektor, oder den Impedanzverläufen (Z(z), Zmeas(ω)) die Parameter der Ersatzschaltung (R, L, C) ermittelt werden können, was für Schaltungen mit niederer Modellordnung besonders leicht ist.
Hierzu ist das eingangs aufgestellte Gleichungssystem entsprechend umzustellen.
[0034] Alternativ kann der Parametervektor aber auch mit der gewichteten rekursiven Methode der Fehlerquadratminimierung bestimmt werden, was eine andere Ausprägung als die Anwendung der direkten Methode der kleinsten Quadrate (OLS) darstellt. Es wird das gleiche Modell, wie bei der OLS-Methode verwendet.
[0035] Diese Form der Wichtung wird auch als Lambda-Wichtung bezeichnet. Auch die direkte Methode der kleinsten Quadrate kann unter Verwendung einer Wichtungsmatrix angewendet werden. Das Verfahren wird in der Fachliteratur auch als Wighted Least Square WLS bezeichnet [15]. In der Rekursionsformel sind E die Einheitsmatrix, 1 ein Vergessensfaktor, P die sog. Präzisionsmatrix. Der Vergessensfaktor und die Präzisionsmatrix müssen geeignet initialisiert werden.
[0036] Eine weitere Ausprägunsvariante des Verfahrens stellt die Anwendung der Total-Least-Square-Methode (TLS) zur Bestimmung des Parametervektors, für das eingangs entwickelte Modell dar. Das Verfahren kann man mit dem untenstehenden Ansatz zur Ermittlung der Parameter des Ersatzkreises implementieren. Es wird das gleiche Modell wie bei der OLS-Methode verwendet. Dazu muss zuerst muss das Gleichungssystem kann nach [1] mit der sog. Singulärwertzerlegung (SVD) rechnerisch in Zwischenergebnisse V umgewandelt werden.
[0037] In einem zweiten Schritt werden aus den Zwischenergebnissen V die Parameter 0 ermittelt, aus welchen sich dann einen Impedanzverlauf generieren lässt.
[0038] Eine weitere alternative Ausprägungsvariante des Verfahrens stellt die Anwendung eines Kaiman-Filters (Kaiman-Filter) zur Bestimmung des Parametervektors für die eingangs beschriebenen Modelle nach Gig. (6, 7, 8, 12, 13) dar.
[0039] Der Parametervektor und die Impedanz kann nach folgender Vorschrift ermittelt werden.
[0040] In dieser Vorschrift ist K das Kalman-Gain. Die Matrizen Q, R und P sind geeignet zu initialisieren. Die mit n indizierten Variablen e und w stellen die entsprechenden Fehlerterme dar.
[0041] Ein weiteres Verfahren löst ein im Zeitbereich aufgestelltes Gleichungssystem, welches den Zusammenhang zwischen Strom und Spannung und den Impedanzparametern herstellt.
[0042] Die Gleichungssysteme für u(t) und i(t) können für jede beliebige Schaltung mit den existierenden Methoden der Netzwerkberechnung aufgestellt werden.
[0043] Mit den Gleichungen für u & i für unterschiedliche Zeiten tkwird ein neues Gleichungssystem mit Hilfe der gleitenden Mittelwerte aufgestellt. Deswegen besitzt es glättende Eigenschaften. Wobei die Mittelwerte in unterschiedlichen Fenstern gebildet werden. In der Vorschrift sind die unterschiedlichen Fenster durch unterschiedliche Indizes (z.B. k, l oder kl...km) dargestellt. Im Folgenden werden die Grössen i(t) und u(t) durch die Variablen x(t) und y(t) substituiert.
[0044] Sind die Parameter (R, L, C oder deren Verknüpfungen), die nun als Θ bezeichnet werden, in Linearkombinationen im Gleichungssystem enthalten, lässt es sich leicht in Matrizenschreibweise aufstellen.
[0045] Die Parameter können nun bestimmt werden mit den folgenden Methode:
[0046] Aus den Parametern kann man die Impedanz mit dem Zusammenhang Z=f(ω, Θ) bestimmen. Da die Verfahren erfindungsgemäss nicht mehr an sinusförmige Signale gebunden sind, können anstelle der sinusförmige Oszillationen erzeugenden Hardware, erfindungsgemäss beliebige Signalformen so auch einfache Recht- oder Dreieck-förmige Signale erzeugende Schaltungen, wie die H-Schaltung, in Kombination mit Kapazitäten verwendet werden. Die Stimulation der Testobjekte kann mit Spannung oder Strom erfolgen. Einfache FET-Stromquellen sind ebenfalls einsetzbar. Das Verfahren stellt keine besonderen Genauigkeitsanforderungen an die Signalformen, was die Stimulationskreise vereinfacht.
Wesentliche Merkmale und Vorteile des Verfahrens und der Vorrichtung sind:
[0047] 1. Das Verfahren ist nicht an spezielle Signalformen gebunden; lediglich die Anzahl der Frequenzkomponenten des Spektrums des Stimulationssignales sollte der Modellordnung angepasst sein. Beispielsweise sollten vorzugsweise mindestens 2 Frequenzen für die Identifikation der Parameter eines Systems 1. Ordnung gewählt werden.
2. Das Verfahren kann sowohl für stationäre als auch für in-stationäre Signalformen angewendet werden, was die Analysezeit erheblich verkürzt.
3. Die Stimulationskreise können ohne Genauigkeitsverluste vereinfacht werden, was insbesondere für die Miniaturisierung und für niedrigere Kosten vorteilhaft ist.
4. Mit relativ einfachen Hardwaremodifikationen lassen sich unterschiedliche Multifrequenz-Signalformen erzeugen, die sich für die schnelle und genaue Breitbandimpedanzspektroskopie eignen.
5. Im Gegensatz zu den bislang existierenden Standardmethoden zur Bestimmung der (Bio)-Impedanz, die als Fehlermodell normalverteiltes Rauschen voraussetzen, können nun spezielle Prozess- bzw. Randbedingungen besser, z.B. im Modellfehler, modelliert bzw. berücksichtigt werden.
6. Die Vorrichtung zeichnet sich durch besondere Einfachheit gegenüber den bekannten Standardvorrichten aus.
[0048] Eine Signalmodell-basierte Methode für die Impedanzberechnung, welche auch bei instationären Prozessen und Signalen korrekte Impedanzmesswerte gibt, wird im Folgenden dargestellt. Sie beruht auf dem Ansatz, die instationären Signalkomponenten explizit in das Signalmodell mit aufzunehmen. Diese Komponenten, wie z.B. Exponentialglieder, werden durch die Polynome einer Taylorreihe bis zur Ordnung K-1 modelliert.
[0049] Für die Berechnung der Impedanz Zmeaswerden die Anteile der Polynome, die mit den Parametern β quantifiziert werden, nicht verwendet, sondern – analog zur Offsetkompensationsmethode – nur die stationären Oszillationsanteile für Strom und Spannung, die mit den Parametern α<<u>>in der Spannung und α<<i><>>im Strom quantifiziert werden.
[0050] Alternativ können diese Parameter auch mit dem Verfahren RLS, WLS, GLS und ELS sowie Kaiman-Filtern ermittelt werden. Die unterschiedlichen Ausprägungen unterscheiden sich im Modell durch a) die Anzahl der Polynome, welche vorzugsweise mit K=2...3 zu wählen ist, und b) die Anzahl der Oszillationen M, welche applikationsabhängig ist.
[0051] Bei Vorhandensein eines sog. BIAS dieser explizit berechnet werden und zur Korrektur verwendet werden, was nur für vorzugsweise nur weisse Rauschsignale sinnvolle Ergebnisse liefert. Das angewendete Verfahren ist als CLS bekannt. [18] Alternativ kann die Berechnung des Parametervektors eines Modells für die Impedanz zur Umgehung des BIAS-Problems aber auch mit das sog. Verfahren der Hilfsvariablen (IV...Instrumental Variables) angewendet werden. Hierzu wird eine transponierte Hilfsvariablenmatrix WIV benutzt, mit welcher der Modellfehler durch-multipliziert wird. Als Hilfsvariablen lass sich vorzugsweise die Abtastwerte der Eingangssignale oder die Schätzwerte der ungestörten Ausgangswerte verwenden. Die nicht-rekursive und rekursive Lösung der so entstandenen Gleichungssysteme ergibt das Verfahren, nach welchem der Parametervektor alternativ bestimmt werden kann.
[0052] Die oben beschriebenen Verfahren nach Ansprüchen 1 bis 10 sind für die Anwendung für die Bestimmung der Parameter von Mehrgrössensystemen prädestiniert, sog. MIMO-Systemen, was insbesondere für Impedanz-Tomographie geeignet ist. Im Prinzip wird für jeden Modellausgang ein separates Teilmodell verwendet. Dazu sind die Mess- und Parametervektoren entsprechen [19, 20] zu gruppieren. Die Lösung des Gleichungssystems wendet die oben genannten Verfahren an.
[0053] Die Vorrichtungen zur Anwendung der Verfahren gemäss Anspruch 10 bestehen grundsätzlich aus am Testobjekt angebrachten 2-, 4- oder Mehr-Elektrodensystemen, an die Spannungen angelegt oder durch die Ströme eingespeist werden, um das Testobjekt zu stimulieren bzw. anzuregen. Die Ströme und Spannungen werden zeitaufgelöst erfasst. Die Ströme und Spannungen können mit Strom- oder Spannungsquellen erzeugt werden. Durch das Verfahren werden können auch präzise Ergebnisse erzielt werden, wenn die Stimulationskreise sehr einfach gemäss Anspruch 11 und 12 ausgeführt sind. Das ermöglicht Platz- und Preisvorteile, was insbesondere für die Miniaturisierung und Einwegprodukte, sog. Disposals, von Vorteil ist. Auch können leicht Single-Chip-Lösungen aufgebaut werden.
[0054] Bei Vorrichtungen Mehrelektrodensystemen, welche insbesondere für EIT-Systeme verwendet werden, sind dann die Elektroden so am Testobjekt anzubringen, dass die Ströme die interessierenden Teile durchströmen, was mit Filz-, Klebe-, Mikronadel- oder auch Textilelektroden erfolgen kann.
Wege zur Ausführung der Erfindung:
[0055] Da das Verfahren und die Vorrichtung zahlreiche Ausführungsalternativen, die in Kombination an- und verwendet werden können,
Vorrichtungsausführungsbeispiel Gesamtsystem:
[0056] Die Vorteile dieser Ausführungsform der Vorrichtung, welche das Verfahren ermöglicht, ergeben sich aus den Ansprüchen 1, 2 und 3.
[0057] Fig. 4 zeigt exemplarische Ersatzschaltungen für die elektrisch zu charakterisierenden Testobjekte (TO, 2). Fig. 1 zeigt wie die Testobjekte, TO (2), durch Elektroden (5a...5d) z.B. nach der 4-Elektroden-Methode zu kontaktieren sind. Nach dieser-Methode müssen die Elektroden (5a) und (5d) so mit dem Testobjekt kontaktiert werden, dass über sie der Strom, der das Testobjekt an den relevanten Regionen durchströmen soll, eingespeist wird. Die Elektroden (5b) und (5c) sind so mit dem Testobjekt zu kontaktieren, dass die Spannung über den relevanten, vom Strom durchflossenen, Regionen, gemessen werden kann. Der das Testobjekt durchströmende Strom, i(t) kann durch die Messgeräte (2a) oder/und (2b) zeitaufgelöst erfasst werden, während die über dem Testobjekt abfallende Spannung, u(t), mit dem Messgerät (3) ebenfalls zeitaufgelöst erfasst wird. Grundsätzlich ist nur eines der beiden Strommessgeräte für die Strommessung erforderlich. Die Stimulations-, Auswerte- und Kommunikationseinheit (1) erzeugt den Stimulationsstrom, i(t), welcher das Testobjekt durchströmt, wenn das Testobjekt mit Strom stimuliert bzw. angeregt wird, oder die Stimulationsspannung u(t), welche am Testobjekt angelegt wird, wenn das Testobjekt mit Spannung stimuliert bzw. angeregt wird. Die Stimulation erfolgt mit Strom- oder Spannungsquellen nach [7, 8]. Das Verfahren gestattet anstelle von spannungsgesteuerten Stromquellen, welche sinusförmige Oszillationen generieren, erfindungsgemäss auch die Verwendung von einfachen FET-Konstantstromquellen aus Fig. 5 als Stimulationskreise, mit denen sich sequenziell Rechteck-förmige Ströme mit jeweils entgegengesetzter Richtung erzeugen lassen, ohne Genauigkeits- oder Zeitverluste bei der Bestimmung der Impedanz im Vergleich zur Standardmethode zu erhalten. Mit Hilfe der Filtereinheiten (7a...7d), die als passive RC-Schaltungen ausgeführt sein können, lassen sich aus den Rechtecksignalen leicht andere Signalformen, wie beispielsweise Sägezahn-ähnliche Signale, generieren, was Flexibilität im spektralen Inhalt der Multifrequenzsignale gibt.
[0058] Die Messung des Stromes kann dabei mit Sensoren (6a, 6b) erfolgen, welche mit unterschiedlichen Technologien, welche stromproportionale Grössen ausgeben wie z.B. Shunts, Transformatoren, Hall-Sensoren etc. erfolgen. Fig. 1 und 2 stellen exemplarisch 2 mögliche Positionen für die Strommessung dar. Die Stromsensoren versorgen die Messgeräte (6a und 6b) mit den zeitkontinuierlichen Eingangsgrössen.
Die Messgeräte 3 und 2a, b übermitteln an die Stimulations- und Auswerteeinheit die zeitaufgelösten Messwerte zur Ermittlung der gesuchten Werte. Die Einheiten (7a...7d) sind mögliche Schutz- und Filtereinheiten, die an den eingezeichneten Orten installiert sein können. Die Schutz- und Filtereinheiten können Komponenten zum Zwecke der elektromagnetischen Verträglichkeit und zur Einhaltung des Abtasttheorems enthalten. Die Komponenten 2, 3, 6, 7 können auch in 1 integriert sein, was insbesondere bei der miniaturisierten Impedanzspektroskop nach Anspruch erforderlich ist.
Vorrichtungsausführungsbeispiel Stimulation mit H-Brückenschaltung:
[0059] Die Vorteile dieser Ausführungsform des Verfahrens und der Vorrichtung ergeben sich aus den Ansprüchen 1, 2 und 3. Fig. 2 zeigt für den Fall, dass die Stimulation mit einer H-Brücke erfolgt, wie das Testobjekt, TO (2), durch Elektroden (5) nach der 2-Elektroden-Methode zu kontaktieren ist, so dass es von einem Strom, i(t), der durch das Messgerät (2) zeitaufgelöst erfasst wird, durchflossen wird, während die über dem Testobjekt abfallende Spannung, u(t), mit dem Messgerät (3) ebenfalls zeitaufgelöst erfasst wird. Mit den Schaltelementen (8), die in der H-Schaltung angeordnet sind, lässt sich der Strom abwechselnd in beide Richtungen in das Testobjekt einspeisen. Die Messung des Stromes kann dabei mit einem Sensor (6) erfolgen.
Vorrichtungsausführungsbeispiel Impedanzspektroskop in einer Tablettenkapsel:
[0060] Die Vorteile dieser Ausführungsform der miniaturisierten Vorrichtung, welche eine In-Vivo-Analyse für Diagnostische Zwecke ergeben sich aus Anspruch 5.
Fig. 3 zeigt schematisch, ein auf Tablettengrösse miniaturisiertes Impedanzspektroskop mit den Komponenten 1...Elektroden, 2...Energiespeicher, 3...drahtlose Kommunikationseinheit, 4...digitale Recheneinheit, 5...Stimulationseinheit mit Messgeräten zur zeitaufgelösten Erfassung der Ströme und Spannungen. Im Unterschied zu der oben dargestellte Ausführung, werden hier die elektrischen Eigenschaften der die Tablette umgebenden Materie bestimmt. Für die die Stimulationseinheit sind vorzugsweise die H-Brücke bzw. die FET-Stromquellen zu verwenden.
Vorrichtungsausführungsbeispiel Impedanzspektroskop für Lowcost-Anwendungen:
[0061] Die Vorteile dieser Ausführungsform der miniaturisierten Vorrichtung, welche eine In-Vivo-Analyse für Diagnostische Zwecke ergeben sich aus Anspruch 5.
Fig. 3 zeigt schematisch, ein auf Tablettengrösse miniaturisiertes Impedanzspektroskop mit den Komponenten 1...Elektroden, 2...Energiespeicher, 3...drahtlose Kommunikationseinheit, 4...digitale Recheneinheit, 5...Stimulationseinheit mit Messgeräten zur zeitaufgelösten Erfassung der Ströme und Spannungen. Im Unterschied zu der oben dargestellte Ausführung, werden hier die elektrischen Eigenschaften der die Tablette umgebenden Materie bestimmt. Für die die Stimulationseinheit sind vorzugsweise die H-Brücke bzw. die FET-Stromquellen zu verwenden.
Verfahrensausführungsbeispiel OLS:
[0062] Die Messgeräte (2) und (3) werden mit einer Abtastrate fA=1/TAbetrieben und liefern eine Reihe von Spannungswerten u0, u1,...uN-1und eine Reihe von Stromwerten i0, i1...iN-1. Aus diesen Werten bestimmt man die Ersatzparameter R und C für den Fall, dass das Testobjekt (4) mit einer einfachen Parallelschaltung von R & C Fig. 4 approximiert werden kann.
Die folgende Differenzialgleichung beschreibt den Zusammenhang zwischen Strom und Spannung bei einer Parallelschaltung von R und C.
[0063] Aus dieser Beziehung kann mit eingangs erwähnten Entwurf zeitdiskreter Darstellungen kontinuierlicher Prozesse via Laplace z-Transformation und bilineare Substitution s=AD(z-1)/(z+1), AD=TA/2 die zeitdiskrete Gleichung abgeleitet werden [16].
Diese Gleichung ist Grundlage für die Ermittlung der gesuchten Werte R und C. Hierzu werden N Messungen des Stromes (i0, i1...iN-1) und der Spannung (u0, u1...uN-1) benötigt, um daraus ein Gleichungssystem mit (N-1) Gleichungen für α0und α1aufzustellen, so dass die Parameter berechnet werden können.
[0064] In Matrizenschreibweise erhält man folgendes Gleichungssystem.
[0065] Dieses Gleichungssystem wird in folgenden alternative Verfahren zur Berechnung des Parametervektors verwendet: OLS, RLS, TLS, GLS, ELS, WLS und Kalmanfilter.
Aus den Parametern können dann über die Lösung des inversen Problems (R, C)=f<-1>(Z(z)) die Kapazität C und der Widerstand R ermittelt werden. Aus dem Parametervektor kann aber auch direkt die Impedanz Zmeasermittelt werden.
[0066] Der Parametervektor kann nun mit verschiedenen Methoden ermittelt werden. Vorzugsweise ist N mindestens grösser als die Modellordnung +2, in diesem Fall N>3. Die Parameterwerte können in diesem Fall mit höherer Genauigkeit ermittelt werden. Beispielsweise können die Parameter mit der Fehlerquadratminimierung nach Gauss bestimmt werden.
Verfahrensausführungsbeispiel RLS:
[0067] Ebenso können die Parameter des Ersatzkreises aber auch mit er rekursiven Methode der Fehlerquadratminimierung (RLS) bestimmt werden.
[0068] Der Vergessensfaktor sollte vorzugsweise zwischen 0,8 und 0.99 liegen und P, die sog. Präzisionsmatrix, zwischen 10<3>E und 10<5>E.
Verfahrensausführungsbeispiel TLS:
[0069] Mit dem untenstehenden Ansatz lassen sich die Parameter des Ersatzkreises ebenso mit der Total Least Square - Methode (TLS)
[0070] Zuerst muss das Gleichungssystem kann nach [1] mit der sog. Singulärwertzerlegung (SVD) rechnerisch in Zwischenergebnisse V umgewandelt werden. In einem zweiten Schritt werden aus den Zwischenergebnissen V die Parameter Θ ermittelt, welche sich dann in einen Impedanzverlauf umwandeln lassen.
Verfahrensausführungsbeispiel Kalman-Filter:
[0071] Oder die Parameter können die Parameter aber auch mit einem Kaiman-Filter bestimmt werden.
Verfahrensausführungsbeispiel Offset:
[0072] Das folgend ausgeführte Ausführungsbeispiel ermöglicht die Hardware von Impedanzspektroskopen zu vereinfachen, da die Kompensation von Offsets oder Driften per Rechnung erfolgt. Dazu wird das Offset explizit durch den eingefügten Parameter ξ, modelliert bzw. ermittelt. In der Ermittlung der Impedanz verwendet man verwendet nur die Parameter α und β nicht den Parameter ξ.
[0073] In der Ermittlung der Impedanz verwendet man verwendet nur die Parameter α und β nicht den Parameter ξ.
Verfahrensausführungsbeispiel Signalmodell Kalmanfilter:
[0074] Die folgende Ausführungsform des Kaimanfilters, welche sog. «cordische» Darstellungsformen für die Ausführung trigonometrischer Funktionen verwendet, soll beschrieben werden. Das Signalmodell für Ström und Spannung besteht aus 1 Frequenzkomponente für eine bestimmte Kreisfrequenz und einem Term, der einen linearen Trend zur Offsetkompensation modelliert. Daraus ergeben sich die unten dargestellten applikationsspezifischen Matrizen.
[0075] Möglich wären aber auch folgende applikationsspezifische Matrizen:
[0076] Die Ausführungsvorschrift, die für u und i ausgeführt werden muss, erhält dann folgendes Aussehen.
[0077] Mit den Parametervektoren für u und i kann dann die Impedanz für die bereits festgelegte Kreisfrequenz ermittelt werden.
Verfahrensausführungsbeispiel Direct Estimation:
[0078] Die RC-Parallelschaltung nach Fig. 4b : Die Strom/Spannungs-Beziehung in diesem RC-Glied stellt die Gleichung unten dar. Diese Gleichung wird für unterschiedliche Zeitpunkte tkaufgestellt.
[0079] Summiert man die Spannungen u(tk) in 2 aufeinanderfolgenden Fenstern, ergeben sich 2 unterschiedliche Gleichungen, in denen die Parameter R & C jeweils in einer Linearkombination enthalten sind. Die Fenster müssen nicht zwingend aufeinanderfolgend sein. Sie dürfen sich auch überlappen. Die Länge der Fenster für die Summation und die relative Lage der Summationsfenster zueinander, bestimmen die Filtereigenschaften des direkten Parameterestimationsalgorithmus.
[0080] Arrangiert man das Gleichungssystem in Matrizenschreibweise, lässt es sich nach dem Vektor, der die Parameter R & C enthält, durch Anwendung der Matrizeninversion umstellen.
[0081] Die Parameter können aber auch nach (43b) aus den integrierten Gleichungen (41, 42) ermittelt werden
[0082] Aus den Parametern kann man nun die Impedanz errechnen.
[0083] RC-Serienschaltung nach Fig. 4c : Das Verfahren ermittelt über eine direkte Estimation die Parametern R & C eines einfachen RC-Gliedes im Zeitbereich. Die Strom/Spannungs-Beziehung in diesem RC-Glied stellt die Gleichung unten dar. Diese Gleichung wird für unterschiedliche Zeitpunkte tkaufgestellt.
[0084] Summiert man die Spannungen u(tk) in 2 aufeinanderfolgenden Fenstern, ergeben sich 2 unterschiedliche Gleichungen, in denen die Parameter R & C jeweils in einer Linearkombination enthalten sind. Die Fenster müssen nicht zwingend aufeinanderfolgend sein. Sie dürfen sich auch überlappen. Die Länge der Integrationsfenster, der Fenster für die Summation und die relative Lage der Summationsfenster zueinander, bestimmen die Filtereigenschaften des direkten Parameterestimationsalgorithmus.
[0085] Arrangiert man das Gleichungssystem in Matrizenschreibweise, lässt es sich nach dem Vektor, der die Parameter R & C enthält, durch Anwendung der Matrizeninversion umstellen.
[0086] Die Parameter können aber auch nach (47b) aus den differenzierten Gleichungen (45, 46) ermittelt werden
[0087] Aus den Parametern kann man nun die Impedanz errechnen.
[0088] Diese beiden Verfahren sind nicht mehr an die Periodendauer der tiefsten Harmonischen gebunden. Ebenso nicht an eine spezielle Signalform. Das Beobachtungsfenster lässt sich an die Applikation anpassen. Mit diesem Algorithmus kann ein robustes Impedanzspektroskop aufgebaut werden, mit kürzerer Analysezeit als die Periodendauer der niedrigsten Oszillation.
Ausführungsbeispiel Impedanzspektroskop für Lowcost-Anwendungen:
[0089] Mit den heute kommerziell verfügbaren Technologien für das sog. «Embedded Signalprocessing and Control» auf einem einzigen Chip nach Fig. 5 , ist es möglich, ein Breitband-Impedanz-Spektroskop mit kommerziellen Mixed-Signal-Microcontrollern aufzubauen. Aus Fig. 7 sind die Vorteile ersichtlich, da funktionell nur ein Chip, eine Batterie, 2 Shunt-Widerstände sowie serielle digitale Kommunikationsverbindungen benötigt werden.
Ausführungsbeispiel Elektrischer Impedanz-Tomograph für die nichtinvasive Diagnose von Brustkrebs oder Lungenödemen:
[0090] Das System besteht vorzugsweise aus einem Mehrelektrodensystem, welches beispielsweise in einem breiten Gürtel, T-Shirt oder BH-Körbchen in Form von Textilelektroden eingearbeitet sein können. Die Stimulationskreise können nach Fig. 2 oder 5 ausgeführt sein. Die Ströme und Spannungen können mit Messgeräten zeitaufgelöst erfasst werden und werden zur Auswerteeinheit übermittelt. Das Ergebnis wird an eine Anzeigeeinheit kommuniziert, wo sie ausgegeben wird, was optisch oder auch akustisch erfolgen kann. Die Auswerteeinheit wendet die vorgängig beschriebenen Verfahren für Mehrgrössensysteme gemäss den Ansprüchen 1 bis 9 an.
Formelzeichen
[0091] u(t), i(t) Zeitwerte von Spannung, Strom un, in Abtastwerte von Spannung, Strom Uk, Ik komplexe Zeiger von Strom und Spannung Zk komplexer Wechselstromwiderstand, Impedanz Ta Abtastzeitschritt N Anzahl der Abtastwerte in einem Beobachtungsfenster Δf Frequenzschritt,- -Differenz, -Unterschied a, b Koeffizienten der Differenzialgleichung (indiziert mit a0, b0etc.) Ku Modellordnung der Spannung (max. Ordnung der Ableitung der Spannung) Ki »Modellordnung des Stromes (max. Ordnung der Ableitung der Stromes u(z), i(z) z-Transformierte von Spannung, Strom α, β Koeffizienten der Differenzengleichung (indiziert mit α0, β0etc.) Mu zeitdiskrete Modellordnung der Spannung Mi zeitdiskrete Modellordnung des Stromes Θ Parametervektor mn Messvektor, Datenvektor yn Ausgangswert, Modellausgang en, εn, wn Modellfehler, Approximationsfehler, Störung ξ Parameter für das Offset κ, γ weitere Koeffizienten des Parametervektors (indiziert mit κ0, γ0etc.) e ́ korreliertes bzw. farbiges Rauschen, z-transformiert e ́(z) ε nichtkorreliertes Rauschen, z-transformiert ε(z), indiziert εn F(z<-1>), z-Transformierte einer Filterfunktion, C(z<-><1>), D(z<-1>) z-Transformierte einer Filterfunktion, f Filterkoeffizienten (indiziert f0, etc.) f() Funktion von f<1>() Umkehrfunktion von w Wichtungskoeffizient (indiziert mit w0, wnetc.) W Wichtungsmatrix Z Impedanz z z-Operator Zmeas zu messende Impedanz des Testobjektes TO Testobjekt R, L, C Widerstand, Kapazität, Induktivität K Kaiman-Gain, P, Q, R zu initialisierende Matrizen im Kaiman-Filter E Einheitsmatrix t Zeit j imaginäre Einheit Indizes n Index für Abtastwert kl, km, l, k Index für Summen in unterschiedlichen Beobachtungsfenstern
Referenzen
[0092] [1] Seoane, F., Lindecrantz, K.: A Novel Approach for Estimation of Electrical Bioimpedance: Total Least Square. Springer, ICEBI 2007
[2] C. Polonschii et.al; Multi Frequency, MultiChannel, Differential Impedance Analyzer for Rapid Assays. Springer, ICEBI 2007
[3] A. Ronk, et.al., 2007. Simultaneous multi-frequency bio-impedance measurement applying synchronised uniform or non-uniform sampling Springer, ICEBI 2007
[4] A.Paterno et.al, Phase/Magnitude Retrieval Algorithms in Electrical Bioimpedance Spectroscopy, Springer, IFMBE, 2009
[5] Vauhkonen, M. Karjalainen, P.A. Kaipio, J.P., IEEE transactions on Biomedical Engineering, Vol.45, Issue 4, 1998, A Kaiman Filter Approach to track Fast Impedance Changes in Electrical Impedance Tomography
[6] Hussein El Dib, Andrew Tizzard and Richard Bayford, Dynamic Electrical Impedance Tomography for Neonate Lung Function Monitoring using Linear Kaiman Filter,
[7] Horowitz, Hill The Art of Electronics Cambridge, 1995
[8] Tietze Schenk, Halbleiterschaltungstechnik, Springer 1991
[9] SEOANE Fernando, BRAGÖS Ramon, LINDECRANTZ Kaj «Current Source for Multifrequency Broadband Electrical Bioimpedance Spectroscopy Systems. A Novel Approach» dans Proceedings of the 28th IEEE, EMBS Annual International Conference New York City, USA, Aug 30-Sept 3, 2006
[10] BERTEMES-FILHO P., BROWN B.H, WILSON A.J. «A comparison of modified Howland circuits as current generators with current mirror type circuits» dans Physiological Measurement. 21 (2000) 1–6
[11] CHENG Kuo-Sheng, CHEN Cheng-Yu, HUANG Min-Wei, CHEN Chien-Hung «A Multi-Frequency Current Source For Bioimpedance Application» dans
[12]: ROSS ALEXANDER S., SAULNIER G J, NEWELL J C, ISAACSON D «Current source design for electrical impedance tomography» dans Physiological Measurement. 24 (2003) 509–516
[13] CHEN Cheng-Yu, LU Yi-Yu, HUANG Wen-Lung, CHENG Kuo-Shang «The Simulation of Current Generator Design for Multi-Frequency Electrical Impedance Tomograph» dans Proceedings of the 28th IEEE EMBS Annual International Conference New York City, USA, Aug 30-Sept 3, 2006
[14]: WANG W., BRIEN M., GU D-W., YANG J. «A Comprehensive Study on Current Source Circuits» dans IFMBE Proceedings 17, pp. 213–216, 2007
[15] Carroll, R.J. Ruppert D., Transformation and Weighting in Regression, Chapman and Hall, New York, 1988
[16] Stearns, Digitale, Verarbeitung analoger Signale, Oldenbourg, 1995
[17] Stanimirov, Schnelle Kurzschlussidentifikation auf Freileitungen, Dissertation Technische Hochschule Zittau 1996
[18] Stoica, Söderström, BIAS correction in least Square identification, Int. J Control 35, 449–457
[19] Hensel, H. Methoden des rechnergestützten Entwurfs und Echtzeiteinsatzes zeitdiskreter Mehrgrössenregelungen und ihre Realsierung in einem CAD-System, Dissertation TH Darmstadt, 1987, VDI-Fortschrittberichte
[20] Schumann, Digitale Mehrgrössenregelung Dissertation TH Darmstadt, 1982, VDI-Fortschrittberichte
[21] Pliquett, U. Elektrische Netzwerkanalyse biologischer Gewebe, Dissertation IHS Mittweida, 1991
Claims (22)
1. Verfahren für eine Impedanzmessfunktion zur Charakterisierung der elektrischen Eigenschaften von Materialien mit Hilfe des elektrischen Wechselstromwiderstandes, der Impedanz, bei welchem der Strom, der das Material durchströmt und die Spannung, die über dem Material abfällt, mit einer Rate 1/Tagemessen werden, um eine Reihe von N Abtastwerten von Strom, i0...iN-1und Spannung u0...uN-1zu erhalten, deren durch den mathematisch-physikalischer Zusammenhang einer elektrischen Ersatzschaltung von beliebig vielen Elementen, wie Induktivitäten, Widerständen oder Kapazitäten etc., durch eine zeitdiskrete Differenzengleichung, der Form
sowie einer anderen abgeleiteten Form dieser Gleichungen beschrieben werden kann dadurch gekennzeichnet, dass N grösser als die Summe der Modellordnungen für Strom & Spannung, (Mi + Mu), und die Abtastwerte in ein Gleichungssystem, der Form
sowie einer anderen abgeleiteten Form dieser Gleichungen, bei der der Modellausgang auch eine Kombination von Abtastwerten sein kann, in diese Gleichung eingesetzt werden kann, und so ein Gleichungssystem überbestimmen, bei dem die Parameter α & β jeweils Funktionen der Abtastfrequenz, einer Referenzfrequenz ω0 und der Parameter der Elemente der Ersatzschaltung, wie Induktivitäten, Widerständen oder Kapazitäten etc., der allgemeinen Form im z-Bereich
sind, und dass aus dem überbestimmten Gleichungssystem die Impedanz Zmeasder Form
sowie die Ersatzparameter beliebig vieler Induktivitäten, Widerstände oder Kapazitäten etc., (R,L,C etc.) aus den inversen funktionalen Zusammenhängen
unter Verwendung der Fehlerausgleichsrechnung, nach weicher der optimale Parametervektor Θ ermittelt wird, insbesondere der verallgemeinerten, gewichteten der Form
quantitativ bestimmt werden kann, wobei sich im Fall wn=l das sog. ungewichtete Verfahren ergibt.
2. Verfahren für eine Impedanzmessfunktion für die Charakterisierung von Materialeigenschaften mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach Anspruch 1, dadurch gekennzeichnet, dass der für die Impedanzmessfunktion zu bestimmende optimale Parametervektor Θ, alternativ iterativ unter Anwendung der rekursiven gewichteten Fehlerausgleiehsrechnung. der Form
gelöst wird, indem für eine Vielzahl von n ein Näherungswert für Θn bzw. Θ berechnet wird, wobei die Stärke der Wichtung, bzw. die Länge des Gedächtnisses, bzw. Vergessens von der Wahl von λ bestimmt wird, welche vorzugsweise zwischen 0,8...0,99 liegen sollte, während bei der Wahl von λ=1 die Ergebnisse der ungewichteten Methode entsprechen, quantitativ bestimmt werden kann.
3. Verfahren für eine Impedanzmessfunktion für die Charakterisierung von Materialeigenschaften mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach Anspruch 1, dadurch gekennzeichnet, dass der für die Impedanzmessfunktion zu bestimmende optimale Parametervektor Θ, alternativ unter Anwendung der totalen Fehlerausgleichsrechnung, der Form quantitativ bestimmt wird.
wobei zuerst muss das Gleichungssystem nach [1] mit der sog. Singulärwertzerlegung (SVD) rechnerisch in Zwischenergebnisse V umgewandelt werden und anschliessend aus den Zwischenergebnissen V die Parameter Θ ermittelt werden, aus welchen sich dann der gewünschte Impedanzverlauf der Impedanzmessfunktion generieren lässt.
4. Verfahren für eine Impedanzmessfunktion für die Charakterisierung von Materialeigenschaften mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach Anspruch 1, dadurch gekennzeichnet, dass der für die Impedanzmessfunktion zu bestimmende optimale Parametervektor Θ, alternativ unter Anwendung von Kaiman-Filtern, der Form
quantitativ bestimmt werden kann, indem das Gleichungssystem für die Modellabweichungen für eine Vielzahl von n mit einem Näherungswert für Θn gelöst wird.
5. Verfahren für eine Impedanzmessfunktion für die Charakterisierung von Materialeigenschaften mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach einem der vorangegangenen Ansprüche, dadurch gekennzeichnet, dass die Ermittlung der Ersatzparameter des Ersatzkreises, wie Induktivität, Kapazitäten und Widerstände etc. offsetkompensiert aus dem optimale Parametervektor Θ nach den Verfahren nach Anspruch 1 bis 4, alternativ unter Anwendung eines mit einem Signalmodell für das Offset (yn=ξ · 1) kombiniertes Systemmodells, der Form
bestimmt werden kann, wobei sowohl die Spannung als auch Strom als Modellausgang und sowohl Strom als auch Spannung als Modelleingang fungieren können, was die aufwendige Offsetkompensation per Hardware in einigen Anwendungen einspart.
6. Verfahren für eine Impedanzmessfunktion, für die Charakterisierung von Materialeigenschaften, mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung, dadurch gekennzeichnet, dass sich aus einem im Zeitbereich aufgestellten Gleichungssystem für einen beliebigen elektrischen Ersatzkreis, der ans Widerständen, Kapazitäten, Induktivitäten etc. nach Fig. 4 bestehen kann, welches auf dem Zusammenhang von Strom und Spannung bzw. deren Abtastwerte sowie über verschiedenen Fenstern gebildete Mittelwerte der Form
die Parameter der Ersatzschaltung sowie die Impedanzverläufe Z=f(ω, Θ) durch Lösen des Gleichungssystems in der Form
bestimmen lassen, wobei für x und y jeweils Abtastwerte von Strom und Spannung bzw. deren (Mehrfach)-Ableitungen oder (Mehrfach)-Integrale entsprechend deren Kopplung durch den elektrischen Ersatzkreis verwendet werden müssen und die einzelnen Elemente des Parametervektors auch math.-physikalisch verknüpfte Basisparameter, wie Widerstand, Induktivität oder Kapazität, enthalten dürfen.
7. Verfahren für eine Impedanzmessfunktion für die Charakterisierung von Materialeigenschaften mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach Anspruch 5, dadurch gekennzeichnet, dass sich mit der Vorschrift der Form
die Parameter der Parallelschaltung einer Kapazität und eines Widerstandes nach Fig. 4b ) sowie deren Impedanzverläufe Z=f(ω, R. C) nach (44) bestimmen lassen, wobei die Summations- und Integrationfenster so zu wählen sind, dass mindestens 2 Abtastwerte darin liegen sollten.
8. Verfahren für eine Impedanzmessfunktion für die Charakterisierung von Materialeigenschaften mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach Anspruch 5, dadurch gekennzeichnet, dass sich mit der Vorschrift der Form
oder
die Parameter der Parallelschaltung einer Kapazität und eines Widerstandes nach Fig. 4c sowie deren Impedanzverläufe Z=f(ω, R, C) nach (48) bestimmen lassen, wobei die Summations- und Integrationfenster so zu wählen sind, dass mindestens 2 Abtastwerte darin liegen sollten.
9. Verfahren für eine Impedanzmessfunktion für die Charakterisierung von Materialeigenschaften mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung, dadurch gekennzeichnet, dass die Impedanz unter Anwendung der Fehlerausgleichsrechnung für Signalmodelle für Strom und Spannung, der Form
unter Anwendung von
oder einem der vorangegangenen Verfahren gemäss den vorangegangenen Ansprüchen zur Ermittlung des optimalen Parametervektor Θ, bestimmt werden kann.
10. Verfahren für eine Impedanzmessfunktion für die Charakterisierung von Materialeigenschaften mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach Anspruch 8, dadurch gekennzeichnet, dass die Impedanz unter Verwendung des Kaiman-Filters für Signalmodelle der Ordnung (M + 1) für Strom und Spannung, vorzugsweise der Formen
unter Verwendung der Vorschrift, die jeweils für u und i anzuwenden ist,
bestimmt werden kann, wobei k im Intervall [1...2 · M] mit M die Anzahl der Frequenzkomponenten grösser als 1 zu wählen ist und die Zeitauflösung der gewünschten Werte durch Anwendung des adaptiven Verfahrens erheblich erhöht werden kann, insbesondere wenn Impedanzwerte bei tiefen Frequenzen ermittelt werden.
11. Verfahren für eine Impedanzmessfunktion für die Charakterisierung von Materialeigenschaften mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach Anspruch 1, dadurch gekennzeichnet, dass der optimale Parametervektor Θ für die Impedanzbestimmung, alternativ unter Anwendung der Verallgemeinerten Methode der kleinsten Quadrate (GLS) nach (12), unter Verwendung des Modells Erweiterten Methode der kleinsten Quadrate (ELS) nach (13) sowie dem Verfahren mit Biaskorrektur (CLS), welche spezielle Umgebungsbedingungen bzw. Randbedingungen zusätzlich zu den Modellen, zur Ermittlung der Impedanz, modellieren, oder der Totalen Methode der kleinsten Quadrate (TLS) nach (18, 19) quantitativ bestimmt werden kann.
12. Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien mittels Messung des Wechselstromwiderstandes der zu charakterisierenden Materialien zur Ausführung der Verfahren nach einem oder mehreren der vorangegangenen Ansprüche umfassend
– Stimulationseinheit, welche aus Strom oder Spannungsquellen bestehen kann, sowie Auswerte-und Kommunikations- und Ausgabeeinheit (Fig. 1 ( 1 )),
– Schutz- und Filtereinheit (Fig. 1 ( 7 )), welche unter anderem die Einhaltung der elektromagnetischen Verträglichkeit sowie des Abtasttheoremes garantieren
– Strom- und Spannungsmessgeräte (Fig. 1 ( 2 , 3 )), welche Strom und Spannung zeitaufgelöst erfassen
– Elektroden (Fig. 1 ( 5 ), 7 – 10 ), welche als Klebe- oder Nadelelektroden sowie mit Mikronadeln besetzte Flächen oder Flächen, Stab-, Zylinder-, Kugel, Halbkugel-oder Plattenelektroden sowie flexiblen Textilelektroden beliebiger Form ausgeführt sein können und als Antennen für die Kommunikation fungieren können, wobei
– die Elektroden (Fig. 1 ( 5 ), 7 – 10 ) mit der Stimulations-, Auswerte- und Kommunikationseinheit (Fig. 1 ( 1 ), 7 – 10 ) derart miteinander verbunden sind, dass von der Stimulationseinheit generierte der Strom zum einen das Testobjekt durchströmen kann und zum anderen vom Messgerät (Fig. 1 ( 2 )), während die über dem Testobjekt abfallende Spannung durch das Messgerät (Fig. 1 ( 3 ) erfasst werden kann,
dadurch gekennzeichnet, dass die Stimulationseinheit, welche im einfachsten Fall nur mit einer FET-Stromquelle zur Strominjektion gemäss Fig. 5 oder einem einfachen Schaltelement, über welches eine Spannung angelegt wird, ausgeführt sein kann, das Testobjekt mit geeigneten, frei wählbaren, der Applikation angepassten, Strömen via Elektroden durchfliessen lässt bzw. Spannungen, die ebenfalls frei wählbar der Applikation anpassbar sind, via Elektroden an das Testobjekt anlegt, welche dann mit den Messgeräten zeitaufgelöst erfasst und der Auswerteeinheit übermittelt, die ein Auswertemittel für die Bestimmung der Impedanz Z=f(un, in), der Impedanzparameter (α, β=f(un, in)) oder der Parameter eines elektrischen Ersatzkreises (R,L,C=f(un, in)) enthält, genutzt werden können und deren Ergebnisse ausgegeben werden können.
13. Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien mittels Messung des Wechselstromwiderstandes der zu charakterisierenden Materialien nach Anspruch 11, umfassend
– Stimulationsschaltung, bestehend aus einer H-Schaltung, bei der die Schaltelemente (Fig. 2 ( 8 )) elektronisch oder elektromechanisch ausgeführt sein können
– Schutz- und Filtereinheit, die als einfache passive RC-Filter (Fig. 1 ( 7 )) ausgeführt sein können, wobei die Messgeräte (Fig. 2 ( 2 , 3 ) ) so angeordnet sind, dass sie den Strom durch das Testobjekt und die Spannung über dem Testobjekt erfassen,
dadurch gekennzeichnet, dass über die H-Schaltung einfach sequenziell entgegengesetzt fliessende Ströme in das Testobjekt eingespeist werden, welche vorzugsweise bereits auch mit binären Ausgängen konventioneller Mikrokontroller (BO) gemäss Fig. 7 ausgeführt werden kann, wobei die Signalform, beispielsweise eine Rechteckform, einfach durch die vorzugsweise passive Filtereinheit, beispielsweise durch eine RC-Schaltung in eine andere Signalform, beispielsweise eine Sägezahnform, umgewandelt werden kann, wobei sich die Stimulationssignale Strom und Spannung einfach mit den Spannungsfällen über Shunts gemäss Fig. 7 erfassen lassen, wobei grundsätzlich nur ein Shunt erforderlich wäre.
14. Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien mittels Messung des Wechselstromwiderstandes der zu charakterisierenden Materialien nach den Ansprüchen 11 und 12 dadurch gekennzeichnet, dass die Vorrichtung auch zusammen mit mehr als 4 Elektroden für elektrische Impedanz-Tomographie (EIT) verwendet werden kann, wobei die Anzahl der Messstellen für Strom und Spannung die Raumauflösung erhöht.
15. Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien mittels Messung des Wechselstromwiderstandes der zu charakterisierenden Materialien nach Anspruch 11 dadurch gekennzeichnet, dass das die Stimulationssignale auch kapazitiv über Kapazitäten bzw. Kondensatoren unterschiedlicher Geometrie oder induktiv, über Spulen unterschiedlicher Geometrie oder Transformatoren, zur Stimulation des Testobjektes eingekoppelt werden können.
16. Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien mittels Messung des Wechselstromwiderstandes der zu charakterisierenden Materialien nach einem oder mehrerer der vorangegangenen Ansprüche dadurch gekennzeichnet, dass das die gemessenen Signale von Strom und Spannung auch kapazitiv, über Kapazitäten bzw. Kondensatoren unterschiedlicher Geometrie, oder induktiv, über Spulen unterschiedlicher Geometrie, Transformatoren, zur Impedanz- bzw. Parameterbestimmung ausgekoppelt werden können, wodurch die gewünschte Sensorfunktion auch ohne direkten galvanischen Kontakt zum Testobjekt, also von «aussen», wenn sich das Testobjekt innerhalb eines elektrisch und magnetisch nicht schirmenden Gefässes oder Rohres befindet, ausführbar ist.
17. Anwendung des Verfahrens für eine Impedanzmessfunktion zur Charakterisierung der elektrischen Eigenschaften von Materialien mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach einem oder mehreren der vorangegangenen Ansprüche
– zur Detektion von Substanzen, wie Eiweisse (Proteine), Schadstoffen, Öl (Lipide), Glukose, Wasser, Kohlenwasserstoffe, sowie deren Gemische, in Materialien, Geweben, Flüssigkeiten, Dispersionen und Emulsionen zur Charakterisierung von biologischen und nicht-biologischen Materialien, Flüssigkeiten, Dispersionen, Emulsionen und Geweben,
– zur Aufnahme des Impedanzspektrums unterschiedlicher biologischer sowie nicht-biologischer Materialien, Geweben, Flüssigkeiten, Schichten (auch Ionenschichten),
– als Mittel zur schnellen Ermittlung des Impedanzspektrums von Zysten, Lipomen, Karzinomen, Sarkomen, Ödemen anhand dem Körper entnommener Proben
– als Mittel zur Durchführung von zytologischen, toxikologischen, ökotoxikologischen und allergologischen Tests,
– als Mittel zur Impedanzmyographie, wobei Stimulation und Messung auch über Textilelektroden, welche an einem oder mehreren Gürteln angebracht sein können,
– zur Beobachtung bzw. Monitoring, Steuerung und Regelung Verfahrens- und biotechnologischer Prozesse, wie der Biotreib- und Rohstoffproduktion, unter Nutzung von Algen, Cyanobakterien sowie anderen Mikroorganismen oder Enzymen sowie zur Bestimmung von Reaktionskonstanten, wie z.B. der Michaelis-Mentenkonstanten in enzymatischen Reaktionen,
– zur Qualitätskontrolle von Lebensmitteln, wie Getränke z.B. Säfte, Wein, Bier, Milch und deren Produkte, Fleisch, Fisch, Honig, Eier, Mehl, Gemüse, wie Bananen, Tomaten, Kirschen, Pilzen Futtermittel, wie Heu oder Silage, Pharmaka sowie Zwischenprodukten, in Produktion, Auslieferung oder Lagerung
– zur Qualitätskontrolle von Lebensmitteln hinsichtlich der Verifikation, Identifikation der Herkunft von Fleisch, Fisch (Delfin, Wal, Hai, Schwein etc.)
– zur Qualitätskontrolle von Spermien, Stammzellen in Produktion, Auslieferung oder Lagerung
– zur Qualitätskontrolle von kühl bzw. cryotechnisch gelagerten Produkten, wie Stammzellen, Spermien, zu transplantierenden Organen oder Blut
– zum Monitoring von Pflanzenwachstum, sowie dem Monitoring Pflanzen-physiologischer Parameter
– als Mittel zur Impedanzmessung für sog. Herz- und Hirnschrittmacher, für die Steuerung, das Monitoring
– als Mittel zur Fruchtwasseruntersuchung,
– als Mittel zur optimalen Platzierung und der Überwachung/Feststellung der aseptischen Lockerung von Implantaten
– zur Korrosionsanalyse, Korrosionsdetektion und Korrosionsmonitoring,
– zur Bestimmung der Zusammensetzung, des Feuchtegehaltes und der Viskosität von Materialien, Dispersionen und Emulsionen
– zum In-Situ-Monitoring der Qualität bzw. Zusammensetzung von Trink- und Abwasser sowie verfahrenstechnischer Produkte in Prozessen,
– zur Beobachtung bzw. Monitoring von zellulären Wachstums- und Differenzierungs- und Transformationsprozessen, Zelle-Zelle-Interaktionen, Zellhaftung, Zellform, Zellgrösse, Zellmenge, Aktivität der Ionenkanäle, Zellmetabolismus, insbesondere in Invivo und Invitro-Systemen,
–als Mittel zur Bestimmung des Zustandes von Filtern oder Membranen in Filtern, (Mikro)-(Bio)-Reaktoren oder künstlichen Organen, wie einer künstlichen Pankreas oder einem Dialysesystem, wobei die Spacer/Abstandshalter/Stabilisatoren welche etwaige die Membranen oder Filter evtl. umgeben multifunktional als Elektroden benutzt werden könnten
– zur Bestimmung der Dicke und Ausdehnung und Einwachsen von biologischen und nichtbiologischen Schichten,
– zur Beobachtung, bzw. dem Monitoring, Steuerung und Regelung von Prozessen in (Mikro)-(Bio-)-Reaktoren, oder disposable (Membran)-Reaktoren
– zur Charakterisierung und Identifikation von Bakterien, Viren und anderen Mikroorganismen wie z.B. Identifikation der unterschiedlicher Stadien des Malaria-Zyklus in Blut
– zur Beobachtung der Aktivität von Ionenkanälen in Zellen oder des Zell-Metabolisums
– zur Charaktersierung und Identifikation von Nanopartikeln
– zur Beobachtung der vulkanischen und seismischen Aktivität in Regionen
– als Mittel zum Monitoring und Feststellung der Wiederherstellung der Barrierefunktion der Haut und zur objektiven Erkennung der Propriorezetion bzw. der Funktion der Propriorezeptoren,
– als Mittel zum Monitoring des Wasserhaushaltes von Menschen, oder Glukose in Schweiss, Blut, Urin, zur Monitoring, Feststellung der Blutgerinnung
– zur Entwicklung intelligenter medizinischer Geräte, wie beispielsweise das intelligente Skalpell, welches vor dem Schnitt das Vorhandensein von Krebszellen verifiziert
dadurch gekennzeichnet dass, die jeweils gewünschte Sensorfunktion in Echtzeit durch Anwendung der Impedanz-Messung «Marker-frei», d.h. dass die Substanzen detektiert werden können, ohne Zugabe spezieller Marker-Substanzen, welche sich nach unterschiedlichen Mechanismen an die zu detektierenden Substanzen heften, sowie sehr schnell ein Breitbandiges Spektrum erzeugt werden kann, wobei die Funktion in einem Bruchteil der Periodendauer der tiefsten Oszillation des Spektrums betragen kann, ein ausführbar ist.
18. Verwendung der Vorrichtung zur Ausführung des Verfahrens für eine Impedanzmessfunktion zur Charakterisierung der elektrischen Eigenschaften von Materialien mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach einem oder mehreren der vorangegangenen Ansprüche nach einem oder mehreren der vorangegangenen Ansprüche als Mittel
– zur Detektion von Substanzen, wie Eiweisse (Proteine), Schadstoffen, Öl (Lipide), Glukose, Wasser, Kohlenwasserstoffe, sowie deren Gemische, in Materialien, Geweben, Flüssigkeiten, Dispersionen und Emulsionen zur Charakterisierung von biologischen und nicht-biologischen Materialien, Flüssigkeiten, Dispersionen, Emulsionen und Geweben,
– zur Aufnahme des Impedanzspektrums unterschiedlicher biologischer sowie nicht-biologischer Materialien, Geweben, Flüssigkeiten, Schichten (auch Ionenschichten),
– als Mittel zur schnellen Ermittlung des Impedanzspektrums von Zysten, Lipomen, Karzinomen, Sarkomen, Ödemen anhand dem Körper entnommener Proben
– als Mittel zur Durchführung von toxikologischen, ecotoxikologischen und allergologischen Tests,
– als Mittel zur Impedanzmyographie, wobei Stimulation und Messung auch über Textilelektroden, welche an einem oder mehreren Gürteln angebracht sein können,
– zur Beobachtung bzw. Monitoring, Steuerung und Regelung Verfahrens- und biotechnologischer Prozesse, wie der Biotreib- und Rohstoffproduktion, unter Nutzung von Algen, Cyanobakterien sowie anderen Mikroorganismen oder Enzymen sowie zur Bestimmung von Reaktionskonstanten, wie z.B. der Michaelis-Mentenkonstanten in enzymatischen Reaktionen,
– zur Qualitätskontrolle von Lebensmitteln, wie Getränke z.B. Säfte, Wein, Bier, Milch und deren Produkte, Fleisch, Fisch, Honig, Eier, Mehl, Gemüse, wie Bananen, Tomaten, Kirschen, Pilzen, Futtermittel wie Heu oder Silage, Pharmaka sowie Zwischenprodukten, in Produktion, Auslieferung oder Lagerung
– zur Qualitätskontrolle von Lebensmitteln hinsichtlich der Verifikation, Identifikation der Herkunft von Fleisch, Fisch (Delfin, Wal, Schwein)
– zur Qualitätskontrolle von Spermien, Stammzellen in Produktion, Auslieferung oder Lagerung
– zur Qualitätskontrolle von kühl bzw. cryotechnisch gelagerten Produkten, wie Stammzellen, Spermien, zu transplantierenden Organen oder Blut
– zum Monitoring von Pflanzenwachstum, sowie dem Monitoring Pflanzen-physiologischer Parameter
– als Mittel zur Impedanzmessung für sog. Herz- und Hirnschrittmacher, für die Steuerung, das Monitoring
– als Mittel zur Fruchtwasseruntersuchung,
– als Mittel zur optimalen Platzierung und der Überwachung/Feststellung der aseptischen Lockerung von Implantaten
– zur Korrosionsanalyse, Korrosionsdetektion und Korrosionsmonitoring,
– zur Bestimmung der Zusammensetzung, des Feuchtegehaltes und der Viskosität von Materialien, Dispersionen und Emulsionen
– zum In-Situ-Monitoring der Qualität bzw. Zusammensetzung von Trink- und Abwasser sowie verfahrenstechnischer Produkte in Prozessen,
– zur Beobachtung bzw. Monitoring von zellulären Wachstums- und Differenzierungs- und Transformationsprozessen in Invitro-Systemen,
– als Mittel zur Bestimmung des Zustandes von Filtern oder Membranen in Filtern, (Mikro)-(Bio)-Reaktoren oder künstlichen Organen, wie einer künstlichen Pankreas oder einem Dialysesystem, wobei die Spacer/Abstandshalter/Stabilisatoren welche etwaige die Membranen oder Filter evtl. umgeben multifunktional als Elektroden benutzt werden könnten
– zur Bestimmung der Dicke und Ausdehnung und Einwachsen von biologischen und nichtbiologischen Schichten,
– zur Beobachtung, bzw. dem Monitoring, Steuerung und Regelung von Prozessen in (Mikro)-(Bio) -oder disposable (Membran)-Reaktoren
– zur Charakterisierung von Bakterien
– zur Charaktersierung und Identifikation von Nanopartikeln
– zur Beobachtung der vulkanischen und seismischen Aktivität in Regionen
– als Mittel zum Monitoring und Feststellung der Wiederherstellung der Barrierefunktion der Haut und zur objektiven Erkennung der Propriorezetion bzw. der Funktion der Propriorezeptoren,
– als Mittel zum Monitoring des Wasserhaushaltes von Menschen, oder Glukose in Schweiss, Blut, Urin, zur Monitoring, Feststellung der Blutgerinnung,
dadurch gekennzeichnet dass, die jeweils gewünschte Sensorfunktion durch Anwendung der Impedanz-Messung «Marker-frei», d.h. dass die Substanzen detektiert werden können, ohne Zugabe spezieller Marker-Substanzen, welche sich nach unterschiedlichen Mechanismen an die zu detektierenden Substanzen heften, sowie sehr schnell ein breitbandiges Spektrum erzeugt werden kann, wobei die Funktion in einem Bruchteil der Periodendauer der tiefsten Oszillation des zu analysierten Spektrums betragen kann, ein ausführbar ist.
19. Anwendung des Verfahrens für eine Impedanzmessfunktion zur Charakterisierung der elektrischen Eigenschaften von Materialien mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach einem oder mehreren der vorangegangenen Ansprüche nach einem oder mehreren der vorangegangenen Ansprüche in zeitlich synchroner Kombination mit optischen und akustischen Analysemethoden, wie Mikroskopie, IR-Spektroskopie, VIS-Spektroskopie, Fluoreszenz-Spektroskopie und/oder (Ultra)-Schall-Analyse bzw. -Spektroskopie zur Bestimmung von Materialeigenschaften, dadurch gekennzeichnet dass, die gewünschte Sensorfunktion durch die verknüpfte Anwendung und Auswertung mehrerer unterschiedlicher Sensorprinzipien aufgrund der Korrelationsfreiheit der jeweiligen Störungen, aufgrund der unterschiedlichen physikalischen Entstehungsmechanismen, sensitiver und auch schneller ausführbar ist.
20. Verwendung der Vorrichtung zur Ausführung des Verfahrens für eine Impedanzmessfunktion zur Charakterisierung der elektrischen Eigenschaften von Materialien mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach einem oder mehreren der vorangegangenen Ansprüche nach einem oder mehreren der vorangegangenen Ansprüche in zeitlich synchroner Kombination mit optischen und akustischen Analysemethoden, wie Mikroskopie, IR-Spektroskopie, VIS-Spektroskopie, Fluoreszenz-Spektroskopie und/oder (Ultra)-Schall-Analyse bzw. -Spektroskopie zur Bestimmung von Materialeigenschaften, dadurch gekennzeichnet dass, die gewünschte Sensorfunktion durch die verknüpfte Anwendung und digitale Auswertung mehrerer unterschiedlicher Sensorprinzipien aufgrund der Korrelationsfreiheit der jeweiligen Störungen, aufgrund der unterschiedlichen physikalischen Entstehungsmechanismen, sensitiver und auch schneller ausführbar ist.
21. Anwendung des Verfahrens für eine Impedanzmessfunktion zur Charakterisierung der elektrischen Eigenschaften von Materialien mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach einem oder mehreren der vorangegangenen Ansprüche innerhalb einer verschluckbaren endoskopischen Kapsel, welche zus. Bild- und Temperaturinformation auch Impedanzmesswerte nach aussen drahtlos übermitteln kann, oder innerhalb einer Pipette, bei der die erforderlichen Elektroden in die Pipettenspitze integriert sind,
– zur Erhöhung der Qualität der Kapselendoskopie durch Impedanzmesswerte
– zur Effizienzerhöhung der Pipettierung
dadurch gekennzeichnet dass, die Impedanzmessfunktion kostengünstig miniaturisiert implementiert werden kann, sowie sehr schnell ein breitbandiges Impedanzspektrum, aus dem die Leitfähigkeit abgeleitet werden kann, erzeugt werden kann, wobei die Funktion in einem Bruchteil der Periodendauer der tiefsten Oszillation des zu analysierten Spektrums betragen kann, ein ausführbar ist.
22. Verwendung der Vorrichtung zur Ausführung des Verfahrens für eine Impedanzmessfunktion zur Charakterisierung der elektrischen Eigenschaften von Materialien mit Ermittlung der Impedanz aus den Abtastwerten von Strom und Spannung nach einem oder mehreren der vorangegangenen Ansprüche innerhalb einer verschluckbaren endoskopischen Kapsel, welche zus. Bild- und Temperaturinformation auch Impedanzmesswerte nach aussen drahtlos übermitteln kann, oder innerhalb einer Pipette, bei der die erforderlichen Elektroden in die Pipettenspitze integriert sind, oder innerhalb einer Pipette, bei der die erforderlichen Elektroden in die Pipettenspitze integriert sind,
– zur Erhöhung der Qualität der Kapselendoskopie durch Impedanzmesswerte
– zur Effizienzerhöhung der Pipettierung
dadurch gekennzeichnet dass, die Impedanzmessfunktion kostengünstig miniaturisiert implementiert werden kann, sowie sehr schnell ein breitbandiges Impedanzspektrum, aus dem die Leitfähigkeit abgeleitet werden kann erzeugt werden kann, wobei die Funktion in einem Bruchteil der Periodendauer der tiefsten Oszillation des zu analysierten Spektrums betragen kann, ein ausführbar ist.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CH14342012A CH706854A2 (de) | 2012-08-21 | 2012-08-21 | Verfahren und Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien. |
| PCT/IB2013/001826 WO2014030054A2 (de) | 2012-08-21 | 2013-08-21 | Verfahren und vorrichtung zur bestimmung der elektrischen eigenschaften von materialien |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CH14342012A CH706854A2 (de) | 2012-08-21 | 2012-08-21 | Verfahren und Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien. |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CH706854A2 true CH706854A2 (de) | 2014-02-28 |
Family
ID=49596333
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CH14342012A CH706854A2 (de) | 2012-08-21 | 2012-08-21 | Verfahren und Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien. |
Country Status (2)
| Country | Link |
|---|---|
| CH (1) | CH706854A2 (de) |
| WO (1) | WO2014030054A2 (de) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019243596A1 (en) * | 2018-06-22 | 2019-12-26 | Université De Rennes 1 | System and method for the estimation of physical parameters of a medium |
| RU2793973C1 (ru) * | 2022-05-25 | 2023-04-11 | Федеральное государственное бюджетное учреждение науки Сибирский федеральный научный центр агробиотехнологий Российской академии наук (СФНЦА РАН) | Способ диагностики грибных болезней земляники садовой |
| CN120908123A (zh) * | 2025-07-18 | 2025-11-07 | 武汉惠强新能源材料科技有限公司 | 一种锂离子电池隔膜喷涂涂层检测系统及检测方法 |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102017011176B4 (de) * | 2017-12-05 | 2022-07-07 | Andreas Güttler | Kapazitives Tomographiegerät und kapazitives Tomographieverfahren |
| CN110702993B (zh) * | 2019-11-19 | 2025-06-24 | 宁波爱得智能科技有限公司 | 一种纺织结构电极的低频阻抗测试装置及方法 |
| CN115868594A (zh) * | 2021-09-26 | 2023-03-31 | 青岛海尔电冰箱有限公司 | 加热装置及用于加热装置的控制方法 |
| CN113866224B (zh) * | 2021-12-01 | 2022-02-15 | 广东利诚检测技术有限公司 | 基于磁电耦合方法检测食品中有害金属元素含量的方法 |
| CN116149186B (zh) * | 2023-02-16 | 2023-08-11 | 大连交通大学 | 一种卫星姿态控制系统故障估计的Kalman滤波方法 |
| CN117538401B (zh) * | 2023-11-20 | 2025-04-18 | 南京航空航天大学 | 一种基于oect与eit结合的成像方法 |
| CN119482692A (zh) * | 2025-01-14 | 2025-02-18 | 济南德明电源设备有限公司 | 一种储能变流器的动态电压控制方法及系统 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2540495B1 (fr) | 1983-02-07 | 1986-02-14 | Roussel Uclaf | Nouveaux derives de o-mercaptopropanamide et de ses homologues, leur procede de preparation, leur application comme medicaments, les compositions les renfermant et les nouveaux intermediaires obtenus |
| US6058325A (en) * | 1996-04-16 | 2000-05-02 | Cardiotronics | Method and apparatus for high current electrode, transthoracic and transmyocardial impedance estimation |
| DE50012309D1 (de) | 2000-11-06 | 2006-04-27 | Abb Research Ltd | Verfahren zum Messen des Widerstands und der Induktivität einer Leitung |
| DE102005034859A1 (de) * | 2005-07-26 | 2007-02-01 | Knorr-Bremse Systeme für Nutzfahrzeuge GmbH | Meßanordnung zur Messung des Induktivitäts- und des Widerstandswertes eines induktiven Sensors |
| WO2011002793A1 (en) * | 2009-06-30 | 2011-01-06 | The University Of Connecticut | Multiple excitation capacitance polling for enhanced electronic capacitance tomography |
-
2012
- 2012-08-21 CH CH14342012A patent/CH706854A2/de not_active Application Discontinuation
-
2013
- 2013-08-21 WO PCT/IB2013/001826 patent/WO2014030054A2/de not_active Ceased
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019243596A1 (en) * | 2018-06-22 | 2019-12-26 | Université De Rennes 1 | System and method for the estimation of physical parameters of a medium |
| RU2793973C1 (ru) * | 2022-05-25 | 2023-04-11 | Федеральное государственное бюджетное учреждение науки Сибирский федеральный научный центр агробиотехнологий Российской академии наук (СФНЦА РАН) | Способ диагностики грибных болезней земляники садовой |
| CN120908123A (zh) * | 2025-07-18 | 2025-11-07 | 武汉惠强新能源材料科技有限公司 | 一种锂离子电池隔膜喷涂涂层检测系统及检测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2014030054A2 (de) | 2014-02-27 |
| WO2014030054A3 (de) | 2014-04-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CH706854A2 (de) | Verfahren und Vorrichtung zur Bestimmung der elektrischen Eigenschaften von Materialien. | |
| Showkat et al. | A review of bio-impedance devices | |
| Aguiar Santos et al. | System description and first application of an FPGA-based simultaneous multi-frequency electrical impedance tomography | |
| Wood | Physiology, biophysics, and biomedical engineering | |
| Critcher et al. | Localized bioimpedance measurements with the max3000x integrated circuit: Characterization and demonstration | |
| Wang et al. | Biomedical sensors and measurement | |
| Varghese et al. | Design, fabrication and evaluation of a stretchable high-density electromyography array | |
| Sirtoli et al. | Design and evaluation of an electrical bioimpedance device based on DIBS for myography during isotonic exercises | |
| Marcôndes et al. | Parasitic effects on electrical bioimpedance systems: Critical review | |
| Arabsalmani et al. | Exploring bio-impedance sensing for intelligent wearable devices | |
| Rendon-Morales et al. | Non-invasive electrocardiogram detection of in vivo zebrafish embryos using electric potential sensors | |
| Heuer et al. | Motion artefact correction for capacitive ECG measurement | |
| Yoshida et al. | The theory of peripheral nerve recording | |
| Raus-Jarzabek et al. | A practical guide to ECG device performance testing according to international standards | |
| Pîslaru-Dănescu et al. | Design and manufacturing of equipment for investigation of low frequency bioimpedance | |
| Bahador et al. | Robust removal of slow artifactual dynamics induced by deep brain stimulation in local field potential recordings using svd-based adaptive filtering | |
| WO2011022851A1 (de) | Verfahren zur eichung einer diagnostischen messvorrichtung | |
| EP3183563B1 (de) | Monitoring der wirkung von substanzen auf in vitro gewebe | |
| Khoma et al. | Advanced computing methods for impedance plethysmography data processing | |
| DE102008032980B4 (de) | Sonde und Vorrichtung für die Messung dielektrischer Materialeigenschaften | |
| DE102015007875A1 (de) | Messanordnung und Verfahren zur in-vivo Bestimmung der Laktatkonzentration in Blut mittels elektrochemischer Impedanzspektroskopie | |
| Pérez-García et al. | Towards bio-impedance based labs: A review | |
| EP3881753B1 (de) | Vorrichtung zur erzeugung eines elektromagnetischen feldes und messung der absorption davon durch ein leitfähiges medium | |
| Nasrollaholhosseini et al. | Impedance spectroscopy of tripolar concentric ring electrodes with Ten20 and TD246 pastes | |
| Mohsen et al. | Modeling woody plant tissue using different fractional-order circuits |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| AZW | Rejection (application) |