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

f = [5 10];          % frekvenser
sigma2_ = [2 2];     % energi för dessa frekvenserna
N = 500;             % Antal simuleringar
dt = 1/(2*max(f))/3; % sampel avstånd
t = 0:dt:20;         % tid
fprintf('Frekvenser:\n')
disp(num2str(f))
fprintf('Energier:\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)
[rayamp]=enkelsumma(f,sigma2_,N,t);
%
% För att använda spekgui
clear data
data.x = rayamp(:,5);
data.dt = dt;