0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014 function [fa,u,F]=LSM_generate_random(Nf,f_min,f_max,par)
0015
0016 bandwidth = par(1);
0017 mean_a = par(2);
0018 sigma_a = par(3);
0019 K = par(4);
0020 dx = par(5);
0021
0022 range_f = (-Nf:Nf)*dx;
0023 [X,Y] = meshgrid(range_f,range_f);
0024 xv = X(:);
0025 yv = Y(:);
0026
0027
0028 [u,F] = generate_F_Gauss_2D(mean_a,sigma_a,bandwidth,K);
0029
0030
0031 fa = signal_from_F_2D(u,F,xv,yv,[[0,0];[0,1];[1,1]]);
0032 fa_max = max(fa(:,1));
0033 fa_min = min(fa(:,1));
0034 fa(:,1) = floor((f_max-f_min) * (fa(:,1)-fa_min)/(fa_max-fa_min) + f_min);
0035 fa(:,2) = fa(:,2)*(f_max-f_min)/(fa_max-fa_min);
0036 return