Fehler: Vektor verkleinert sich

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Long

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 21.06.2009, 15:02 Titel: Fehler: Vektor verkleinert sich

Hallo liebe Leute,

ich habe einen Fehler, den ich irgendwie nicht erklären kann.
Der folgende Code verursacht den Fehler

??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==> gallo_bb at 104
if x*P*x' > lb

Das liegt daran, dass x irgendwann auf einmal nur noch 99 Komponenten statt 100 hat. Aber dieses verstehe ich nicht, denn ich verkleinere doch nirgendwo im Code den Vektor x.

Woran könnte es liegen?

Code:

function [x,ub] = gallo_bb(c,w,P)

n = length(w);

U = []; % fixed auf 1
Z = []; % fixed auf 0
F = 1:n; % free
M = []; % organisiert Backtracking
test0 = zeros(1,n);
test1 = zeros(1,n);

c_start = c;

% 1. obere Schranke
[ub,x,v] = gallo_ub(c,w,P);
% 1. untere Schranke
lb = dualgreedy(c,w,P);

while (1~=0)
ub
if ub > lb % noch nicht ausgelotet
% zu verzweigende Variable r bestimmen
F1 = find(x==1);
eta = zeros(1,n);
for j=F1
eta(j) = v(j)/w(j);
end

[maxiumum,r] = max(eta);

else
if test1(r)==1 %x_r = 1 auch untersucht: backtracking
r = M(length(M));
% i = 0;
while (mark(r)==1)
test0(r) = 0;
test1(r) = 0;
M = M(M~=r);
F = [F,r]; % Variablen wieder freimachen
%r
% i = i+1;
r = M(length(M));
if (r == M(1)) && (mark(r)==1)
return
end
end
end
end
r

if test0(r)==1 %Verzweigen an x_r = 1
x(r) = 1;
mark(r) = 1;
test1(r) = 1;
U = [U,r];
F = F(F~=r);

% neues c
for i=U
c = c_start-w(i);
end

% neue v
v=zeros(1,n);
for j=F
spalte_j = P(:,j);
[weightsum,v(j),x] = ckp(c,w(F),spalte_j(F)');
v(j) = v(j)+2*sum(spalte_j(U))+2*P(r,j);
end

% neue ub mit Zusatz x_r=1
[weightsum,ub,neue_x] = ckp(c,w(F),v(F));
x(F) = neue_x;
% x
% if x*P*x' > lb
% lb = x*P*x';
% end
disp('branching 1')

else % Verzweigen an x_r = 0
x(r) = 0;
M = [M,r];
Z = [Z,r];
F = F(F~=r);
test0(r) = 1;

% neues c
for i=U
c = c_start-w(i);
end

% neue v
v=zeros(1,n);
for j=F
spalte_j = P(:,j);
[weightsum,v(j)] = ckp(c,w(F),spalte_j(F)');
v(j) = v(j)+2*sum(spalte_j(U));
end

% neue ub mit Zusatz x_r=0
[weightsum,ub,neue_x] = ckp(c,w(F),v(F));
x(F) = neue_x;
x
if x*P*x' > lb
lb = x*P*x';
end
disp('branching 0')
% ub = x*v' (alle x nicht aus neue_x sind 0)

end

end

Funktion ohne Link?

Long

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 21.06.2009, 15:31 Titel:

Fehler gefunden.


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.