All pastes #754264 Raw Edit

Unnamed

public text v1 · immutable
#754264 ·published 2007-10-29 21:50 UTC
rendered paste body
x = [0.0671 0.7232 0.7298 0.7222 0.4260 0.2549 0.4499 
0.0631 0.7001 0.6392 0.8440 0.3598 0.3269 0.5214 
0.1659 0.6273 0.7240 0.8380 0.3161 0.3847 0.6055 
0.1700 0.7735 0.5819 0.7884 0.4989 0.2999 0.5350 
0.1304 0.5894 0.7385 0.8334 0.4520 0.3191 0.5563 
0.1293 0.6487 0.7133 0.6759 0.3717 0.2506 0.4956 
0.0740 0.7776 0.7098 0.7052 0.3220 0.2323 0.5309 
0.1561 0.7501 0.6696 0.7038 0.3423 0.3615 0.6053 
0.1705 0.7748 0.6041 0.7073 0.4207 0.2638 0.4583 
0.0507 0.7191 0.5527 0.6994 0.4469 0.2933 0.4774 
0.1515 0.7796 0.5683 0.6672 0.3385 0.2644 0.5983 
0.0346 0.7121 0.5959 0.6977 0.4560 0.3410 0.6043 
0.0733 0.7296 0.7308 0.7698 0.4895 0.3811 0.4920 
0.1524 0.6213 0.7334 0.6928 0.4687 0.3538 0.4485 
0.0135 0.6465 0.7363 0.8436 0.3590 0.4229 0.5520 
0.0183 0.7087 0.6207 0.6966 0.3499 0.2465 0.4175 
0.1651 0.5923 0.7177 0.7240 0.4579 0.4248 0.4472 
0.0399 0.6514 0.5480 0.7250 0.3189 0.3899 0.6159 
0.0294 0.7133 0.5686 0.7155 0.3404 0.2758 0.4954 
0.0344 0.6370 0.7276 0.7317 0.4157 0.2497 0.4204 
0.0392 0.7576 0.7179 0.6909 0.5117 0.3736 0.5903 
0.1033 0.6879 0.5996 0.6657 0.3659 0.2684 0.5995 
0.0320 0.5914 0.6161 0.7933 0.4913 0.3499 0.4474 
0.0707 0.6354 0.5438 0.7329 0.4434 0.2482 0.5698 
0.2013 0.7317 0.6985 0.7010 0.3526 0.2838 0.5588 
0.2028 0.7162 0.6513 0.8231 0.4739 0.3648 0.5503 
0.1715 0.7828 0.6395 0.8170 0.4417 0.2924 0.5483 
0.0789 0.7764 0.6397 0.6717 0.4482 0.3341 0.4604 
0.1523 0.6024 0.7349 0.6951 0.3457 0.3588 0.4161 
0.0745 0.6438 0.5510 0.6723 0.3954 0.3945 0.5902 
0.0118 0.7342 0.7117 0.7992 0.5045 0.3624 0.5645 
0.0763 0.7655 0.6116 0.6654 0.4571 0.4228 0.4474 
0.0906 0.6075 0.6321 0.8274 0.3918 0.2817 0.4869 
0.1559 0.7189 0.7301 0.8201 0.4074 0.3571 0.4277 
0.1158 0.6799 0.7147 0.8246 0.4079 0.3885 0.5473 
0.0073 0.6150 0.6411 0.6612 0.3589 0.2967 0.5960 
0.0701 0.6391 0.6287 0.8220 0.3503 0.3328 0.5065 
0.0724 0.6648 0.7194 0.8058 0.4902 0.3226 0.5759 
0.0341 0.6018 0.6171 0.7282 0.4102 0.4251 0.4845 
0.0579 0.7057 0.6469 0.6862 0.4107 0.3304 0.5847 
0.1686 0.7603 0.6641 0.7667 0.4358 0.2518 0.4477 
0.0901 0.7008 0.5959 0.8104 0.3910 0.2374 0.5332 
0.1190 0.6289 0.5596 0.8402 0.3653 0.2720 0.4259 
0.1286 0.6980 0.5613 0.7809 0.4020 0.2444 0.4909 
0.0571 0.7737 0.6680 0.8292 0.4418 0.3908 0.5031 
0.2036 0.6079 0.6476 0.7627 0.3703 0.3053 0.4290 
0.1163 0.7560 0.5712 0.6878 0.3270 0.3960 0.4429 
0.1843 0.6917 0.7348 0.6776 0.3231 0.3072 0.4986 
0.0876 0.6729 0.6175 0.8350 0.4475 0.4235 0.4487 
0.1571 0.6636 0.6330 0.6612 0.4259 0.3056 0.5746 
0.9481 0.4583 0.4329 0.3796 1.0415 0.6633 0.4715 
0.8133 0.5300 0.2537 0.4636 1.1592 0.5593 0.5146 
0.9322 0.5786 0.2935 0.4616 1.1581 0.5319 0.4713 
0.9225 0.4631 0.3501 0.4930 0.9930 0.4884 0.4425 
0.9339 0.5701 0.3483 0.3687 1.0646 0.6689 0.4419 
0.8609 0.5375 0.3122 0.3785 1.1060 0.5326 0.4599 
0.9640 0.5592 0.3928 0.4941 1.1627 0.5602 0.3946 
0.9402 0.5628 0.2869 0.3515 1.0636 0.4986 0.4983 
0.9857 0.4745 0.3602 0.3562 1.0666 0.5633 0.4268 
0.8532 0.5337 0.2721 0.3583 1.1120 0.6328 0.4999 
0.9573 0.4578 0.3140 0.4171 1.1092 0.5135 0.5311 
0.8885 0.5274 0.4518 0.4154 1.1140 0.6097 0.5046 
0.7968 0.5137 0.4418 0.4656 1.1244 0.4804 0.5277 
0.9236 0.4186 0.3587 0.3867 1.0774 0.4915 0.4858 
0.8785 0.4202 0.2656 0.3579 1.0298 0.6686 0.5645 
0.9063 0.5806 0.3911 0.3913 1.1298 0.6667 0.3941 
0.9326 0.4424 0.3822 0.4357 1.1524 0.5282 0.5004 
0.8617 0.4783 0.3465 0.4841 1.0158 0.6232 0.5414 
0.8774 0.4679 0.4436 0.5192 1.1075 0.5409 0.5342 
0.9894 0.4429 0.4524 0.4311 0.9690 0.5382 0.5067 
0.8938 0.3996 0.2919 0.5384 1.0410 0.5077 0.3713 
0.9234 0.5305 0.2581 0.3634 1.1465 0.6326 0.4849 
0.9323 0.5763 0.3629 0.4984 1.0899 0.6775 0.3780 
0.9923 0.4239 0.4353 0.4760 1.0603 0.4895 0.4604 
0.9552 0.4408 0.3918 0.3948 1.1105 0.6239 0.5150 
0.9197 0.5809 0.2833 0.3836 1.1317 0.4915 0.4760 
0.8333 0.5028 0.3112 0.3612 1.0692 0.6566 0.4914 
0.8045 0.5733 0.2866 0.5083 1.0901 0.6493 0.3912 
0.9251 0.4059 0.3475 0.4796 0.9945 0.5726 0.3826 
0.9115 0.4961 0.3038 0.3906 0.9710 0.5568 0.5198 
0.9103 0.4645 0.2919 0.4454 1.1105 0.6533 0.5167 
0.9942 0.5563 0.3500 0.5236 0.9719 0.6072 0.5571 
0.8034 0.4400 0.3882 0.4602 0.9808 0.6557 0.4659 
0.9735 0.5354 0.3839 0.5372 1.0020 0.5070 0.5194 
0.8579 0.4756 0.3207 0.3449 1.0346 0.5292 0.5671 
0.9763 0.4529 0.3084 0.3615 1.0159 0.6601 0.4332 
0.9098 0.4199 0.4004 0.3653 1.0446 0.6523 0.5397 
0.9591 0.4321 0.4323 0.4638 0.9958 0.5121 0.5450 
0.9909 0.5032 0.3252 0.3892 1.0578 0.6084 0.4174 
0.8888 0.4824 0.2678 0.4962 0.9769 0.5128 0.5238 
0.8369 0.4333 0.3321 0.4391 1.0697 0.6048 0.4934 
0.9935 0.4548 0.3080 0.4786 1.0696 0.5352 0.4841 
0.9598 0.5712 0.3417 0.4276 1.1670 0.6026 0.5571 
0.9776 0.5459 0.3907 0.4047 1.0794 0.6072 0.5221 
0.8060 0.4059 0.4220 0.3774 0.9752 0.5423 0.5318 
0.8422 0.5004 0.4376 0.4089 1.0089 0.5527 0.5635 
0.9386 0.5116 0.4093 0.4362 1.0143 0.5434 0.5261 
0.9221 0.5149 0.3603 0.5268 1.1239 0.5452 0.5161 
0.8828 0.5013 0.2757 0.4068 1.0922 0.6587 0.5610 
0.7992 0.4050 0.3334 0.3624 1.0264 0.6155 0.5643 
0.5145 0.8482 0.3052 0.9107 0.3303 0.1535 0.3836 
0.5960 1.0140 0.3213 0.9352 0.2991 0.0687 0.3646 
0.5651 0.9496 0.1961 0.8502 0.3038 0.2159 0.3912 
0.4357 0.9039 0.3051 0.9699 0.2438 0.1480 0.3918 
0.4299 0.8827 0.1323 0.9324 0.3680 0.2317 0.4826 
0.4320 0.9313 0.2366 1.0288 0.3211 0.0761 0.5244 
0.4021 0.9313 0.3185 1.0246 0.3563 0.2071 0.5193 
0.4899 0.8562 0.2550 0.9486 0.3056 0.2218 0.4971 
0.4631 0.9401 0.2880 0.8516 0.4041 0.2037 0.4918 
0.5830 0.8985 0.2260 0.8625 0.4177 0.2029 0.4597 
0.5986 0.9208 0.2176 1.0107 0.2561 0.1869 0.3757 
0.4792 0.8811 0.2583 1.0339 0.2339 0.2305 0.4696 
0.5956 0.8348 0.1406 0.9929 0.3966 0.1193 0.4723 
0.4587 0.8373 0.3067 0.8604 0.3809 0.0922 0.4202 
0.5179 0.9482 0.2382 0.8866 0.3675 0.1155 0.4756 
0.4794 0.8997 0.2139 0.8748 0.4084 0.0715 0.5162 
0.5260 0.8343 0.2582 0.9956 0.3867 0.2288 0.4401 
0.4965 0.9091 0.2569 0.8604 0.3415 0.1090 0.3855 
0.4578 0.9853 0.2139 0.9503 0.4168 0.1430 0.4592 
0.4342 0.9780 0.2047 1.0003 0.3012 0.2220 0.3610 
0.4291 0.9487 0.1503 1.0306 0.3298 0.0998 0.3867 
0.4256 0.8255 0.1779 1.0139 0.3780 0.2061 0.4081 
0.4954 1.0143 0.2883 0.9268 0.4242 0.1066 0.4373 
0.5003 0.9176 0.2071 0.9488 0.3959 0.1089 0.4092 
0.4087 0.8818 0.2944 0.9813 0.4092 0.0835 0.4304 
0.4855 0.8900 0.3144 0.8769 0.3103 0.2606 0.3448 
0.4628 0.8462 0.1548 0.8778 0.2515 0.1736 0.4839 
0.4829 0.9880 0.2926 0.8467 0.3617 0.0896 0.3746 
0.5813 0.8498 0.1731 1.0150 0.3753 0.0832 0.4952 
0.4692 0.8534 0.2396 1.0296 0.3202 0.1752 0.4711 
0.3998 0.9882 0.3095 0.8488 0.2524 0.1420 0.5249 
0.5234 0.9549 0.2229 1.0125 0.3973 0.1756 0.3499 
0.4585 0.9413 0.1682 0.9915 0.3425 0.0984 0.5165 
0.4399 0.8226 0.2694 0.9251 0.3086 0.2478 0.5201 
0.4942 0.8244 0.2956 0.8564 0.2508 0.2065 0.3609 
0.4109 0.8336 0.2762 0.8516 0.3719 0.1114 0.4565 
0.5021 0.9360 0.1891 1.0234 0.2512 0.1257 0.4958 
0.5982 0.8344 0.3262 1.0294 0.4245 0.1798 0.5216 
0.4631 0.8741 0.1438 0.9114 0.2912 0.1477 0.3383 
0.4569 0.8648 0.2448 1.0189 0.2837 0.1855 0.4360 
0.4253 0.8216 0.1968 0.9170 0.4315 0.2551 0.4183 
0.5865 0.9203 0.3258 1.0364 0.3769 0.1967 0.5040 
0.4056 0.8941 0.1504 0.9118 0.2342 0.2303 0.3365 
0.4843 0.8276 0.2670 0.9105 0.3447 0.1521 0.3615 
0.4123 0.9442 0.2574 0.8802 0.3143 0.2084 0.4174 
0.4416 0.9391 0.2991 1.0225 0.4098 0.2597 0.4854 
0.4730 0.8345 0.1534 1.0323 0.3349 0.1094 0.4157 
0.5653 0.8494 0.1959 1.0085 0.2914 0.0831 0.5171 
0.4328 0.8982 0.2288 1.0071 0.3139 0.2262 0.3353 
0.5908 0.9919 0.3289 0.8950 0.3292 0.1321 0.3429 
0.3431 0.7221 1.0572 0.1524 0.9868 0.2813 1.0502 
0.3390 0.7864 0.9952 0.1734 0.8834 0.3176 1.1009 
0.4286 0.7146 1.0069 0.2255 0.8686 0.2492 1.0822 
0.4179 0.7124 1.0328 0.1780 0.8868 0.4083 1.1511 
0.3712 0.7153 1.0322 0.0712 0.8345 0.2477 1.0056 
0.4640 0.7589 1.1105 0.2036 0.9007 0.3888 1.0672 
0.3304 0.7781 0.9941 0.1563 1.0175 0.3170 1.1601 
0.3776 0.7243 1.0791 0.1137 0.9079 0.2840 1.0861 
0.4291 0.7785 1.1504 0.0771 0.9334 0.3210 1.0792 
0.3895 0.8698 1.1176 0.0659 0.9901 0.3815 1.0687 
0.4387 0.6880 1.0639 0.1968 0.8498 0.3352 1.0408 
0.4631 0.7275 1.0077 0.0979 0.9686 0.3230 1.0106 
0.3292 0.8210 1.1469 0.1812 0.9819 0.2649 1.1542 
0.4655 0.8376 1.1035 0.2212 0.9058 0.3553 1.0000 
0.4874 0.7585 1.0217 0.1739 0.9104 0.3916 1.1056 
0.3424 0.7073 1.0637 0.2325 1.0119 0.3706 1.0442 
0.3903 0.6775 1.0154 0.1880 0.9797 0.3037 1.0797 
0.3413 0.6945 1.0726 0.2068 0.8803 0.3167 1.0088 
0.3081 0.7787 1.0952 0.0474 1.0175 0.3231 1.1398 
0.4516 0.7794 1.0527 0.1912 0.8692 0.4122 1.0174 
0.3996 0.8568 1.1244 0.1666 0.8741 0.2994 1.0672 
0.3539 0.8655 1.0592 0.1807 0.9966 0.4209 1.1768 
0.3719 0.8015 1.0292 0.2098 0.9270 0.3851 1.0784 
0.4833 0.6821 1.1429 0.1831 0.9509 0.3130 1.0346 
0.4580 0.8701 0.9791 0.0905 1.0036 0.3589 1.1243 
0.3934 0.6768 0.9975 0.2360 0.9939 0.2265 1.0969 
0.3742 0.8307 1.0498 0.1448 0.9780 0.3419 1.0250 
0.3129 0.8087 1.1516 0.1095 0.8568 0.2372 1.0616 
0.3579 0.7513 1.1426 0.1841 0.9775 0.3472 1.1295 
0.4422 0.7056 1.0720 0.2200 0.8560 0.3874 1.0136 
0.4821 0.7395 0.9532 0.1923 0.8834 0.3667 1.0159 
0.4601 0.8198 1.1262 0.2147 0.9739 0.3575 1.1801 
0.4553 0.8426 1.0801 0.0715 0.9783 0.3699 1.1126 
0.4356 0.7240 1.0501 0.1250 0.9808 0.3763 1.0475 
0.4180 0.7224 1.0753 0.0940 1.0199 0.3757 1.0961 
0.4800 0.7482 0.9690 0.2300 0.9598 0.3467 1.1787 
0.2935 0.7363 1.1452 0.0607 0.9408 0.3174 1.1432 
0.4627 0.7132 1.0446 0.2209 0.8260 0.3757 1.1589 
0.3537 0.7232 1.0510 0.1773 1.0105 0.3115 1.0992 
0.3107 0.7870 1.0393 0.0765 0.8524 0.4104 1.0270 
0.4907 0.8142 1.1444 0.1954 0.8218 0.2887 1.1085 
0.3523 0.8106 1.0213 0.2364 0.8924 0.2275 1.1622 
0.4875 0.7384 1.0279 0.1899 0.9102 0.2436 1.1982 
0.3982 0.7030 0.9784 0.2131 0.9633 0.4233 1.1628 
0.4474 0.8288 1.1331 0.0609 0.8640 0.3938 1.1483 
0.3720 0.7371 1.0470 0.0692 0.9205 0.3296 1.1017 
0.3365 0.8631 1.0557 0.1321 0.9770 0.3704 1.1202 
0.3024 0.6905 1.0793 0.2155 0.9071 0.3653 1.0137 
0.2948 0.7777 1.1260 0.2000 0.9815 0.2628 1.1426 
0.3340 0.8334 0.9866 0.0815 0.9145 0.2385 1.1190 
];

[m,n] = size(x);
mu = sum(x)/m;
sigma = zeros(n,n);
for i = 1:m
        sigma = sigma + ((x(i,:)-mu)')*(x(i,:)-mu);
end
sigma = sigma/m-1

[w,d] = eig(sigma)
w1 = w(:,n);
w2 = w(:, n-1);
w3 = w(:, n-2);
close all
hold on

for i=1:m
    plot3(x(i,:)*w1, x(i,:)*w2, x(i,:)*w3, 'or')
end