מאמר שיטה

שיטות נסיוניות ללימוד בקרת הפוסטאורל האנושית

DOI:

10.3791/60078

⸱

11 בספטמבר 2019

 ,  , 
1Department of Biomedical Engineering, McGill University

במאמר זה

סיכום

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

מאמר זה מציג מסגרת ניסיונית/אנליטית לחקר שליטה בגוף האדם. הפרוטוקול מספק הליכים צעד-אחר-צעד לביצוע ניסויים עומדים, מדידת הקינמטיקה ואותות קינטיקה, וניתוח התוצאות כדי לספק תובנות לגבי המנגנון שבבסיס השליטה האנושית.

תקציר

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

רכיבים רבים של מערכות העצבים ושריר השלד פועלים בהופעה כדי להשיג את היציבה האנושית היציבה והזקופה. ניסויים מבוקרים המלווים בשיטות מתמטיות מתאימות נדרשים כדי להבין את התפקיד של מערכות המשנה השונות המעורבות בשליטה אנושית. מאמר זה מתאר פרוטוקול לביצוע ניסויים עומדים ומודאג, רכישת נתונים ניסיוניים, וביצוע הניתוח המתמטי הבא, במטרה להבין את התפקיד של מערכת שריר-שלד ושליטה מרכזית באדם יציבה זקופה. התוצאות שנוצרו על ידי שיטות אלה חשובות, כי הם מספקים תובנה לתוך בקרת איזון בריא, ליצור את הבסיס להבנת האטיולוגיה של איזון לקוי אצל חולים וקשישים, וסיוע בתכנון של התערבויות לשפר בקרה ויציבות. ניתן להשתמש בשיטות אלה כדי ללמוד את התפקיד של מערכת המגע, הקשיחות הפנימית של מפרק הקרסול, ומערכת חזותית בבקרת הפוסט-שליטה, וניתן גם להרחיב אותה לחקור את תפקידה של מערכת שיווי המידה. השיטות הן לשימוש במקרה של אסטרטגיית קרסול, שבה הגוף נע בעיקר על מפרק הקרסול ונחשב למטוטלת הפוכה של קישור יחיד.

מבוא

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

השליטה האנושית ממומשת באמצעות אינטראקציות מורכבות בין מערכת העצבים המרכזית לבין מערכות שריר-שלד1. הגוף האנושי בעמידה אינו יציב באופן מיסודו, בכפוף למגוון פנימי (כגון: נשימה, פעימות לב) וחיצוני (למשל, הכבידה) רטבאליות. היציבות מושגת על ידי בקר מבוזר עם רכיבים מרכזיים, רפלקס ופנימיים (איור 1).

השליטה הפוסט מושגת על ידי: בקר פעיל, בתיווך על ידי מערכת העצבים המרכזית (CN) ואת חוט השדרה, אשר משנה הפעלת שריר; בקר נוקשות פנימית המתנגד לתנועה משותפת ללא שינוי בהפעלת שרירים (איור 1). הבקר המרכזי משתמש במידע חושי כדי ליצור פקודות בסדר יורד המייצרים כוחות שרירים מתקנות כדי לייצב את הגוף. מידע חושי מושפע ממערכות הראייה, השמיעה והמגע. באופן ספציפי, מערכת המגע מפיקה מידע בנוגע למשטח התמיכה ולזוויות משותפות; ראיה מספקת מידע בנוגע לסביבה; ומערכת השמיעה מפיקה מידע בנוגע למהירות זוויתית מראש, האצת לינארית ואוריינטציה ביחס לכוח הכבידה. הבקר המרכזי, הלולאה הסגורה, פועל עם עיכובים ארוכים שעלולים להיות בעלי יציבות של2. המרכיב השני של הבקר הפעיל הוא קשיות רפלקס, אשר מייצרת פעילות שרירים עם השהיה קצרה ומייצרת התנגדות תנועה משותפת.

קיימת השהיה המשויכת לשני הרכיבים של בקר פעיל; כתוצאה מכך, נוקשות פנימית משותפת, הפועלת ללא דיחוי, ממלאת תפקיד חשוב בבקרת הפוסט-שלוש. נוקשות פנימית נוצרת על ידי ויסקו פסיבי מאפיינים גמישים של שרירי הקבלנות, רקמות רכות ותכונות אינרציאליות של הגפיים, אשר מייצרת torques התנגדות מיידי בתגובה לכל תנועה משותפת4. התפקיד של הקשיחות המשותפת (נוקשות פנימית ורפלקס) בבקרת הפוסט אינו מובן בבירור, שכן הוא משתנה עם תנאי הפעולה, המוגדר על ידי הפעלת שרירים4,5,6 ומיקום משותף ד , מיכל סבן , 8, שניהם משתנים עם הגוף מתנדנד, הטבועה לעמוד.

זיהוי התפקידים של הבקר המרכזי ונוקשות משותפת בשליטה פוסט-אוראל חשוב, כפי שהוא מספק את הבסיס: אבחון האטיולוגיה של ליקויי איזון; תכנון התערבויות ייעודיות למטופלים; הערכת הסיכון לנפילות; פיתוח אסטרטגיות למניעת נפילות בקשישים; ואת העיצוב של התקנים מסייעת כגון orthotics ו תותבות. עם זאת, זה קשה, כי מערכות משנה שונות לפעול יחד ורק את הגוף הכולל כתוצאה מכך, torques משותף, ו אלקטרומגנטים שרירים ניתן למדוד.

לכן, חיוני לפתח שיטות נסיוניות ואנליטיות המשתמשות במשתני הפוסט-מדידים כדי להעריך את תרומתו של כל מערכת משנה. בעיה טכנית היא שמדידה של משתני הפוסטאורל נעשית בלולאה סגורה. כתוצאה מכך, התשומות והתפוקות (סיבה ותוצאה) הם קשורים. כתוצאה מכך, יש להחיל את הרטבאליות החיצוניים (כתשומות) לעורר תגובות פוסט-אוראטיות בתגובות (כפלטי מדפסת), ו-b) להעסיק שיטות מתמטיות מיוחדות כדי לזהות מודלים של מערכת ולגרום לתוצאה ולאפקט9.

