function [out]=sigexpand(d,M)
% 將輸入的序列擴展成間隔為N-1個0的序列
%------------------------輸入參數
%M:在這里指系統采樣頻率，即在一個碼元寬度內采多少個樣點，設計時要考慮載波頻率與基帶信號采樣率。M為二者的乘積。
%d:輸入的二進制代碼
%----------------------輸出(返回)參數
%out：為1*(M*length(d))維的矩陣
%例如d=[1 1],M=3,則out=[1 0 0 1 0 0]
N=length(d);
out=zeros(M,N);
out(1,:)=d;
out=reshape(out,1,M*N);