Fs = 1000；％采样频率首先使用上面的代码生成一个1Hz，1.5s的采样信号。然后，我们添加振幅为0.7，频率为120 S的正常正弦波信号= 0.7 * sin（2 * pi * 50 * t）+ sin（2 * pi * 120 * t）；向信号添加零平均白噪声，方差为4X = S + 2 * randn（size（t））;在时域图中绘制噪声信号（1000 * t（1:50），X（1:50））标题（＆＃39;零均值随机噪声破坏的信号）xlabel（＆amp; ＃39; t（毫秒）＆y;＃39; ylabel（＆＃39; X（t）＆＃39;）计算信号的FFTY = fft（X）;计算两侧频谱P2，再计算一侧频谱P1，平均信号长度LP2 = abs（Y / L）； P1 = P2（1：L / 2 + 1）; P1（2：end-1）= 2 * P1（2：end-1）;定义时域f，并画出P1f = Fs *（0：（L / 2））/ L;图（f，P1） X（t）的TItle（＆＃39;单面幅度谱）xlabel（＆amp;＃39; f（Hz）＆amp;＃39;）ylabel（＆amp;＃39; | P1（ f）|＆＃39;）原始信号的傅立叶变换Y = fft（S）; P2 = abs（Y / L）; P1 = P2（1：L / 2 + 1）; P1（2：end- 1）= 2 * P1（2：end-1）;图（f，P1）TItle（＆amp;＃39; S（t）＆amp;＃39;）xlabel（＆＃39; f（Hz）＆＃39;）ylabel（＆amp;＃39; | P1（f）|＆amp;＃39;）