מתנד הרמוני

מַתְנֵד הרמוני (או אוסצילטור הרמוני, מאנגלית: Harmonic oscillator) הוא מערכת מכנית שבה פועל על גוף נתון כוח מתכונתי (יחסי) להעתק הגוף ובכיוון מנוגד לו: $\ {\vec {F}}=-K{\vec {\Delta x}}$, כאשר $\ {\vec {F}}$ הוא הכוח, $\ K$ הוא קבוע המתנד ו- $\ {\vec {\Delta x}}$ הוא ההעתק. כוח זה נקרא גם כוח מחזיר שמשמעותו היא כוח שתמיד פועל לכיוון נקודת שיווי המשקל.

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

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

מתנד הרמוני פשוט

מתנד הרמוני פשוט, הוא מערכת בה פועל כח מחזיר בלבד, שכן כשמו הוא מנסה להחזיר את המערכת לנקודת שיווי המשקל שלה, לאחר שזו סטתה ממנו. באופן כללי משוואת המתנד ההרמוני היא משוואה דיפרנציאלית מהצורה ${\ddot {\xi }}+\omega ^{2}\xi =0$ כאשר $\xi$הוא הגודל המחזורי, ו${\ddot {\xi }}$ הוא הנגזרת השנייה בזמן של הגודל המחזורי.

כאן נסקור את בעיית המתנד ההרמוני במכניקה הקלאסית. כח מחזיר, הוא כח מהצורה $F=-kx$ כאשר $k$ הוא קבוע המאפיין את המערכת כמו קבוע הקפיץ בחוק הוק, ו$x$ הוא העתק המסה מנקודת שיווי המשקל.

ניתן לרשום את החוק השני של ניוטון עבור המערכת הבנויה ממסה המחוברת לקפיץ, תוך כדי כתיבת התאוצה בדרך מעט שונה $a={\ddot {x}}$:

$-kx=m{\ddot {x}}$. נשים לב, כי המסה תסומן באות $m$. כעת, נחלק את שני אגפי המשוואה ב$m$ ,נגדיר $\omega _{0}^{2}={\frac {k}{m}}$ ונעביר אגפים כדי לקבל את המשוואה ${\ddot {x}}+\omega _{0}^{2}x=0$.

נשים לב לנקודות הבאות:

• הגודל $\omega _{0}$ חיובי, היות שהמסה חיובית תמיד, וכך גם קבוע הקפיץ
• הגודל $\omega _{0}$ הוא גודל קבוע, היות שהמסה או קבוע הקפיץ אינם משתנים (הערה: אם המסה משתנה, משוואת הכוחות שכתבנו קודם לכן, אינה מתאימה ויש להיעזר בפיתוח המשוואה למסה משתנה)

משוואה זו היא משוואה דיפרנציאלית ליניארית הומוגנית מסדר שני, שכדי לפתור אותה, נכתוב את הפולינום האופייני של המערכת: $\lambda ^{2}+\omega _{0}^{2}=0$. משוואת הפולינום האופייני תיפתר על ידי נוסחת השורשים, ושורשיה נתונים על ידי $\lambda =\pm \omega _{0}i$ כאשר $i$ הוא היחידה המדומה. פתרון המשוואה הדיפרנציאלית, יהיה מהצורה $x(t)=\alpha e^{\lambda _{1}t}+\beta e^{\lambda _{2}t}$. על ידי הצבת הפתרונות של הפולינום האופייני נקבל $x(t)=\alpha e^{\omega _{0}it}+\beta e^{-\omega _{0}it}$. וכעת ישנה בעיה חדשה. הפתרון אינו פונקציה ממשית במלואה. כלומר, חלק מהפתרון מרוכב. העתק המסה הוא גודל ממשי לחלוטין, ולכן עלינו לקחת רק את החלק הממשי של הפתרון. לשם כך, נעזר בנוסחת אוילר כדי לפשט את הפתרון:

