Verfasst am: 05.02.2013, 16:18
Titel: Zeitvektoren überschneiden sich - wie ist das zu ermitteln?
Moin,
ich hab mal wieder ein Brett vorm Kopf. Folgendes Problem
Es existieren 3 Zeitvektoren, alle 3 sind unterschiedlich lang. Teilweise sind die Werte in diesen Vektoren auch nicht sortiert.
Die meisten Einträge in diesen Zeitvektoren sind nur einmal vorhanden, andere Einträge sind in zwei der drei Vektoren vorhanden und ganz wenige Einträge sogar in allen drei.
Nun gilt es herauszufinden an welcher Stelle in diesen Vektoren Einträge existieren, die 3mal, 2mal oder nur 1mal vorhanden sind. Ansich sollte das machbar sein, aber ich sitze nur wie doof vorm PC.
Was ich schon habe sind die Werte, die 3mal vorhanden sind:
Code:
%Generiere Zelle mit allen 3 Zeitvektoren
cell_zeit{1}=zeit1;
cell_zeit{2}=zeit2;
cell_zeit{3}=zeit3;
dreimal = cell_zeit{1};
for i=2:length(cell_zeit)
dreimal = intersect(cell_zeit{i}, dreimal);
end;
%Finde die dazugehörigen Einträge
for i=1:1:length(isc)
ind1(i)=find(zeit1==dreimal(i));
ind2(i)=find(zeit2==dreimal(i));
ind3(i)=find(zeit3==dreimal(i));
end
die Einträge pro Vektor sind aber eindeutig?
intersect liefert auch gleich den Index zurück, sodass die find-Schleife
unnötig ist (was ist im Übrigen isc?).
Wenn es egal ist, wie oft ein Wert mehrfach vorkommt, sonder nur wichtig
ist, ob oder nicht, kannst Du auch sowas versuchen:
das ist ein guter Ansatz, danke! Mit biegen und brechen habe ich es auch geproggt bekommen, dass mir bekannt ist, ob ein Wert mehrfach vorkommt. Die Frage ist halt auch wie häufig und vorallem in welchem Vektor.
Nach einiger Zeit bin ich dann auf die Mengenlehre gekommen.
3 Kreise (=Mengen) die sich überschneiden repräsentieren die 3 Vektoren. Mit einfachen Verknüpfungen ist es mir dann gelungen die Einträge zu finden, die mir wichtig sind. Etwa so: (A,B,C sind hier die Zeitvektoren)
Code:
%------------------------------------------------------------------------
%1. A\(B U C); B\(A U C); C\(A U B)
%Generiert drei Ergebnisse, wo es keine Überschneidungen gibt, ergo keine
%Addition von Zielgröße Y
zeit_alleine1=setdiff(zeit1,union(zeit2,zeit3));
zeit_alleine2=setdiff(zeit2,union(zeit1,zeit3));
zeit_alleine3=setdiff(zeit3,union(zeit1,zeit2));
%Finde die dazugehörigen Einträge
for i=1:1:length(zeit_alleine1)
I_zeit_alleine1(i)=find(zeit1==zeit_alleine1(i));
end for i=1:1:length(zeit_alleine2)
I_zeit_alleine2(i)=find(zeit2==zeit_alleine2(i));
end for i=1:1:length(zeit_alleine3)
I_zeit_alleine3(i)=find(zeit3==zeit_alleine3(i));
end
%Addition vom Zielwert::
Y=[Y1(I_zeit_alleine1,:); Y2(I_zeit_alleine2,:); Y3(I_zeit_alleine3,:)];
zeit=[zeit_alleine1 zeit_alleine2 zeit_alleine3];
%-----------------------------------------------------------------------
%2. A n B n C
%Generiert ein Ergebnisse, wo sich alle 3 Überschneiden
%Generiere Zelle mit allen 3 Zeitvektoren
cell_zeit{1}=zeit1;
cell_zeit{2}=zeit2;
cell_zeit{3}=zeit3;
zeit_all = cell_zeit{1};
for i=2:length(cell_zeit)
zeit_all = intersect(cell_zeit{i}, zeit_all);
end;
%Finde die dazugehörigen Einträge
for i=1:1:length(zeit_all)
zeit_all1_I(i)=find(zeit1==zeit_all(i));
zeit_all2_I(i)=find(zeit2==zeit_all(i));
zeit_all3_I(i)=find(zeit3==zeit_all(i));
end
%Addition vom Zielwert::
Y_alle3=Y1(zeit_all1_I,:)+Y2(zeit_all2_I,:)+Y3(zeit_all3_I,:);
zeit_all;
Y=[Y; Y_alle3];
zeit=[zeit zeit_all];
%-------------------------------------------------------------------------
%3. A n B \ (A n B n C) A n C \ (A n B n C) C n B \ (A n B n C) % zeit_all = (A n B n C)
%Generiert drei Ergebnisse, jeweils wo es eine Schnittmenge von 2
%Zeitvektoren gibt
zeit_dualA=setdiff(intersect(zeit1,zeit2),zeit_all); %Vec 1 + 2
zeit_dualB=setdiff(intersect(zeit2,zeit3),zeit_all); %Vec 2 + 3
zeit_dualC=setdiff(intersect(zeit3,zeit1),zeit_all); %Vec 3 + 1
%Finde die dazugehörigen Einträge Wert Y(ACHTUNG! Hier müssen
%nun die 2 Vektorenmengen betrachtet werden, mit den passenden Indicies)
%Für Zeit_dualA:
for i=1:1:length(zeit_dualA)
I_zeit_dual1(i)=find(zeit1==zeit_dualA(i));
I_zeit_dual2(i)=find(zeit2==zeit_dualA(i));
end
Y_dualA=Y1(I_zeit_dual1,:)+Y2(I_zeit_dual2,:);
clear I_zeit_dual1; clear I_zeit_dual2;
%Für Zeit_dualB:
for i=1:1:length(zeit_dualB)
I_zeit_dual2(i)=find(zeit2==zeit_dualB(i));
I_zeit_dual3(i)=find(zeit3==zeit_dualB(i));
end
Y_dualB=Y2(I_zeit_dual2,:)+Y3(I_zeit_dual3,:);
clear I_zeit_dual3; clear I_zeit_dual2;
%Für Zeit_dualC:
for i=1:1:length(zeit_dualC)
I_zeit_dual3(i)=find(zeit3==zeit_dualC(i));
I_zeit_dual1(i)=find(zeit1==zeit_dualC(i));
end
Y_dualC=Y3(I_zeit_dual3,:)+Y1(I_zeit_dual1,:);
clear I_zeit_dual3; clear I_zeit_dual1;
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
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.
RSS