Сравнение чисел с порогом

clc, clear all, close all
 
n = 1000; 
x = rand(1,n).*100; %наш рендомный массив
 
[var0 idx0] = find(x<=50); %находим элементы (var0) и их номер (idx0) которые <=50
[var1 idx1] = find(x>50); %находим элементы (var1) и их номер (idx1) которые >50
k = size(idx0); %находим размер первого массива чисел, состоящего из индексов 
k(2) %и выводим его на экран
l = size(idx1); l(2) %аналогично с idx0

clear;
clc;
c=50;
for k=1:100
    x=rand(1,100);
        if x>c
        %подсчитать сколько раз больше
    elseif x<c
        %подсчитать сколько раз меньше
    end
end
%и соответственно вывести на экран
disp(n1);
disp(n2);

clc, clear all, close all
 
c = 50;
n = 10;
x = rand(1,n).*100;
m = x; %чтобы запомнить наш массив, для проверки :-)
for k = 1:n
        if x(k) > c
        x(k) = 1;
        else
        x(k) = 0;
    end
end
r = find(x == 1);
r = size(r);
disp(r(2))
s = find(x == 0);
s = size(s);
disp(s(2))

size(x(x<=50),2)
size(x(x>50),2)

clear;
clc;
N=20;%колличество отсчетов
M=100;%объем выборки
 
for k=1:N
    X1=randn(1,k)*1;
    X2=randn(1,k)*1.1;
end
Dx1=var(X1);%дисперсия процесса X1
Dx2=var(X2);%дисперсия процесса X2
disp(Dx1);
disp(Dx2);
C=log(Dx1/Dx2);%порог
disp(C);
n1=0;
n2=0;
n11=0;
n22=0;
ST1=0;
ST2=0;
 
for K=1:M
    for i=1:N
        X(i)=randn(1,1);
        ST1=ST1+sum(X(i).^2)*((1/(2*Dx1^2))-(1/(2*Dx2^2)));
        
    end
    for ii=1:N
        X(ii)=randn(1,1);
        ST2=ST2+sum(X(ii).^2)*((1/(2*Dx1^2))-(1/(2*Dx2^2)));
        
    end
    if ST1<C
        n1=n1+1;%верно
    elseif ST1>C
        n2=n2+1;%ошибка
    end
    if ST2>C
        n11=n11+1;%верно
    elseif ST2<C
        n22=n22+1;%ошибка
    end
 
end
 
% disp(ST1);
% disp(ST2);
% disp(n1);
% disp(n2);
% disp(n11);
% disp(n22);