Лабы МОИИ / 4 / 8587 / 8587 / 6лаба / 6

for j=1:100
t=70:250;
y=square(2*pi/140*t);
y1=(y+1)/2;
t1=t-70;
plot(t1,y1);
n=0.3*randn(1,181);
for i=1:181
f1(i)=y1(i)+n(i);
end
hold on;
plot(t1,f1);
Up=1.12;
plot(t1,Up,'r')
A = f1
B=find(A > Up)
t1=B(1)
t2=max(B)
T=t2-t1
AM=A(B)
AMsr=mean(AM)
Pogr(j)=AMsr-1
Mpogr=mean(Pogr)
Dpogr=var(Pogr)
SKOpogr=std(Pogr)
end