$$x(t)=\alpha e^{\omega _{0}it}+\beta e^{-\omega _{0}it}=\alpha (\cos(\omega _{0}t)+i\sin(\omega _{0}t))+\beta (\cos(\omega _{0}t)-i\sin(\omega _{0}t))=\alpha \cos(\omega _{0}t)+\beta \cos(\omega _{0}t)+\alpha i\sin(\omega _{0}t)-\beta i\sin(\omega _{0}t)=\cos(\omega _{0}t)(\alpha +\beta )+i\sin(\omega _{0}t)(\alpha -\beta )$$

כעת ניתן להפריד את החלק הממשי והחלק המדומה. נשים לב כי ${\mathfrak {R}}{\mathfrak {e}}\{x(t)\}=\cos(\omega _{0}t)(\alpha +\beta )$ ו${\mathfrak {I}}{\mathfrak {m}}\{x(t)\}=\sin(\omega _{0}t)(\alpha -\beta )$. כפי שאמרנו, ניקח רק את החלק הממשי. נגדיר $\alpha +\beta =A$ ונכתוב את הפתרון אליו הגענו: $x(t)=A\cos(\omega _{0}t)$. כדי להגיע לפתרון המלא אשר יתאר באופן יותר כללי ונרחב את המתנד ההרמוני, נוסיף הזזה אפשרית בתוך הפונקציה, ונאמר שהפתרון הכללי הוא $x(t)=A\cos(\omega _{0}t+\phi )$.

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

הפאזה שתסומן באות $\phi$, תימצא גם היא מתנאי ההתחלה. גודל זה מאפיין את ההזזה של פונקציית הקוסינוס. מבחינה פיזיקלית, הפאזה תופיע, אם מדידת הזמן החלה לא בדיוק ברגע בו שוחררה המסה ממנוחה.

כעת נבחן מקרה פשוט יותר של המתנד, בו אין פאזה. במקרה כזה, העתק המסה נתון על ידי המשוואה $x(t)=A\cos(\omega _{0}t)$. היות שהפונקציה המתארת את העתק המסה היא פונקציה טריגונומטרית כלשהי, ניתן להסיק כי התנועה מחזורית. אם כך, קיים זמן $T$ מינימלי, שדרוש כדי שהמערכת תשלים מחזור אחד (תבצע תנועה אחת שלמה) ותחזור לאותה הנקודה ממנה התחילה. המחזור של פונקציית הקוסינוס הוא $2\pi$, ולכן נדרוש שכאשר נציב את $T$ לתוך הפונקציה, הארגומנט של הפונקציה, יהיה שווה ל$2\pi$. כלומר, נדרוש שיתקיים $2\pi =\omega _{0}T$.

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

מכאן נקבל את הזמן המינימלי הדרוש למערכת כדי לבצע תנועה אחת שלמה ולחזור לנקודה ממנה החלה התנועה, ומעתה הוא יקרא בשם זמן מחזור: $T={\frac {2\pi }{\omega _{0}}}$.

וכן, התדירות של המערכת, לפי הגדרתה $f={\frac {1}{T}}$ תוצג כ$f={\frac {\omega _{0}}{2\pi }}$. בגלל הקשר המיוחד בין התדירות לגודל $\omega _{0}$, הוא יקרא מעתה והלאה בשם תדירות זוויתית. התדירות הזוויתית, וכך גם התדירות נקבעות אך ורק על ידי מאפייני המערכת, ואינם תלויות בזמן או בגדלים משתנים אחרים).

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

$x(t)=A\cos(\omega _{0}t+\phi )$

$v(t)=-A\omega _{0}\sin(\omega _{0}t+\phi )$

$a(t)=-A\omega _{0}^{2}\cos(\omega _{0}t+\phi )$

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

$$E(t)={\frac {1}{2}}mv^{2}+{\frac {1}{2}}kx^{2}={\frac {1}{2}}m(-A\omega _{0}\sin(\omega _{0}t+\phi ))^{2}+{\frac {1}{2}}k(A\cos(\omega _{0}t+\phi ))^{2}={\frac {1}{2}}mA^{2}\omega _{0}^{2}\sin ^{2}(\omega _{0}t+\phi )+{\frac {1}{2}}kA^{2}\cos ^{2}(\omega _{0}t+\phi )={\frac {1}{2}}A^{2}m{\frac {k}{m}}\sin ^{2}(\omega _{0}t+\phi )+{\frac {1}{2}}kA^{2}\cos ^{2}(\omega _{0}t+\phi )={\frac {1}{2}}kA^{2}(\sin ^{2}(\omega _{0}t+\phi )+\cos ^{2}(\omega _{0}t+\phi ))={\frac {1}{2}}kA^{2}$$

