IONEX писал(а): По мере написания кода возник еще один вопрос: как предусмотреть пропуск данных (в векторе данных содержится 0), весовых матриц нет, как например это есть при МНК оценке?

что вы имеете ввиду под пропуском данных?

Хорошо БПФ и Герцель сошлись, это доказывает правильность работы. Непойму заначение результата который они вывели (около 1000), чему соответсвует полученное знание (амплитуда гармоники на частоте f0 ведь равна 1)? Если говорить строго то амплитуда спектра синусоиды равна 0.5 на частоте синусоиды ес...

это в октаве программа? и есть ли источник как получилось это представление окна Дольфа-Чебышева в частотной области? Да в октаве если надо в матлабе замените символы комментария # на % и все будет работать. Если есть октава то наберите в консоли type chebwin и вы увидите ту же функцию (я ее оттуда...

Так записывается окно Чебышева в частотной области. Реально расчет окна Чебышева может быть представлено так: n = 64; at = 80; # beta calculation gamma = 10^(-at/20); beta = cosh(1/(n-1) * acosh(1/gamma)); # freq. scale k = (0:n-1); x = beta*cos(pi*k/n); # Chebyshev window (freq. domain p = cheb(n-1...

Подскажите в чем может быть ошибка? Вот код для матлаба для алгоритма Герцеля думаю вам поможет Fs = 1; f0 = 1/400; N = 3000; t = (0:N-1)/Fs; k = fix(N*f0/Fs); s = cos(2*pi*f0*t); alpha = 2*cos(2*pi*k/N); W = exp(i*2*pi*k/N); a = [1 -alpha 1]; S = filter(1,a,s); S0 = fft(s); SQ = abs(S0(k+1)) SR = ...

Вот эта строчка мне не понятна. grath+=arctan((h[i]*sin(x*i))/(h[i]*cos(x*i))) Если всё привести то получается вот это grath+=x*i Да правильно сначала надо накопить потом арктангенс вычислять. Вот поэтому у меня и было предположение, что осыпаться АЧХ начинает именно из-за недостаточной точности вы...

eiyawii писал(а): верно ли?

не совсем думаю верно если так (без модулей):

на которой чётко видна жирная зеленая линия, которая на всей оси частот равна нулю - это мнимая ЧХ(!), то есть мнимая часть равна нулю. Синяя кривая - это действительная ЧХ, точнее её модуль. Красная кривая - АЧХ. Теперь всё-таки если мнимая часть равна нулю, то она не должна вносить вклад в АЧХ, н...

вы представляете себе что такое 500 дБ? Зачем вам такое подавление? Я думаю что 500 дБ на сегодняшний день это не реально

В вашем случае это неверно, так как период частотной характеристики у вас 2\cdot\pi , в моем случае период чх равен \pi , поэтому мнимая частотная характеристика равна нулю. :o понимаете у вас ких фильтр задан импульсной характеристикой и его чх периодическая с периодом 2*pi. Вы никак не сделаете ч...