Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   
Bücher:

Fachkräfte:
Promotion Roboterführung auf Basis von Texturmerkmalen (m/w)
Entwicklung von Verfahren zur kamerabasierten Lokalisation
Robert Bosch GmbH - Hildesheim

System-Entwickler (m/w)
Analyse und Spezifikation von Systemanforderungen & Abstimmung von Anforderungen
MicroNova AG - Friedrichshafen

Elektromobilität: E-Drive Entwicklungsingenieur HV-Energiemanagement und Betriebsstrategie (m/w)
Entwicklung von Teilfunktionen für Energiebordnetze
MBtech Group GmbH & Co. KGaA - Sindelfingen

Systemingenieur (m/w) Funktionsentwicklung Automotive
Konzeption und Spezifikation von spezifischen Funktionen in elektronischen Steuergeräten
DRÄXLMAIER Group - Vilsbiburg bei Landshut

Ingenieur (w/m) modellbasierte Entwicklung / Embedded Programmierung
Organisation und Durchführung von Projekten zur Integration von TargetLink in aktuelle Entwicklungsprozesse
dSPACE GmbH - Pfaffenhofen

weitere Angebote

Partner:


Vermarktungspartner


Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Probleme bei den constraints von fmincon

 

Neuling89
Forum-Newbie

Forum-Newbie


Beiträge: 1
Anmeldedatum: 15.12.17
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 15.12.2017, 11:53     Titel: Probleme bei den constraints von fmincon
  Antworten mit Zitat      
Hey!

Ich bin gerade an der Optimalsteuerung einer Pendelbewegung (Doppelpendel) dran und habe Probleme bei den Bedingungen, die ich stellen will.

Das Ziel ist es, den Winkelbereich einzuschränken, sodass das 1. Pendel sich nur zwischen -pi/4 und pi/2 und das 2. Pendel nur zwischen -pi/2 und pi/2 schwingen kann.

Ich habe verstanden, dass c(X) < 0 erfüllt sein muss.

X besteht einerseits aus den Lagevariablen q und andererseits aus den Aktuierungsmomenten tau

Code:

function [c,ceq,gc,gceq] = constraints(x,J,parameter)  % mit J = 40

Nq = parameter.Nq;  % Dimension q = 2 (aus structure parameter)
Ntau = parameter.Ntau;   % Dimension tau = 2

% Ueberpruefe die Dimension von x
if length(x)~=((Nq+Ntau)*(J+1))    %x muss die gleiche Dimension haben wie [q0,q1,...,qN,tau0,tau1,...,tauN]
    error('beschr: x hat falsche Dimension');
end
% Extrahiere Lagekoordinaten
Q = x(1:Nq*(J+1));      %alle [q0,q1,...,qN]

% Extrahiere Aktuierungsmomente
Tau = x(Nq*(J+1)+1:end);    %alle [tau0,tau1,...,tauN]

q0 = Q(1:Nq);
%ceq = .....
 


ceq(X) = 0 funktioniert, das habe ich daher nicht aufgeführt.

Ich habe nun erst einmal versucht, q nach unten zu begrenzen

Code:

%nichtlineare Randbedingung
c = zeros((J+1+2)*Nq,1);              % Anzahl Gleichungsbeschraenkungen
gc = zeros((J+1+2)*Nq,(J+1)*(Nq+Ntau));    % Zugehoeriger Gradient
gctau = zeros((J+1+2)*Nq,(J+1)*Ntau);       % Zugehoeriger Gradient von tau, spielt aber eigentlich keine Rolle hier(?), da ja nur Lagekoordinaten begrenzt werden

phi1min = -pi/4;
phi2min = -pi/2;

q0 = Q(1:Nq);
c(1:Nq) = q0 + [phi1min;phi2min];  %soll ja <0 sein...
gc(1:Nq,1:Nq) =  eye(Nq);               %Gradient


q1 = Q(Nq+1:2*Nq);
c(Nq+1:2*Nq) = q1 + [phi1min;phi2min];
gc(Nq+1:2*Nq,Nq+1:2*Nq) =  eye(Nq);

j = 1;
while j<J
q = Q;
c((j+1)*Nq+1:(j+2)*Nq) = q(j:j+1) + [phi1min;phi2min];
gc((j+1)*Nq+1:(j+2)*Nq,(j+1)*Nq+1:(j+2)*Nq) =  eye(Nq);
j = j+1;
end

gc(:,(J+1)*Nq+1:end) = gctau;
gc = gc.';


end
 


Matlab führt alles ohne Fehlermeldung durch, wenn ich das aber anschließend plotte, bewegt sich das Pendel leider gar nicht..

Bin für jede Hilfe dankbar!
Liebe Grüße
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
.


goMatlab ist ein Teil des goForen-Labels
goForen.de goMATLAB.de goLaTeX.de goPCB.de


 Impressum  | Nutzungsbedingungen  | Datenschutz  | Werbung/Mediadaten | Studentenversion | FAQ | goMatlab RSS Button RSS


Copyright © 2007 - 2018 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
Partner: LabVIEWforum.de

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.