Differenzieren

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

mike793

Forum-Century


	Beiträge: 119

	Anmeldedatum: 11.08.08

	Wohnort: ---

	Version: ---

Verfasst am: 26.08.2008, 18:04 Titel: Differenzieren

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);

Funktion ohne Link?

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

Danke

nschlange

Ehrenmitglied


	Beiträge: 1.318

	Anmeldedatum: 06.09.07

	Wohnort: NRW

	Version: R2007b

Verfasst am: 26.08.2008, 21:06 Titel:

Hi,

das geht mit

Code:

diff

Funktion ohne Link?

, 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);

Funktion ohne Link?

_________________

Viele Grüße
nschlange

"Chuck Norris ejakuliert fluessigen Stahl!"

mike793

Themenstarter

Forum-Century


	Beiträge: 119

	Anmeldedatum: 11.08.08

	Wohnort: ---

	Version: ---

Verfasst am: 27.08.2008, 03:56 Titel: zum Differenzieren

Ich habe diesen Befehl bei meinem Problem angewendet, jedoch bekomme ich jedesmal eine Fehlermeldung. beim plot befehl!

Code:

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

Funktion ohne Link?

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

steve

Ehrenmitglied


	Beiträge: 2.022

	Anmeldedatum: 03.09.07

	Wohnort: Wien

	Version: R2023b

Verfasst am: 27.08.2008, 07:32 Titel:

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?

mike793

Themenstarter

Forum-Century


	Beiträge: 119

	Anmeldedatum: 11.08.08

	Wohnort: ---

	Version: ---

Verfasst am: 27.08.2008, 08:15 Titel:

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)

Funktion ohne Link?

Code:

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

Funktion ohne Link?

Ich ihr könnt damit etwas anfangen

steve

Ehrenmitglied


	Beiträge: 2.022

	Anmeldedatum: 03.09.07

	Wohnort: Wien

	Version: R2023b

Verfasst am: 27.08.2008, 08:32 Titel:

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

Funktion ohne Link?

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?

mike793

Themenstarter

Forum-Century


	Beiträge: 119

	Anmeldedatum: 11.08.08

	Wohnort: ---

	Version: ---

Verfasst am: 27.08.2008, 08:47 Titel: Hallo

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);

Funktion ohne Link?

das ' nach dem y?

und

Code:

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

Funktion ohne Link?

den 1 vor dem end-1?
Danke nochmals

steve

Ehrenmitglied


	Beiträge: 2.022

	Anmeldedatum: 03.09.07

	Wohnort: Wien

	Version: R2023b

Verfasst am: 27.08.2008, 09:15 Titel:

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?


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 |

RSS

Hosted by:

Copyright © 2007 - 2024 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.