%%?艦船輻射噪聲仿真
%?連續譜?最高頻率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