Jump to content

    

C2000

Свой
  • Content Count

    515
  • Joined

  • Last visited

Community Reputation

0 Обычный

About C2000

  • Rank
    Знающий
  • Birthday 01/12/1984

Контакты

  • Сайт
    Array
  • ICQ
    Array

Информация

  • Город
    Array

Recent Profile Visitors

3510 profile views
  1. Добрый день! Существуют ли библиотеки аналогичные библиотекам питона и матлаба? Например для синтеза фильтров, расчёта эллиптических функций якоби, лагорифмов, работы с комплексными числами, numpy.convolve и т.д.
  2. Спасибо еще раз. Разобрался всё сходиться!!!
  3. Спасибо! Теперь все вроде бы ясно, только вот результат не сходится ((
  4. Можно подробнее, пожалуйста. В каких именно скобках?
  5. Добрый день! Есть формула: С вычислением B0 и B1 все более менее понятно сумма комплексных чисел, где модуль комплексного числа вычисляется как биномиальный коэффициент Но дальше непонятно что такое B0* и D* (сопряженный числа)? ЧТО ОБОЗНАЧАЮТ ФИГУРНЫЕ СКОБКИ ln{D/D*}???? Что значит логарифм комплексного числа? И как из комплексных чисел в итоге получается угол, который совсем не комплексный? Автор приводит пример: для wp = 0.3pi и N = 6 получается Фи0 = 1.536422 Но как это значение высчитывается подробно не описано. Прошу гуру математики помочь разобраться с этими формулами. Есть ли у кого какие предположения как рассчитать B0* и D* и затем Фи0? Дополню: Вот эта формула: Получается как то из вот этой: Возможно это даст какое-то понимание)
  6. Добрый день! Есть методика расчёта, в ней содержится такая вот формула: Подскажите, пожалуйста, как это рассчитать в Phyton или в Matlab или еще где-то
  7. А после этого можно этот комплексный сигнал повернуть на нужный угол и восстановить, т.е. на выходе тот же массив выборок АЦП но сдвинутый на определенный угол во всём диапазоне частот? Как бы аналог преобразования Гильберта.
  8. Как из питоновского скрипта запустить компиляцию Си файла (запуск gcc + заранее созданый make файл) для создания dll, которая в дальнейшем будет использована в этом же скрипте?
  9. Можно подробнее, если возможно с примером или ссылкой на подробное описание
  10. А если подать единичный сигнал, получить импульсную характеристику и из нее уже ачх и фчх?
  11. Подскажите, как можно получить (в matlab) АЧХ и ФЧХ фильтра реализованного в Си функции (как подключить Си функцию в скрипт известно). Си функция принимает выходные выборки фильтрует и на выходе соответственно отфильтрованный сигнал. Первая мысль в лоб, это прогонять через фильтр сигналы на всех частотах и анализировать что на выходе, но возможно есть более правильный способ?
  12. Добрый день! Подскажите, есть ли аналог ftool для phyton, для построения АЧХ, ФЧХ, импульсных характеристик? И что предпочтительнее в качестве бесплатной замены для Matlab, python или octave?
  13. В принципе сумма в приведенном примере сходиться уже после 5 элементов, но результат почему-то некорректный. Пробую рассчитать Elliptic Halfband Filter, по методике указанной в этой статье (страница 3, начиная с пункта (i)): IEEE_Elliptic filter design for a class of generalized halfband filters.pdf Беру за пример для расчёта Example II (стр.4) из этой же статьи, но коэффициенты никак не хотят совпадать с теми что автор рассчитал, и как следствие желаемые характеристики фильтра не получаются. С коэффициентами из статьи все гуд. Прошу помощи у тех кто имеет больший опыт с матлабом и опыт в дизайне фильтров. Возможно в статье ошибка? Например в формуле (36) по-моему должна быть сумма для m от 1 до inf, а не для i. Вот мой скрипт, гляньте пожалуйста, что не так: % (i) We are given input data Ap = 0.001; % is the maximum permissible passband magnitude deviation As = 0.001; % is the maximum permissible stopband magnitude deviation Wp = 0.48*pi; % frequency pass Ws = pi - Wp; % frequency stop L = 6; % % (ii) Compute A as in (20). A = min(As, (2*Ap - Ap^2)^0.5); % the required stopband magnitude deviation for the analog filter %(iii) Compute k k = tan(Wp/2)^2; k_ = (1-k^2)^0.5; %(iv) If p = 0.5 * (1 - k_^0.5) / (1 + k_^0.5); q = p + 2*(p^5) + 15*(p^9) + 150*(p^13); % (iv) From (27) and (31), obtain N = 2L + 1; N = 2*L + 1; % (vi) Compute by approximating the expansion of (23): ai = zeros(L,1); for i = 1:L syms m f = ((-1)^m) * ( q^(m*(m+1)) ) * sin( ((2*m + 1)*pi*i) / N ); sum_1 = symsum(f, m, 0, 50); syms m f = ((-1)^m)* (q^(m^2)) * cos(2*m*pi*i/N); sum_2 = symsum(f, m, 1, 50); Om = (2*(q^0.25) * sum_1) / (1 + sum_2); % (vii) Compute cos using (25) and (26) ri = abs( ( (1 - k * (Om^2))*(1 - (Om^2)/k) ) )^0.5; cosO = ( ((-1)^(i+1)) * ri ) / (1 + Om^2); ai_ = (1 - cosO) / (1 + cosO); format long %выводить много цифр после запятой if (ai_ <= 1) ai(i) = ai_; else ai(i) = 1/ai_; end; end ai
  14. Подскажите как в Matlab посчитать вот это: При использовании symsum с аргументов inf Матлаб виснет.