# how to contaminate ecg signal with noise

1 view (last 30 days)

Show older comments

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

##### 11 Comments

Daniel M
on 20 Nov 2019

### Answers (2)

Daniel M
on 20 Nov 2019

Edited: Daniel M
on 21 Nov 2019

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 Comments

Daniel M
on 23 Nov 2019

Did you not see my response that phi must be << than 1? Did you even try it? It works.

### See Also

### Categories

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!