цос 2сем / pr1 / prct1

clc
clear
close all

Nb =input("Вставьте номер бригады: "); 		%% Номер бригады
N =  64                                   		%% Период (длина) последовательности 
Fs = 2000*(mod(Nb,5)+1)                		%% Частота дискретизации
T=1/Fs;
A1 = 1+0.01*Nb                    			%% Амплитуды дискретных гармоник 
A2 = 2*A1									    %% Амплитуды дискретных гармоник 
f1 = Fs/8										%% Частоты дискретных гармоник 
f2 = f1*(mod(Nb,2)+2)									%% Частоты дискретных гармоник

%%
e1 = 1e-7;
n = 0:(N-1); % ДИСКРЕТНОЕ НОРМИРОВАННОЕ ВРЕМЯ
k = 0:(N-1); % ДИСКРЕТНАЯ НОРМИРОВАННАЯ ЧАСТОТА
w1 = 2*pi*f1/Fs; w2 = 2*pi*f2/Fs;
x = A1*cos(w1*n+pi/4)+A2*cos(w2*n+pi/16);
X = fft(x);
MOD = (2/N)*abs(X); % АМПЛИТУДНЫЙ СПЕКТР ПЕРИОДИЧЕСКОЙ ПОСЛЕДОВАТЕЛЬНОСТИ
MOD(1) = (1/N)*abs(X(1));
[MODm,m] = fft_e1(MOD,e1);
A1 = MODm(1) 
A2 = MODm(2)
k1 = m(1); k2 = m(2); % ДИСКРЕТНЫЕ НОРМИРОВАННЫЕ ЧАСТОТЫ
f1 = k1*Fs/N 
f2 = k2*Fs/N
PHASE = angle(X); % ФАЗОВЫЙ СПЕКТР ПЕРИОДИЧЕСКОЙ ПОСЛЕДОВАТЕЛЬНОСТИ
for i = 1:N 
    if (abs(X(i)) < 1e-4) 
    PHASE(i)=0; 
    end
end
figure('Name','Amplitude Spectrum','NumberTitle', 'off') 
subplot(2,1,1), stem(k,MOD,'MarkerSize',3,'Linewidth',2), grid 
xlabel('k'), ylabel('1/N|X(k)|') 
title(strcat(['Amplitude Spectrum of the Periodic Sequence N = ',num2str(N)])) 
subplot(2,1,2), stem(k*(Fs/N),MOD,'MarkerSize',3,'Linewidth',2),grid 
xlabel('f (Hz)'), ylabel('1/N|X(f)|') 
title(strcat(['Amplitude Spectrum of the Periodic Sequence N = ',num2str(N)])) 
%%%%%%
figure('Name','Phase Spectrum','NumberTitle', 'off') 
subplot(2,1,1), stem(k, PHASE,'MarkerSize',3,'Linewidth',2), grid 
xlabel('k'), ylabel('arg{X(k)} (rad)') 
title(strcat(['Phase Spectrum of the Periodic Sequence N = ',num2str(N)])) 
subplot(2,1,2), stem(k*(Fs/N),PHASE,'MarkerSize',3,'Linewidth',2) 
grid, xlabel('f (Hz)'), ylabel('arg{X(f)} (rad)')
title(strcat(['Phase Spectrum of the Periodic Sequence N = ',num2str(N)])) 
%%%%
MOD_K = abs(fft(x)); % МОДУЛЬ ДПФ КОНЕЧНОЙ ПОСЛЕДОВАТЕЛЬНОСТИ
figure(100);
subplot(2,1,1), stem(k,MOD_K,'MarkerSize',3,'Linewidth',2), grid 
xlabel('k'), ylabel('|X(k)|') 
title('DFT Modulus of the Finite Sequence') 
subplot(2,1,2), stem(k,PHASE,'MarkerSize',3,'Linewidth',2), grid 
xlabel('k'), ylabel('arg{X}') 
title('DFT Arguments of the Finite Sequence') 
%%%%
x2=A1*cos(w1*n);
x3=A2*cos(w2*n);
for i=1:N
if i<N/2-1
    x1(i)=x2(i);
elseif N/2<i
        x1(i)=x3(i);
end
end
figure(200);
subplot(2,1,1), stem(n,x,'MarkerSize',3,'Linewidth',2), grid 
xlabel('n'), ylabel('x') 
title('График конечной последовательности x(n)') 
subplot(2,1,2), stem(n,x1,'MarkerSize',3,'Linewidth',2), grid 
xlabel('n'), ylabel('x1') 
title('График конечной последовательности x1(n)') 
%%%
j=zeros(1,N);
j(1)=1;
x4=abs(fft(j));
figure(300);
subplot(2,1,1),stem(n,j,'MarkerSize',3,'Linewidth',2), grid 
xlabel('n'), ylabel('j') 
title('График цифрового едичного импульса') 
subplot(2,1,2),stem(n,x4,'MarkerSize',3,'Linewidth',2), grid 
xlabel('n'), ylabel('x4') 
title('График ффт') 
%%%
C=1;
v=pi/8;
Fe=0;
w=pi/2;
fe=0;
m=0.5;
x5=C.*(1+m.*cos(v*2.*n+Fe)).*cos(w*2.*n+fe);
K=fft(x5);
MOD5 = (2/N)*abs(K); % АМПЛИТУДНЫЙ СПЕКТР ПЕРИОДИЧЕСКОЙ ПОСЛЕДОВАТЕЛЬНОСТИ
MOD5(1) = (1/N)*abs(K(1));
figure(400);
subplot(2,1,1), stem(k,MOD5,'MarkerSize',3,'Linewidth',2), grid 
xlabel('k'), ylabel('Mod') 
title('ГРафик амплитудного спектра x5') 
subplot(2,1,2), stem(n,x5,'MarkerSize',3,'Linewidth',2), grid 
xlabel('n'), ylabel('x5') 
title('График конечной последовательности x5(n)') 
%%%
Fs=2000;
T=1/Fs;
nT=-500*(N-1)*T:T:500*T*(N-1);
x6=sinc(pi.*nT);
L=abs(fft(x6));
figure(500);
subplot(2,1,1), stem(nT,x6,'MarkerSize',3,'Linewidth',2), grid 
xlabel('nT'), ylabel('x6') 
title('ГРафик последовательности x6') 
subplot(2,1,2), stem(nT,L,'MarkerSize',3,'Linewidth',2), grid 
xlabel('k'), ylabel('L') 
title('График модуля ДПФ') 
%%%
a=0.0005;
om1=pi/12;
n=-3*(N-1):3*(N-1);
x7=exp(-a.*n.^2).*cos(om1.*n);
L1=abs(fft(x7));
figure(501);
subplot(2,1,1), stem(n,x7,'MarkerSize',3,'Linewidth',2), grid 
xlabel('n'), ylabel('x7') 
title('ГРафик последовательности x7') 
subplot(2,1,2), stem(n,L1,'MarkerSize',3,'Linewidth',2), grid 
xlabel('n'), ylabel('L') 
title('График модуля ДПФ')