Hi! I have this assignment due in Matlab, and I just have no idea how to implement the model in matlab ... My main issue at this point is that after I have defined a 50X50 grid to hold the state of each region (where the state can be 0,1,2 or 3 and depends on things like time elapsed in the current state, current state of the neighboring regions, etc), how do I define a separate 50X50 grid to store the time elapsed for the current state in each region in the original grid, where in this grid each region should be reset to 0 every time the corresponding region in the original grid changes state. In case that didn't make much sense, I'll post the actual assignment wording below - At this point I'm having problems with the part about defining a grid called S_duration ... if anyone has any tips or suggestions that would be greatly appreciated!(adsbygoogle = window.adsbygoogle || []).push({});

Assignment:

"Consider a 50×50 square grid, each square representing a small 2×2 mm2 region of

electrical activity in the ventricles of the heart. The grid is wrapped around into a

cylinder, so that the left and right edges are assumed to be in contact. Every region has

exactly eight neighbours, with the exception of those squares on the top and bottom edges

alone.

Each square in the grid can be in one of four states:

0 - quiescent the region is excited on the next time step if at least one if

its eight neighbours is currently excited.

3 - excited the region is excited, and can excite neighbouring regions

on the next time step.

2 - absolute refractory the region cannot be excited, nor can it excite any of its

neighbours.

1 - relative refractory the region is excited on the next time step if more than one

if its eight neighbours is currently excited. The number of

excited neighbours required is dependent on the time

elapsed since entering this state:

Time elapsed (in ms)/Number of neighbours required

≤ 2/8

≤ 4/7

≤ 6/6

≤ 8/5

≤ 12/4

≤ 20/3

≤ 50/2

Once a region is excited, it will move through states 3 → 2 → 1 → 0 in that order, unless

it is prematurely excited again.

The duration of the excited state is fixed at Tex and the duration of the relative refractory

period is fixed at 50 ms. To account for the inhomogeneous refractory properties of the

tissue, the duration of the total refractory period (absolute + relative) is a random variable

for each region, drawn from a normal distribution of mean Trp and standard deviation σrp.

The ventricles are excited by a periodic impulse from the AV node, of period T, located

in the small square on the top border of the grid one quarter of the way in from the left

edge. Whenever applied, this impulse can only excite the small square if its current state

is not absolute refractory.

All model parameters are given below:

Parameter Description/Value/Units

Tex Duration of excited state/70/ms

Trp Mean total refractory period/250/ms

σrp Standard deviation of total refractory period/100/ms

T AV node pacing period/0.8/s

dT Time step increment/2/ms

Your tasks for this assignment are:

1. Implement this model in Matlab, solving the heart's electrical activity for a

duration of 1.7 s.

2. Decrease the pacing period to 0.2 seconds, and solve the model again. What do

you observe?

HINTS AND SUGGESTIONS:

a) Define a 50×50 array H_state to store the current state of each region in the

grid. Update the values of the array during each time increment.

b) Similarly, define a 50×50 array S_duration to store the time elapsed for the

current state in each region in the grid. Whenever a region changes state, its elapsed

time should be reset to zero. Whilst remaining in its current state, the time elapsed

should be incremented by dT on every time step.

c) To plot the grid states, use the command

pcolor(H_state), caxis([0 3]), axis('square');

d) To generate a square array of random total refractory periods, use

RP = MRP + SD*randn(size(H_state));

e) To determine the total number of excited neighbours at any instant, use the code:

padded_exc = [zeros(1,52); (H_state(:,50) == 3), ...

(H_state == 3), (H_state(:,1) == 3); zeros(1,52)];

EN = padded_exc(1:50,1:50) + padded_exc(1:50,2:51) + ...

padded_exc(1:50,3:52) + padded_exc(2:51,1:50) + ...

padded_exc(2:51,3:52) + padded_exc(3:52,1:50) + ...

padded_exc(3:52,2:51) + padded_exc(3:52,3:52);

This provides an extra "padded" row and column around all edges of H_state,

allowing the direct summation of the "eight" neighbours.

Please hand in your Matlab code listing and plots of cardiac electrical activity at the

instant t = 1.7 s for both cases above. "

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Homework Help: MATLAB model of ventricular electrical activity

Can you offer guidance or do you also need help?

Draft saved
Draft deleted

**Physics Forums | Science Articles, Homework Help, Discussion**