|
Revision 1105, 0.7 kB
(checked in by zimamiro, 15 years ago)
|
|
|
| Line | |
|---|
| 1 | function losses=mc_study2(system,range,n)
|
|---|
| 2 |
|
|---|
| 3 | load 'eta.txt';
|
|---|
| 4 | load 'Beta.txt';
|
|---|
| 5 | load 'ny.txt';
|
|---|
| 6 | H1(:,:,1)=eta;
|
|---|
| 7 | H1(:,:,2)=Beta;
|
|---|
| 8 | u1=ny;
|
|---|
| 9 |
|
|---|
| 10 | load 'eta2.txt';
|
|---|
| 11 | load 'Beta2.txt';
|
|---|
| 12 | load 'ny2.txt';
|
|---|
| 13 | H2(:,:,1)=eta2;
|
|---|
| 14 | H2(:,:,2)=Beta2;
|
|---|
| 15 | u2=ny2;
|
|---|
| 16 | for i=1:n
|
|---|
| 17 | apriori(3)=range.P0;
|
|---|
| 18 | system.b=range.b0 + apriori(3)*randn;
|
|---|
| 19 | apriori(2)=range.b0 + apriori(3)*randn;
|
|---|
| 20 | apriori(1)=range.y0 + range.y0_range*randn;
|
|---|
| 21 | while (system.b==0)
|
|---|
| 22 | system.b=range.b0 + apriori(3)*randn;
|
|---|
| 23 | end
|
|---|
| 24 | while (apriori(2)==0)
|
|---|
| 25 | apriori(2)=range.b0 + apriori(3)*randn;
|
|---|
| 26 | end
|
|---|
| 27 | seed=randn(1,size(H1,1)+2);
|
|---|
| 28 | losses(1,i)=rizeni(H1,u1,system,apriori,seed);
|
|---|
| 29 | % losses(2,i)=rizeni(H2,u2,system,apriori,seed);
|
|---|
| 30 | end
|
|---|
| 31 | end |
|---|
