噪声谱特征提取

%% 舰船辐射噪声仿真
% 连续谱 最高频率200Hz 总20480Hz 
% 线谱 基频10Hz，4根 10 20 30 40 
% N为时域采样点数（信号长度），NFFT为做FFT的点数，N=NFFT=fs ND为做CPSD的点数
% 变量说明：fs采样频率；f1低频线谱基频；T为1s；t为时间轴；ff:0:1:N-1；
%          sig_linea低频线谱；Sig_linea低频线谱的频谱；
%          noise生成的N点高斯白噪声；noise1经过带通滤波器（[-4096:4096]）的高斯白噪声；NOISE为noise1的频域；
%          Sig_cont连续谱的频谱；sig_cont连续谱的时域；signal低频线谱与连续谱的时域叠加；Signal为signal的频域；
%          mod包络函数（时域）；sig_cont_mod为在时域调制后的连续谱 signal_final最终的时域谱
% HILBERT变换
clc;close all;

%% 线谱
% 线谱时域2048个点，做2048点FFT，每点间隔频率1Hz
fs=10240;
f1=10; %线谱基频10Hz
n=4;
T=1;N=fs*T;NFFT = fs;
t = 0:1/fs:T-1/fs;
% sig_linea = 0;
% Fai_n = 2*pi*rand（14）;
% An=10e11*rand（14）; % 随机幅值和相位 
% for i=1:1:n
%     sig_li=An（i）*cos（i*（2*pi*f1*t）+Fai_n（i））;
%     sig_linea = sig_linea + sig_li;
% end
sig_linea = 0.3*10e12