WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Numerik-Differentialgleichungen

 

maschinenbaustudent

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 05.01.2014, 14:23     Titel: Numerik-Differentialgleichungen
  Antworten mit Zitat      
Hallo,
ich habe ein riesengroßes Problem. Und zwar habe ich von meinem Mathematik Professor 3 Aufgaben bekommen, welche ich mit Hilfe von Matlab lösen soll. Folgendes Problem... Ich besitze zwar matlab, habe allerdings noch nie damit gearbeitet.
Die Aufgaben lauten:

Aufgabe 1: Bakterienkultur
Bakterien in einer Lösung vermehren sich mit eine Reproduktionsrate r = 0.01 pro Stunde.
In der geschlossenen Kultur sei nur Platz für maximal N = 100000 Bakterien, die
Vermehrung geschieht entsprechend der Differenzialgleichung
y' = r(1 -(y/N))*Y

Ist y0 = 1000 die Anzahl der Bakterien zum Zeitpunkt 0, wieviele Bakterien gibt es in
der Lösung nach 30 Tagen? Wieviele gäbe es ohne die Beschränkung nach oben?

Aufgabe 2: PISA Test
Ein Baseball falle am Turm von Pisa aus einer Höhe von 50 m in genau 4 sec auf den
Boden. Wie gross ist sein Luftwiderstandswert r, wenn die folgende Differenzialgleichung
den freien Fall beschreibt (y(t) die Höhe über dem Boden):
y'' = -g + r *(y')², y(0) = 50
Bestimme r auf etwa drei Nachkommastellen.

Aufgabe 3: Erzwungene Schwingung
Ein Pendel ist frei beweglich an einem horizontalen Stab aufgehängt. Das Pendel ist zum
Zeitpunkt t0 = 0 in Ruhe, wenn eine oszillierende Bewegung y(t)=Y *sin(w* t) auf die
Aufhängung einwirkt. Der Winkel z der Auslenkung folgt dann einer Differenzialgleichung
z'' = -(g/L)*sin(z)+(w²/L)*Y*cos(z)*sin(wt)

mit g = 9.81 m/s², L = 1.0 m, Y = 0.25 m und w = 2.5 rad/s.
Löse die Differenzialgleichung und stelle die Auslenkung des Pendels über einen Zeitraum
von 60 Sekunden graphisch dar, Graphik mit Beschriftung, Grid und Legende.


Da ich leider überhaupt keine Ahnung hab und nichtmal weiß wie ich anfangen soll, wäre ich für jede Hilfe, eventuell auch Musterlösung dankbar.

Aufgaben2013 Numerik.pdf
 Beschreibung:

Download
 Dateiname:  Aufgaben2013 Numerik.pdf
 Dateigröße:  140.71 KB
 Heruntergeladen:  515 mal


Harald
Forum-Meister

Forum-Meister


Beiträge: 24.501
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 05.01.2014, 15:11     Titel:
  Antworten mit Zitat      
Hallo,

mit ode45 kannst du Differentialgleichungen numerisch lösen.

Schau dir doch mal die Hilfe dieser Funktion an
Code:


Differentialgleichungen 2. Ordnung musst du zunächst in Systeme 1. Ordnung umwandeln.

Für 2. könntest du beispielsweise fzero verwenden, um die Nullstelle der Funktion "Höhe nach 4s" zu bestimmen.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
maschinenbaustudent

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 05.01.2014, 15:30     Titel:
  Antworten mit Zitat      
Habe mit Hilfe eines anderen Studenten probiert Aufgabe 3 zu lösen..
Könnte der Ansatz so stimmen?
Ist einfach mal ein Versuch...

Code:
function dy = pendel1(t,y)

g = 9.81;   L = 1;  Y = 0.25;   w = 2.5;

y1 = y(1);
y2 = y(2);

dy1 = y1;
dy2 = -(g / L) * sin(z1) + ((w^2) / L) * Y * cos(z1) * sin(w *t);
dz = [dz1;dz2];
end


Schon mal vielen Dank Smile
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.501
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 05.01.2014, 15:43     Titel:
  Antworten mit Zitat      
Hallo,

du solltest dich auf y oder z einigen.

Es sollte dir zu denken geben, dass y2 auf der rechten Seite nicht vorkommt - da ist also noch ein Fehler drin.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
.





 Impressum  | Nutzungsbedingungen  | Datenschutz | FAQ | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.