function [out]=sigexpand(d,M)
%將輸入序列擴展為間隔為N－1個0的序列
% d為輸入信號 M為采樣點數 
% d=[1 1 1 0 1 1 0 1];
% M=20;
N=length(d);
out=zeros(M,N);
out(1,:)=d;
out=reshape(out,1,M*N);