Advertisement
Guest User

Untitled

a guest
May 13th, 2019
161
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.00 KB | None | 0 0
  1. testw = [0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
  2. 0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
  3. 0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
  4. 0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
  5. 0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
  6. 0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
  7. 0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
  8. 0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
  9. 0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ...
  10. 0.3568 0.452 0.494 0.2982 0.3892 0.675 0.576 0.489 0.745 0.225 ];
  11. %testw = breite;
  12.  
  13. %Datenvorbereitung
  14. gerundet = round(testw,1);
  15.  
  16. maximum = max(gerundet),1;
  17. anzahl = ceil(maximum / 0.1);
  18.  
  19. werte = zeros([1 anzahl]);
  20. sz = size(werte);
  21. haeufigkeit = zeros(sz);
  22.  
  23. for i= 1:anzahl
  24. werte(i) = 0.1*i;
  25. end
  26.  
  27. wertebereich = zeros(sz);
  28.  
  29. for i = 1:anzahl
  30. wertebereich(i) = 0.1;
  31. end
  32.  
  33. %Häufigkeitsverteilung
  34. for i = 1:numel(gerundet)
  35. if gerundet(i) <= 0.2
  36. haeufigkeit(1) = haeufigkeit(1) + 1
  37. else
  38. for j = 2:anzahl
  39. if gerundet(i) >= werte(j-1) && gerundet(i) < werte(j)
  40. haeufigkeit(j) = haeufigkeit(j) + 1
  41. end
  42. end
  43. end
  44. end
  45.  
  46. %Grauwerte
  47. prozentual = zeros(sz);
  48.  
  49. for i = 1:anzahl
  50. prozentual(i) = haeufigkeit(i)/sum(haeufigkeit)
  51. end
  52.  
  53. grauwerte = zeros(sz);
  54.  
  55. for i = 1:anzahl
  56. grauwerte(i) = round(1-1*prozentual(i),3);
  57. end
  58.  
  59. colormap gray
  60. b = bar([wertebereich; nan(size(wertebereich))],'stacked');
  61. % colorbar('Direction','reverse')
  62. colorbar('Ticks',[0, 0.2, 0.4, 0.6 ,0.8, 1],...
  63. 'TickLabels',{'100%','80%','60%','40%','20%', '0%'})
  64.  
  65. for i = 1:anzahl
  66. b(1,i).FaceColor = [grauwerte(i) grauwerte(i) grauwerte(i)];
  67. b(1,i).EdgeColor = 'none';
  68. end
  69.  
  70. hold on
  71. med = median(gerundet);
  72. plot([.5 1.5], [med med], '-', 'color', [.2 .6 1], 'linewidth', 2)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement