dmakessfigure13.m
来自「beamlet变化的工具箱」· M 代码 · 共 15 行 · 第 1/2 页
M
15 行
function MakeSSFigure13(varargin)
% MakeSSFigure13.m -- Figure 13 of Slant Stack Paper
% Usage: MAKESSFIGURE13 will start a new figure window and plot Figure 13 of Slant Stack Paper.
% MAKESSFIGURE13(n) for any argument type n will plot the figure on the current figure
% window and plot Figure 13 of Slant Stack Paper.
if nargin == 0
figure;
else
gcf;
end
% MakehnkFig -- Process output of RefinedBiDensityPartEvals
% for n=16,32,64,128,256
% *** n = 16 ***
%
% k= lambdamin lambdamax ratio Norm hnk^ext
hn16 = [ ...
1 3.785 5.658 1.495 0.806 ...
2 3.867 5.657 1.463 0.772 ...
3 3.937 5.657 1.437 0.742 ...
4 3.966 5.657 1.426 0.727 ...
5 3.982 5.657 1.420 0.702 ...
6 3.997 5.657 1.415 0.695 ...
7 4.000 5.657 1.414 0.669 ...
8 4.000 5.657 1.414 0.667 ...
9 4.000 5.657 1.414 0.635 ...
10 4.000 5.657 1.414 0.639 ...
11 3.999 5.657 1.414 0.594 ...
12 3.999 5.657 1.415 0.605 ...
13 3.983 5.654 1.420 0.526 ...
14 3.966 5.633 1.420 0.556 ...
15 3.968 5.528 1.393 0.000 ...
];
Q16 = hn16( 5:5:(15*5) );
A16 = (1:15)./16;
%figure
plot(A16,Q16,'r');
hold on;
axis([0 1 0 1]);
xlabel('\alpha = k/n')
ylabel('||H_{n,k}^{ext}||')
plot([.5 .6],[.85 .85],'-r');
th = text(.65,.85,'n=16')
set(th,'Color',[1 0 0]); % red
set(th,'FontSize',12);
%
% *** n = 32 ***
%
% k= lambdamin lambdamax ratio Norm hnk^ext
hn32 = [ ...
1 5.310 8.002 1.507 0.819 ...
2 5.389 8.000 1.484 0.790 ...
3 5.458 8.000 1.466 0.771 ...
4 5.512 8.000 1.451 0.759 ...
5 5.552 8.000 1.441 0.746 ...
6 5.568 8.000 1.437 0.738 ...
7 5.585 8.000 1.432 0.727 ...
8 5.602 8.000 1.428 0.721 ...
9 5.620 8.000 1.424 0.710 ...
10 5.635 8.000 1.420 0.705 ...
11 5.647 8.000 1.417 0.695 ...
12 5.652 8.000 1.415 0.692 ...
13 5.655 8.000 1.415 0.681 ...
14 5.657 8.000 1.414 0.679 ...
15 5.657 8.000 1.414 0.667 ...
16 5.657 8.000 1.414 0.667 ...
17 5.657 8.000 1.414 0.654 ...
18 5.657 8.000 1.414 0.655 ...
19 5.657 8.000 1.414 0.640 ...
20 5.657 8.000 1.414 0.642 ...
21 5.657 8.000 1.414 0.625 ...
22 5.657 8.000 1.414 0.629 ...
23 5.656 8.000 1.414 0.607 ...
24 5.655 8.000 1.415 0.614 ...
25 5.654 8.000 1.415 0.585 ...
26 5.654 8.000 1.415 0.595 ...
27 5.642 8.000 1.418 0.553 ...
28 5.621 8.000 1.423 0.569 ...
29 5.602 7.997 1.427 0.496 ...
30 5.596 7.981 1.426 0.532 ...
31 5.618 7.903 1.407 0.000 ...
];
Q32 = hn32( 5:5:(31*5) );
A32 = (1:31)./32;
plot(A32,Q32,'b');
height = .80
plot([.5 .6],[height height],'-b');
th = text(.65,height,'n=32')
set(th,'Color',[0 0 1]); % blue
set(th,'FontSize',12);
%
% *** n = 64 ***
%
% k= lambdamin lambdamax ratio Norm hnk^ext
hn64 = [ ...
1 7.478 11.316 1.513 0.825 ...
2 7.561 11.314 1.496 0.800 ...
3 7.632 11.314 1.482 0.786 ...
4 7.686 11.314 1.472 0.776 ...
5 7.729 11.314 1.464 0.768 ...
6 7.763 11.314 1.457 0.761 ...
7 7.789 11.314 1.453 0.755 ...
8 7.802 11.314 1.450 0.750 ...
9 7.817 11.314 1.447 0.744 ...
10 7.832 11.314 1.445 0.740 ...
11 7.848 11.314 1.442 0.735 ...
12 7.864 11.314 1.439 0.731 ...
13 7.880 11.314 1.436 0.726 ...
14 7.896 11.314 1.433 0.723 ...
15 7.912 11.314 1.430 0.718 ...
16 7.926 11.314 1.427 0.716 ...
17 7.940 11.314 1.425 0.711 ...
18 7.953 11.314 1.423 0.708 ...
19 7.964 11.314 1.421 0.704 ...
20 7.972 11.314 1.419 0.702 ...
21 7.978 11.314 1.418 0.697 ...
22 7.984 11.314 1.417 0.695 ...
23 7.989 11.314 1.416 0.691 ...
24 7.993 11.314 1.415 0.689 ...
25 7.996 11.314 1.415 0.684 ...
26 7.998 11.314 1.415 0.683 ...
27 7.999 11.314 1.414 0.678 ...
28 8.000 11.314 1.414 0.678 ...
29 8.000 11.314 1.414 0.673 ...
30 8.000 11.314 1.414 0.672 ...
31 8.000 11.314 1.414 0.667 ...
32 8.000 11.314 1.414 0.667 ...
33 8.000 11.314 1.414 0.661 ...
34 8.000 11.314 1.414 0.661 ...
35 8.000 11.314 1.414 0.656 ...
36 8.000 11.314 1.414 0.656 ...
37 8.000 11.314 1.414 0.650 ...
38 8.000 11.314 1.414 0.651 ...
39 8.000 11.314 1.414 0.644 ...
40 8.000 11.314 1.414 0.645 ...
41 8.000 11.314 1.414 0.638 ...
42 8.000 11.314 1.414 0.640 ...
43 8.000 11.314 1.414 0.631 ...
44 8.000 11.314 1.414 0.634 ...
45 8.000 11.314 1.414 0.625 ...
46 7.999 11.314 1.414 0.627 ...
47 7.999 11.314 1.414 0.617 ...
48 7.999 11.314 1.414 0.621 ...
49 7.999 11.314 1.414 0.609 ...
50 7.999 11.314 1.414 0.613 ...
51 7.998 11.314 1.415 0.599 ...
52 7.997 11.314 1.415 0.604 ...
53 7.996 11.314 1.415 0.587 ...
54 7.995 11.314 1.415 0.594 ...
55 7.995 11.314 1.415 0.575 ...
56 7.993 11.314 1.415 0.581 ...
57 7.974 11.314 1.419 0.567 ...
58 7.951 11.314 1.423 0.564 ...
59 7.930 11.314 1.427 0.555 ...
60 7.913 11.314 1.430 0.543 ...
61 7.907 11.311 1.431 0.519 ...
62 7.919 11.300 1.427 0.517 ...
63 7.959 11.243 1.413 0.000 ...
];
Q64 = hn64( 5:5:(63*5) );
A64 = (1:63)./64;
plot(A64,Q64,'c');
height = .75
plot([.5 .6],[height height],'-c');
th = text(.65,height,'n=64')
set(th,'Color',[0 1 1]); % yellow
set(th,'FontSize',12);
% *** n = 128 ***
%
% k= lambdamin lambdamax ratio Norm hnk^ext ...
hn128 = [ ...
1 10.552 16.004 1.517 0.828 ...
2 10.649 16.000 1.503 0.805 ...
3 10.729 16.000 1.491 0.793 ...
4 10.787 16.000 1.483 0.785 ...
5 10.830 16.000 1.477 0.779 ...
6 10.865 16.000 1.473 0.774 ...
7 10.893 16.000 1.469 0.770 ...
8 10.917 16.000 1.466 0.766 ...
9 10.934 16.000 1.463 0.762 ...
10 10.944 16.000 1.462 0.759 ...
11 10.955 16.000 1.461 0.756 ...
12 10.967 16.000 1.459 0.753 ...
13 10.979 16.000 1.457 0.750 ...
14 10.992 16.000 1.456 0.748 ...
15 11.005 16.000 1.454 0.745 ...
16 11.018 16.000 1.452 0.743 ...
17 11.032 16.000 1.450 0.741 ...
18 11.045 16.000 1.449 0.739 ...
19 11.059 16.000 1.447 0.736 ...
20 11.072 16.000 1.445 0.734 ...
21 11.086 16.000 1.443 0.732 ...
22 11.099 16.000 1.442 0.730 ...
23 11.112 16.000 1.440 0.728 ...
24 11.125 16.000 1.438 0.726 ...
25 11.137 16.000 1.437 0.724 ...
26 11.149 16.000 1.435 0.722 ...
27 11.161 16.000 1.434 0.720 ...
28 11.173 16.000 1.432 0.719 ...
29 11.184 16.000 1.431 0.717 ...
30 11.195 16.000 1.429 0.715 ...
31 11.206 16.000 1.428 0.713 ...
32 11.216 16.000 1.427 0.712 ...
33 11.225 16.000 1.425 0.710 ...
34 11.232 16.000 1.424 0.708 ...
35 11.240 16.000 1.424 0.706 ...
36 11.247 16.000 1.423 0.705 ...
37 11.253 16.000 1.422 0.703 ...
38 11.260 16.000 1.421 0.702 ...
39 11.266 16.000 1.420 0.700 ...
40 11.272 16.000 1.419 0.699 ...
41 11.277 16.000 1.419 0.697 ...
42 11.283 16.000 1.418 0.696 ...
43 11.287 16.000 1.418 0.694 ...
44 11.292 16.000 1.417 0.693 ...
45 11.296 16.000 1.416 0.691 ...
46 11.299 16.000 1.416 0.690 ...
47 11.302 16.000 1.416 0.688 ...
48 11.304 16.000 1.415 0.687 ...
49 11.306 16.000 1.415 0.685 ...
50 11.308 16.000 1.415 0.685 ...
51 11.310 16.000 1.415 0.682 ...
52 11.311 16.000 1.415 0.682 ...
53 11.312 16.000 1.414 0.680 ...
54 11.313 16.000 1.414 0.679 ...
55 11.313 16.000 1.414 0.677 ...
56 11.313 16.000 1.414 0.677 ...
57 11.314 16.000 1.414 0.674 ...
58 11.314 16.000 1.414 0.674 ...
59 11.314 16.000 1.414 0.672 ...
60 11.314 16.000 1.414 0.672 ...
61 11.314 16.000 1.414 0.669 ...
62 11.314 16.000 1.414 0.669 ...
63 11.314 16.000 1.414 0.667 ...
64 11.314 16.000 1.414 0.667 ...
65 11.314 16.000 1.414 0.664 ...
66 11.314 16.000 1.414 0.664 ...
67 11.314 16.000 1.414 0.662 ...
68 11.314 16.000 1.414 0.662 ...
69 11.314 16.000 1.414 0.659 ...
70 11.314 16.000 1.414 0.659 ...
71 11.314 16.000 1.414 0.657 ...
72 11.314 16.000 1.414 0.657 ...
73 11.314 16.000 1.414 0.654 ...
74 11.314 16.000 1.414 0.655 ...
75 11.314 16.000 1.414 0.652 ...
76 11.314 16.000 1.414 0.652 ...
77 11.314 16.000 1.414 0.649 ...
78 11.314 16.000 1.414 0.650 ...
79 11.314 16.000 1.414 0.647 ...
80 11.314 16.000 1.414 0.648 ...
81 11.314 16.000 1.414 0.644 ...
82 11.314 16.000 1.414 0.645 ...
83 11.314 16.000 1.414 0.642 ...
84 11.314 16.000 1.414 0.643 ...
85 11.314 16.000 1.414 0.639 ...
86 11.313 16.000 1.414 0.640 ...
87 11.313 16.000 1.414 0.636 ...
88 11.313 16.000 1.414 0.638 ...
89 11.313 16.000 1.414 0.634 ...
90 11.313 16.000 1.414 0.635 ...
91 11.313 16.000 1.414 0.631 ...
92 11.313 16.000 1.414 0.632 ...
93 11.313 16.000 1.414 0.628 ...
94 11.313 16.000 1.414 0.629 ...
95 11.313 16.000 1.414 0.625 ...
96 11.313 16.000 1.414 0.627 ...
97 11.313 16.000 1.414 0.621 ...
98 11.312 16.000 1.414 0.623 ...
99 11.312 16.000 1.414 0.618 ...
100 11.312 16.000 1.414 0.620 ...
101 11.312 16.000 1.414 0.616 ...
102 11.311 16.000 1.415 0.616 ...
103 11.311 16.000 1.415 0.613 ...
104 11.311 16.000 1.415 0.613 ...
105 11.311 16.000 1.415 0.611 ...
106 11.310 16.000 1.415 0.608 ...
107 11.310 16.000 1.415 0.608 ...
108 11.310 16.000 1.415 0.604 ...
109 11.309 16.000 1.415 0.606 ...
110 11.308 16.000 1.415 0.598 ...
111 11.307 16.000 1.415 0.604 ...
112 11.306 16.000 1.415 0.593 ...
113 11.305 16.000 1.415 0.602 ...
114 11.305 16.000 1.415 0.588 ...
115 11.305 16.000 1.415 0.600 ...
116 11.306 16.000 1.415 0.584 ...
117 11.298 16.000 1.416 0.598 ...
118 11.277 16.000 1.419 0.577 ...
119 11.254 16.000 1.422 0.596 ...
120 11.231 16.000 1.425 0.567 ...
121 11.209 16.000 1.427 0.594 ...
122 11.191 16.000 1.430 0.547 ...
123 11.180 16.000 1.431 0.584 ...
124 11.178 16.000 1.431 0.525 ...
125 11.191 15.997 1.429 0.542 ...
126 11.224 15.990 1.425 0.509 ...
127 11.275 15.949 1.414 0.000 ];
Q128 = hn128( 5:5:(127*5) );
A128 = (1:127)./128;
plot(A128,Q128,'g');
height = .70
plot([.5 .6],[height height],'-g');
th = text(.65,height,'n=128')
set(th,'Color',[0 1 0]); % green
set(th,'FontSize',12);
% hn256 = [ ...
% 1 5.602 16.732 2.987 0.413 ...
% 2 5.953 16.918 2.842 0.465 ...
% 3 6.124 17.016 2.779 0.485 ...
% 4 6.237 17.079 2.739 0.494 ...
% 5 6.322 17.125 2.709 0.497 ...
% 6 6.393 17.160 2.684 0.499 ...
% 7 6.453 17.188 2.664 0.498 ...
% 8 6.507 17.212 2.645 0.496 ...
% 9 6.556 17.231 2.628 0.494 ...
% 10 6.601 17.248 2.613 0.492 ...
% 11 6.644 17.263 2.598 0.489 ...
% 12 6.684 17.276 2.585 0.485 ...
% 13 6.722 17.287 2.572 0.482 ...
% 14 6.759 17.298 2.559 0.478 ...
% 15 6.794 17.307 2.547 0.475 ...
% 16 6.828 17.316 2.536 0.471 ...
% 17 6.861 17.324 2.525 0.467 ...
% 18 6.894 17.331 2.514 0.464 ...
% 19 6.925 17.338 2.504 0.460 ...
% 20 6.956 17.344 2.493 0.456 ...
% 21 6.986 17.349 2.483 0.452 ...
% 22 7.016 17.355 2.474 0.448 ...
% 23 7.045 17.360 2.464 0.445 ...
% 24 7.074 17.364 2.455 0.441 ...
% 25 7.102 17.369 2.446 0.437 ...
% 26 7.130 17.373 2.437 0.434 ...
% 27 7.157 17.376 2.428 0.430 ...
% 28 7.184 17.380 2.419 0.426 ...
% 29 7.211 17.383 2.411 0.423 ...
% 30 7.237 17.386 2.402 0.419 ...
% 31 7.263 17.389 2.394 0.416 ...
% 32 7.289 17.392 2.386 0.412 ...
% 33 7.314 17.395 2.378 0.409 ...
% 34 7.339 17.397 2.370 0.405 ...
% 35 7.364 17.399 2.363 0.402 ...
% 36 7.388 17.401 2.355 0.398 ...
% 37 7.413 17.403 2.348 0.395 ...
% 38 7.437 17.405 2.340 0.391 ...
% 39 7.461 17.407 2.333 0.388 ...
% 40 7.484 17.408 2.326 0.385 ...
% 41 7.508 17.410 2.319 0.382 ...
% 42 7.531 17.411 2.312 0.378 ...
% 43 7.554 17.413 2.305 0.375 ...
% 44 NaN NaN NaN 0.372 ...
% 45 7.600 17.415 2.292 0.369 ...
% 46 7.622 17.416 2.285 0.366 ...
% 47 7.644 17.417 2.278 0.362 ...
% 48 7.666 17.418 2.272 0.359 ...
% 49 7.688 17.418 2.266 0.356 ...
% 50 7.710 17.419 2.259 0.353 ...
% 51 7.732 17.419 2.253 0.350 ...
% 52 7.753 17.420 2.247 0.347 ...
% 53 7.774 17.420 2.241 0.344 ...
% 54 7.796 17.420 2.235 0.341 ...
% 55 7.817 17.420 2.229 0.338 ...
% 56 7.837 17.421 2.223 0.335 ...
% 57 7.858 17.421 2.217 0.333 ...
% 58 7.879 17.420 2.211 0.330 ...
% 59 7.899 17.420 2.205 0.327 ...
% 60 7.920 17.420 2.200 0.324 ...
% 61 7.940 17.420 2.194 0.321 ...
% 62 7.960 17.419 2.188 0.318 ...
% 63 7.980 17.419 2.183 0.316 ...
% 64 8.000 17.418 2.177 0.313 ...
% 65 8.020 17.417 2.172 0.310 ...
% 66 8.040 17.417 2.166 0.307 ...
% 67 8.059 17.416 2.161 0.305 ...
% 68 8.079 17.415 2.156 0.302 ...
% 69 8.098 17.414 2.150 0.299 ...
% 70 8.118 17.413 2.145 0.297 ...
% 71 8.137 17.411 2.140 0.294 ...
% 72 8.157 17.410 2.134 0.291 ...
% 73 8.176 17.409 2.129 0.289 ...
% 74 8.196 17.407 2.124 0.286 ...
% 75 8.215 17.405 2.119 0.284 ...
% 76 8.234 17.404 2.114 0.281 ...
% 77 8.254 17.402 2.108 0.278 ...
% 78 8.273 17.400 2.103 0.276 ...
% 79 8.292 17.397 2.098 0.273 ...
% 80 8.312 17.395 2.093 0.271 ...
% 81 8.331 17.393 2.088 0.268 ...
% 82 8.351 17.390 2.082 0.266 ...
% 83 8.370 17.387 2.077 0.263 ...
% 84 8.390 17.384 2.072 0.261 ...
% 85 8.410 17.381 2.067 0.258 ...
% 86 8.430 17.378 2.061 0.255 ...
% 87 8.450 17.375 2.056 0.253 ...
% 88 8.471 17.371 2.051 0.250 ...
% 89 8.491 17.367 2.045 0.248 ...
% 90 8.512 17.363 2.040 0.245 ...
% 91 8.533 17.359 2.034 0.242 ...
% 92 8.555 17.354 2.029 0.240 ...
% 93 8.577 17.349 2.023 0.237 ...
% 94 8.599 17.344 2.017 0.234 ...
% 95 8.621 17.339 2.011 0.232 ...
% 96 8.644 17.333 2.005 0.229 ...
% 97 8.668 17.327 1.999 0.226 ...
% 98 8.692 17.320 1.993 0.223 ...
% 99 8.717 17.313 1.986 0.220 ...
% 100 8.742 17.305 1.980 0.217 ...
% 101 8.769 17.297 1.973 0.214 ...
% 102 8.796 17.289 1.966 0.211 ...
% 103 8.824 17.280 1.958 0.208 ...
% 104 8.853 17.270 1.951 0.205 ...
% 105 8.884 17.259 1.943 0.202 ...
% 106 8.916 17.247 1.934 0.198 ...
% 107 8.950 17.234 1.926 0.194 ...
% 108 8.985 17.221 1.917 0.191 ...
% 109 9.023 17.205 1.907 0.187 ...
% 110 9.063 17.189 1.897 0.183 ...
% 111 9.106 17.170 1.886 0.179 ...
% 112 9.152 17.150 1.874 0.174 ...
% 113 9.202 17.127 1.861 0.169 ...
% 114 9.257 17.102 1.847 0.164 ...
% 115 9.317 17.073 1.832 0.159 ...
% 116 9.383 17.039 1.816 0.153 ...
% 117 9.457 17.001 1.798 0.147 ...
% 118 9.540 16.957 1.777 0.140 ...
% 119 9.635 16.904 1.754 0.132 ...
% 120 9.745 16.842 1.728 0.124 ...
% 121 9.872 16.767 1.698 0.114 ...
% 122 10.022 16.676 1.664 0.102 ...
% 123 10.199 16.567 1.624 0.089 ...
% 124 10.408 16.438 1.579 0.073 ...
% 125 10.650 16.292 1.530 0.053 ...
% 126 10.916 16.144 1.479 0.029 ...
% 127 11.171 16.029 1.435 0.000 ];
% Q256 = hn256( 5:5:(127*5) );
% A256 = (1:127)./256;
% plot(2*A256,Q256,'k');
% height = .65
% plot([.5 .6],[height height],'-k');
% th = text(.65,height,'n=256')
% set(th,'Color',[0 0 0]); % black
% set(th,'FontSize',12);
%
%% Part of BeamLab Version:200% Built:Friday,23-Aug-2002 00:00:00% This is Copyrighted Material% For Copying permissions see COPYING.m% Comments? e-mail beamlab@stat.stanford.edu%⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?