Location: Tong_2011_V1 @ c34dd0cd0ede / code / Distribution-Moment-Approximation-Physiome-paper-1.0 / ZahM / lengthvar.m

Author:
Leyla <noroozbabaee@gmail.com>
Date:
2022-06-13 21:13:50+12:00
Desc:
Final Version Figure 1.
Permanent Source URI:
https://staging.physiomeproject.org/workspace/85c/rawfile/c34dd0cd0ede67cbef05bcb02934b699325457b3/code/Distribution-Moment-Approximation-Physiome-paper-1.0/ZahM/lengthvar.m

function [alpha,V] =lengthvar(t,f)

L=1*(t<180)+(1+(0.01*sin(2*pi*f*t)))*(t>=180 && t<=245)+(1+(0.04*sin(2*pi*f*t)))*(t>245);
gamma=25;
alpha=-gamma*(L-1);

%This
%is the V to use for DM

V=0*(t<180)+(-2*pi*gamma*0.01*f*cos(2*pi*f*t))*(t>=180 && t<=245)+(-2*pi*gamma*0.04*f*cos(2*pi*f*t))*(t>245);