怎么在这段MATLAB程序导入的心电信号上添加一段50HZ的噪声？

t = (0:0.001:1)

drift= sin(2*pi*50*t)

这样是以采样率为1000Hz采样得到的正弦信号，叠加在心电信号上即可，如果需要的话进行转置。

在matlab中无论是wgn还是awgn函数，实质都是由randn函数产生的噪声。即，wgn函数中调用了randn函数，而awgn函数中调用了wgn函数。

根据awgn的实现代码可以知道“向已知信号添加某个信噪比（SNR）的高斯白噪声”,

即：awgn(x,snr,’measured’,'linear’)，命令的作用是对原信号x添加信噪比（比值）为SNR的噪声，在添加之前先估计信号x的强度。

直接对原始信号添加噪声：

y=x+rand(length(x),1)

y=x+randn(length(x),1)）

1、首先启动MATLAB软件。

2、首先设定好波形的基本参数，采样点数，采样频率，采样间隔，时间间隔，最高采样频率等，注意要符合采样定理才能保证信号不失真。

3、然后定义声音信号的波形，可以是普通的余弦波，也可以复杂的波形，然后用subplot语句将波形绘制出来。

4、然后创建噪声的波形，这里是调用了取值在0~1的随机函数，用来产生高斯噪声，再绘制出它的波形。

5、然后根据傅里叶变换的原则，将噪声的时域表达转换为频域表达，再绘制出它的波形。

6、通过将原始信号和噪声信号相加，即可得到混合噪声后的信号，再调用绘图语句画出它的波形。

---