цос 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('График модуля ДПФ')