% sin_wav.m: 子程序，產(chǎn)生正弦波。
% 參數(shù)：周期T1，諧波次數(shù)Kh，初相位phs，采樣點(diǎn)數(shù),N。

function x=sin_wav(T1,Kh,phs,N)
pi=3.1415926535897932;
h=T1/N;
w=Kh*h*2*pi/T1;
fai=phs-w;
for i=1:N
    fai=fai+w;
    x(i)=sin(fai);
end

% 能量歸一化
A=sqrt(2/T1);
x=A*x;