קיבלנו כי האנרגיה לא באמת תלויה בזמן, אלא בכל נקודה בזמן היא שווה לגודל קבוע. מכאן, שהאנרגיה נשמרת.

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

מתוך משוואות שימור האנרגיה, ניתן לראות שכאשר כל האנרגיה הפוטנציאלית מתאפסת, האנרגיה הקינטית מקסימלית. למעשה, הנקודה בה האנרגיה הפוטנציאלית מתאפסת, היא בדיוק נקודת שיווי המשקל (באנלוגיה של הקפיץ, האנרגיה הפוטנציאלית מתאפסת כאשר הקפיץ רפוי, כלומר כאשר המערכת בשיווי משקל). לכן, המהירות מקסימלית, בנקודת שיווי המשקל.

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

מתנד הרמוני מרוסן

מערכת של מסה, קפיץ ובוכנה המפעילה כוח חיכוך מתכונתי למהירות התנועה

בתיאור הקודם, הנחנו כי אין חיכוך במערכת. נראה איך מתנהג המתנד במקרה בו יש חיכוך.

חיכוך הוא כוח מרסן נפוץ בטבע - חיכוך עם האוויר או עם נוזל, התנגדות של מעגל חשמלי, ועוד. כוח משכך זה תלוי ביחס ישר במהירות (נקרא גם כוח סטוקס). נדמה כי כוח החיכוך נגרם על ידי בוכנת אוויר, כבאיור, ונניח כי הכוח מתכונתי למהירות הגוף ונתון על ידי $\ F_{friction}=-bv=-b{\dot {x}}$, כאשר b הוא מקדם הפרופורציה.

כעת משוואת הכוחות נראית כך: $F=-kx-b{\dot {x}}=m{\ddot {x}}$. על ידי סידור אלגברי מקבלים $\ m{\ddot {x}}+b{\dot {x}}+kx=0$ או

$\ {\ddot {x}}+{\frac {b}{m}}{\dot {x}}+\omega _{0}^{2}x=0$, כאשר $\ \omega _{0}={\sqrt {\frac {k}{m}}}$, כמקודם.

מתנד הרמוני פשוט הוא מקרה פרטי של מתנד הרמוני מרוסן, בו $\ b=0$.

בהנחה שהחיכוך קטן מספיק, הפתרון הכללי למשוואה זו הוא:

$\ x(t)=X_{0}e^{-t/\tau }\cos(\omega t+\phi )$, כאשר:

• $\ \tau ={\frac {2m}{b}}$, הוא הזמן האופייני של המערכת. זהו הזמן בו משרעת התנודה דועכת בפקטור e.
• $\omega ={\sqrt {\omega _{0}^{2}-{\frac {1}{\tau ^{2}}}}}$, היא תדירות התנועה החדשה.
• הביטוי $\ X_{0}e^{-{\frac {t}{\tau }}}$, הוא משרעת התנודה, הדועכת באופן מעריכי עם הזמן.

פתרון מפורט מופיע בהמשך הערך. זוהי תנודה הרמונית, המתבצעת בתדירות נמוכה מעט מהתדירות הטבעית (תמיד מתקיים $\ \omega _{0}>\omega$). במערכת זו, בניגוד לתנודה ההרמונית הפשוטה, יש איבוד אנרגיה מתמיד לחיכוך.

אם החיכוך גדול מדי, התנועה תדעך מהר מאוד, באופן מעריכי, כך שלא יסתיים אפילו מחזור תנודה יחיד. קיים ערך מקסימלי של שיכוך, שבערכי שיכוך גדולים ממנו לא יתרחשו כלל תנודות. ערך זה נקרא השיכוך הקריטי, והוא נתון על ידי $\ b_{critical}=2{\sqrt {km}}$.

