-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathplotCluster.m
63 lines (45 loc) · 1.23 KB
/
plotCluster.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
function plotCluster
for j = 2:13
s = 2;
strg=sprintf('%%.%dd',s);
indxStr=sprintf(strg,j);
v = textread(['C:\Alex\R\R-2.3.1\avgV_786OVHL30_NaCl',indxStr,'.txt'],'%n');
dx = sort(v);
le = length(v);
c = textread(['C:\Alex\R\R-2.3.1\avgV_786OVHL30_NaCl',indxStr,'_distribution_3G.txt'],'%n');
n = c(1);
p =zeros(le,1);
figure,histogram(v,'smooth')
hold on
for i = 1:n
g = le * c(i+1+2*n) * normpdf(dx,c(i+1),sqrt(c(i+1+n)));
p = p + g;
plot(dx,g,'r--')
clear g
hold on
end
plot(dx,p,'r-')
hold off
end
for j = 1:12
s = 2;
strg=sprintf('%%.%dd',s);
indxStr=sprintf(strg,j);
v = textread(['C:\Alex\R\R-2.3.1\avgV_786Opar_NaCl',indxStr,'.txt'],'%n');
dx = sort(v);
le = length(v);
c = textread(['C:\Alex\R\R-2.3.1\avgV_786Opar_NaCl',indxStr,'_distribution_3G.txt'],'%n');
n = c(1);
p =zeros(le,1);
figure,histogram(v,'smooth')
hold on
for i = 1:n
g = le * c(i+1+2*n) * normpdf(dx,c(i+1),sqrt(c(i+1+n)));
p = p + g;
plot(dx,g,'r--')
clear g
hold on
end
plot(dx,p,'r-')
hold off
end