image thumbnail
from RJ-MCMC algorithm for sinusoids parameter estimation by Sebastien PARIS
RJ-MCMC algorithm for sinusoids estimation

generate_y_sinu(N , param_sinu , sig)
function y  = generate_y_sinu(N , param_sinu , sig)

% y  = generate_y_sinu(N , param_sinu , sig);


if( nargin <3 | isempty(sig) )
    
    sig     = 1;
    
end

t           = (1 : N)';

[d , k]     = size(param_sinu);


if( (d ~= 2) & (k ~= 0) )
    
    error('param_sinu must be (2 x k)');
    
end

y           = zeros(N , 1);

if (k > 0)
    
    tp_wk       = param_sinu(2 , :);
    
    y           = cos(tp_wk(ones(N , 1) , :).*t(: , ones(1 , k)))*(sqrt(2)*sig*10.^(param_sinu(1 , :)/20))';
    
end

Contact us at files@mathworks.com