בנוגע לטיפול במתנד הרמוני מאולץ, ראו: תהודה.

בנוגע לטיפול במתנד הרמוני לפי מכניקת הקוונטים ראו: מתנד הרמוני קוונטי.

דוגמאות למערכות שמקיימות תנודות הרמוניות

מטוטלת

ערך מורחב – מטוטלת מתמטית

מטוטלת מתמטית פשוטה היא מערכת שמקיימת תנועה הרמונית רק בקירוב. תנועתה המלאה היא מחזורית וסוטה במקצת מתנועה הרמונית, אך בקירוב זוויות קטנות אפשר להתייחס לתנועתה כתנועה הרמונית.

נניח שמסת המטוטלת היא m והיא מחוברת לתקרה בחוט אידיאלי חסר משקל באורך $\ l$ בתוך שדה כבידה אחיד g. נסמן ב $\theta$ את הזווית בין המטוטלת לאנך.

הכוח המחזיר יהיה הרכיב של כוח הכובד בכיוון המשיק לתנועה הרדיאלית של המטוטלת, כלומר: $\ F=-mg\sin {(\theta )}$. לכן, משוואת התנועה תהיה:

$m{\ddot {l\theta }}=-mg\sin \theta$

בקירוב של זוויות קטנות $\sin {\theta }\approx \theta$, ומאחר ש-$\ l$ נשאר קבוע, מקבלים את משוואת התנועה:

${\ddot {\theta }}=-{\frac {g}{l}}\theta$

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

מעגל LC

מעגל LC אידיאלי (ללא התנגדות) - זהו מעגל המכיל קבל וסליל, ובו האנרגיה האלקטרומגנטית עוברת בין הקבל לסליל באופן מחזורי הרמוני.

מכך שהמתח הכללי על המעגל הוא אפס (על פי חוקי קירכהוף) ניתן להסיק כי:

$\ L{dI \over dt}+{q \over C}=0$

כאשר:

• q הוא המטען שאגור בקבל.
• C הוא הקיבול של הקבל.
• I הוא הזרם במעגל.
• L הוא ההשראות של הסליל.

ומאחר ש $\ I={\frac {dq}{dt}}$, נובע ש

$\ L{\frac {d^{2}q}{dt^{2}}}+{\frac {q}{C}}=0$, ולכן $\ {\frac {d^{2}q}{dt^{2}}}+{q \over LC}=0$ וכאן תדירות התנועה היא $\omega ={\sqrt {1 \over LC}}$.

מעגל RLC, הוא אותו מעגל, בתוספת נגד. מעגל שכזה הוא מתנד הרמוני משוכך, כיוון שהוא מקיים את המשוואה:

$L{\frac {d^{2}q}{dt^{2}}}+R{\frac {dq}{dt}}+{\frac {q}{c}}=0$.

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

מתנד מכני	מתנד חשמלי
מסה - m	השראות הסליל- L
קבוע הקפיץ - k	1 חלקי הקיבול - ${\frac {1}{C}}$
העתק - x	מטען חשמלי - q
מהירות - v	זרם - I
שיכוך - b	התנגדות - R
$\ \omega _{0}={\sqrt {\frac {k}{m}}}$	$\ \omega _{0}={\sqrt {\frac {1}{LC}}}$
$\tau ={\frac {2m}{b}}$	$\tau ={\frac {2L}{R}}$

פתרון המשוואות

הניתוח המובא פה הוא למקרה המתואר למעלה: מסה $\ m$ מחוברת לקפיץ בעל קבוע $\ K$, עם מקדם חיכוך $\ b$ (כאמור, מתנד הרמוני פשוט הוא מקרה פרטי שבו מתקיים $\ b=0$) . המשוואה המתארת את המצב היא:

$\ {\ddot {x}}+{\frac {b}{m}}{\dot {x}}+{\frac {K}{m}}x=0$

