| ax(1) =gca; |
| set(gcf, 'Position', [30, 30, 700, 720]); |
| set(gca, 'Position', [.115.10.78.8]); |
| f1=0:0.0001:1 |
| for i=1:7 |
| A = [0.1, 0.2, 0.5, 1, 2, 5, 10]; |
| B = [10, 5, 2, 1, 0.5, 0.2, 0.1]; |
| r1 = A(i); |
| r2 = B(i); |
| F1 =(r1*f1.^2+f1-f1.^2)./(r1*f1.^2+2*(f1-f1.^2)+r2*(1-f1).^2); |
| plot(f1, F1, 'linewidth', 2); |
| hold on |
| end |
| xlabel('{\itf}_1'); |
| ylabel('{\itF}_1'); |
| hold on |
| ax(2) = axes('position', get(ax(1), 'position')); |
| f2=1-f1 |
| for i=1:7 |
| A = [0.1, 0.2, 0.5, 1, 2, 5, 10]; |
| B = [10, 5, 2, 1, 0.5, 0.2, 0.1]; |
| r1 = A(i); |
| r2 = B(i); |
| F2=(r2*f2.^2+f2-f2.^2)./(r2*f2.^2+2*(f2-f2.^2)+r1*(1-f2).^2); |
| plot(f2, F2, '--', 'linewidth', 2); |
| hold on |
| end |
| set(ax(1), 'box', 'off'); |
| set(ax(1), 'ylim', [0, 1], 'yTick', [0:0.1:1]); |
| set(ax(2), 'XAxisLocation', 'top', 'YAxisLocation', 'right', 'box', 'off', 'color', 'none', 'XDir', 'reverse'); |
| set(ax(2), 'ylim', [0, 1], 'yTick', [0:0.1:1]); |
| xlabel('{\itf}_2'); |
| ylabel('{\itF}_2'); |
)
