当前位置:编程学习 > Matlab >>

怎么在这段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 清除图形

CopyRight © 2012 站长网 编程知识问答 www.zzzyk.com All Rights Reserved
部份技术文章来自网络,