function [y, t, v, P] =  CreateNPSK(NPSK, fc, amp, fs, RB, M)
% 產生PSK信號，其中NPSK為進制數，典型的為2、4、8等，fc為載波頻率，amp為載波基準幅度， fs為采樣頻率, RB為碼元速率
% 載波為 c = amp*sin(2*pi*fc*t);
% y返回已調信號，t為時間矢量，v為碼元矢量, P為采樣信號的平均功率
% [y, t, v] = Create2ASK(fc, fs, N), 其中M為信息碼元數，隨機產生N個碼元；
% [y, t, v] = Create2ASK(fc, fs, V), 其中V為給定的碼元序列，為行矢量；
if (length(M) == 1)
    N = M;
    V = randint(1, N, NPSK);
else
    V = M;
    N = length(V);
end;

v = V;

TotalTime = 1/RB * N;       %總時間

ts = 1 / fs;                %采樣周期
index = 1;
curTime = 0;
while(1)
    symbIndex = floor(curTime * RB) + 1;    %碼元序號
    if (symbIndex > N)
        break;
    end
    symb = V(symbIndex);                    %碼元符號
    y(index) = amp * cos(2*pi*fc*curTime + 2*pi/NPSK * symb + pi/NPSK);         %調制過程，產生已調信號瞬時幅度
         
    t(index) = curTime;                    
    curTime = curTime + ts;
    index = index + 1;
    if (curTime > TotalTime)
        break;
    end;
end;

P = mean(abs(y).^2);