Excitation cost allowable error and error center from xml

rcnl-org · Jun 14, 2023 · 5cac387 · 5cac387
1 parent 89844ce
commit 5cac387
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 3 deletions.
diff --git a/src/MuscleTendonPersonalization/Optimization/CostTerms/calcMuscleExcitationPenaltyCost.m b/src/MuscleTendonPersonalization/Optimization/CostTerms/calcMuscleExcitationPenaltyCost.m
@@ -27,13 +27,15 @@
 
 function cost = calcMuscleExcitationPenaltyCost(modeledValues, ...
     experimentalData, costTerm)
+errorCenter = valueOrAlternate(costTerm, "errorCenter", 0.5);
+maximumAllowableError = valueOrAlternate(costTerm, "maxAllowableError", 0.25);
 muscleExcitationsConstraint = modeledValues.muscleExcitationsNoTDelay(: , ...
     setdiff(1 : size(modeledValues.muscleExcitationsNoTDelay, 2), ...
     [experimentalData.synergyExtrapolation.missingEmgChannelGroups{:}]), ...
     experimentalData.numPaddingFrames + 1 : ...
     size(modeledValues.muscleExcitationsNoTDelay, 3) - ...
     experimentalData.numPaddingFrames);
-cost = 120 * (muscleExcitationsConstraint - 0.5) .^ 8;
+cost = 30 / maximumAllowableError * (muscleExcitationsConstraint - errorCenter) .^ 8;
 cost(isnan(cost))=0;
-cost = sum([sqrt(0.1) .* cost ].^ 2, 'all');
+cost = sum((sqrt(0.1) .* cost).^ 2, 'all');
 end
diff --git a/src/MuscleTendonPersonalization/Optimization/calcMtpCost.m b/src/MuscleTendonPersonalization/Optimization/calcMtpCost.m
@@ -76,7 +76,7 @@
                     synxModeledValues, modeledValues, experimentalData, costTerm);
             case "muscle_excitation_penalty"
                 cost = calcMuscleExcitationPenaltyCost( ...
-                    synxModeledValues, experimentalData);%, costTerm);
+                    synxModeledValues, experimentalData, costTerm);
             otherwise
                 throw(MException('', 'Cost term %s is not valid for MTP', ...
                     costTerm.type))