power_spec_signal.m

来自「高斯白噪声输入」· M 代码 · 共 19 行

function Power_Input_X=Power_Spec_Signal(Input_X,flag)

%--Input_X为输入的时域信号
%--Power_Input_X为输入信号的离散功率谱密度
%--时域离散，对应 频域连续
%--N=512表示进行512点的傅立叶变换
%flag
N=512;

Spect_Input_X = fft(Input_X,N);

Power_Input_X = Spect_Input_X.* conj(Spect_Input_X) / 512;

if lower(flag)=='half'
    
    Power_Input_X=Power_Input_X(1:N/2);
    
end;