Revision 1351, 511 bytes
(checked in by zimamiro, 13 years ago)
|
|
Line | |
---|
1 | function mc_study(system,apriori,n)
|
---|
2 |
|
---|
3 | for i=1:n
|
---|
4 | system.x=apriori.H(:,1)+sqrt(diag(apriori.H(:,2:5))).*randn(4,1);
|
---|
5 | [losses(1,i) losses(2,i) losses(3,i)]=srovnej(system,apriori);
|
---|
6 | end
|
---|
7 |
|
---|
8 | %losses(1,:)=losses(1,:)./losses(3,:);
|
---|
9 | losses(2,:)=losses(2,:)./losses(1,:);
|
---|
10 | %losses(3,:)=losses(3,:)./losses(3,:);
|
---|
11 |
|
---|
12 | %hist(losses(1,:));
|
---|
13 | %title([mean(losses(1,:)) median(losses(1,:))])
|
---|
14 | %AXIS([0 5 ylim])
|
---|
15 |
|
---|
16 | figure
|
---|
17 | hist(losses(2,:));
|
---|
18 | title([mean(losses(2,:)) median(losses(2,:))])
|
---|
19 | AXIS([0.5 1.5 ylim])
|
---|
20 |
|
---|
21 | end |
---|