Есть синусоидальный сигнал из суммы двух синусоид. Необходимо построить фильтр нижних и верхних частот и их АЧХ. Думаю использовать фильтр Баттерворда. Помогите, пожалуйста, очень мало опыта работы в данной сфере
((Вот такой код у меня:
Код: Выделить всё
close all
Fs = 1000; %количество точек общая
dt=1/Fs; % шаг общая
L = 2000; % длина сигнала
t = (0:L-1)*dt; %длительность общая
Fc1 = 4; %частота первой синусоиды
a1 = 1; %амплитуда первой синусоиды
sin1 = a1*sin(2*pi*Fc1*t); %график первой синусоиды
Fc2 = 16; %частота второй синусоиды
a2 = 2; %амплитуда второй синусоиды
sin2 = a2*sin(2*pi*Fc2*t); %график второй синусоиды
signal = (sin1+sin2); % сумма синусоид
Y = fft(signal); % преобразование Фурье (сигнала)
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1); % спектр сигнала (односторонний)
f = Fs*(0:(L/2))/L; % частотный диапазон