geht soweit, erzeugt aber noch Murks

1 files changed, 18 insertions, 1 deletions

diff --git a/werteunit.pas b/werteunit.pas
index 745ed68..d26b65a 100644
--- a/werteunit.pas
+++ b/werteunit.pas
@@ -73,6 +73,7 @@ type
     procedure kantenFilter(betraege: tLLWerte; xFak,yFak: extended; filterTyp: tKantenFilterTyp; einseitig: boolean; out maxPos: tIntPoint); overload;
     procedure fenstereWerte(xMi,xMa,tMi,tMa: int64; xFen,tFen: tFenster; hg: extended);
     procedure verschiebe(richtung: tIntPoint; xV,xB,tV,tB: longint; komplex: boolean);
+    procedure ermittlePhasenWinkel(xMi,xMa: longint);
   end;
   tLLWerteSingle = specialize tLLWerte<single>;
   tLLWerteDouble = specialize tLLWerte<double>;
@@ -1961,7 +1962,6 @@ begin
     end;
   end;
   setLength(maxima,mCnt);
-  writeln(length(maxima),' (von ',betraege.params.xSteps*betraege.params.tSiz,')');
 
   betraege.sortiereMaxima(maxima);
   maxPos:=maxima[1];
@@ -2084,6 +2084,23 @@ begin
       end;
 end;
 
+procedure tLLWerte.ermittlePhasenWinkel(xMi,xMa: longint);
+var
+  i,j,ts2: longint;
+begin
+  ts2:=params.tSiz div 2;
+  for j:=0 to ts2-1 do
+    for i:=xMi to xMa do
+      if (werte[i+(j+ts2)*params.xSteps]=0) and (werte[i+j*params.xSteps]=0) then
+        werte[i+j*params.xSteps]:=0
+      else
+        werte[i+j*params.xSteps]:=
+          arctan2(
+            werte[i+(j+ts2)*params.xSteps],
+            werte[i+j*params.xSteps]
+          );
+end;
+
 // tWavelet ********************************************************************
 
 function tWavelet.setzeTyp(s: string): boolean;