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

function codeAudioGen(code,audioFileName) 

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

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

codeWave=[];

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

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

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

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