doc_num = 3;
data1=[1651 1 2349 1 3754 1 4142 1 4727 3 7191 2 7221 1 7584 2 7603 1 7702 1 8326 19 8328 2 8675 2 8948 1 9161 1 9615 10 10087 3 10156 1 10365 2 10632 23 10704 1 10771 1 10790 1 11053 1 11220 1 11656 1 12655 5 12871 1 13458 1 13646 1 13921 1 13922 1 13927 1 14147 1 14762 1 15759 1 15851 1 15951 1 16155 1 16210 1 16539 1 16598 1 17071 1 17700 1 18524 3 18561 1 19078 2 19848 1 19864 1 20037 1 20260 1 20494 1 21806 1 21818 1 22463 1 22693 1 22783 1 22789 1 22791 1 22943 8 23170 1 23509 1 23722 1 24350 1 24819 1 24968 1 25772 1 25917 1 25944 1 26052 1 26229 1 26231 1 26284 1 26344 2 26389 1 26616 2 26681 1 26852 2 27087 3 27107 3 27261 1 28014 3 29732 1 30061 3 30403 3 30835 1 31565 1 31956 1 32326 2 32471 4 32475 1 33061 1 35677 1 35691 1 36007 1 36055 1 36232 4 36951 5 37562 1 37965 1 38136 2 39192 2 39514 2 39522 1 40333 1 40452 1 40466 1 40953 1 41036 1 41365 1 42113 1 42379 1 42736 1 43153 1 43313 2 43469 1 43486 1 43515 1 43742 2 43747 1 44463 1 44857 1 45098 3 45115 1 45587 2 45698 1 46332 1 46816 1 46840 2 46950 3 46981 1 47331 3 47476 1 47719 1 47722 1 47732 1 48197 1 49691 1 49924 1 49976 2 50050 1 50903 1 51329 1 51522 1 51740 1 51759 3 51829 1 51920 1 52285 1 52329 1 52347 5 52559 1 52572 1 52859 1 53922 3 53928 1 53955 1 54147 1 54645 3 55041 1 55982 1 56193 1 56257 1 56318 1 56686 1 56792 1 57014 1 57069 2 57083 1 57220 1 57252 2 58193 1 58845 1 59301 1 59375 5 59394 1 60012 1 60851 2 60945 3 60955 1 61107 3 61254 4 61311 1 61519 1 61940 1 61977 1 62297 1 62520 1 62727 1 62825 1 63348 2 63406 1 63576 1 63797 2 63926 1 64078 1 64186 2 64547 1 64823 4 65072 1 66417 1 66439 2 66440 1 67146 1 67231 2 67352 1 67426 1 67679 3 68753 1 68976 1 69040 1 69109 2 69277 1 69297 1 69346 1 69475 1 69509 1 69611 1 69662 2 70079 1 70084 1 70246 1 ];
data2=[1070 1 1084 2 1611 2 1651 2 1730 1 2162 1 2194 1 2242 1 2294 1 2349 1 2539 2 2540 5 3120 1 3167 5 3350 1 3395 1 3580 1 3608 1 3744 1 3897 2 4142 4 4513 1 4541 6 4775 2 5122 1 5141 1 5250 1 5275 2 5285 1 5291 1 5613 2 6022 3 6027 1 6030 1 6223 1 6291 3 6627 1 7198 1 7204 1 7221 1 7279 1 7356 3 7403 1 7424 3 7425 1 7576 1 7758 1 7779 1 7809 1 7943 1 7974 2 8078 1 8104 1 8148 1 8324 1 8326 6 8328 4 8592 3 8620 1 8639 2 8675 6 8680 2 8702 1 8788 1 8987 1 9161 1 9208 4 9430 1 9454 1 9605 2 9824 1 9872 5 9919 2 9978 1 10087 5 10137 1 10183 1 10188 1 10365 1 10377 1 10588 1 10597 1 10618 1 10632 74 10704 6 10712 1 10771 4 10985 1 11146 1 11315 1 11418 1 12089 1 12251 1 12331 1 12351 2 12530 1 12558 1 12655 3 12759 3 12771 2 12927 4 12929 1 13034 6 13039 2 13040 1 13131 3 13297 1 13309 1 13458 5 13534 1 13921 6 13957 1 13961 18 14253 2 14326 1 14358 1 14446 1 14610 2 14633 1 14958 6 15271 1 15675 1 15851 5 15935 1 16296 2 16506 1 16577 3 16928 4 17018 1 17034 1 17071 3 17088 1 17089 2 17200 1 17324 1 17339 1 17443 4 17620 1 17638 1 17871 1 17937 3 18034 2 18263 1 18463 1 18500 1 18568 3 19078 1 19145 1 19307 6 19761 1 19848 1 19853 2 19917 1 20037 8 20178 1 20236 1 20494 1 20495 3 20707 1 20910 2 20948 1 20950 1 21777 2 21813 1 21814 1 22321 1 22323 1 22376 1 22781 6 22796 2 23018 1 23066 2 23180 1 23406 1 23410 2 23447 2 23475 1 23570 3 23576 1 23589 1 23728 2 23872 1 23977 1 24232 1 24445 1 24478 1 24567 1 24600 1 24731 2 24819 1 24992 2 25090 2 25353 1 25391 3 25394 1 25816 1 25883 2 25896 1 25917 3 25933 1 25944 6 25979 1 26161 1 26230 1 26284 5 26677 1 26863 1 27087 3 27177 1 27261 3 27299 1 27307 1 27422 1 28052 1 28200 3 28516 1 28526 1 28636 1 28648 6 28722 1 29031 2 29417 1 29720 1 29731 1 29738 2 30061 4 30099 2 30760 1 31807 2 31847 1 31857 1 31999 1 32021 1 32053 1 32307 1 32326 3 32471 11 32475 2 32975 1 33003 1 33045 1 33479 1 33520 1 33527 1 33948 1 35305 1 35479 1 36007 1 36030 1 36181 2 36630 2 36901 1 37378 2 37562 1 37965 1 38082 3 38187 4 38197 1 38378 4 38553 2 39359 1 40659 1 40781 1 41023 24 41099 2 41162 1 41331 1 41636 1 41898 1 42110 1 42342 2 42379 1 42672 1 42939 1 43010 1 43203 2 43266 1 43469 1 43600 3 43610 1 43613 62 43729 1 43742 1 44244 2 45005 1 45036 1 45148 8 45196 2 45574 1 45623 25 45649 1 45663 2 46066 5 46076 1 46647 1 46888 1 46947 1 47011 1 47129 1 47147 2 47151 2 47182 1 47258 1 47331 14 47350 4 47357 1 47358 1 47363 1 47394 5 47399 1 47690 1 47771 1 47791 1 48362 1 48622 4 49128 2 49547 1 49639 1 49647 1 49648 1 49942 3 49945 1 49958 1 50073 1 50223 1 50684 1 50708 1 50901 2 50903 1 51329 1 51584 6 51631 1 51645 1 51921 1 51946 2 51963 6 52012 1 52347 1 53518 6 53922 8 53968 1 54027 1 54080 1 54267 1 54361 1 54560 1 54607 1 54645 4 54687 1 54924 1 55393 5 55566 3 56257 1 56447 3 56462 1 56483 1 56780 4 56793 3 56797 2 56861 1 57044 7 57069 2 57163 1 57252 2 57404 1 57618 1 57622 2 57941 3 58277 1 58315 2 58338 2 58477 1 58547 1 58842 4 58843 1 58845 1 59625 1 60026 1 60851 1 60871 1 60945 3 60950 1 60965 2 61107 6 61146 2 61197 1 61218 2 61313 1 61519 1 61608 1 61940 2 61962 5 61978 2 61995 2 62136 1 62278 3 62281 3 62449 7 62461 1 62743 1 62790 1 62825 1 62872 1 63005 1 63025 1 63103 1 63179 3 63348 7 63376 1 63435 1 63457 4 63461 1 63547 2 63576 1 63597 1 63797 5 63926 1 64078 5 64186 1 64684 3 64737 5 64791 1 64797 2 64823 10 64930 3 65017 1 65108 1 65427 1 65504 3 65559 1 65654 4 66412 1 66441 1 66444 2 66528 1 66552 1 67109 8 67187 1 67421 2 67428 1 67433 1 67436 1 67510 1 67679 3 67911 5 68107 1 68439 1 68708 1 68820 1 68852 2 69056 1 69075 1 69106 1 69124 1 69145 1 69180 1 69320 3 69495 1 69496 1 69508 1 69509 9 69631 1 69650 1 69796 1 69940 1 70079 1 70083 1 70084 9 70085 44 70087 2 70230 2 70257 1 70424 1 72502 1];
data3=[20228 1 20236 2 20260 1 20494 5 1651 4 20496 1 21314 1 21604 1 21806 3 21814 1 21818 1 22322 1 22463 4 22647 1 22802 1 22885 1 22934 1 23086 1 23164 2 23385 3 23486 2 23509 1 23727 1 23872 1 23977 2 24281 1 24567 1 24697 1 24703 2 24708 2 24819 1 24953 6 25399 1 25453 1 25587 1 26154 1 26275 1 26392 1 26568 6 26677 2 26862 1 27087 4 27101 3 27167 1 27232 2 27422 1 27457 1 28371 1 28396 1 28404 1 28722 4 29029 1 29086 1 29245 1 29249 1 29731 1 29738 1 30061 4 30240 1 30306 1 30326 1 30663 1 30998 1 31565 1 31847 1 31868 1 31891 1 31931 3 32046 2 32326 1 32471 9 33002 1 33045 1 33048 1 33105 4 33479 1 33488 4 33680 1 34736 1 35195 1 35497 1 35677 2 35984 1 36007 2];
%关键求出d1中每个元素在所有文件中出现的数目并得出最后总和,求p值时候会用到
%只要找到一个,包含它的数量就加1,总数会在求q时候用到
% m = 1;
% eval(['d',num2str(m)]) = zeros(max(size(eval(['data',num2str(m)]))/2),1);
d1 = zeros(max(size(data1)/2),1); %初始化d1矩阵
d1_doc = zeros(max(size(data1)/2),1); %初始化d1包含文档总数的矩阵
for j=1:2:max(size(data1))
for i=1:doc_num
n=find(eval(['data',num2str(i)])-data1(j)==0);
if n %说明找到了
d1((j+1)/2) = d1((j+1)/2) + eval(['data',num2str(i),'(n+1)']);
d1_doc((j+1)/2) = 1 + d1_doc((j+1)/2);
end
end
d1_ben((j+1)/2) = data1(j+1); %d1在自己文档中的矩阵列
end
p1 = d1_ben / sum(d1); %p1的值
q1 = log((1+doc_num)./d1_ben);
p_q1 = p1.*q1;