Revision 1351, 0.8 kB
(checked in by zimamiro, 13 years ago)
|
|
Line | |
---|
1 | function rid_y0(sidp_parameters, compare_parameters,system,apriori,N) |
---|
2 | m=3; |
---|
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 | subplot(2,1,1); |
---|
21 | plot((0:n)*y_max/n,quality(1,:),['',col(j)]); hold on |
---|
22 | [H C]=load_H; |
---|
23 | subplot(2,1,2); |
---|
24 | plot(H,C,['',col(j)]); hold on; a=max(a,max(C(:))); b=min(b,min(C(:)));end |
---|
25 | |
---|
26 | end |
---|
27 | xlabel('\fontsize{18} y_0'); ylabel('\fontsize{18} relativn�tr�'); |
---|
28 | legend('0)','1)','2)','3)'); |
---|
29 | set( get(gcf, 'Children'), 'FontSize', 18); |
---|
30 | end |
---|