|
Revision 1105, 0.6 kB
(checked in by zimamiro, 14 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+randn*range.P0_range; |
|---|
| 18 | system.b=apriori(3)*randn; |
|---|
| 19 | apriori(2)=apriori(3)*randn; |
|---|
| 20 | apriori(1)=range.y0_range*randn; |
|---|
| 21 | while (system.b==0) |
|---|
| 22 | system.b=apriori(3)*randn; |
|---|
| 23 | end |
|---|
| 24 | while (apriori(2)==0) |
|---|
| 25 | apriori(2)=apriori(3)*randn; |
|---|
| 26 | end |
|---|
| 27 | losses(1,i)=rizeni(H1,u1,system,apriori); |
|---|
| 28 | losses(2,i)=rizeni(H2,u2,system,apriori); |
|---|
| 29 | end |
|---|
| 30 | end |
|---|