המאמר הנוכחי מתמקד בשליטה משותפת כאשר משתמשים באסטרטגיית קרסול, כלומר, כאשר התנועות מתרחשות בעיקר על מפרק הקרסול. במצב זה, הגוף העליון והגפיים התחתונות לנוע יחד, כתוצאה מכך, הגוף יכול להיות מעוצב כמטוטלת הפוכה הקישור יחיד במישור משונן10. אסטרטגיית הקרסול משמשת כאשר משטח התמיכה מוצק והרטבאליות הם קטנים1,11.

מנגנון עמידה המסוגל ליישם מכני מתאים (בקינטיביות) והחושים החזותיים הרטבאליות והקלטת את הגוף הקינמטיקה, הקינטיקה ופעילות השרירים פותחה במעבדה שלנו12. המכשיר מספק את הסביבה הניסיונית הנחוצה כדי ללמוד את התפקיד של קשיות הקרסול, מנגנוני הבקרה המרכזיים, והאינטראקציות שלהם על ידי יצירת תגובות באמצעות שימוש בגירויים חזותיים או/והמגע הבין-חושי. ניתן גם להאריך את המכשיר כדי ללמוד את התפקיד של מערכת שיווי המידה על ידי יישום של גירוי חשמלי ישיר לתהליכי מסטואיד, שיכולים ליצור תחושה של מהירות ראש ולעורר תגובות פוסט12,13 .

אחרים פיתחו גם התקנים דומים כדי ללמוד בקרת הדואר האנושי, שם piezo ליניארי מפעילים חשמליים11, מנועים חשמליים רוטרי14,15, ומנועי חשמל ליניארי16,17 , 18 שימשו להחלת רטבאליות מכני לקרסול בעמידה. התקנים מורכבים יותר פותחו גם כדי ללמוד שליטה מרובת מקטעים בקרה, שם ניתן להחיל רטבאליות מרובים על המפרקים קרסול והירך בו19,20.

מתקנים עומדים

שני סרוו מבוקרת אלקטרו הידראולי מפעילים רוטרי להעביר שני דוושות להחיל רטבאליות נשלט של מיקום הקרסול. הפעילים יכולים לייצר טורקים גדולים (> 500 ננומטר) הדרושים לבקרת הפוסט-שליטה; זה חשוב במיוחד במקרים כגון להישען קדימה, שם מרכז המסה של הגוף הוא רחוק (קדמי) מציר הקרסול של סיבוב, וכתוצאה מכך ערכים גדולים של מומנט הקרסול עבור שליטה הפוסט.

כל מפעיל רוטרי נשלט על-ידי שסתום סרוו פרופורציונלי נפרד, באמצעות משוב פדאל מיקום, שנמדד על ידי מדידת פוטנציאל ביצועים גבוהים על פיר המפעיל (טבלת חומרים). הבקר מיושם באמצעות מערכת עיבוד אותות דיגיטליים בזמן אמת MATLAB מבוסס xPC. למפעיל/סרוו-שסתום יחד יש רוחב פס של יותר מ-40 Hz, הרבה יותר גדול מרוחב הפס של מערכת הבקרה הכוללת של המערכת, נוקשות במפרק הקרסול והבקר המרכזי21.

התקן וסביבה של מציאות וירטואלית

