1. clear all
2. close all
3. clear clc
4.  
5. %cwiczenie1
6.  
7. M=8;
8. A3=1.099;
9.  
10. N=3;
11. mi=10;
12. sigma=1;
13.  
14. Samples=sigma*randn(N,M)+mi;
15.  
16.  
17. subplot(3,1,1)
18.  
19. for i=1:1:N
20.     plot(Samples(i,:),'.');
21.     xlabel('numer próbki')
22.     ylabel('wartość')
23.     hold on
24. end
25.  
26.  
27.  
28. Average=zeros(M);
29.  
30. for m=1:1:M
31.     for n=1:1:N
32.         Average(m)=Average(m)+Samples(n,m);
33.     end
34. end
35.  
36. for m=1:1:M
37.         Average(m)=Average(m)/N;
38. end
39.  
40. subplot(3,1,2)
41. plot(Average(:,1),'.')
42. xlabel('numer próbki')
43. ylabel('wartość średnia')
44.  
45. Standard=zeros(N,M);
46. Dev=zeros(M);
47. for m=1:1:M
48.     for n=1:1:N
49.         Standard(n,m)=Samples(n,m)-Average(m);
50.     end
51. end
52.  
53. for m=1:1:M
54.     for n=1:1:N
55.         Dev(m)=Dev(m)+(Standard(n,m)*Standard(n,m));
56.     end
57. end
58.  
59. for m=1:1:M
60.         Dev(m)=sqrt(Dev(m));
61. end
62.  
63. subplot(3,1,3)
64. plot(Dev(:,1),'.')
65. xlabel('numer próbki')
66. ylabel('wartość odchylenia')
67.  
68. figure(2)
69. %ćwiczenie 2
70. boxplot(Samples)
71.  
72. %ćwiczenie 3
73. AveAverage=0;
74. for a=1:1:M
75. AveAverage=AveAverage+Average(a);
76. end
77. AveAverage=AveAverage/M;
78.  
79. S=mean(std(Samples));
80.  
81. GLK=AveAverage+A3*S;
82. DLK=AveAverage-A3*S;
83.  
84. %ćwiczenie 4
85.  
86. B3=0;
87. B4=2.568;
88.  
89. DLK2=B3*S;
90. GLK2=B4*S;
91.  
92.  
93.  
94. figure(3)
95. subplot(3,1,1)
96.  
97. for i=1:1:N
98.     plot(Samples(i,:),'.');
99.     xlabel('numer próbki')
100.     ylabel('wartość')
101.     hold on
102. end
103.  
104.  
105.  
106. Average=zeros(M);
107.  
108. for m=1:1:M
109.     for n=1:1:N
110.         Average(m)=Average(m)+Samples(n,m);
111.     end
112. end
113.  
114. for m=1:1:M
115.         Average(m)=Average(m)/N;
116. end
117.  
118. subplot(3,1,2)
119. plot(Average(:,1),'.')
120. xlabel('numer próbki')
121. ylabel('wartość średnia')
122. refline([0 GLK])
123. refline([0 DLK])
124.  
125. Standard=zeros(N,M);
126. Dev=zeros(M);
127. for m=1:1:M
128.     for n=1:1:N
129.         Standard(n,m)=Samples(n,m)-Average(m);
130.     end
131. end
132.  
133. for m=1:1:M
134.     for n=1:1:N
135.         Dev(m)=Dev(m)+(Standard(n,m)*Standard(n,m));
136.     end
137. end
138.  
139. for m=1:1:M
140.         Dev(m)=sqrt(Dev(m));
141. end
142.  
143. subplot(3,1,3)
144. plot(Dev(:,1),'.')
145. xlabel('numer próbki')
146. ylabel('wartość odchylenia')
147. refline([0 GLK2])
148. refline([0 DLK2])
149.  
150.