From 15462f85106ed9ebde3d78ab6fe665c88bf8b32e Mon Sep 17 00:00:00 2001
From: JimPotvin <jpotvin1@cogeco.ca>
Date: Sun, 28 May 2017 22:56:37 -0400
Subject: [PATCH] removed recovery components

---
 ...e_Potvin_Fuglevand_2017_Muscle_Fatigue_Model.m | 15 ++++-----------
 1 file changed, 4 insertions(+), 11 deletions(-)

diff --git a/MatLab_Code_Potvin_Fuglevand_2017_Muscle_Fatigue_Model.m b/MatLab_Code_Potvin_Fuglevand_2017_Muscle_Fatigue_Model.m
index b84a471..5c52f15 100644
--- a/MatLab_Code_Potvin_Fuglevand_2017_Muscle_Fatigue_Model.m
+++ b/MatLab_Code_Potvin_Fuglevand_2017_Muscle_Fatigue_Model.m
@@ -1,10 +1,9 @@
 % Motor Unit Based Muscle Fatigue Model by Jim Potvin & Andrew Fuglevand
 % front end (rested size-principle) based on Fuglevand, Winter & Patla (1993)
-% last updated 2017-05-21 by Jim Potvin
+% last updated 2017-05-28 by Jim Potvin
 
 %clear all
 clc
-clf('reset')   %clears all graphics
 
 %% Model input parameters
     nu = 120;           % number of neurons (ie. motor units) in the modeled pool ("n") 
@@ -36,7 +35,7 @@
 
 %     % Create isotonic data -----------------------------------
 
-        fthscale = 0.50             % sets %MVC level for the trial duration (100% MVC is 1.00)
+        fthscale = 0.5             % sets %MVC level for the trial duration (100% MVC is 1.00)
         con = '0.50';               % for output file names
         fthtime = 100;              % duration to run trial (seconds)
         
@@ -178,16 +177,11 @@
     b2 = log(fat)/(nu-1);       
     mufatrate = exp(b2 * (n-1));  
 
-    b3 = log(rec)/(nu-1);       
-    murecrate = exp(b3 * (n-1));  
-
     fatigue = zeros(1,nu);
-    recovery = zeros(1,nu);
     for mu = 1:nu
          fatigue(mu) = mufatrate(mu) * (FatFac / fat) * P(mu);    
             % the full fatigue rate is mufatrate(mu) * [FatFac / fat] * Pr(mu,act) * P(mu)
             % the only variable is the relative force: Pr(mu,act), so this part is calculated once here
-         recovery(mu) = 0;   %set to zero for now
     end
     
 % Establishing the rested excitation required for each target load level (if 0.1% resolution, then 0.1% to 100%)    
@@ -337,7 +331,7 @@
         % Calculating the fatigue (force loss) for each motor unit
 
         for mu = 1:nu
-            if  mufrFAT(mu, i) >= recminfr                                          % Force loss of each MU for each interval
+            if  mufrFAT(mu, i) >= 0                                                 % Force loss of each MU for each interval
                 Pchange(mu,i) = -1 * (fatigue(mu) / samprate) * PrFAT(mu, i);       % based on % of MU fusion force
             elseif mufrFAT(mu, i) < recminfr                                        
                 Pchange(mu,i) = recovery(mu) / samprate;   
@@ -377,7 +371,6 @@
         end        
     end
     
-clf('reset')   %clears all graphics
                 
 %% Output
                 
@@ -385,7 +378,7 @@
 
 endurtime
             
-for mu = 0:mujump:nu
+for mu = 0:nu
     if mu == 0
         mu = 1;
     end