אוזניית מציאות וירטואלית (VR) (טבלת חומרים) משמשת לperturb החזון. האוזנייה האלחוטית מכילה מסך LCD (כפול AMOLED 3.6 ' מסך עם רזולוציה של 1080 x 1200 פיקסלים לעין) המספקת למשתמש תצוגה סטריאוסקופי של המדיה הנשלחת למכשיר, המציעה תפיסת עומק תלת-ממדית. קצב הרענון הוא 90 Hz, מספיק כדי לספק תחושה וירטואלית מוצקה למשתמשים22. שדה התצוגה של המסך הוא 110 °, מספיק כדי ליצור רטבאליות חזותיים דומה למצבים בעולם האמיתי.

האוזנייה עוקבת אחר סיבוב ראש המשתמש ומשנה את התצוגה הווירטואלית בהתאם, כך שהמשתמש שקוע באופן מלא בסביבה הווירטואלית; לכן, היא יכולה לספק את המשוב החזותי הרגיל; והוא יכול גם לperturb חזון על-ידי סיבוב השדה החזותי במישור משונן.

מדידות קינטית

כוח תגובה אנכי נמדד על ידי ארבעה תאים לטעון, דחוקה בין שתי צלחות מתחת לרגל (טבלת חומרים). מומנט הקרסול נמדד ישירות על ידי התמרה מומנט עם קיבולת של 565 ננומטר ונוקשות טורמיעית של 104 kNm/rad; כמו כן ניתן למדוד אותו בעקיפין מהכוחות האנכיים העומדים בפני תאי העומס, תוך שימוש במרחקים שלהם לציר הקרסול של סיבוב23, בהנחה שהכוחות האופקיים החלים על הרגליים בעמידההם קטנים,24. מרכז הלחץ (COP) נמדד במישור משונן על ידי חלוקת מומנט הקרסול על ידי כוח אנכי הכולל, נמדד על ידי העומס התאים23.

מדידות קימטיות

זווית כף הרגל זהה לזווית הדוושה, כי כאשר משתמשים באסטרטגיית קרסול, הרגל של הנבדק זזה עם הדוושה. זווית שאנק ביחס האנכי מושגת בעקיפין מן העקירה הליניארית של shank, נמדד על ידי מוצא לייזר טווח (טבלת חומרים) עם רזולוציה של 50 יקרומטר ו רוחב פס של 750 Hz25. זווית הקרסול היא הסכום של הרגל ואת הזוויות שאנק. זווית הגוף עם כבוד האנכי מתקבל בעקיפין מן העקירה הליניארית של אמצע הנקודה בין שמאל וימין האחורי מעולה השדרה (PSIS), נמדד באמצעות טווח לייזר מאתר (טבלת חומרים) עם רזולוציה של 100 יקרומטר ורוחב פס של 750 הרץ23. מיקום וסיבוב הראש נמדדים ביחס למערכת הקואורדינטות הגלובלית של סביבת VR על-ידי תחנות הבסיס של מערכת VR שפולטות פולסים של אינפרא-אדום מתוזמן (IR) ב-60 פולסים לשנייה, הנאסף על-ידי חיישני אינפרא-אדום עם האוזנייה האלחוטית עם תת-מילימטר דיוק.

רכישת נתונים

כל האותות מסוננים עם מסנן נגד החלקה עם תדירות פינתית של 486.3 ולאחר מכן שנדגמו ב-1000 הרץ עם ביצועים גבוהים 24-bit/8-ערוץ, דגימה בו, כרטיסי רכישת אות דינאמי (לוח חומרים) עם דינאמי טווח של 20 וולט.

מנגנוני בטיחות

שישה מנגנוני בטיחות שולבו במנגנון עמידה כדי למנוע פציעות לנושאים; הדוושות נשלטות בנפרד ולעולם לא יפריעו זה לזה. (1) פיר המפעיל יש מצלמת, אשר מפעילה בצורה מכנית שסתום שינתק לחץ הידראולי אם סיבוב הפיר חורג ± 20 ° ממיקומה האופקי. (2) שתי עצירות מכניות מתכווננות מגבילות את טווח התנועה של המפעיל; אלה מוגדרים לטווח התנועה של כל נושא לפני כל ניסוי. (3) הנושא והנסבית מחזיקים כפתור בהלה; לחיצה על לחצן מתנתק כוח הידראולי מן הפעילים וגורם להם להיות רופף, כך שהם יכולים להיות מועברים באופן ידני. (4) מעקות הממוקמים משני צדי הנושא זמינים לספק תמיכה במקרה של אי יציבות. (5) הנושא לובש רתמה גוף מלא (טבלת חומרים), המצורפת מרובי נוקשה בתקרה כדי לתמוך בהם במקרה של נפילה. הרתמה היא מרווח ואינה מפריעה לעמידה הרגילה, אלא אם כן הנושא הופך לבלתי יציב, כאשר הרתמה מונעת מהנושא ליפול. במקרה של נפילת, תנועות הדוושה תופסק באופן ידני או על ידי הנושא, באמצעות לחצן הפניקה או על ידי הנסניסה. (6) מנגנוני הסרוו עוצרים את סבב הפעילים באמצעות מנגנונים אל-כשל במקרה של הפרעה באספקת החשמל.

הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.

פרוטוקול

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

כל השיטות הנסיוניות אושרו על ידי מועצת האתיקה מחקר של אוניברסיטת מקגיל והנושאים לחתום על הסכמה מושכלת לפני השתתפות.

1. ניסויים

הערה: כל ניסוי כרוך בשלבים הבאים.

  1. טרום בדיקה
    1. הכן מיתאר מובהק של כל המבחנים שיש לבצע ולבצע רשימת פעולות לביצוע עבור איסוף נתונים.
    2. ספק לנושא טופס הסכמה עם כל המידע הדרוש, בקש מהם לקרוא אותו ביסודיות, לענות על שאלות ולאחר מכן לחתום על הטופס.
    3. רשום את המשקל, הגובה והגיל של הנושא.
  2. הכנה לנושא
    1. מדידת אלקטרומגרפיה
      1. השתמש באלקטרודות דיפרנציאליות בודדות (טבלה של חומרים) עם מרחק בין אלקטרודה של 1 ס מ למדידה של אלקטרויוגרפיה (emg) של שרירי הקרסול.
      2. השתמש במגבר (טבלת חומרים) עם רווח כולל של 1000 ורוחב פס של 20-2000 Hz.
      3. כדי להבטיח האות גבוה ליחס הרעש (SNR) ומינימלית לדבר הצלב, לאתר ולסמן את הקבצים המצורפים האלקטרודה לפי ההנחיות המסופקות על ידי פרויקט Seniam26, כמו להלן: (1) עבור הסקופ האמצעי (MG), הבליטה הבולטת ביותר של שרירים (2) עבור הגקינממיוס הצדדיים (LG), 1/3 של הקו בין ראש הפיבולה לעקב; (3) לסוליה (SOL), 2/3 של הקו בין קונדיקלס המדיאלי של עצם הירך והקרסולית המדיאלי; (4) עבור הקרסולית הקדמית (TA), 1/3 של הקו בין קצה הפיבולה והקצה של האמצעי.
      4. לגלח את האזורים מסומנים עם תער ולנקות את העור עם אלכוהול. הניחו לעור להתייבש ביסודיות.
      5. לגלח את האזור הגרמי על צלחית עבור האלקטרודות התייחסות, ונקי עם אלכוהול.
      6. הנושא משקר בתנוחה רגועה.
      7. מניחים את האלקטרודה ההפניה על האזור המגולח של patella.
      8. לצרף את האלקטרודות אחד אחד לאזורים מגולח של השרירים, באמצעות קלטת כפולה, לטפל כדי להבטיח כי אלקטרודות מקובעים את העור בצורה מאובטחת.
      9. לאחר הצבת כל אלקטרודה, בקשו מהנושא לבצע התכווצות באמצעות כיווץ כנגד התנגדות ולבחון את צורות ה-גל על האולוסקופ כדי לוודא שהאות EMG יש SNR גבוהה. אם האות SNR הוא עני, הזז את האלקטרודות עד שיימצא מיקום בעל SNR גבוהה.
      10. ודא שתנועות הנושא אינן מעכבת את כבלי EMG.
    2. מדידות קימטיות
      1. הצמד סמן רפלקטיבי לסכין עם רצועה, לשמש למדידת זווית שאנק.
        הערה: מניחים את סמן הסכין גבוה ככל האפשר על הלהב כדי ליצור את העקירה הקווית הגדולה ביותר האפשרית עבור סיבוב נתון, ולכן, שיפור רזולוציה זוויתית.
      2. . שהנושא ישים על הרתמה של הגוף
      3. לצרף סמן משקפים את המותניים של הנושא עם רצועה, לשמש מדידה זווית הגוף העליון. ודא כי סמן מותניים רפלקטיבי ממוקם באמצע הנקודה בין השמאל והימני, ושהבגדים של הנושא אינו מכסה את המשטח המשקף את המותניים.
      4. . שהנושא יהיה על המנגנון העומד
      5. להתאים את מיקום הרגל של הנושא כדי ליישר את הmalleoli לרוחב והמדיאלי של כל רגל לציר של הדוושה של סיבוב.
      6. התווה את עמדות הרגל של הנושא באמצעות סמן והנחה אותם לשמור על רגליהם באותם מיקומים במהלך הניסויים. זה מבטיח את צירי סיבוב של קרסוליים ומפעילים נשארים מיושרים לאורך כל הניסויים.
      7. כוונן את המיקום האנכי של מפלסי טווח הלייזר כדי להצביע על מרכז סמנים רפלקטיבי. כוונן את המרחק האופקי בין מאתר טווח הלייזר לבין סמנים רפלקטיבי, כך שהטווח מאתר את העבודה בטווח הביניים ואינו מבצע רוויה במהלך העמידה השקטה.
      8. יש את הנושא להישען קדימה ואחורה על הקרסול ולהבטיח כי לייזרים להישאר בתוך טווח העבודה שלהם.
      9. למדוד את הגובה של טווח לייזר מפלסי ביחס לציר הקרסול של סיבוב.
        הערה: גבהים אלה משמשים להמרת displacements ליניארי לזוויות.
    3. פרוטוקולים ניסיוניים
      1. ליידע את הנושא של מה לצפות עבור כל תנאי ניסיון.
      2. הורה לנושא לעמוד בשקט עם ידיים לצדו תוך שהוא מסתכל קדימה, ולשמור על האיזון שלהם כפי שהם עושים, כאשר מתמודדים מול רטבאליות העולם האמיתי.
      3. לקבלת מבחני דאגה, התחילו בהפרעות והניחו לנושא להסתגל אליו.
      4. הפעל את רכישת הנתונים לאחר שהנושא יצר התנהגות יציבה.
      5. לספק את הנושא עם תקופת מנוחה מספקת לאחר כל משפט כדי למנוע עייפות. תקשר איתם כדי לראות. אם הם זקוקים לזמן נוסף
      6. בצע את המבחנים הבאים.
        1. לבדיקת מכשירים, בצע בדיקה של 2 דקות כדי לבדוק את נתוני החיישן 2 h לפני הגעתו של הנבדק. חפש רעשים או היסטים גדולים ביותר בנתוני החיישן המוקלט. אם קיימות בעיות, פתור אותן לפני שהנושא מגיע.
        2. לעמידה שקט, לבצע משפט שקטה 2 דקות לעמוד ללא רטבאליות.
          הערה: משפט זה מספק הפניה, הצורך לקבוע אם/כיצד משתנים שלאחר השינוי בתגובה לרטבאליות.
        3. לניסויים משרכים, הפעל את הפרטורציה ורכוש נתונים עבור 2-3 דקות. החלת הדוושה רטבאליות אם המטרה היא לחקור את התפקיד של מערכת המגע נוקשות/קרסול בעמידה. החלת רטבאליות חזותית אם המטרה היא לבחון את התפקיד של החזון בבקרת הפוסט-אוראל. החלת רטבאליות חזותיים ופדלים בו אם המטרה היא לבחון את האינטראקציה בין שתי המערכות בבקרת הדואר.
          הערה: רטבאליות פדאלים מוחלים כסיבוב של דוושות ההתקן העומד. באופן דומה, רטבאליות חזותיים מוחלים על-ידי סיבוב השדה החזותי הווירטואלי, באמצעות אוזניית VR. הזווית של שדה הדוושה/חזותי מלווה באות, הנבחרת בהתאם ליעדי המחקר. סעיף הדיון מספק פרטים בנוגע לסוגי ההפרעות, המשמשים לחקר השליטה הפוסט-שיתית והיתרונות של כל הפרטורציה.
      7. בצע מינימום של 3 מבחנים עבור כל הפרטורציה מסוימת.
        הערה: מספר ניסויים מתבצע כדי להבטיח אמינות של המודלים בעת ביצוע הניתוח על הנתונים שנאספו; לדוגמה, ניתן לחצות את אימות המודלים.
      8. לבצע את המבחנים בסדר אקראי כדי להבטיח את הנושאים לא ללמוד להגיב לתגובה מסוימת; פעולה זו גם מאפשרת לבדוק התנהגות שונות בזמן.
      9. בדוק את הנתונים ויזואלית לאחר כל משפט כדי לוודא שהאותות הנרכשים הם באיכות גבוהה.

2. זיהוי בקרת השליטה האנושית

  1. זיהוי לא פרמטרי של היחס הדינאמי של זווית הגוף לרטבאליות חזותי
    1. ניסוי
      1. לרכוש מבחנים באופן חזותי עבור 2 דקות בהתאם לשלבים בסעיפים 1.1 ו 1.2.
      2. השתמש באות טרפז (בסדר) עם משרעת שיא לשיא של 0.087 rad ומהירות של 0.105 rad/s.
      3. החזק את מיקום הדוושה קבוע בזווית האפס.
    2. ניתוח
      הערה: ניתוח נתונים במקטעים 2.1.2 ו-2.2.2 מבוצע באמצעות MATLAB.
      1. להשמיד את זווית הגוף הגולמי ואותות ההפרעות החזותיות (כך שהתדירות הגבוהה ביותר היא 10 Hz), תוך שימוש בפקודות הבאות:
        figure-protocol-1
        figure-protocol-2
        איפה
        figure-protocol-3
        figure-protocol-4
        figure-protocol-5
        הערה: עבור קצב דגימה של kHz 1, יחס הפענוח חייב להיות 50 בתדר הגבוה ביותר של 10 Hz.
      2. בחר את תדירות הריבית הנמוכה ביותר, שתקבע את אורך החלון עבור שערוך צריכת חשמל.
        הערה: כאן, תדירות מינימלית של 0.1 Hz נבחרה, כך שאורך החלון עבור שערוך החשמל הוא 1/0.1 Hz = 10 ס מ. רזולוציית התדר זהה לתדר המינימלי, ולכן החישובים נעשים עבור 0.1, 0.2, 0.3,..., 10 Hz.
      3. בחר את סוג החלון ומידת החפיפה כדי למצוא את ספקטרום העוצמה.
        הערה: לאורך ניסיון של 120 s, 10 s הנינג חלונות עם 50% התוצאות חופפים בממוצע של 23 מקטעים עבור אומדן ספקטרום הכוח. מאז שאנחנו הושמדה את הנתונים ל 20 הרץ, חלון של 10 יש אורך של 200 דגימות.
      4. השתמש figure-protocol-6 בפונקציה כדי לאתר את תגובת התדירות (FR) של המערכת:
        figure-protocol-7
        איפה
        figure-protocol-8
        figure-protocol-9
        figure-protocol-10
        figure-protocol-11
        הערה: הפונקציה המוצגת figure-protocol-12 מחשבת את הרוחב בין הפרטורציה VR הושמדה ואת זווית הגוף בתדרים שצוין על-ידי figure-protocol-13 , באמצעות חלון hanning עם האורך שצוין על-ידי figure-protocol-14 ומספר חופפים שווה ל figure-protocol-15 (כלומר, 50% חופפים). באופן דומה, הוא מחשב את הספקטרום האוטומטי של קלט VR. לאחר מכן, באמצעות הספקטרום הצולב והספקטרום האוטומטי, הוא מחשב את FR של המערכת.
      5. מצא את הרווח והשלב של FR המשוער בשלב 2.1.2.4, באמצעות הפקודות הבאות:
        figure-protocol-16
        figure-protocol-17
        איפה
        figure-protocol-18
        figure-protocol-19
      6. חשב את הפונקציה קוהרנטית באמצעות הפקודה הבאה:
        figure-protocol-20
        איפה
        figure-protocol-21
        הערה figure-protocol-22 : הפונקציה מלווה בהליך דומה figure-protocol-23 כדי למצוא figure-protocol-24 את figure-protocol-25 הקוהרנטיות בין ו.
      7. מתווה את הרווח, השלב והקוהרנטיות כפונקציה של תדירות.
        figure-protocol-26
        figure-protocol-27
        figure-protocol-28
        הערה: ניתן להרחיב את השיטה המוצגת למקרה שבו הן רטבאליות חזותיות והן מכניות מוחלות, כאשר יש להשתמש בשיטת זיהוי מרובה-קלט(mimo) באמצעות מספר רב של משתמשים. הזיהוי יכול להיעשות גם באמצעות שיטה תת חלל (אשר ביסודו עוסק מערכות MIMO)27 או באמצעות שיטות העברה פרמטרית פונקציות כמו Mimo Box-ג'נקינס28. הן תת-חלל והן קופסא-ג'נקינס (ושיטות אחרות) מיושמות בארגז הכלים לזיהוי מערכת MATLAB.
  2. זיהוי פרמטרי של קשיות הקרסול הפנימי בעמידה
    1. ניסוי
      1. ביצוע משפטי מכני מכנית עבור 2 דקות. השתמש ברצפי מדומה אקראית (PRBS) עם משרעת שיא לשיא של 0.02 rad ומרווח מיתוג של 200 ms. ודא שזווית ממוצע הדוושה היא אפס.
    2. ניתוח
      1. להבדיל את האות ברגל פעם אחת כדיfigure-protocol-29להשיג מהירות הרגל (figure-protocol-30 , פעמיים כדי להשיג תאוצה ברגלfigure-protocol-31 (ושלוש פעמים כדי לקבל את הנתירה שלה (בדומה להבדיל את מומנט הכוח כדי להשיג את מהירות ההאצה שלה, באמצעות הבאים פקודה
        figure-protocol-32
        איפה
        figure-protocol-33
        figure-protocol-34
        figure-protocol-35
      2. חשב את המיקום של קיצון המקומי מקסימה ומקומית של מהירות כף הרגל כדי לאתר פולסים, באמצעות הפקודה הבאה:
        figure-protocol-36
        figure-protocol-37
        איפה
        figure-protocol-38
        figure-protocol-39
        figure-protocol-40
        figure-protocol-41
        הערה figure-protocol-42 : הפונקציה מוצאת את כל המקסימה המקומית (מהירות הרגל החיובית) ואת מיקומם. כדי למצוא את minima המקומי, הפונקציה זהה משמש, אבל את הסימן של מהירות זווית הרגל חייב להיות הפוך.
      3. עיצוב שלהזמנה 8 באטורת לעבור מסנן נמוך עם תדר פינתי של 50 Hz, באמצעות הפקודה הבאה:
        figure-protocol-43
        figure-protocol-44
        figure-protocol-45
        figure-protocol-46
        figure-protocol-47
      4. סנן את כל האותות באמצעות משמרת אפס-שלב באמצעות המסנן ' בדלורת ':
        figure-protocol-48
        figure-protocol-49
        figure-protocol-50
        הערה:הפונקציה filtfiltאינה גורמת לשינוי באות המסונן. אין להשתמש בפונקציה "filter" , מכיוון שהיא מפיקה משמרת.
      5. העלילה מהירות הרגל, ויזואלית למצוא הערכה של תקופת הזמן בין האקסטרמה של מהירות הרגל ואת ההתחלה של הדופק (שהיא הנקודה הראשונה עם מהירות הרגל אפס לפני מהירות שיא). עבור הפרטורציה במחקר זה, הנקודה הזאת התרחשה 25 ms לפני מהירות אקסטרמה נמצא בשלב 2.2.2.2.
      6. עבור כל פעימה, לחשב את מומנט הרקע של הקרסול כמו הממוצע של מומנט הקרסול של 25 ms לפני תחילת הדופק, כלומר, את הממוצע של מומנט בקטע החל 50 ms עד 25 ms לפני אקסטרמה מהירות. עשה זאת עבור הפולסה -k עם מהירות חיובית באמצעות הפקודה הבאה:
        figure-protocol-51
        figure-protocol-52
        figure-protocol-53
        הערה: הדבר נעשה עבור המהירויות המרביות והמינימליות (מהירות הרגל שלילית) נמצא בשלב 2.2.2.2.
      7. מצא את המינימום והמקסימום של כל ממדי הרקע עבור כל הפולסים, באמצעות הפקודה הבאה:
        figure-protocol-54
        figure-protocol-55
      8. עבור כל פעימה, לחלץ את נתוני מומנט של 65 ms לאחר הדופק להתחיל (כמו קטע מומנט הפנימי), באמצעות הפקודה הבאה:
        figure-protocol-56
        figure-protocol-57
        הערה: זה נעשה גם עבור הנגזרת הראשונה והשנייה של מומנט הקרסול (כדי לספק את הפעם הראשונה והשנייה של הפיתול הפנימי), כמו גם, זווית כף הרגל, מהירות הרגל, האצת כף הרגל, ואת הרגל נתירה.
      9. חשב את השינוי במקטע הפיתולהפנימי ה -k מהערך ההתחלתי שלו, באמצעות הפקודה הבאה:
        figure-protocol-58
        הערה: פעולה זו מתבצעת באופן דומה לזווית הרגל כדי להשיגfigure-protocol-59.
      10. לחלק את טווח מומנט (מושגת בשלב 2.2.2.7) לתוך 3 לסלים רחב ננומטר ולמצוא את הפולסים עם מומנט הרקע של כל סל.
        הערה: פעולה זו מתבצעת באמצעות הפונקציה "find" ויצירת אינדקס. מניחים כי הנוקשות הפנימית היא קבועה בכל סל, מאז מומנט הרקע של הקרסול אינו משתנה באופן משמעותי.
      11. העריכו את הפרמטרים של הנוקשות הפנימית של המודל הפנימי המורחב (היים)29,עבור ה-j סלfigure-protocol-60באמצעות פולסים בקבוצה j ().
        1. לשרשר את כל התגובות מומנט הפנימישל ה-j סל figure-protocol-61 כדי ליצור את הווקטור:
          figure-protocol-62
          איפה figure-protocol-63 הוא ה-i (figure-protocol-64) מומנט התגובה הפנימית בקבוצה j.
          הערה: באופן דומה, לשרשר זווית הרגל, מהירות, והאצת, ואת הנגזרות הראשון והשני של מומנט הפנימי שלהקבוצה j לשמש בשלב 2.2.2.11.2.
        2. מניחים את זווית כף הרגל, מהירות, האצה, נתירה, כמו גם את הנגזרת הראשונה והשנייה של מומנט של הקבוצה j יחד כדי ליצור את המטריצה regressor:
          figure-protocol-65
        3. מצא את הפרמטרים הפנימיים של הקשיות עבורהקבוצה j באמצעות האופרטור קו נטוי הלוכסן (\):
          figure-protocol-66
        4. לחלץ את האלמנט figure-protocol-67 הרביעי כמו קשיותfigure-protocol-68נמוכה בתדר הפנימי.
      12. בצע שלבים במקטע 2.2.2.11 עבור כל הקבוצות (סלים) והערכת הקשיחות הפנימית המתאימה בתדר נמוך.
      13. חלק את כל ערכי הקשיחות בתדר נמוך לפי הנוקשות הקריטית של הנבדק:
        figure-protocol-69
        כאשר m הוא מסה של הנושא, g היא האצת כבידה, והוא figure-protocol-70 הגובה של מרכז הגוף של המסה מעל ציר הקרסול של סיבוב, נגזר נתונים אנתרופומטרים30. זה נותן את הקשיחות המנורמלת (figure-protocol-71).
      14. להמיר את הקרסול רקע מומנט לקרסול רקעfigure-protocol-72שוטר העמדה () על ידי חלוקת הרקע הקרסול torques עם כוחות אנכיים נמדד המקביל.
      15. מגרש figure-protocol-73 כפונקציה של מרכז הלחץ.
        figure-protocol-74
        איפה
        figure-protocol-75
        figure-protocol-76

הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.

תוצאות

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

אותות מדומים ברצף (PRT) ושידורי לטרעון

איור 2A מציג אות prt, אשר מופק על ידי שילוב של פרופיל מהירות אקראית מדומה. עבור כל אחת מהדוגמאות figure-results-1 , מהירות האות עשויה להיות שווה לאפס, או לרכוש ערך חיובי או figure-results-2 שלילי מוגדר מראש. על-

הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.

דיון

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

מספר שלבים הם קריטיים בביצוע ניסויים אלה כדי ללמוד שליטה בפני האדם. שלבים אלה משויכים למדידה הנכונה של האותות וכוללים: 1. יישור נכון של ציר הקרסול השאנק של הרוטציה לדוושות, למדידת המדידה הנכונה של הקרסול. 2) הגדרת מפלסי הטווח הנכונים כדי להבטיח שהם פועלים בטווח שלהם ואינם רוויים במהלך הניסויים. 3) מדידה של EMG עם איכות טובה ודיבור צולב מינימלי. 4) יישום של רטבאליות מתאים, שעוררו תגובות מספיקות, אך לא משבשים את השליטה הרגילה. 5) מבחר של אורך ניסיון מתאים, המבוסס על הניתוח המתוכנן, תוך הימנעות משינוי הגוף והעייפות. בנוסף ל...

הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.

גילויים

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

. למחברים אין מה לגלות

תודות

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

מאמר זה התאפשר על ידי מענק NPRP #6-463-2-189 מן המחקר הלאומי קטאר מענק מגב #81280 מכונים הקנדי של מחקר הבריאות.

הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.

חומרים

רשימת החומרים שנעשה בהם שימוש במאמר זה
שםחברהמספר קטלוגהערות
פוטנציומטר 5KMaurey112P19502מודד זווית פיר מפעיל
8 ערוצים משטח Bagnoli מגברי EMG ואלקטרודותDelsysמודד את ה-EMG של שרירי הקרסול
מחשב נייד AlienWareDell Inc.P69F001-כומר A02מחשב נייד מוכן ל-VR
כרטיס רכישת נתוניםמכשירים לאומיים4472דוגם את האותות האנלוגיים מהחיישנים
שסתום כיווניREXROTH4WMR10C3Xעוקף את הזרימה אם זווית פיר המפעיל חורגת מ-± 20°
רתמת גוף מלאהJelco740הגן על הנבדקים מפני נפילה
מד טווח לייזרמיקרו-אפסילון 1302-1001507307מודד תזוזה ליניארית של השוק
מד טווח לייזרמיקרו-אפסילון 1302-2001509074מודד תזוזה ליניארית של הגוף
תא עומסאומגהLC302-100מודד כוחות
תגובה אנכיתשסתום סרוו פרופורציונליMOOGD681-4718שולט בזרימה ההידראולית למפעילים הסיבוביים
מפעיל סיבוביRotac26R21VDEISFTFLGMTGמיישם הפרעות מכניות
מתמר מומנטLebow2110-5kמודד מומנט קרסול
סביבה וירטואלית עוקבי תנועהHTC inc.1551984681עוקב אחר תנועת הראש
אוזניות מציאות מדומהHTC inc.1551984681מספק הפרעות חזותיות