זוהי משוואה דיפרנציאלית הומוגנית וליניארית, מסדר שני. במשוואות מסוג זה, מרחב הפונקציות המקיימות את המשוואה הוא מרחב וקטורי מממד שתיים. לכן מספיק למצוא שתי פונקציות (בלתי תלויות) שמקיימות את המשוואה - הן יהיו בסיס שלם למרחב כל הפתרונות האפשריים. מנחשים פתרון מהסוג $\ x(t)=Ae^{\lambda t}$, (המקדם $\lambda$ יכול להיות מרוכב). כאשר מציבים את הפונקציה הזו במשוואה, ולפי תכונת האקספוננט, $\ {\frac {d}{dt}}e^{\lambda t}=\lambda e^{\lambda t}$, מתקבל פולינום ממעלה שנייה שפתרונו פשוט: $\ \lambda ^{2}+{\frac {b}{m}}\lambda +{\frac {K}{m}}=0$.

זוהי משוואה ריבועית שפתרונותיה הם $\ \lambda _{1,2}=-\alpha \pm {\sqrt {\alpha ^{2}-\omega _{0}^{2}}}$. כאן $\ \omega _{0}={\sqrt {\frac {K}{m}}},\ \alpha ={\frac {b}{2m}}$. מכיוון שהפתרונות יכולים להיות מרוכבים, יש תמיד פתרון למשוואה. אם אין חיכוך, $\ b=0$ ואז $\ \lambda =i\omega _{0}$.

כעת ישנם שלושה מקרים מובחנים:

• אם $\alpha ^{2}-\omega _{0}^{2}=0$, ישנו פתרון ממשי שלילי עבור $\ \lambda$. הפתרון למקרה הכללי הוא חיבור של שתי פונקציות אקספוננציאליות דועכות:

$(\alpha ,\omega _{0}>0),\ \ X(t)=Ae^{-\alpha t}+Bte^{-\alpha t}=Ae^{-\omega _{0}t}+Bte^{-\omega _{0}t}$

זהו מצב של שיכוך-קריטי, שבו לא יתרחשו תנודות, אלא התנועה תדעך לאפס. התנאי לשיכוך קריטי הוא שהדיסקרימיננטה תהיה אפס, משמע $\ {\frac {b^{2}}{4m^{2}}}-{\frac {K}{m}}=0$, כלומר $\ b^{2}=4Km$.

• אם $\alpha ^{2}-\omega _{0}^{2}>0$, ישנם שני פתרונות ממשיים שליליים עבור $\ \lambda$. הפתרון למקרה הכללי הוא חיבור של שני פונקציות אקספוננציאליות דועכות:

$(\lambda _{1},\lambda _{2}<0),\ \ X(t)=Ae^{\lambda _{1}t}+Be^{\lambda _{2}t}$

זהו מצב של שיכוך-יתר, שבו לא יתרחשו תנודות, אלא התנועה תדעך לאפס. התנאי לשיכוך יתר הוא שהדיסקרימיננטה תהיה חיובית, משמע $\ {\frac {b^{2}}{4m^{2}}}-{\frac {K}{m}}>0$, כלומר $\ b^{2}>4Km$.

• אם $\alpha ^{2}-\omega _{0}^{2}<0$, מתקבלים שני פתרונות מרוכבים צמודים. זהו המקרה המעניין יותר, וממילא גם השכיח יותר, שבו השיכוך קטן ולכן מתרחשות תנודות הרמוניות דועכות. בהצבת פתרונות המשוואה הריבועית, ושימוש בנוסחת אוילר, מקבלים את שני הפתרונות למשוואה הדיפרנציאלית:

$\ y_{1}(t)=e^{(-\alpha +i{\sqrt {\omega _{0}^{2}-\alpha ^{2}}})t}=e^{-\alpha t}(\cos({\sqrt {\omega _{0}^{2}-\alpha ^{2}}}\cdot t)+i\sin({\sqrt {\omega _{0}^{2}-\alpha ^{2}}}\cdot t))$

$\ y_{2}(t)=e^{(-\alpha -i{\sqrt {\omega _{0}^{2}-\alpha ^{2}}})t}=e^{-\alpha t}(\cos({\sqrt {\omega _{0}^{2}-\alpha ^{2}}}\cdot t)-i\sin({\sqrt {\omega _{0}^{2}-\alpha ^{2}}}\cdot t))$

