singen6_6.m

用dsp解压mp3程序的算法

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% singen6_6 - This routine generates the sinewave corrupted by white
%             noise for experiments in Section 6.6
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Fs = 10000;     % sampling rate
Ts = 1/Fs;      % total samples
F = 2000;       % sinewave frequency
sinewave = sin(2*pi*F*(0:Ts:1));              % generate sinewave
noise = sqrt(0.01).*randn(1,length(0:Ts:1));  % var = 0.01
in = sinewave+noise;             % combine sinewave with noise
save in.mat in                   % save in in.mat
in = in';
in = in(1:300);                  % using first 300 samples only
save in.dat in -ascii ;          % save in file in.dat
in_int = round(32768*in./max(abs(in))); % Convert to 16-bit integer
fid = fopen('in_int.dat','w');   % save signal to file in_int.dat
fprintf(fid,'%4.0f\n',in_int);   % in integer format
fclose(fid);