Ableitung in Scilab

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

sirrbouk

Forum-Anfänger


	Beiträge: 15

	Anmeldedatum: 08.10.15

	Wohnort: ---

	Version: ---

Verfasst am: 28.10.2015, 07:07 Titel: Ableitung in Scilab

Hallo

ich mochte im Scilab erste und zweite Ableitung plotten, welche Befehle brauche ich???

ich habe das geschrieben aber die erste Ableitung sieht kommisch aus!!!!

Code:

t=0:0.005:0.05;
koeff=[0 0 0 400000 -12000000 96000000];
p=poly(koeff,"t","coeff");
horner(p,t);
plot(t,horner(p,t));xgrid;

p1=diff(p);
horner(p1,t);
plot(t,horner(p1,t),'r');xgrid

Funktion ohne Link?

[EDITED, Jan, Bitte code-Umgebung verwenden - Danke!]

Winkow

Moderator


	Beiträge: 3.842

	Anmeldedatum: 04.11.11

	Wohnort: Dresden

	Version: R2014a 2015a

Verfasst am: 28.10.2015, 08:02 Titel:

kenn mich mit scilab nicht aus aber bei matlab muss man auch noch diff(t) mit berücksichtigen.

Code:

t=linspace(-10,10,1000)
x=t.^2
dx=diff(x);
dt=diff(t);
plot(t,x)
hold on
plot(t(2:end),dx)
plot(t(2:end),dx./dt)

Funktion ohne Link?

_________________

richtig Fragen

sirrbouk

Themenstarter

Forum-Anfänger


	Beiträge: 15

	Anmeldedatum: 08.10.15

	Wohnort: ---

	Version: ---

Verfasst am: 29.10.2015, 04:45 Titel:

Hallo vielen Dank für die Antworte, leider klappt bei diff nicht.ich habe auch mit derivat probiert ich kann immer nur erste Gleichung (Weg) plotten:

Code:

t=0:0.005:0.05;

koeff=[0 0 0 400000 -12000000 96000000];
p=poly(koeff,"t","coeff");
horner(p,t);

plot(t,horner(p,t));xgrid;
label('Time[sec]'),ylabel('Pos[rad]');

p1=derivat(p);
plot(t,horner(p1,t));xgrid;
xlabel('Time[sec]'),ylabel('Speed[rad/s]');

p2=derivat(p1);
plot(t,horner(p2,t));xgrid;
xlabel('Time[sec]'),ylabel('Acc[1e3rad/s^2]');

Funktion ohne Link?

hat jemand vielleicht noch eine Idee.

Lg

[EDITED, Bitte Code-Umgebung verwenden - Danke!]

Jan S

Moderator


	Beiträge: 11.057

	Anmeldedatum: 08.07.10

	Wohnort: Heidelberg

	Version: 2009a, 2016b

Verfasst am: 30.10.2015, 09:00 Titel:

Hallo sirrbouk,

Bitte verwende die Code-Umgebung für Code.

Was genau bedeutet "klappt bei diff nicht"? Je genauer Du das erklärst, desto leichter ist es eine Antwort zu finden.

Gruß, Jan

sirrbouk

Themenstarter

Forum-Anfänger


	Beiträge: 15

	Anmeldedatum: 08.10.15

	Wohnort: ---

	Version: ---

Verfasst am: 30.10.2015, 09:13 Titel:

Hallo,

ich habe gemeint :

mit dem Befehl "diff" ich habe keine Ableitung gefunden um die Geschwindigkeit zu plotten.

LG
Sirrbouk

helmat

Forum-Fortgeschrittener


	Beiträge: 92

	Anmeldedatum: 18.12.11

	Wohnort: ---

	Version: R2020a

Verfasst am: 05.11.2015, 13:12 Titel:

Hallo sirrbouk,

Versuch es mal so:

Code:

xdel();
t=linspace(-0.01,0.07,100);
koeff=[0 0 0 400000 -12000000 96000000];
p=poly(koeff,"tau","coeff");
y=horner(p,t);
p1=derivat(p);
y1=horner(p1,t);
p2=derivat(p1);
y2=horner(p2,t);
p3=derivat(p2);
y3=horner(p3,t);
//Bitte Skalierungen beachten
plot(t,20*y,'r',t,y1,'b',t,0.01*y2,'g',t,0.0001*y3,'k');xgrid

Funktion ohne Link?

_________________

Herzliche Grüße
helmat

p.s.
Hard work beats talent if talent doesn't work hard.

sirrbouk

Themenstarter

Forum-Anfänger


	Beiträge: 15

	Anmeldedatum: 08.10.15

	Wohnort: ---

	Version: ---

Verfasst am: 08.11.2015, 05:00 Titel: Bild verbessern

Hallo

Erstmal besten Dank für den Vorschlag Smile

ich habe Folgendes geschrieben:

Code:

t=0:0.005:0.05;
koeff=[0 0 0 400000 -12000000 96000000];
p=poly(koeff,"tau","coeff");
y=horner(p,t);
p1=derivat(p);
y1=horner(p1,t);
p2=derivat(p1);
y2=horner(p2,t);
//p3=derivat(p2);
//y3=horner(p3,t);
//Position
subplot(1,3,1);
plot2d(t,y);xgrid
xlabel('Time [s]');
ylabel('Pos[rad]');
//Speed
subplot(1,3,2);
plot2d(t,1.1*y1);xgrid
xlabel('Time [s]');
ylabel('Speed[rad/s]');
//Acc
subplot(1,3,3);
plot2d(t,0.0012*y2);xgrid
xlabel('Time[s]');
ylabel('Acc[1e3 rad/s^2]');
Title('Bewegungsprofil mit Polynomen 5. Ordnung')

Funktion ohne Link?

was raus gekommen ,wie ich mir wünsche, leider bei Speed und Acc habe ich Spitzen und das sollt fein rund, (siehe Anhäng ) Wie kann ich das verbessern.

Vielen im voraus
Sirrbouk

Pos,Speed,Acc.pdf

Beschreibung:

Download

Dateiname:

Pos,Speed,Acc.pdf

Dateigröße:

35.98 KB

Heruntergeladen:

616 mal

helmat

Forum-Fortgeschrittener


	Beiträge: 92

	Anmeldedatum: 18.12.11

	Wohnort: ---

	Version: R2020a

Verfasst am: 08.11.2015, 08:43 Titel:

Einfach den t-Schritt verkleinern - und title statt Title.
_________________

Herzliche Grüße
helmat

p.s.
Hard work beats talent if talent doesn't work hard.


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.