tsym = 0.16
fs=16000
fmid=1500
fsym = 1/tsym
nstart=1
erg = 0
% Phasensprung: 0 bis 0.5, letzteres entspricht pi
%sprung = 0.5
sprung = 0.0
for i = 0:7
  f = fmid+i*fsym
  perioden = floor (f*tsym)
  t=0:1/fs:(perioden+sprung)/f;
  y = sin(2*pi*f*t);
  erg = [erg, y];
end

N=length(erg);
nn=1:N;
T=N/fs;
t=0:1/fs:T-1/fs;
ff=(mod(nn-1+N/2,N)-N/2)/T;
U=fft(erg)*2/N;
data1 = abs(U);
plot(t, erg);
semilogy(ff(1:3850)*1e-3,data1(1:3850));
axis tight
ylim([0.0001 0.2]);

audiowrite('fsktest.wav', erg, fs);