Составить математическую модель и построить график сигнала

clear all, clc
 
T=10; % длительность
dt = 0.01; %шаг времени
t = 0:dt:T;
df = 1/T;
F = 1/dt;
f = -F:2*df:F;
s = (4+t).*(t>=0 & t<1) + (6-t).*(t>=1 & t<=6) + 0*(t>6 & t<=10);
 
subplot(3,1,1)
plot(t,s,'-r')
ylim([-1 6])
xlabel('t')
title('s(t)')
 
S = fft(s);
 
subplot(3,1,2)
stem(f, fftshift(abs(S)))
xlim([-5 5])
xlabel('f')
title('|S(j \omega)|')
 
subplot(3,1,3)
plot(f, fftshift(angle(S)))
xlim([-5 5])
xlabel('f')
title('arg(S(j \omega))')

'|S(j\omega)|'

'|S(jf)|'

clear all;
clc;
close all;
T    = 100;         % длительность
dt   = 0.01;        %шаг времени
t    = -T:dt:T;
df   = 1/T;
F    = 1/dt;
f    = -F:df:F;
Um   = 5;
tau1 = 1;
tau2 = 6;
s    = 4*heaviside(t)+t./tau1.*heaviside(t)-...
      (t-tau1)./tau1.*heaviside(t-tau1)-...
      5*(t-tau1)./(tau2-tau1).*heaviside(t-tau1)+...
      6*(t-tau2)./tau2.*heaviside(t-tau2);
subplot(3,1,1)
plot(t,s,'-r')
ylim([-1 6])
xlim([-1 7])
xlabel('t')
title('s(t)')
grid on 
S = fft(s);
subplot(3,1,2)
plot(f, fftshift(abs(S)))
xlim([-5 5])
xlabel('f')
title('|S(j \omega)|')
grid on 
subplot(3,1,3)
plot(f, fftshift(angle(S)))
xlim([-0.5 0.5])
xlabel('f')
title('arg(S(j \omega))')
grid on

clear, clc
 
per = 0.05;
T  = 5;
t = 0:per:T-per;
F = 1/per;
df = 1/T;
f = -F+2*df:2*df:F;
w=2*pi/T;
% 
for i=1:1:length(t)
    if t(i) >= 0 && t(i) <= 1
           s(i) = 1.2*t(i)+3.8;
    else if t(i) > 1 && t(i) <= 6
           s(i) = 6 - t(i);
        else if t(i) > 6
                s(i) = 0;
            end
        end
    end
end
N = length(s);
i = sqrt(-1);
% прямое дискретное преобразование фурье
for k = 0:N-1
    S(k+1) = sum( s.*exp(-w*i*(0:N-1)/N*k) );
end
Sa = abs(S);
Sf = angle(S);
figure(1)
plot(f,Sa)
title('Амплитудный спектр')
figure(2)
plot(f,Sf)
title('Фазовый спектр')
grid on
 P = Sa.*exp(i*Sf); 
% обратное преобразование фурье
for k = 0:N-1
    p(k+1) =1/N*sum( P.*exp(w*i*(0:N-1)/N*k) );
end
 
figure(3)
plot(t,s)
grid on
title('Исходный сигнал')
 
figure(4)
plot(t,real(p))
grid on
title('ситезированный сигнал')

clear, clc
 
per = 0.05;
T  = 10;
t = 0:per:T-per;
F = 1/per;
df = 1/T;
f = -F+2*df:2*df:F;
w=2*pi/T;
% 
for i=1:1:length(t)
    if t(i) >= 0 && t(i) <= 1
           s(i) = 1.2*t(i)+3.8;
    else if t(i) > 1 && t(i) <= 6
           s(i) = 6 - t(i);
        else if t(i) > 6
                s(i) = 0;
            end
        end
    end
end
N = length(s);
i = sqrt(-1);
% прямое ПФ
for k = 0:N-1
    S(k+1) = sum( s.*exp(-2*pi*i*(0:N-1)/N*k) );
end
Sa = abs(S);
Sf = angle(S);
figure(1)
plot(f,Sa)
title('Амплитудный спектр')
figure(2)
plot(f,Sf)
title('Фазовый спектр')
grid on
 P = Sa.*exp(i*Sf); 
% обратное преобразование фурье
for k = 0:N-1
    p(k+1) = 1/N*sum( P.*exp(2*pi*i*(0:N-1)/N*k) );
end
 
figure(3)
plot(t,s)
grid on
title('Исходный сигнал')
 
figure(4)
plot(t,real(p))
grid on
title('ситезированный сигнал')
plot(t,s)