Simulating 100000 AMP data
clear
A = xlsread('200sample', 'W1:AL16
                                       );#A is the state transition probability matrix TR
B = zeros(1,100000);
B(1) = 16;
i = 16;
s = zeros(1,16);
n = 2;
for n = 2:100000
r = rand(1);
i = B(n-1);
s(1) = A(i,1);
for j = 2:16
s(j) = s(j-1)+A(i,j);
if r >= 0&&r <= s(1)
B(n) = 1;
elseif r >= s(j-1)&&r <= s(j)
B(n) = j;
end
end
end
B
for i = 1:100000
r = rand(1);
B(i) = 1310+B(i)*5-5*r;
end
```
#### **References**

B'


© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Energies* Editorial Office E-mail: energies@mdpi.com www.mdpi.com/journal/energies

MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel: +41 61 683 77 34

www.mdpi.com

ISBN 978-3-0365-5712-0