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

Differenzieren
 

mike793
Forum-Century

Forum-Century

Beiträge: 119
Anmeldedatum: 11.08.08
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 26.08.2008, 17:04     Titel: Differenzieren
  Antworten mit Zitat      
Hallo

Ich möchte gerne eine funktion differenzieren, der dazu benötigte befehl sollte die gleichen Eigenschaften haben wie der für das Integrieren (cumtrapz).
Code:

z.B.:

 Z = cumtrapz(y);
 

Gibt es diesen Befehl und wenn ja wie weißt er?

Danke
Private Nachricht senden Benutzer-Profile anzeigen E-Mail senden


nschlange
Ehrenmitglied

Ehrenmitglied


Beiträge: 1.320
Anmeldedatum: 06.09.07
Wohnort: NRW
Version: R2007b
     Beitrag Verfasst am: 26.08.2008, 20:06     Titel:
  Antworten mit Zitat      
Hi,

das geht mit
Code:
, allerdings ist der Ergebnisvektor um ein Element gegenüber der Eingabe verkürzt:
Code:
clear all;
clc;

x=-5:0.02:5;
y=-0.5*x.^3+3*x-5;
dy=diff(y)./diff(x);
plot(x,y);
hold on;
plot(x(1:end-1),dy,'r');
hl=legend('f(x)','f''(x)');
set(hl,'FontSize',20);

_________________

Viele Grüße
nschlange

"Chuck Norris ejakuliert fluessigen Stahl!"
Private Nachricht senden Benutzer-Profile anzeigen E-Mail senden
 
mike793
Themenstarter

 Forum-Century

Forum-Century

Beiträge: 119
Anmeldedatum: 11.08.08
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 27.08.2008, 02:56     Titel: zum Differenzieren
  Antworten mit Zitat      
Ich habe diesen Befehl bei meinem Problem angewendet, jedoch bekomme ich jedesmal eine Fehlermeldung. beim plot befehl!

Code:

??? Error using ==> ./
Matrix dimensions must agree.
 


Deshalb wäre es für mich wichtig, zu wissen, ob es eine andere Methode gibt eine Funktion, die mit dem Befehl ODE erstellt wurde dann über die Zeit differenziert anzuzeigen.

Ich habe eine 1x501 double array funktion zu ploten,
weis jemand wie man das macht, bzw welchen befehl man dabei verwenden muss?

Danke
Private Nachricht senden Benutzer-Profile anzeigen E-Mail senden
 
steve
Ehrenmitglied

Ehrenmitglied


Beiträge: 2.029
Anmeldedatum: 03.09.07
Wohnort: Wien
Version: R2024a
     Beitrag Verfasst am: 27.08.2008, 06:32     Titel:
  Antworten mit Zitat      
Moin,

dann zeig doch mal den Code mit Beispieldaten bei dem dieser Fehler reproduzierbar ist.

Gruß
Alex
_________________

>> I told me to.

____________________________________
Matlab Cheat Sheet
goMatlab-Knigge - dran gehalten?!
Schon in den FAQ gesucht?
Ist vielleicht bei den Skripten oder den Tutorials was für dich dabei?
Private Nachricht senden Benutzer-Profile anzeigen
 
mike793
Themenstarter

 Forum-Century

Forum-Century

Beiträge: 119
Anmeldedatum: 11.08.08
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 27.08.2008, 07:15     Titel:
  Antworten mit Zitat      
Also mein Code Lautet


Code:

t0=0;
te=3;
x0=[0;0];
u=1;
a0=1;
a1=15;
a2=10;
b0=1;
b1=23;
b2=22;
[t,x]=ode45('gg',t0,te,x0);
y=[(b0-a0*b2)(b1-a1*b2)]*x'+b2*u;
dya=diff(y)./diff(t);
plot(t(t0:end-1),dya)
 


Code:


function dx = gg(t,x)
dx=[x(2);-16*x(1)-10*x(2)+1]
 


Ich ihr könnt damit etwas anfangen
Private Nachricht senden Benutzer-Profile anzeigen E-Mail senden
 
steve
Ehrenmitglied

Ehrenmitglied


Beiträge: 2.029
Anmeldedatum: 03.09.07
Wohnort: Wien
Version: R2024a
     Beitrag Verfasst am: 27.08.2008, 07:32     Titel:
  Antworten mit Zitat      
Hmmm,

ich denke du solltest dir nochmal anschauen, wie die odesolver-Aufrufe erfolgen und wie man auf einen Bereich eines Vektors zugreift.

Dieser Code läuft bei mir jetzt zumindest fehlerfrei durch!
Code:

function odeaufruf

t0=0;
te=3;
x0=[0;0];
u=1;
a0=1;
a1=15;
a2=10;
b0=1;
b1=23;
b2=22;
[t,x]=ode45(@gg,[t0 te],x0);
y=[(b0-a0*b2) (b1-a1*b2)]*x'+b2*u;
dya=diff(y')./diff(t);
plot(t(1:end-1),dya)


function dx = gg(t,x)
    dx=[x(2);-16*x(1)-10*x(2)+1];
end

end


Gruß
Alex
_________________

>> I told me to.

____________________________________
Matlab Cheat Sheet
goMatlab-Knigge - dran gehalten?!
Schon in den FAQ gesucht?
Ist vielleicht bei den Skripten oder den Tutorials was für dich dabei?
Private Nachricht senden Benutzer-Profile anzeigen
 
mike793
Themenstarter

 Forum-Century

Forum-Century

Beiträge: 119
Anmeldedatum: 11.08.08
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 27.08.2008, 07:47     Titel: Hallo
  Antworten mit Zitat      
Ich habe zu deinem Code einige Fragen. Ich hoffe ich belässtige dich damit nicht. aber was bedeutet bei dieser Zeile

Code:

dya=diff(y')/diff(t);
 

das ' nach dem y?

und

Code:

plot(t(1:end-1),dya)
 


den 1 vor dem end-1?
Danke nochmals
Private Nachricht senden Benutzer-Profile anzeigen E-Mail senden
 
steve
Ehrenmitglied

Ehrenmitglied


Beiträge: 2.029
Anmeldedatum: 03.09.07
Wohnort: Wien
Version: R2024a
     Beitrag Verfasst am: 27.08.2008, 08:15     Titel:
  Antworten mit Zitat      
Das ' macht nichts weiter als den Vektor zu transponieren, also aus einem 1x5-Vektor macht es einen 5x1-Vektor.

Zu deiner zweiten Frage:
Genau das meinte ich mit: "Anschauen wie man auf Bereiche von Vektoren zugreift":
Dein Vektor "t" besteht aus z Elementen. Du möchtest jetzt auf alle Elemente zugreifen, nur das letzte Element soll nicht berücksichtigt werden. Also willst du die Elemente 1 bis end-1 haben.
Dein Code war an der Stelle insofern fehlerhaft, dass du mit y0=0 versucht hast auf das 0te Element des Vektors zuzugreifen und das ist nicht möglich, da es kein 0tes Element gibt Wink

Ich hoffe ich konnte helfen!

Gruß
Alex
_________________

>> I told me to.

____________________________________
Matlab Cheat Sheet
goMatlab-Knigge - dran gehalten?!
Schon in den FAQ gesucht?
Ist vielleicht bei den Skripten oder den Tutorials was für dich dabei?
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.