Fachhochschule Nordwestschweiz (FHNW) Hochschule für Technik Institut für Geistes- und Naturwissenschaft Lösung Serie 1 (Differentialgleichungen 1-ter Ordnung) Dozent: Roger Burkhardt Klasse: Studiengang ST Büro: 4.613 Semester: 2 (VZ) Modul: MDS Datum: FS2009 1. Aufgabe Gegeben seien die nachfolgenden gewöhnlichen Differentialgleichungen erster Ordnung (Anfangswertprobleme). Welche der Differentialgleichungen sind separierbar und welche sind linear (a) t2 x0 + tx = e− 2 , x (0) = 1 Lösung: t2 Lineare DGL mit p(t) = t und s(t) = e− 2 . (b) x0 x = t, x (0) = 3 Lösung: Separierbare DGL: t dx = ⇒ xdx = tdt dt x (c) x0 − sin (t) x = cos (t) , x (0) = 1 Lösung: Lineare DGL mit p(t) = − sin(t) und s(t) = cos(t). (d) x0 + x = e−t , x (0) = 0 Lösung: Lineare DGL mit p(t) = 1 und s(t) = e−t (lineare DGL mit konstanten Koeffizienten). (e) x0 + tan (t) x = sin (t) , x (0) = 0 Lösung: Lineare DGL mit p(t) = tan(t) und s(t) = sin(t). (f) x0 = t2 , x2 x (0) = 1 Lösung: Separierbare DGL: dx t2 = 2 ⇒ x2 dx = t2 dt dt x MDS Lösung Serie 1 (Differentialgleichungen 1-ter Ordnung) FS 2009 (g) x0 − 3x = t2 , x (0) = − 2 27 Lösung: Lineare DGL mit p(t) = −3 und s(t) = t2 (lineare DGL mit konstanten Koeffizienten). (h) x0 − √ t= √ tx2 , x (0) = 0 Lösung: Separierbare DGL: √ 1 dx √ = t 1 + x2 ⇒ dx = tdt dt 1 + x2 2. Aufgabe (a) Bestimme die allgemeine und die partikulären Lösungen der separierbaren Differentialgleichungen. Lösung: • zu 1b: – Allgemeine Lösung: Z Z √ t x2 t2 dx = ⇒ xdx = tdt ⇒ = + C ⇒ x = ± C + t2 dt x 2 2 – Partikuläre Lösung: x(0) = 3 ⇒ 3 = √ √ C + 02 ⇒ C = 9 ⇒ x = 9 + t2 • zu 1f: – Allgemeine Lösung: Z Z √ dx t2 x3 t3 3 2 = 2 ⇒ x dx = t2 dt ⇒ = + C ⇒ x = C + t3 dt x 3 3 – Partikuläre Lösung: x(0) = 1 ⇒ 1 = √ √ C + 03 ⇒ C = 1 ⇒ x = 1 + t3 • zu 1h: – Allgemeine Lösung: dx √ = t 1 + x2 ⇒ dt Z √ 1 2√ 3 dx = tdt ⇒ arctan(x) = t +C 1 + x2 3 √ 2 3 ⇒ x = tan t +C 3 Z Seite 2 / 8 MDS Lösung Serie 1 (Differentialgleichungen 1-ter Ordnung) FS 2009 – Partikuläre Lösung: √ √ 2 3 2 3 x(0) = 0 ⇒ 0 = tan 0 + C ⇒ C = 0 ⇒ x = tan t 3 3 (b) Bestimme die allgemeine und die partikulären Lösungen der linearen Differentialgleichungen. Lösung: • zu 1a: – Homogene DGL: Z Z t2 1 t2 0 dx = −tdt ⇒ ln |x| = − + C ⇒ xhom = Ce− 2 x + tx = 0 ⇒ x 2 – Inhomogene DGL (Variation der Konstanten): Ansatz: t2 t2 t2 x = C(t)e− 2 ⇒ x0 = C 0 (t)e− 2 − tC(t)e− 2 Einsetzen: t2 t2 t2 t2 C 0 (t)e− 2 − tC(t)e− 2 + tC(t)e− 2 = e− 2 Z ⇒ C(t) = 1dt = t + C t2 t2 t2 ⇒ x(t) = (t + C) e− 2 = te− 2 + Ce− 2 – Partikuläre Lösung: 02 02 t2 x(0) = 1 ⇒ 1 = 0e− 2 + Ce− 2 ⇒ C = 1 ⇒ x = (t + 1) e− 2 • zu 1c: – Homogene DGL: Z 0 x − sin (t) x = 0 ⇒ 1 dx = x Z sin(t)dt ⇒ ln |x| = − cos(t) + C ⇒ xhom = Ce− cos(t) – Inhomogene DGL (Variation der Konstanten): Ansatz: x = C(t)e− cos(t) ⇒ x0 = C 0 (t)e− cos(t) + sin(t)C(t)e− cos(t) Einsetzen: C 0 (t)e− cos(t) + sin(t)C(t)e− cos(t) − sin(t)C(t)e− cos(t) = cos(t) Z ⇒ C(t) = cos(t)ecos(t) dt =??? Es gibt keine exakte Lösung! • zu 1d: Seite 3 / 8 MDS Lösung Serie 1 (Differentialgleichungen 1-ter Ordnung) FS 2009 – Homogene DGL: Z Z 1 0 dx = −1dt ⇒ ln |x| = −t + C ⇒ xhom = Ce−t x +x=0⇒ x – Inhomogene DGL (Variation der Konstanten): Ansatz: x = C(t)e−t ⇒ x0 = C 0 (t)e−t − C(t)e−t Einsetzen: C 0 (t)e−t − C(t)e−t + C(t)e−t = e−t Z ⇒ C(t) = 1dt = t + C ⇒ x(t) = (t + C) e−t = te−t + Ce−t – Partikuläre Lösung: x(0) = 0 ⇒ 0 = 0e−0 + Ce−0 ⇒ C = 0 ⇒ x = te−t • zu 1e: – Homogene DGL: Z 0 x + tan(t)x = 0 ⇒ 1 dx = x Z − tan(t)dt ⇒ ln |x| = ln(cos(t)) + C ⇒ xhom = C cos(t) – Inhomogene DGL (Variation der Konstanten): Ansatz: x = C(t) cos(t) ⇒ x0 = C 0 (t) cos(t) − C(t) sin(t) Einsetzen: C 0 (t) cos(t) − C(t) sin(t) + tan(t)C(t) cos(t) = sin(t) Z ⇒ C(t) = tan(t)dt = − ln(cos(t)) + C ⇒ x(t) = (− ln(cos(t)) + C) cos(t) – Partikuläre Lösung: x(0) = 0 ⇒ 0 = (− ln(cos(0)) + C) cos(0) ⇒ C = 0 ⇒ x = − cos(t) ln(cos(t)) • zu 1g: – Homogene DGL: 0 x − 3x = 0 ⇒ Z 1 dx = x Z 3dt ⇒ ln |x| = 3t + C ⇒ xhom = Ce3t Seite 4 / 8 MDS Lösung Serie 1 (Differentialgleichungen 1-ter Ordnung) FS 2009 – Inhomogene DGL (Variation der Konstanten): Ansatz: x = C(t)e3t ⇒ x0 = C 0 (t)e3t + 3C(t)e3t Einsetzen: C 0 (t)e3t + 3C(t)e3t − 3 C(t)e3t = t2 Z 1 ⇒ C(t) = t2 e−3t dt = − 2 + 6t + 9t2 e−3t + C 27 −3t 1 2 ⇒ x(t) = − 2 + 6t + 9t e + C e3t 27 – Partikuläre Lösung: 2 2 x(0) = − ⇒ − = 27 27 −3∗0 1 2 + C e3∗0 − 2+6∗0+9∗0 e 27 ⇒C=0⇒x=− 1 2 + 6t + 9t2 27 3. Aufgabe Erstelle zu den Differentialgleichungen 1b, 1d und 1g ein ausagekräftiges Richtungsfeld und skizziere damit die gesuchten Lösungskurven der Anfangswertprobleme. Lösung: • zu 1b: Mit MATLAB: >> richtungsfeld(t/x,[t,x],[-6,6,-6,6],[-2,-1,-0.5,0,0.5,1,2],20) >> hold on >> ezplot(sqrt(9+t^2),[-6,6]) Seite 5 / 8 MDS Lösung Serie 1 (Differentialgleichungen 1-ter Ordnung) FS 2009 • zu 1d: Mit MATLAB: >> richtungsfeld(exp(-t)-x,[t,x],[-6,6,-6,6],[-5,-2,-1,-0.5,0,0.5,1,2,5],20) >> hold on >> ezplot(t*exp(-t),[-6,6]) • zu 1g: Mit MATLAB: >> richtungsfeld(t^2+3*x,[t,x],[-2,2,-2,2],[-5,-2,-1,-0.5,0,0.5,1,2,5],20) >> hold on >> ezplot(-(2+6*t+9*t^2)/27,[-2,2]) Seite 6 / 8 MDS Lösung Serie 1 (Differentialgleichungen 1-ter Ordnung) FS 2009 4. Aufgabe Modelliere die DGL der 1. Aufgabe in SIMULINK und vergleiche die Lösungen mit den berechneten exakten Lösungen. Lösung: Modelle: MATLAB-Befehle für die Graphen: >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> tt=simout.time; x1=simout.signals.values(:,1); x2=simout.signals.values(:,2); x3=simout.signals.values(:,3); x4=simout.signals.values(:,4); x5=simout.signals.values(:,5); x6=simout.signals.values(:,6); x7=simout.signals.values(:,7); x8=simout.signals.values(:,8); subplot(2,4,1) plot(tt,x1,’r*’) hold on ezplot((t+1)*exp(-t^2/2),[0,2]) subplot(2,4,2) plot(tt,x2,’r*’) hold on ezplot(sqrt(9+t^2),[0,2]) Seite 7 / 8 MDS Lösung Serie 1 (Differentialgleichungen 1-ter Ordnung) >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> >> subplot(2,4,3) plot(tt,x3,’r*’) subplot(2,4,4) plot(tt,x4,’r*’) hold on ezplot(t*exp(-t),[0,2]) subplot(2,4,5) plot(tt,x5,’r*’) hold on ezplot(-cos(t)*log(cos(t)),[0,2]) subplot(2,4,6) plot(tt,x6,’r*’) hold on ezplot(sqrt(1+t^3),[0,2]) subplot(2,4,7) plot(tt,x7,’r*’) hold on ezplot(-1/27*(2+6*t+9*t^2),[0,2]) subplot(2,4,8) plot(tt,x8,’r*’) hold on ezplot(tan(2/3*sqrt(t^3)),[0,2]) Graphen: Seite 8 / 8 FS 2009
© Copyright 2024 ExpyDoc