source: orbit/octave/GenerateAudioCommand/finalCodeAudioGen.m @ 6aa182f

function finalCodeAudioGen(code,audioFileName,sps,repeatTime) 

% code: array of 1s and 0s. 1 means bit "1", 0 means bit "0".
% audioFileName: like 'audio.wav'.

bps = 16;       % bits per sample
%sps = 96000;     % sample rate [samples/s]

codeWave=[halfSineGen('u',800,0,1,sps);halfSineGen('d',400,0,1,sps);halfSineGen('u',800,0,1,sps);halfSineGen('d',400,0,1,sps)];

for i=1:size(code,2)
        codeWave = [codeWave; bitGen(code(i),1,sps)];
endfor

fullCodeWave=[];
for i=1:repeatTime
        fullCodeWave=[fullCodeWave; codeWave; zeros(sps/20,1)];
endfor

%special wave form for carrier
carrierWave=[halfPeriodGen('u',30*1.215,-0.353,0.723,sps);
halfPeriodGen('d',22*1.215,-0.015,0.384,sps);
halfPeriodGen('u',31*1.215,0.261,0.661,sps);
halfPeriodGen('d',28*1.215,0.392,0.530,sps);
halfPeriodGen('u',25*1.215,0.246,0.384,sps);
halfPeriodGen('d',32*1.215,-0.092,0.723,sps);
halfPeriodGen('u',25*1.215,-0.430,0.384,sps);
halfPeriodGen('d',26*1.215,-0.515,0.469,sps)];

fullCodeWave=[fullCodeWave fullCodeWave].*0.5;

sizeCW=size(carrierWave,1);
cWaveRepeat=[];
n=floor(size(fullCodeWave,1)/sizeCW);
for i=1:n
        cWaveRepeat=[cWaveRepeat; carrierWave];
endfor

fullCodeWave(1:n*sizeCW,2)+=cWaveRepeat;
wavwrite(fullCodeWave, sps, bps, audioFileName);