久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

標(biāo)題: 一個對帕里托定律的復(fù)現(xiàn)的C語言小實驗 [打印本頁]
作者: wwwwr    時間: 2020-2-16 17:01
標(biāo)題: 一個對帕里托定律的復(fù)現(xiàn)的C語言小實驗
“帕累托定律”（Pareto law），即只有很少一部分人能獲得最多的資本，并且發(fā)現(xiàn)最成功的人其天賦中庸，但他們卻非常幸運。而這樣的結(jié)果也與現(xiàn)實世界驚人的相似。對該定律做了一個小實驗，進行了復(fù)現(xiàn)。在本文設(shè)計的模型中，假設(shè)共有1000個個體，他們的天賦服從均值為0.6、方差為0.1的正態(tài)分布，并且為每個個體賦予了10的初始資本。將這1000個體隨機的分布在一個正方形區(qū)域內(nèi)，他們被一些事件所包圍。這些事件分為幸運事件和不幸事件，在初始時，它們被隨機地分配到正方形區(qū)域內(nèi)。本實驗對這1000個個體進行40年的模擬，每六個月，這些事件將進行一次移動，并且與個體發(fā)生交互。當(dāng)這些事件與個體相遇時，個體的資本便產(chǎn)生相應(yīng)的變化。其變化規(guī)則為：若是個體與幸運事件相遇，則產(chǎn)生一個（0,1）的隨機數(shù)，若這個隨機數(shù)小于該個體的天賦則將個體的資本翻一倍，否則資本保持不變；若是遇到不幸事件，則將該個體的資本減半。

