|
|
|
DGL Ergebnis in Funktion einbinden |
|
| Jojo92 |
Forum-Newbie
|
 |
Beiträge: 3
|
|
|
|
Anmeldedatum: 17.05.12
|
|
|
|
Wohnort: Magdeburg
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 17.05.2012, 18:01
Titel: DGL Ergebnis in Funktion einbinden
|
 |
Hallo,
ich habe folgendes Problem mit Matlab.
Für eine kleine Hausarbeit, soll ich eine Fahrzeugsimulation in Matlab durchführen. Die DGL's habe ich auch schon aufgestellt und gelöst. Mein Problem ist jetzt, dass ich, um die Steigung für die Hangabtriebskraft zu berechen, den Winkel "alpha" brauche. Diesen will ich mit Hilfe des arctan von meiner Steigungsfunktion berechnen, um jedoch die aktuelle Steigung zu berechnen, muss ich in diese Funktion den aktuellen Wert des Weges, den ich aus der DGL erhalte, einbinden. Nun zu meiner Frage: "Wie geht das?"
Die DGL's habe ich mit ode45 gelöst und dadurch erhalte ich ja einen Ergebnisvektor, wenn ich diesen aber in die Steigungsfunktion einbinden will, erhalte ich nur eine Error.
Ich hoffe es kann mir jemand helfen. Wenn ihr die m-files braucht, sagt bescheid, kann die auch noch hochladen.
PS: Ich arbeite erst knapp 5 Wochen mit Matlab.
|
|
|
|
|
|
| Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
|
Verfasst am: 17.05.2012, 18:18
Titel:
|
|
Hallo,
zumindest solltest du sagen, welchen Fehler du erhältst und bei welchem Befehl. Beispielcode ist immer hilfreich, am besten ein möglichst überschaubares Beispiel.
Grüße,
Harald
|
|
|
|
| Jojo92 |
Themenstarter
Forum-Newbie
|
|
Beiträge: 3
|
|
|
|
Anmeldedatum: 17.05.12
|
|
|
|
Wohnort: Magdeburg
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 17.05.2012, 18:48
Titel:
|
|
| |
 |
 |
Hallo Harald,
der Fehler lautet:
Undefined function 's' for input arguments of type 'double'.
Error in ..... (line 20)
alpha(i)=atan(Steigung(s(i)));
Das .... steht für den Ordnernamen.
Dieser Fehler taucht auf, wenn ich das Skript ausführen will. Ich schreib es einfach mal hier rein, weil mir kein Beispiel im Moment einfällt.
Hauptskript:
close all
clear all
% Konstanten
FA=3000;
m=1500;
d=20;
cw=0.3;
rho=1.293;
vw=4.167;
g=9.81;
A=2;
t=[0:0.1:600];
n=length(t);
for i=0:n %Berechnung von alpha in
alpha(i)=atan(Steigung(s(i))); %Abhängigkeit des zurück-
if s(i)<=3000 %gelegten Weges.
alpha(i)=0;
else
alpha(i)=atan(Steigung(s(i)));
end
[t,s]=DGL2_Loesung(FA,m,d,rho,cw,A,g,vw,0,0,i);
[t,v]=DGL1_Loesung(FA,m,rho,cw,A,vw,g,d,0,i);
end
figure;
hold on;
subplot(2,1,1) % 1.Subplot (Weg-Zeit-Diagramm)
plot(t,s,'k');
xlabel('t in s');
ylabel('s in m');
subplot(2,1,2) % 2.Subplot (Geschwindigkeit-Zeit-Diagramm)
plot(t,v,'k');
xlabel('t in s');
ylabel('v in m/s');
1.DGL
% Funktion zur Lösung der DGL 1.Ordnung, um die Geschwindigkeit % zu berechnen.
function [t,vLoesung]=DGL1_Loesung(FA,m,rho,cw,A,vw,g,d,alpha,v0,tend)
[t,vLoesung]=ode45(@DGL1,[0 tend],v0);
function dvdt=DGL1(t,v)
dvdt=FA/m-d/m*v-(rho*cw*A*(v-vw)^2)/(2*m)-g*sin(alpha);
end
end
2.DGL
% Funktion zur Lösung der DGL 2.Ordnung, um den Weg zu berechnen. Zerlegt in DGL-System 1.Ordnung.
function [t,sLoesung]=DGL2_Loesung(FA,m,d,rho,cw,A,g,vw,alpha,s10,s20,tend)
[t,sLoesung]=ode45(@DGL2,[0,tend],[s10,s20]);
function dsdt=DGL2(t,s)
dsdt=zeros(2,1); % Erstellt einen Spaltenvektor.
dsdt(1)=s(2);
dsdt(2)=FA/m-d/m*s(2)-(rho*cw*A*(s(2)-vw)^2)/(2*m)-g*sin(alpha);
end
end
Alpha-Berechnung:
function Winkel=alpha(s)
Winkel=atan(Steigung)
function mg=Steigung
mg=1/10000*(s-3000);
if s<=3000
mg=0;
else
mg=1/10000*(s-3000);
end
end
end
Ich hoffe, dass das hilft.
MfG
|
|
|
|
| Jan S |
Moderator
|
|
Beiträge: 11.057
|
|
|
|
Anmeldedatum: 08.07.10
|
|
|
|
Wohnort: Heidelberg
|
|
|
|
Version: 2009a, 2016b
|
|
|
|
|
|
|
Verfasst am: 17.05.2012, 23:43
Titel:
|
|
Hallo Jojo92,
Bitte benutze die Code-Umgebung für Programm-Code.
Nach dem "clear all" is die Variable "s" nicht mehr definiert. Was ist "s"?
Nebenbei möchte ich erwähnen, das "clear all" zwar infaltionär eingesetzt wird, es ist aber trotzdem nicht nützlich. Neben dem Löschen aller Variablen (und das ist schon von zweifelhaftem Nutzen), werden auch alle in den Speicher geladenen Funktionen gelöscht, so dass ein zeitaufwändiges Nachladen von der Festplatte erforderlich wird. Dies ist dann sinnvoll, wenn man alle bereits gelesenen M-files während des Programmablaufs verändert hat - was aber wiederum eine wirklich schlecht Idee wäre...
Gruß, Jan
|
|
|
|
| Jojo92 |
Themenstarter
Forum-Newbie
|
|
Beiträge: 3
|
|
|
|
Anmeldedatum: 17.05.12
|
|
|
|
Wohnort: Magdeburg
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 18.05.2012, 14:23
Titel:
|
|
Hallo Jan,
danke für den Hinweis. "s" ist der Weg, der mir durch die 2.DGL berechnet wird. Nur leider gibt mir die DGL den Weg ja in der Form [t,s] aus und ich weis jetzt nicht, wie ich ihn in der Funktion zur Berechnung von alpha aufrufen soll.
|
|
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
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.
|
|