מקורות

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. Horak, F. B. Postural orientation and equilibrium: what do we need to know about neural control of balance to prevent falls? Age and Ageing. 35, 7-11 (2006).
  2. Morasso, P. G., Schieppati, M. Can muscle stiffness alone stabilize upright standing? Journal of Neurophysiology. 82 (3), 1622-1626 (1999).
  3. Kearney, R. E., Hunter, I. W. System identification of human joint dynamics. Critical Reviews in Biomedical Engineering. 18 (1), 55-87 (1990).
  4. Mirbagheri, M. M., Barbeau, H., Kearney, R. E. Intrinsic and reflex contributions to human ankle stiffness: variation with activation level and position. Experimental Brain Research. 135 (4), 423-436 (2000).
  5. Weiss, P. L., Hunter, I. W., Kearney, R. E. Human ankle joint stiffness over the full range of muscle activation levels. Journal of Biomechanics. 21 (7), 539-544 (1988).
  6. Golkar, M. A., Sobhani Tehrani, E., Kearney, R. E. Linear Parameter Varying Identification of Dynamic Joint Stiffness during Time-Varying Voluntary Contractions. Frontiers in Computational Neuroscience. 11, 35(2017).
  7. Weiss, P. L., Kearney, R. E., Hunter, I. W. Position dependence of ankle joint dynamics--I. Passive mechanics. Journal of Biomechanics. 19 (9), 727-735 (1986).
  8. Weiss, P. L., Kearney, R. E., Hunter, I. W. Position dependence of ankle joint dynamics--II. Active mechanics. Journal of Biomechanics. 19 (9), 737-751 (1986).
  9. Engelhart, D., Boonstra, T. A., Aarts, R. G. K. M., Schouten, A. C., van der Kooij, H. Comparison of closed-loop system identification techniques to quantify multi-joint human balance control. Annual Reviews in Control. 41, 58-70 (2016).
  10. Kiemel, T., Elahi, A. J., Jeka, J. J. Identification of the plant for upright stance in humans: multiple movement patterns from a single neural strategy. Journal of Neurophysiology. 100 (6), 3394-3406 (2008).
  11. Loram, I. D., Lakie, M. Direct measurement of human ankle stiffness during quiet standing: the intrinsic mechanical stiffness is insufficient for stability. Journal of Physiology-London. 545 (3), 1041-1053 (2002).
  12. Fitzpatrick, R., Burke, D., Gandevia, S. C. Loop gain of reflexes controlling human standing measured with the use of postural and vestibular disturbances. Journal of Neurophysiology. 76 (6), 3994-4008 (1996).
  13. Dakin, C. J., Son, G. M. L., Inglis, J. T., Blouin, J. S. Frequency response of human vestibular reflexes characterized by stochastic stimuli. The Journal of Physiology. 583 (3), 1117-1127 (2007).
  14. Vlutters, M., Boonstra, T. A., Schouten, A. C., vander Kooij, H. Direct measurement of the intrinsic ankle stiffness during standing. Journal of Biomechanics. 48 (7), 1258-1263 (2015).
  15. Casadio, M., Morasso, P. G., Sanguineti, V. Direct measurement of ankle stiffness during quiet standing: implications for control modelling and clinical application. Gait and Posture. 21 (4), 410-424 (2005).
  16. Sakanaka, T. E. Causes of Variation in Intrinsic Ankle Stiffness and the Consequences for Standing. , University of Birmingham. Doctoral dissertation (2017).
  17. Sakanaka, T. E., Lakie, M., Reynolds, R. F. Sway-dependent changes in standing ankle stiffness caused by muscle thixotropy. Journal of Physiology. 594 (3), 781-793 (2016).
  18. Peterka, R. J., Murchison, C. F., Parrington, L., Fino, P. C., King, L. A. Implementation of a Central Sensorimotor Integration Test for Characterization of Human Balance Control During Stance. Frontiers in Neurology. 9, 1045(2018).
  19. Engelhart, D., Schouten, A. C., Aarts, R. G., van der Kooij, H. Assessment of Multi-Joint Coordination and Adaptation in Standing Balance: A Novel Device and System Identification Technique. IEEE Transactions on Neural Systems and Rehabilitation Engineering. 23 (6), 973-982 (2015).
  20. Boonstra, T. A., Schouten, A. C., van der Kooij, H. Identification of the contribution of the ankle and hip joints to multi-segmental balance control. Journal of Neuroengineering and Rehabilitation. 10, 23(2013).
  21. Forster, S. M., Wagner, R., Kearney, R. E. A bilateral electro-hydraulic actuator system to measure dynamic ankle joint stiffness during upright human stance. Proceedings of the 25th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. , Cancun, Mexico. (2003).
  22. Davis, J., Hsieh, Y. -H., Lee, H. -C. Humans perceive flicker artifacts at 500 Hz. Scientific Reports. 5, 7861(2015).
  23. Amiri, P., Kearney, R. E. Ankle intrinsic stiffness changes with postural sway. Journal of Biomechanics. 85, 50-58 (2019).
  24. van der Kooij, H., van Asseldonk, E., van der Helm, F. C. Comparison of different methods to identify and quantify balance control. Journal of Neuroscience Methods. 145 (1-2), 175-203 (2005).
  25. Amiri, P., MacLean, L. J., Kearney, R. E. Measurement of shank angle during stance using laser range finders. International Conference of the IEEE Engineering in Medicine and Biology. , Orlando, FL. (2016).
  26. The SENIAM project. , Available from: http://www.seniam.org/ (2019).
  27. Jalaleddini, K., Tehrani, E. S., Kearney, R. E. A Subspace Approach to the Structural Decomposition and Identification of Ankle Joint Dynamic Stiffness. IEEE Transactions on Biomedical Engineering. 64 (6), 1357-1368 (2017).
  28. Amiri, P., Kearney, R. E. A Closed-loop Method to Identify EMG-Ankle Torque Dynamic Relation in Human Balance Control. Conference Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society. , Berlin, Germany. (2019).
  29. Sobhani Tehrani, E., Jalaleddini, K., Kearney, R. E. Ankle Joint Intrinsic Dynamics is More Complex than a Mass-Spring-Damper Model. IEEE Transactions on Neural Systems and Rehabilitation Engineering. 25 (9), 1568-1580 (2017).
  30. NASA. Anthropometry and biomechanics. , Available from: http://msis.jsc.nasa.gov/sections/section03.htm (1995).
  31. Peterka, R. J. Sensorimotor integration in human postural control. Journal of Neurophysiology. 88 (3), 1097-1118 (2002).
  32. Amiri, P., Kearney, R. E. Ankle intrinsic stiffness is modulated by postural sway. Conference Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society. , Seogwipo, South Korea. (2017).
  33. Jeka, J. J., Allison, L. K., Kiemel, T. The dynamics of visual reweighting in healthy and fall-prone older adults. Journal of Motor Behavior. 42 (4), 197-208 (2010).
  34. Jilk, D. J., Safavynia, S. A., Ting, L. H. Contribution of vision to postural behaviors during continuous support-surface translations. Experimental Brain Research. 232 (1), 169-180 (2014).
  35. Winter, D. A., Patla, A. E., Prince, F., Ishac, M., Gielo-Perczak, K. Stiffness control of balance in quiet standing. Journal of Neurophysiology. 80 (3), 1211-1221 (1998).
  36. Pasma, J. H., Boonstra, T. A., van Kordelaar, J., Spyropoulou, V. V., Schouten, A. C. A Sensitivity Analysis of an Inverted Pendulum Balance Control Model. Frontiers in Computational Neuroscience. 11, 99(2017).
  37. Pasma, J. H., et al. Changes in sensory reweighting of proprioceptive information during standing balance with age and disease. Journal of Neurophysiology. 114 (6), 3220-3233 (2015).
  38. Pasma, J. H., et al. Impaired standing balance: The clinical need for closing the loop. Neuroscience. , 157-165 (2014).
  39. Engelhart, D., et al. Impaired Standing Balance in Elderly: A New Engineering Method Helps to Unravel Causes and Effects. Journal of the American Medical Directors Association. 15 (3), (2014).
  40. Pasma, J. H., Boonstra, T. A., Campfens, S. F., Schouten, A. C., Van der Kooij, H. Sensory reweighting of proprioceptive information of the left and right leg during human balance control. Journal of Neurophysiology. 108 (4), 1138-1148 (2012).
  41. Goodworth, A. D., Peterka, R. J. Sensorimotor integration for multisegmental frontal plane balance control in humans. Journal of Neurophysiology. 107 (1), 12-28 (2012).
  42. Kiemel, T., Zhang, Y., Jeka, J. J. Identification of neural feedback for upright stance in humans: stabilization rather than sway minimization. Journal of Neuroscience. 31 (42), 15144-15153 (2011).
  43. van der Kooij, H., van Asseldonk, E. H. F., Geelen, J., van Vugt, J. P. P., Bloem, B. R. Detecting asymmetries in balance control with system identification: first experimental results from Parkinson patients. Journal of Neural Transmission. 114 (10), 1333(2007).
  44. Fujisawa, N., et al. Human standing posture control system depending on adopted strategies. Medical and Biological Engineering and Computing. 43 (1), 107-114 (2005).
  45. Johansson, R., Magnusson, M., Fransson, P. A., Karlberg, M. Multi-stimulus multi-response posturography. Mathematical Biosciences. 174 (1), 41-59 (2001).
  46. Jeka, J., Oie, K., Schöner, G., Dijkstra, T., Henson, E. Position and Velocity Coupling of Postural Sway to Somatosensory Drive. Journal of Neurophysiology. 79 (4), 1661-1674 (1998).
  47. Peterka, R. J., Benolken, M. S. Role of somatosensory and vestibular cues in attenuating visually induced human postural sway. Experimental Brain Research. 105 (1), 101-110 (1995).
  48. Maki, B. E., Fernie, G. R. A system identification approach to balance testing. Progress in Brain Research. 76, 297-306 (1988).
  49. Johansson, R., Magnusson, M., Akesson, M. Identification of human postural dynamics. IEEE Transactions on Biomedical Engineering. 35 (10), 858-869 (1988).
  50. Maki, B. E., Holliday, P. J., Fernie, G. R. A Posture Control Model and Balance Test for the Prediction of Relative Postural Stability. IEEE Transactions on Biomedical Engineering. BME-34. 10 (10), 797-810 (1987).
  51. Werness, S. A., Anderson, D. J. Parametric analysis of dynamic postural responses. Biological Cybernetics. 51 (3), 155-168 (1984).
  52. Hwang, S., Agada, P., Kiemel, T., Jeka, J. J. Identification of the Unstable Human Postural Control System. Frontiers in Systems Neuroscience. 10, 22(2016).
  53. Ishida, A., Imai, S., Fukuoka, Y. Analysis of the posture control system under fixed and sway-referenced support conditions. IEEE Transactions on Biomedical Engineering. 44 (5), 331-336 (1997).
  54. Ishida, A., Miyazaki, S. Maximum likelihood identification of a posture control system. IEEE Transactions on Biomedical Engineering. 34 (1), 1-5 (1987).
  55. Ljung, L. System Identification: Theory for the User. , Prentice-Hall, Inc. Upper Saddle River, NJ. (1986).
  56. Forssell, U., Ljung, L. Closed-loop identification revisited. Automatica. 35 (7), 1215-1241 (1999).
  57. Horak, F. B., Nashner, L. M. Central programming of postural movements: adaptation to altered support-surface configurations. Journal of Neurophysiology. 55 (6), 1369-1381 (1986).

הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.

הדפסות חוזרות והרשאות

בקש הרשאה לשימוש חוזר בטקסט או באיורים של מאמר JoVE זה

בקש הרשאה

תגיות

מאמרים קשורים