Xlabel("Частота");

ylabel("Модуль Фурье-образа сигнала");

grid on;

figure;

R1=fft(imp);

R2=phase(imp);

subplot(2,1,1),area(t,abs(R1));

title("Амплитудный спектр прямоугольного сигнала");

Xlabel("Частота");

ylabel("Модуль Фурье-образа сигнала");

grid on;

subplot(2,1,2),plot(t,R2);

title("Фазовый спектр прямоугольного сигнала");

Xlabel("Частота");

ylabel("Фаза");

% % % Задание 6

A = [1 1;1 -1];

B = kron(A,A);

C = kron(B,B);

D = kron(C,B);

E = kron(D,A);

R1=fwht(imp);

walshSpectr= fwht(imp);

FurieSpectr=abs(fft(imp));

walshSpectr=walshSpectr(:,1:101);

subplot(2,1,2),area(t,FurieSpectr),grid on;

title("FurieSpectr");

subplot(2,1,1),area(t,walshSpectr), grid on;

title("WalshSpectr");

% % % Задание 7

fs = 1024;

t=0:1/fs:1-1/fs;

for i = 1:100

signal = cos(2*pi*(10+i/10)*t);

sigHann = signal.*hann(1024)';

[ch1(i),ch2(i),ch3(i)]=freqchr(signal);

[chHann1(i),chHann2(i),chHann3(i)]=freqchr(sigHann);

end

figure

hold on

plot(ch1,'k')

plot(ch2,'r')

title('Частотная характеристика 2-х соседних каналов');

Xlabel("Частота");

ylabel("Амплитуда");

grid on;

hold off

figure

plot(chHann1,'k')

hold on

plot(chHann2,'r')

grid on;

title('Частотная характеристика 2-х соседних каналов (окно Ханнинга)')

Xlabel("Частота");

ylabel("Амплитуда");

%%Задание 8

fs = 1024;

t=0:1/fs:1-1/fs;

for i = 1:1000

signal = cos(2*pi*(10+i/100)*t);

sigHann = signal.*hann(1024)';

sigBlackman = signal.*blackman(1024)';

sigHamming = signal.*hamming(1024)';

[ch1(i),ch2(i),ch3(i)]=freqchr(signal);

[chHann1(i),chHann2(i),chHann3(i)]=freqchr(sigHann);

[chBlack1(i),chBlack2(i),chBlack3(i)]=freqchr(sigBlackman);

[chHamming1(i),chHamming2(i),chHamming3(i)]=freqchr(sigHamming);

end

figure

plot(ch1,'k')

hold on

plot(ch2,'r')

plot(ch3)

grid on;

title('Частотная характеристика 3-х соседних каналов');

xlabel("Частота");

ylabel("Амплитуда");

%%Задание 9

figure

subplot(3,1,1)

hold on

plot(chHann1,'k')

plot(chHann2,'r')

plot(chHann3)

title('Частотная характеристика 3-х соседних каналов (окно Ханнинга)');

xlabel("Частота");

ylabel("Амплитуда");

grid on;

hold off

subplot(3,1,2)

hold on

plot(chHamming1,'k')

plot(chHamming2,'r')

plot(chHamming3)

title('Частотная характеристика 3-х соседних каналов (окно Хэмминга)');

xlabel("Частота");

ylabel("Амплитуда");

grid on;

hold off

subplot(3,1,3)

hold on

plot(chBlack1,'k')

plot(chBlack2,'r')

plot(chBlack3)

title('Частотная характеристика 3-х соседних каналов (окно Блэкмена)');

xlabel("Частота");

ylabel("Амплитуда");

grid on;

hold off

function [chan14, chan15, chan16] = freqchr(signal)

f=abs(fft(signal));

chan14=f(14)/512;

chan15=f(15)/512;

chan16=f(16)/512;

end

28