% Simulerar en enkel normalprocess på formen
%
% X(t) = A0 + sum Ak cos(2 pi fk t + phik)
%
% så som finns beskrivet i kompendiet kap 4.3.2

f = [16];       % frekvenser
sigma2_ = [6];   % energi för dessa frekvenserna
N = 100;             % Antal simuleringar
% sampla processen för sällan
dt = 1/(1.5*16) ;        % sampel avstånd
t = 0:dt:20;           % tid
dtok = 1/(10*16);
tok = 0:dtok:20;
fprintf('Frekvenserna är\n')
disp(num2str(f))
fprintf('Energierna är\n')
disp(num2str(sigma2_))
fprintf('Maximala frekvensen i samplad signal %0.5g Hz\n',1/(2*dt))
fprintf('Samplingsavstånd %0.5g sek\n',dt)

[rayamp2]=enkelsumma(f,sigma2_,N,t,0);
clear data2
data2.x = rayamp2(:,5);
data2.dt = dt;
[rayamp2]=enkelsumma(f,sigma2_,N,tok);