Формирование матрицы

T=[0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20]*10^-3;
%амплитуда тока:
I1=4;
%амплитуда и угол напряжения своего конца u1(t)
U1=20; fiu1=pi/2;
%создаем массив тока с изменениями по частоте и углу:
for fi=0:90
    for f=40:60
        for k=1:length(T)
        i_t=I1*sin(2*pi*f*T+fi*pi/180);
        I_t(1,f-39,k)=f;
        I_t(fi+1,1,k)=fi;
        I_t(fi+2,f-38,k)=i_t;
        end
    end
end
I_t

T=[0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20]*10^-3;
%амплитуда тока:
I1=4;
%амплитуда и угол напряжения своего конца u1(t)
U1=20; fiu1=pi/2;
%создаем массив тока с изменениями по частоте и углу:
for fi=0:90
    for f=40:60
        for k=1:length(T)
        i_t=I1*sin(2*pi*f*T+fi*pi/180);
        I_t(1,f-39,1)=f;
        I_t(fi+1,1,1)=fi;
        I_t(1,1,k)=i_t;
        end
    end
end
I_t

T=[0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20]*10^-3;
%амплитуда тока:
I1=4;
%амплитуда и угол напряжения своего конца u1(t)
U1=20; fiu1=pi/2;
%создаем массив тока с изменениями по частоте и углу:
for fi=0:90
    for f=40:60
        for k=1:length(T)
        i_t=I1*sin(2*pi*f*T+fi*pi/180);
        I_t(1,f-39,1)=f;
        I_t(fi+1,1,1)=fi;
        I_t(1,1,k)=T;
        I_t(fi+2,f-38,k+1)=i_t;
        end
    end
end
I_t

clc; close all; clear all;
T=(0:20)*10^-3;
I1=4; % амплитуда тока
U1=20; fiu1=pi/2; % амплитуда и угол напряжения своего конца u1(t)
% Cоздаем массив тока с изменениями по частоте и углу:
fi = 0:90;
f = 40:60;
I_t = cell(length(fi),length(f));
for i = 1:length(fi)
    for j = 1:length(f)
        i_t(i, j) = { I1*sin(2*pi*f(j)*T+fi(i)*pi/180) };
    end
end
plot(T, I_t{10, 5})

clear all
clc
 
t = linspace(0,10,1000); % массив времени
y0 = sin(2*pi*t);
T = 2;
for tau = 0:0.05:10-T
    yb = [-1 1 1 -1 -1];
    xb = [tau,tau,tau+T,tau+T,tau];
    y1 = y0.*((t>tau)&(t<tau+T));
    plot(t,y0,':g',t,y1,'r',xb,yb,'--b')
    axis([0 10 -5 5])
    pause(0.01)
end

M=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 X=[0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20]*10^-3; 
fd=X(2)-X(1); %частота дискретизации 
Td=1/fd; %период дискретизации 
I1=4; %амплитуда тока 
U1=20; 
fiu1=pi/2; %амплитуда и угол напряжения своего конца u1(t) 
t=1; %создаем массив тока с изменениями по частоте и углу: 
for fi=0:90     
  for f=40:60         
    for k=3:(length(X)-2)             
      i_X=I1*sin(2*pi*f*X+fi*pi/180);                 
      Y=i_X((k-2):(k+2));     
      diff_eq=Y(1)+Y(2)+Y(4)+Y(5);             
      C(1,t)=fi;             
      C(2,t)=f;             
      C(3,t)=diff_eq;             
      t=t+1;         
   end     
  end 
end

t = 0:0.01:5;
i = 1;
for fi=0:90
  for f = 20:40
     Y(i,:) = sin(2*pi*f*t + fi*pi/180);
     i = i+1;
  end
end

clear, clc
tmax = 10;
t = 0:0.01:tmax;
f0 = 20; % начальная частота
df = 20; % размах частоты
fi0 = 0; % начальная фаза
dfi = pi/2; 
y = sin( 2*pi*(f0+df/tmax*t).*t + (fi0+dfi/tmax*t) );

% первым делом формируем матрицу Y!!!
f = inline('sqrt(sum(x.^2))'); % функция обработки отсчетов
% в данном случае считается корень суммы квадратов
% важно, чтобы результат был числом!
n = 7; % размер окна
Z = zeros(size(Y)); % инициализируем массив
for i = 1:size(Y,1)
  j = 0;
  while j+n <= size(Y,2) % двигаем, пока окно не достигнет края
    x = Y(i, j+1:j+n); % сохраняем n отсчетов окна
    Z(i, j+fix(n/2)) = f(x); % считаем функцию и сохраняем на позиции центрального отсчета
    % можно сразу:
    % Z(i, j+fix(n/2)) = f( Y(i, j+1:j+n) );
    j = j+1; % сдвигаем окошко
  end
end
% рисуем графики
plot(t,Y(1,:),'b', t, Z(1,:),'r')