clear all;
clc;
hold off;

data=load('W_500_3.txt');
n_an=500;           %Drehzahl 1/min

x_null=-21.4057;         %x-Wert relativ zu Koordinatenurpsrung(Drehachse von Lamelle)) in mm
y_null=43.5653;          %y-Wert relativ zu Koordinatenurpsrung(Drehachse von Lamelle)) in mm


x_PIV=data(:,1);    %x-Wert im Davis Koordinatensystem in mm
y_PIV=data(:,2);    %y-Wert im Davis Koordinatensystem in mm
x_korr=x_PIV+x_null;    %x-Wert im Koordinatensystem mit Ursprung Drehachse Lamelle in mm
y_korr=y_PIV+y_null;    %x-Wert im Koordinatensystem mit Ursprung Drehachse Lamelle in mm

%Berechnung Tangentialgeschwindigkeit durch Lamellendrehung
Z=zeros(1,length(x_korr))';
R=[x_korr y_korr Z];
o=ones(size(x_korr));
omega_z=-2*pi*n_an/60;
omega=[Z Z o*omega_z];    %1/s
v_tang=cross(omega,R); %Tangentialgeschwidnigkeit durch Lamellendrehung in mm/s
v_tang_x=v_tang(:,1)/1000;
v_tang_y=v_tang(:,2)/1000;

for i=1:1:length(x_korr)
    if data(i,3)==0
        v_x(i,1)=0;
        v_y(i,1)=0;
    else
        v_x(i,1)=data(i,3)-v_tang_x(i);
        v_y(i,1)=data(i,4)-v_tang_y(i);
    end;
end;
        

PIV_korr=[x_korr y_korr v_x v_y]; %Matrix im neuen "Drehachsenkoordinatensystem%

quiver(PIV_korr(:,1),PIV_korr(:,2),PIV_korr(:,3),PIV_korr(:,4),'-.or');