File talk:Judd's (r,g) diagram.svg
Jump to navigation
Jump to search
Source
[edit]r=[0 0 0 0 0 0.001 0.005 0.023 0.095 0.136 0.177 0.238 0.287 0.33 0.368 0.399 0.438 0.487 0.55 0.639 0.719 0.801 0.862 0.906 0.951 0.975 0.985 1;
0.001 0.002 0.002 0.002 0.003 0.005 0.014 0.037 0.116 0.157 0.196 0.251 0.294 0.331 0.363 0.389 0.423 0.465 0.52 0.6 0.669 0.743 0.799 0.839 0.882 0.904 0.912 0.928;
0.003 0.004 0.004 0.004 0.006 0.011 0.023 0.053 0.137 0.178 0.214 0.263 0.3 0.331 0.358 0.381 0.409 0.446 0.493 0.561 0.623 0.689 0.739 0.776 0.813 0.834 0.843 0.856;
0.005 0.006 0.006 0.007 0.01 0.017 0.035 0.072 0.159 0.198 0.231 0.274 0.305 0.331 0.354 0.373 0.397 0.427 0.468 0.526 0.58 0.637 0.681 0.713 0.747 0.766 0.774 0.784;
0.008 0.009 0.01 0.011 0.015 0.026 0.049 0.092 0.182 0.218 0.248 0.284 0.31 0.331 0.351 0.366 0.386 0.411 0.444 0.493 0.539 0.588 0.626 0.654 0.683 0.699 0.706 0.716;
0.012 0.013 0.014 0.016 0.022 0.037 0.067 0.117 0.205 0.238 0.264 0.294 0.315 0.332 0.347 0.359 0.375 0.395 0.422 0.462 0.5 0.54 0.572 0.596 0.621 0.635 0.641 0.65;
0.018 0.02 0.021 0.024 0.032 0.052 0.09 0.145 0.229 0.258 0.279 0.303 0.319 0.332 0.344 0.353 0.366 0.381 0.402 0.433 0.463 0.495 0.521 0.54 0.56 0.572 0.576 0.584;
0.028 0.03 0.032 0.036 0.047 0.073 0.121 0.179 0.254 0.277 0.293 0.311 0.323 0.333 0.341 0.348 0.357 0.368 0.383 0.406 0.428 0.452 0.472 0.486 0.501 0.51 0.513 0.52;
0.045 0.048 0.051 0.057 0.073 0.109 0.163 0.22 0.279 0.296 0.307 0.319 0.327 0.333 0.338 0.343 0.348 0.356 0.365 0.381 0.395 0.411 0.424 0.434 0.444 0.45 0.452 0.456;
0.087 0.092 0.096 0.106 0.129 0.174 0.227 0.27 0.306 0.315 0.321 0.326 0.33 0.333 0.336 0.338 0.34 0.344 0.349 0.356 0.363 0.371 0.378 0.382 0.388 0.391 0.392 0.394;
0.142 0.148 0.154 0.165 0.19 0.232 0.273 0.3 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32 0.32];
g=[0 0.002 0.003 0.007 0.017 0.044 0.105 0.203 0.397 0.492 0.562 0.625 0.638 0.624 0.604 0.584 0.551 0.506 0.446 0.359 0.28 0.199 0.138 0.094 0.049 0.025 0.015 0;
0.001 0.004 0.005 0.009 0.02 0.047 0.111 0.209 0.391 0.475 0.534 0.585 0.594 0.582 0.564 0.548 0.52 0.482 0.431 0.355 0.287 0.216 0.161 0.122 0.081 0.059 0.051 0.036;
0.003 0.006 0.007 0.011 0.023 0.052 0.117 0.216 0.386 0.458 0.508 0.548 0.555 0.544 0.53 0.515 0.492 0.459 0.416 0.352 0.293 0.231 0.183 0.148 0.113 0.093 0.084 0.072;
0.005 0.008 0.009 0.014 0.026 0.058 0.125 0.223 0.38 0.442 0.483 0.514 0.519 0.509 0.497 0.485 0.466 0.439 0.403 0.349 0.299 0.246 0.205 0.175 0.143 0.125 0.118 0.108;
0.008 0.011 0.013 0.017 0.031 0.066 0.136 0.232 0.374 0.426 0.458 0.483 0.486 0.478 0.467 0.458 0.442 0.42 0.391 0.347 0.305 0.26 0.225 0.199 0.172 0.158 0.151 0.142;
0.012 0.015 0.017 0.023 0.037 0.075 0.148 0.242 0.368 0.41 0.435 0.453 0.455 0.449 0.44 0.433 0.42 0.403 0.38 0.344 0.31 0.274 0.245 0.224 0.201 0.188 0.183 0.175;
0.018 0.021 0.024 0.031 0.047 0.088 0.164 0.254 0.361 0.394 0.413 0.426 0.427 0.422 0.415 0.41 0.4 0.387 0.369 0.342 0.315 0.287 0.264 0.247 0.229 0.219 0.215 0.208;
0.028 0.032 0.034 0.042 0.061 0.107 0.185 0.268 0.355 0.379 0.392 0.401 0.401 0.397 0.393 0.388 0.381 0.372 0.359 0.339 0.32 0.299 0.282 0.27 0.256 0.249 0.246 0.24;
0.045 0.05 0.053 0.063 0.086 0.138 0.215 0.285 0.348 0.363 0.371 0.377 0.377 0.374 0.372 0.369 0.364 0.358 0.35 0.337 0.325 0.311 0.3 0.292 0.282 0.277 0.275 0.272;
0.087 0.093 0.098 0.111 0.139 0.194 0.26 0.307 0.341 0.348 0.352 0.355 0.354 0.353 0.352 0.35 0.348 0.345 0.341 0.335 0.329 0.322 0.317 0.313 0.308 0.306 0.305 0.303;
0.142 0.149 0.155 0.169 0.197 0.245 0.292 0.319 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337 0.337];
b=[1 0.998 0.997 0.993 0.983 0.955 0.89 0.774 0.508 0.372 0.261 0.137 0.075 0.046 0.028 0.017 0.011 0.007 0.004 0.002 0.001 0 0 0 0 0 0 0;
0.998 0.994 0.993 0.989 0.978 0.948 0.875 0.754 0.493 0.368 0.27 0.164 0.112 0.087 0.073 0.063 0.057 0.053 0.049 0.045 0.044 0.041 0.04 0.039 0.037 0.037 0.037 0.036;
0.994 0.99 0.989 0.985 0.971 0.937 0.86 0.731 0.477 0.364 0.278 0.189 0.145 0.125 0.112 0.104 0.099 0.095 0.091 0.087 0.084 0.08 0.078 0.076 0.074 0.073 0.073 0.072;
0.99 0.986 0.985 0.979 0.964 0.925 0.84 0.705 0.461 0.36 0.286 0.212 0.176 0.16 0.149 0.142 0.137 0.134 0.129 0.125 0.121 0.117 0.114 0.112 0.11 0.109 0.108 0.108;
0.984 0.98 0.977 0.972 0.954 0.908 0.815 0.676 0.444 0.356 0.294 0.233 0.204 0.191 0.182 0.176 0.172 0.169 0.165 0.16 0.156 0.152 0.149 0.147 0.145 0.143 0.143 0.142
0.976 0.972 0.969 0.961 0.941 0.888 0.785 0.641 0.427 0.352 0.301 0.253 0.23 0.219 0.213 0.208 0.205 0.202 0.198 0.194 0.19 0.186 0.183 0.18 0.178 0.177 0.176 0.175;
0.964 0.959 0.955 0.945 0.921 0.86 0.746 0.601 0.41 0.348 0.308 0.271 0.254 0.246 0.241 0.237 0.234 0.232 0.229 0.225 0.222 0.218 0.215 0.213 0.211 0.209 0.209 0.208;
0.944 0.938 0.934 0.922 0.892 0.82 0.694 0.553 0.391 0.344 0.315 0.288 0.276 0.27 0.266 0.264 0.262 0.26 0.258 0.255 0.252 0.249 0.246 0.244 0.243 0.241 0.241 0.24;
0.91 0.902 0.896 0.88 0.841 0.753 0.622 0.495 0.373 0.341 0.322 0.304 0.296 0.293 0.29 0.288 0.288 0.286 0.285 0.282 0.28 0.278 0.276 0.275 0.274 0.273 0.273 0.272;
0.826 0.815 0.806 0.783 0.732 0.632 0.513 0.423 0.353 0.337 0.327 0.319 0.316 0.314 0.312 0.312 0.312 0.311 0.31 0.309 0.308 0.306 0.305 0.305 0.304 0.303 0.303 0.303;
0.716 0.703 0.691 0.666 0.613 0.525 0.435 0.381 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343 0.343];
pr=[0.38 0.37 0.358 0.346 0.339 0.334 0.329 0.325 0.32 0.314 0.312;
0.426 0.405 0.38 0.357 0.344 0.334 0.325 0.318 0.309 0.296 0.294;
0.474 0.439 0.401 0.367 0.349 0.334 0.322 0.312 0.299 0.281 0.277;
0.522 0.472 0.42 0.376 0.353 0.335 0.319 0.306 0.29 0.267 0.262;
0.57 0.505 0.438 0.385 0.357 0.335 0.316 0.301 0.281 0.255 0.249;
0.694 0.582 0.478 0.402 0.364 0.336 0.311 0.291 0.264 0.228 0.221;
0.826 0.654 0.513 0.417 0.37 0.336 0.306 0.282 0.25 0.207 0.198;
0.96 0.723 0.543 0.429 0.375 0.336 0.303 0.275 0.239 0.19 0.179;
1.106 0.788 0.569 0.439 0.379 0.336 0.3 0.27 0.229 0.175 0.163;
0.386 0.386 0.386 0.386 0.386 0.337 0.295 0.261 0.215 0.152 0.138;
0.391 0.391 0.391 0.391 0.391 0.337 0.292 0.254 0.203 0.135 0.12;
0.245 0.245 0.245 0.245 0.245 0.245 0.245 0.245 0.188 0.111 0.094;
0.174 0.174 0.174 0.174 0.174 0.174 0.174 0.174 0.174 0.088 0.069;
0.153 0.153 0.153 0.153 0.153 0.153 0.153 0.153 0.153 0.056 0.033;
0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.035 0.012;
0.024 0.024 0.024 0.024 0.024 0.024 0.024 0.024 0.024 0.024 0];
pg=[0.31 0.304 0.298 0.295 0.295 0.297 0.3 0.303 0.306 0.311 0.312;
0.287 0.276 0.265 0.26 0.262 0.266 0.271 0.276 0.282 0.292 0.294;
0.263 0.248 0.235 0.229 0.232 0.238 0.245 0.252 0.261 0.275 0.277;
0.239 0.222 0.207 0.202 0.206 0.213 0.222 0.23 0.242 0.259 0.262;
0.215 0.195 0.18 0.176 0.181 0.191 0.202 0.211 0.225 0.245 0.249;
0.153 0.133 0.121 0.122 0.131 0.144 0.158 0.171 0.189 0.216 0.221;
0.087 0.075 0.07 0.078 0.091 0.107 0.124 0.138 0.161 0.192 0.198;
0.02 0.02 0.027 0.041 0.058 0.077 0.096 0.113 0.137 0.172 0.179;
-0.053 -0.032 -0.012 0.01 0.03 0.051 0.072 0.091 0.117 0.155 0.163;;
-0.014 -0.014 -0.014 -0.014 -0.014 0.011 0.035 0.057 0.087 0.13 0.138;
-0.046 -0.046 -0.046 -0.046 -0.046 -0.018 0.008 0.031 0.063 0.111 0.12;
-0.004 -0.004 -0.004 -0.004 -0.004 -0.004 -0.004 -0.004 0.032 0.084 0.094;
0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.057 0.069;
-0.041 -0.041 -0.041 -0.041 -0.041 -0.041 -0.041 -0.041 -0.041 0.021 0.033;
-0.002 -0.002 -0.002 -0.002 -0.002 -0.002 -0.002 -0.002 -0.002 -0.002 0.012;
-0.014 -0.014 -0.014 -0.014 -0.014 -0.014 -0.014 -0.014 -0.014 -0.014 0];
pb=[0.31 0.326 0.344 0.359 0.366 0.369 0.371 0.372 0.374 0.375 0.376;
0.287 0.319 0.355 0.383 0.394 0.4 0.404 0.406 0.409 0.412 0.412;
0.263 0.313 0.364 0.404 0.419 0.428 0.433 0.436 0.44 0.444 0.446;
0.239 0.306 0.373 0.422 0.441 0.452 0.459 0.464 0.468 0.474 0.476;
0.215 0.3 0.382 0.439 0.462 0.474 0.482 0.488 0.494 0.5 0.502;
0.153 0.285 0.401 0.476 0.505 0.52 0.158 0.171 0.547 0.556 0.558;
0.087 0.271 0.417 0.505 0.539 0.557 0.57 0.58 0.589 0.601 0.604;
0.02 0.257 0.43 0.53 0.567 0.587 0.601 0.612 0.624 0.638 0.642;
-0.053 0.244 0.443 0.551 0.591 0.613 0.628 0.639 0.654 0.67 0.674;
0.628 0.628 0.628 0.628 0.628 0.652 0.67 0.682 0.698 0.718 0.724;
0.655 0.655 0.655 0.655 0.655 0.681 0.7 0.715 0.734 0.754 0.76;
0.759 0.759 0.759 0.759 0.759 0.759 0.759 0.759 0.78 0.805 0.812;
0.825 0.825 0.825 0.825 0.825 0.825 0.825 0.825 0.825 0.855 0.862;
0.888 0.888 0.888 0.888 0.888 0.888 0.888 0.888 0.888 0.923 0.934;
0.967 0.967 0.967 0.967 0.967 0.967 0.967 0.967 0.967 0.967 0.976;
0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 1];
L1=[410:10:490 495 500:10:670];
L2=[492.5 495 500:10:570 571.8];
plot(g',r','k',pg',pr','m'); hold;
plot([repmat(1/3,1,length(L1)); g(1,:)],[repmat(1/3,1,length(L1)); r(1,:)],':','Color',[0.5 0.5 0.5]);
intersection=(pr(1,:)-pg(1,:))./(1-3*pg(1,:));
plot(repmat([1/3; 0],1,length(L2)),[repmat(1/3,1,length(L2)); intersection],':','Color',[0.75 0.75 0.75]);
text(g(1,:),r(1,:),num2str(L1'));
axis('equal',[0 0.7 0 1]);