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

s=load('08215_2.txt');
I=s(1001:2000);
drift = 100*sin(linspace(0,2*pi,1000));
I=I'+drift;
原程序使在导入的信号上添加了一段正弦信号的漂移
如果要改成添加50HZ的噪声应该怎么加？

追问：我的采样率是250Hz 应该怎么转啊？转完以后是什么样呢？s=load('08215_2.txt');
I=s(1:1000);
t=(0:0.004:1);
drift = 50*sign(linspace(0,2*pi,1000));
I=I'+drift;

C=length(I);
%I=load('ECG250Hz.dat');
%C=length(I);
STEP=100;
COUNT=floor(C/STEP);
n=0;

这样改完了以后运行出来的结果是
??? Error using ==> unknown
Matrix dimensions must agree.

Error in ==> f_detection1 at 16
I=I'+drift;

这个是因为数组大小不一致导致的么？不好意思刚才贴错了。。。
s=load('08215_2.txt');
I=s(1:1000);
t=(0:0.004:1);
drift= sin(2*pi*50*t);
I=I'+drift;

C=length(I);
%I=load('ECG250Hz.dat');
%C=length(I);
STEP=100;
COUNT=floor(C/STEP);
n=0;
运行出来也是那个结果 这个怎么改啊应该？运行出来还是一样的问题T T还是不行T_____T
高人 我能加你QQ问你么？

答案：t = (0:0.001:1);
drift= sin(2*pi*50*t);
这样是以采样率为1000Hz采样得到的正弦信号，叠加在心电信号上即可，如果需要的话进行转置。

其他：，怎么实现二分频呢？例如 sin（100t） 把它分频sin（50t）各位高手怎么实现呢？怎么设置参数？ t=[0:0.1:n] , simulation step=0.1 然后设

上一个：matlab 用后diff 数据长度变小 l
下一个：matlab gui 清除图形