xnoisesig.m

来自「可综合FFT的VHDL实现 已经过测试」· M 代码 · 共 28 行

% generate test data for fixed-point DA FFT
n=16;					% n samples in test vector
ar=0.98;				% peak amplitude of test data - real
ai=0.9;				% peak amplitude of test data - imag
state = 1.0e+009 * [0.4851
                    3.7261];
    
                 
randn('state',state);
                 
% first generate a gaussian sequence

rr=randn(1,n);
ri=randn(1,n);

mx = max(abs(rr));
rr=ar/mx*rr;

mx = max(abs(ri));
ri=ai/mx*ri;

fp=fopen('xnoise.dat','w');
for k=1:n
   fprintf(fp,'%f\n',rr(k));
   fprintf(fp,'%f\n',ri(k));
end
fclose(fp);