לשם נוחות, נסמן $\omega ={\sqrt {\omega _{0}^{2}-\alpha ^{2}}}$.

פתרון המשוואה הכללי הוא מהצורה $\ y(t)=Ay_{1}(t)+By_{2}(t)$. כאן $\ A$ ו-$\ B$ הם קבועים כלשהם, התלויים בתנאי ההתחלה של הבעיה. כיוון שרק פתרון ממשי הוא פיזיקלי, נדרש $\ Ay_{1}={\overline {By_{2}}}$ (הקו העליון מסמן צמוד מרוכב). מאחר שמתקיים ממילא $\ y_{1}={\overline {y_{2}}}$, חייב להתקיים $A={\overline {B}}$. כך מתקבל פתרון ממשי טהור מהצורה

$\ y(t)=(A+{\overline {A}})e^{-\alpha t}\sin(\omega t)+i(A-{\overline {A}})e^{-\alpha t}\cos(\omega t)$.

כל המקדמים בביטוי הם ממשיים, ולכן ניתן לפשט את הסכום לביטוי יחיד המכיל מופע:

$\ y(t)=A_{0}e^{-\alpha t}\sin(\omega t+\phi )$.

כאן $\ A_{0},\phi$ הן המופע והמשרעת ההתחלתיים, הנקבעים על ידי תנאי ההתחלה של הבעיה. פתרון זה מתאר תנודה הרמונית בתדירות זוויתית

$\omega ={\sqrt {\omega _{0}^{2}-\alpha ^{2}}}$

עם משרעת דועכת בזמן. מכיוון שלמקדם $\ \alpha$ יש יחידות של ${\frac {1}{\mbox{sec}}}$, נהוג להגדיר את הזמן האופייני $\ \tau ={\frac {1}{\alpha }}={\frac {2m}{b}}$.

כאשר אין חיכוך, $\ b=0$ ולכן $\ \alpha =0$ ומקבלים את הביטוי המוכר של תנודה הרמונית פשוטה: $\ y(t)=A_{0}\sin(\omega _{0}t+\phi )$.

הכללות של מתנד הרמוני

המתנד ההרמוני מצטיין בתכונה הבסיסית שזמן המחזור של חלקיק שתונד הרמונית אינו תלוי באמפליטודה. לא משנה כמה רחוק החלקיק ישוחרר ממרכז הכוח של המתנד, הוא תמיד ישלים תנודה באותו זמן. תכונה זו אינה נכונה כאשר דנים בחוקי כוח אחרים בהם הכוח לא יחסי למרחק r מהמרכז - משפט ידוע קובע שאם התלות של זמן המחזור באמפליטודה היא מהצורה $T=CA^{n}$ (כאשר A היא האמפליטודה ו-C הוא קבוע פרופורציה) אז הכוח המרכזי הוא מצורה $F=kr^{1-2n}$. למשל אם n = 0 (זמן המחזור לא תלוי באמפליטודה) נקבל: $F=kr$, כלומר נקבל מתנד הרמוני. יותר מכך, אם נציב $n=3/2$ נקבל $F=kr^{-2}$, כלומר נקבל כוח היפוך-ריבוע - במקרה פרטי זה קיבלנו את החוק השלישי של קפלר. חישוב קבוע הפרופורציה C בביטוי לזמן המחזור לפי האמפליטודה לרוב כרוך בחישוב אינטגרל מסובך ולעיתים אין לו פתרון אנליטי.

ראו גם

• תהודה
• תורת התנודות
• אנליזת פורייה: טור פורייה והתמרת פורייה
• מתנד הרמוני קוונטי
• אופני תנודה עצמיים

קישורים חיצוניים

• פתרון המד"ר
• מתנד הרמוני, באתר אנציקלופדיה למתמטיקה (באנגלית)

הערך באדיבות ויקיפדיה העברית, קרדיט,
רשימת התורמים
רישיון cc-by-sa 3.0

34015401מתנד הרמוני