C語言源程序如下:
  1. ub=1;lb=-1;i=1000;
  2. hold on
  3. for i=1:1000
  4. x=randn(1);y=randn(1);
  5. while x>ub || x<lb || y>ub || y<lb
  6. x=randn(1);y=randn(1);
  7. end
  8. pos(i,:)=[x y];
  9. end
  10. title('隨機人員圖');
  11. figure(1);
  12. plot(pos(:,1),pos(:,2),'o');

  14. ub=1;lb=-1;i=250;
  15. hold on
  16. for i=1:250
  17. a=randn(1);b=randn(1);c=rand
  18. while a>ub || a<lb || b>ub || b<lb
  19. a=randn(1);b=randn(1);
  20. end
  21. pos1(i,:)=[a b];
  22. end
  23. plot(pos1(:,1),pos1(:,2),'*r')

  25. ub=1;lb=-1;i=250;
  26. hold on
  27. for i=1:250
  28. a=randn(1);b=randn(1);c=rand
  29. while a>ub || a<lb || b>ub || b<lb
  30. a=randn(1);b=randn(1);
  31. end
  32. pos2(i,:)=[a b];
  33. end
  34. plot(pos2(:,1),pos2(:,2),'*g')

  36. legend('人員','幸運事件','不幸事件');

  38. for i=1:1:1000
  39.     GN(i,1)=0;
  40.     BN(i,1)=0;
  41. end

  43. T=normrnd(0.6,0.1,1000);

  45. for i=1:1:1000
  46.     w(i,1)=10;
  47. end

  49. for i=1:1:80
  50.         for j=1:1:1000
  51.             for k=1:1:250
  52.             if sqrt((pos(j,1)-pos1(k,1))*(pos(j,1)-pos1(k,1))+(pos(j,2)-pos1(k,2))*(pos(j,2)-pos1(k,2)))<1/201
  53.                 if rand(1)<T(j,1)
  54.                 w(j,1)=2*w(j,1);
  55.                 GN(j,1)=GN(j,1)+1;
  56.                 end
  57.             end
  58.             end
  59.             for k=1:1:250
  60.             if sqrt((pos(j,1)-pos2(k,1))*(pos(j,1)-pos2(k,1))+(pos(j,2)-pos2(k,2))*(pos(j,2)-pos2(k,2)))<1/201
  61.                 w(j,1)=0.5*w(j,1);
  62.                 BN(j,1)=BN(j,1)+1;
  63.             end
  64.             end
  65.         end
  66.     for m=1:1:250
  67.         [pos1(m,1),pos1(m,2)]=move(pos1(m,1),pos1(m,2));
  68.     end
  69.     for m=1:1:250
  70.         [pos2(m,1),pos2(m,2)]=move(pos2(m,1),pos2(m,2));
  71.     end
  72. end

  74. figure(2);
  75. i=1:1:1000
  76. plot(i,w(i,1),'b')
  77. title('財富分布圖');
  78. xlabel('人員');
  79. ylabel('財富');

  81. figure(3);
  82. q=1:1:1000
  83. plot(q,T(q,1),'b')
  84. title('天賦分布圖');
  85. xlabel('人員');
  86. ylabel('天賦');

  88. [b,in]=sort(w);
  89. Totalw=sum(w);
  90. Partw=sum(b(801:1000,:),1);
  91. P=Partw/Totalw;
  92. disp('******');
  93. disp(P);
  94. disp('******');

  96. disp(w(in(1000),1)),disp(T(in(1000),1)),disp(GN(in(1000),1)),disp(BN(in(1000),1));disp('******');
  97. disp(w(in(999),1));disp(T(in(999),1));disp(GN(in(999),1));disp(BN(in(999),1));disp('******');
  98. disp(w(in(998),1));disp(T(in(998),1));disp(GN(in(998),1));disp(BN(in(998),1));disp('******');
  99. disp(w(in(997),1));disp(T(in(997),1));disp(GN(in(997),1));disp(BN(in(997),1));disp('******');
  100. disp(w(in(996),1));disp(T(in(996),1));disp(GN(in(996),1));disp(BN(in(996),1));


  103. function  [X,Y]=move(x,y)
  104. r=1.5/201;
  105. theta=randn(1)*pi;
  106. while theta>pi || theta<-pi
  107. theta=randn(1)*pi;
  108. end
  109. X=x+r*cos(theta);
  110. if X>1
  111.     X=x-r*cos(theta);
  112. end
  113. if X<-1
  114.     X=-r*cos(theta)-x;
  115. end
  116. Y=y+r*sin(theta);
  117. if Y>1
  118.     Y=y-r*sin(theta);
  119. end
  120. if Y<-1
  121.     Y=-r*cos(theta)-y;
  122. end
  123. end
復(fù)制代碼

全部資料51hei下載地址：
帕累托定律.zip (1.3 KB, 下載次數(shù): 5)





歡迎光臨 (http://www.zg4o1577.cn/bbs/) Powered by Discuz! X3.1
主站蜘蛛池模板： 国产性色视频 | 成av人电影在线 | 在线视频成人 | 成人免费视屏 | 国产高清视频在线观看播放 | 草草视频在线播放 | 欧美精品一区二区三区在线 | 国产一二三视频在线观看 | 国产激情第一页 | 综合色久 | 精品不卡 | 久久久久国产精品 | 日日干日日操 | 国产精品视频一区二区三区四蜜臂 | 亚洲综合色丁香婷婷六月图片 | 国产欧美一区二区三区免费 | 成年人在线观看视频 | 国产高清视频在线观看 | 天堂综合 | 综合久久久 | 男女那个视频 | 欧美一级欧美三级在线观看 | 天堂中文资源在线 | 最近免费日本视频在线 | av中文字幕在线观看 | 日韩视频免费看 | 成人av免费| 久久看精品 | 人人人人人爽 | 亚洲国产成人av好男人在线观看 | 在线观看不卡av | 欧美精品免费观看二区 | 日韩中文字幕在线播放 | 久久99精品久久久久久秒播九色 | 午夜精品久久久久久久久久久久久 | www.久久精品 | 成人免费一区二区三区牛牛 | 日韩精品免费在线 | 日本三级做a全过程在线观看 | 国产亚洲一区二区三区在线观看 | 免费a大片|