double_half_sine

function y = double_half_sine(tarr, td0, td1, td2, amp1, amp2)
    % double half sine

    y = zeros(size(tarr));

    % 각 구간에 맞는 인덱스를 계산
    idx0 = tarr < td0; % 1/3 시작점
    idx1 = (tarr >= td0) & (tarr < td1); % 2/3 중간점 증분
    idx2 = (tarr >= td1) & (tarr < td2); % 3/3 종료점 증분

    % 각각의 구간에 대한 y 값 계산
    y(idx1) = amp1 * sin(pi / (td1 - td0) * (tarr(idx1) - td0)); % 첫 번째 구간
    y(idx2) = amp2 * sin(pi / (td2 - td1) * (tarr(idx2) - td1)); % 두 번째 구간

end