Revision 1351, 501 bytes
(checked in by zimamiro, 13 years ago)
|
|
Rev | Line | |
---|
[1351] | 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(3,:); |
---|
| 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 ylim]) |
---|
| 20 | |
---|
| 21 | end |
---|