clear;clc;close all
% 初始化
m = [-6,-2,0,2,4,6];
sigma = [1,1,0.5,0.25,0.6,2];
h = [1,2,3,2,2.13,3.14];
x=-10:0.01:10;
sumy = zeros(1,length(x));
% 绘制各分支高斯信号
figure;
for i = 1:length(m)
    y=(h(i)/((sqrt(2*pi))*sigma(i)))*exp(-((x-m(i)).^2)/(2*sigma(i).^2));
    y=awgn(y,30,\'measured\');
    sumy = sumy+y;
    hold on;
    plot(x,y,\'LineWidth\',1);
end
title(\'Gaussian Signal add AWGN\')
xlabel(\'x\');ylabel(\'f(x)\')
% 绘制叠加高斯信号
figure;
plot(x,sumy,\'b\',\'LineWidth\',1)
title(\'Gaussian signal add AWGN\')
xlabel(\'x\');ylabel(\'f(x)\')
% mspeaks寻峰，返回的是峰值点坐标
p=mspeaks(x,sumy,\'SHOWPLOT\',true);