|
Revision 1185, 0.6 kB
(checked in by zimamiro, 14 years ago)
|
|
|
| Line | |
|---|
| 1 | function rid_y0(sidp_parameters, compare_parameters,system,apriori,N)
|
|---|
| 2 | m=4;
|
|---|
| 3 | col=['b' 'g' 'r' 'c' 'm'];
|
|---|
| 4 | n=10;
|
|---|
| 5 | quality=zeros(1,n+1);
|
|---|
| 6 | apriori.b0=0;
|
|---|
| 7 | apriori.P0=10;
|
|---|
| 8 | y_max=apriori.y0;
|
|---|
| 9 | figure
|
|---|
| 10 | for j=1:m;
|
|---|
| 11 | q=0;
|
|---|
| 12 | apriori.Hc=j-1;
|
|---|
| 13 | sidp3(sidp_parameters, compare_parameters,system,apriori);
|
|---|
| 14 | for i=0:n
|
|---|
| 15 | q=q+1;
|
|---|
| 16 | apriori.y0=i*y_max/n;
|
|---|
| 17 | los=mc_study(system,apriori,N);
|
|---|
| 18 | quality(q)=los(2)/los(3);
|
|---|
| 19 | end
|
|---|
| 20 | plot((0:n)*y_max/n,quality(1,:),['',col(j)]); hold on
|
|---|
| 21 | end
|
|---|
| 22 | xlabel('\fontsize{18} y_0'); ylabel('\fontsize{18} relativn�tr�');
|
|---|
| 23 | legend('0)','1)','2)','3)');
|
|---|
| 24 | set( get(gcf, 'Children'), 'FontSize', 18);
|
|---|
| 25 | end |
|---|
