%%% 時間軸:
% 8k[sampling/sec], 0.5秒の時間軸
time = linspace(0, 0.5, 8000);

%%% 基本波形データ作成
wave1 = sin(time*2*pi*800)'; %800Hzの振幅1の波形
wave2 = sin(time*2*pi*1000)'; %1000Hzの振幅1の波形


%%% 各種操作
% 交互につなぐ
tmp1 = [wave1; wave2; wave1; wave2];

% 2つの波形を合成 (振幅が2になるので元の振幅に調整)
tmp2 = (wave1+wave2)/2;

% tmp1とtmp2をつなぐ
wave_all = [tmp1; tmp2];


%%% ファイルに出力
wavwrite(wave_all, 8000, 16, 'test.wav')

%%% Windowsのシステムを使ってでファイルを鳴らす
system('test.wav')