%
% tplot: 波形データの時間軸プロット (11/12/22 改訂： YH)
%       tplot(y, Fs, t1, t2)
%           y : 音響データベクトル
%           Fs: サンプリングレート（デフォールト 44100）
%           t1: 開始時刻（デフォールト： 0）
%           t2: 終了時刻（デフォールト： y の時間長さ）
%               t2 省略時には、t1（第３パラメタ）が終了時刻となる。
%

function [y tt] = tplot(y, Fs, t1, t2, ax)

if (nargin < 4)
    t2 = 0;
end;
if (nargin == 3)
    t2 = t1; t1 = 0;
elseif (nargin < 3)
    t1 = 0;
end;
if (nargin < 2)
    Fs = 44100;
end;
ylen = length(y);
tlen = length(y)/Fs;
if (t2 == 0)
    t2 = tlen;
end;
if (nargin < 5)
    ax = [];
end;

i1 = round(t1 * Fs)+1;
i2 = round(t2 * Fs);
if (i2 > ylen); i2 = ylen; end;

t = (0:(length(y)-1))/Fs;
ymax = max(abs(y(i1:i2)));
if (isempty(ax))
    ax = [t1 t2 -ymax*1.2 ymax*1.2];
end;

clf;
hold on;
axis(ax);
plot(t(i1:i2), y(i1:i2), 'b');
plot([t1 t2], [0 0], 'k');
hold off;
shg;
tt = t(i1:i2);