how to contaminate ecg signal with noise
2 次查看(过去 30 天)
显示 更早的评论
I am working on the project namely "FALSE ALARM REDUCTION IN ATRIAL FIBRILLATION DETCTION USING DEEP BELIEF NETWORKS".
In this i need to download two databases, namely 07910m.mat{MIT_BIH ATRIAL FIBRILLATION DATABASE} which is attached here,and also the 'em' recording {MIT_BIH NOISE STRESS TEST DATABASE} is attached.
The afib database has Fs=250Hz and 10h in length,and 'em' recording is 30 min and has Fs=360Hz
now i need to make the 'em' recording 10h in length
in the paper its told that i can be done using auto regressive model,iam not knowing how to do it
can anyone suggest me the coding for this![screenshot.png](https://www.mathworks.com/matlabcentral/answers/uploaded_files/248616/screenshot.png)
![screenshot.png](https://www.mathworks.com/matlabcentral/answers/uploaded_files/248616/screenshot.png)
11 个评论
Daniel M
2019-11-20
Well you know there will be a loop involved. Can you show me the code you currently have? I can help you get there.
回答(2 个)
Daniel M
2019-11-20
编辑:Daniel M
2019-11-21
I obviously cannot run this code because I don't have the data nor any of the associated functions. By the way, you are getting an error because you are trying to access z at negative index locations.
Here is an example of how you would plot the equation. Keep in mind the assumptions I've made, which I've tried to comment.
clearvars
clc
sz = [1 1000];
e = randn(sz); % epsilon
X = nan(sz); % preallocate
X(1) = 0; % set first value. Don't know what this should be.
c = 1; % a constant, don't know what this should be
p = 5; % set the order of the model. Don't know what this should be.
% note that p cannot be greater than the iterating variable of the loop
% otherwise, it would try to access negative index locations
% So, either initialize your loop index at a value greater than 2 (and have the starting values of X)
% or do it this way, using increasing values of p, up to n-1
phi = rand(1,p); % your values of phi
for n = 2:sz(2)
ii = 1:min(p,n-1); % see the note above
sigma = dot(phi(ii), X(n-ii)); % sum of the products of the elements
X(n) = c + sigma + e(n);
end
7 个评论
Daniel M
2019-11-23
Did you not see my response that phi must be << than 1? Did you even try it? It works.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Measurements and Feature Extraction 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!