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

Numerische Integration

 

whitecrane19
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 27
Anmeldedatum: 25.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 28.10.2013, 10:14     Titel:
  Antworten mit Zitat      
Ok aber wozu die numerische Ableitung. Meine beiden Formeln fx und gx muss ich ja integrieren. Dann nochmal ableitungen?
Private Nachricht senden Benutzer-Profile anzeigen


Harald
Forum-Meister

Forum-Meister


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

in deiner Formel von 25.10.2013, 18:50 treten Ableitungen auf. Wenn man also Werte für f und g gegeben hat, muss man erst diese Ableitungen bestimmen.

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

Forum-Anfänger

Forum-Anfänger


Beiträge: 27
Anmeldedatum: 25.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 28.10.2013, 15:31     Titel:
  Antworten mit Zitat      
Ok, die Ableitungen kann ich auf alle Fälle per Hand machen. Ist am Anfangf übersichtlicher.

Deine Formeln müssten ja die Länge L berechnen. Aber die Vorgehensweise ist mir nicht ganz klar. Für mich sieht es so aus, als hast du der einfachheit halber alles getrennt eingegeben, damit es vor allem übersichtlicher erscheint oder? ich verstehe eine Teil davon nämlich nicht:

df = (fx(3:end)-fx(1:end-2))./(x(3:end) - x(1:end-2));
dg = (gx(3:end)-gx(1:end-2))./(x(3:end) - x(1:end-2));

Kann ich da als Grenzen z.B. das Intervall, was für fx und gx vorher vorgegeben war, einsetzen?
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister


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

die Grenzen werden an sich automatisch berücksichtigt.

Die von dir zitierten Formeln berechnen die zentralen Differenzenquotienten für die inneren Datenpunkte, also
f'(x) ~ (f(xr) - f(xl))/(xr-xl)
wobei xr und xl die rechten bzw. linken Nachbarn der Datenpunkte sind.
Danach hatte ich dann noch einseitige Differenzenquotienten für die äußeren Datenpunkte angefügt.

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

Forum-Anfänger

Forum-Anfänger


Beiträge: 27
Anmeldedatum: 25.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 28.10.2013, 15:56     Titel:
  Antworten mit Zitat      
Hmm, dann ist das sozusagen die Rechenformel für die Gleichung, wobei ich lediglich die Funktionen fx und gx bzw. deren Ableitung angeben muss und die Grenzen z.B. von 0 bis 2*pi?
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister


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

die Grenzen für x gibst du ganz am Anfang im linspace-Befehl an.

Dann die Funktionsauswertung und die der Ableitungen. Wenn du die Ableitung analytisch bestimmen kannst, brauchst du sie natürlich nicht numerisch anzunähern.

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

Forum-Anfänger

Forum-Anfänger


Beiträge: 27
Anmeldedatum: 25.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 28.10.2013, 16:34     Titel:
  Antworten mit Zitat      
In meinem Fall dürfte sinx und cosx mehr als einfachabzuleiten sein, steht ja auch in fast jeden Tabelle ....
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister


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

in diesem Fall: ja.
Es kann ja aber durchaus sein, dass du mal nicht so triviale Funktionen hast, wo der Integrand immer 1 ist.

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

Forum-Anfänger

Forum-Anfänger


Beiträge: 27
Anmeldedatum: 25.10.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 28.10.2013, 23:07     Titel:
  Antworten mit Zitat      
Hab noch eine Fehler:

Error using diff
Difference order N must be a positive integer scalar.

Error in Aufg1 (line 57)
df = diff(fx, x);


Hier das komplette Beispiel:

Code:
% Beispiel 1
t=linspace(0, 2*pi, 100);
f = inline ('cos(t)')
g = inline ('sin(t)')

% Funktionsplot
plot (t,f(t), t,g(t))
grid on
axis equal


%Integration%
df = diff(f, t);
dg = diff(g, t);

x = linspace(0, 2*pi, 1000);
f = @sin;
g = @cos;
fx = f(x);
gx = g(x);
% Ableitungen über zentrale bzw. am Rand einseitige Differenzen
df = (fx(3:end)-fx(1:end-2))./(x(3:end) - x(1:end-2));
dg = (gx(3:end)-gx(1:end-2))./(x(3:end) - x(1:end-2));
df = [(fx(2)-fx(1))/(x(2)-x(1)), df, (fx(end)-fx(end-1))/(x(end)-x(end-1))];
dg = [(gx(2)-gx(1))/(x(2)-x(1)), dg, (gx(end)-gx(end-1))/(x(end)-x(end-1))];
integrand = sqrt(df.^2 + dg.^2);
L=trapz(x, integrand)


Außerdem möchte ich fragen, ob

f = @sin;
g = @cos;

meine bereits abgeleiteten Funktionen sind?
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen

Gehe zu Seite Zurück  1, 2

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.