function [sigr erg]=fitsinus_pony1(sig);
% Das geht jetzt auch mit sehr stark ueberabgetasteten Signalen
% Vorsicht: Das Signal muss vorher DC frei sein

sig=sig(:).';

%load
%n=2345678;
%w=8*2*pi/n
%sig = 1.34*sin(8.123*2*pi*(0:n-1)/n);
%davor =davor-mean(davor);
spd=fft(sig);
[y,ind]=max(abs(spd(1:end/2)));
ind=ind-1;% Wieviele Wellen ?
nn=length(sig)/ind; % Abtastwerte/Welle
dd=round(nn/4); % Das ist der Abstand für 4 Werte/Welle

spd=1; % loeschen

n=length(sig);
s1=sig(1:(n-2*dd))+sig((1+2*dd):n);
s2=sig((1+dd):(n-dd));

%whos

% Regression s1 auf s2
MM= [s2.' ones(length(s1),1)];
coff=inv(MM'*MM)*MM'*s1.';
wn=acos(coff(1)/2)
wn=wn/dd;

%plot(s2,s1,'b.',s2,(MM*coff).','r.')
%plot(s2,s1,'b.',s2,s1,'r.')
M=[cos((0:n-1)*wn).' sin((0:n-1)*wn).'];
coff=inv(M'*M)*M'*sig.';

cc=coff(2)+j*coff(1);
amp=abs(cc);
ph=angle(cc);

sig=5;

sigr=amp*sin(wn*(0:n-1)+ph);

erg=[amp wn ph];
%plot(1:n,sig,1:n,sigr);
return

