function?[TA_out]=add_noise（T_inA_inmeansigma）
global?n
	%ampitude=0.04;
????%A_in=sqrt（Power）.*exp（i*Phase）;
????num_point=length（T_in）;
????dt=T_in（2）-T_in（1）;
????dw=2*pi/dt/num_point;
????W=（-（num_point/2）:（num_point/2-1））*dw;
????Fft_A_in=fftshift（fft（A_in））;
????%Fft_A_out=Fft_A_in+normrnd（meansqrt（sigma）1num_point）;
????nn=1:num_point;
????a=-0.05/（（num_point）^2）;
????b=num_point/2;
????c=sigma;
????enve=（a*（nn-b）.^2+c）/10;
%?????figure
%?????plot（nnenve）
????rand（‘state‘sum（100*clock））;
????Fft_noise=enve.*（normrnd（meansqrt（sigma）1num_point））;
????Fft_A_out=Fft_A_in+Fft_noise;
????
????T=T_in;
????
????figure
%?????subplot（211）
%?????plot（Wabs（Fft_A_out）.^2*dt‘r‘）;
%?????hold?on
%?????plot（Wabs（Fft_A_in）.^2*dt‘b‘）;
????
????A_out=ifft（ifftshift（Fft_A_out））;
???
%?????subplot（212）
????plot（T_inabs（A_out）‘r‘‘LineWidth‘3）
????hold?on
????plot（T_inabs（A_in）‘-b‘‘LineWidth‘3）
????axis（[max（T_in）-n*64*dt?max（T_in）?-0.02?0.15]）