реклама на сайте
подробности

 
 
> Как задать значения ф-ций, Начал изучать матлаб и тут понеслось...
Грендайзер
сообщение Aug 31 2016, 12:20
Сообщение #1


Местный
***

Группа: Участник
Сообщений: 368
Регистрация: 18-04-11
Из: Город-герой Москва
Пользователь №: 64 451



Здравствуйте. Начал изучать матлаб и появилось несколько вопросов на которые никак не найду ответов. Вопросы глупые конечно, но практики маловато так что...
1) Имею значения коэффициентов передаточной ф-ции (допустим фильтра), ну и хочу посмотреть его АЧХ, для этого мне надо построить его дискретное преобразование фурье, т.е.
h(n) = h(n0) + h(n1) + ... + h(N-1);
H(jw) = h(0) + h(n1)*exp(-jw) + h(n2)*exp(-2jw) +...+ h(N-1)*exp(-(N-1)jw)
Собственно если коэффициентов h не много всё нормално, но если их количество растёт, то прописывать умножение каждого из них на экспоненту становится нудной задачкой... Вообщем можно это дело как то оптимизировать? Никак не допру... Матлаб ругается, т.к. ему нужно совпадение размера матриц, а круговую частоту w я обычно задаю так w = [1:500]*pi/500 или что то в этом духе, количество же коэффициентов значительно меньше.
2) И ещё, как задать скажем такую ф-цию:
h(n) = sin(n)/n, n /= 0;
h(n) = 1/pi, n = 0;
Пробовал при помощи if, но что то не выходит... 05.gif
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
Грендайзер
сообщение Sep 15 2016, 13:35
Сообщение #2


Местный
***

Группа: Участник
Сообщений: 368
Регистрация: 18-04-11
Из: Город-герой Москва
Пользователь №: 64 451



Я пробовал использовать эту ф-цию... но всё стало ещё непонятней...
Go to the top of the page
 
+Quote Post
mdmitry
сообщение Sep 15 2016, 14:54
Сообщение #3


Начинающий профессионал
*****

Группа: Свой
Сообщений: 1 215
Регистрация: 25-10-06
Из: СПб
Пользователь №: 21 648



Цитата(Грендайзер @ Sep 15 2016, 17:35) *
Я пробовал использовать эту ф-цию... но всё стало ещё непонятней...

Цитата
Написал кодик для того, что бы построить АЧХ и ФЧХ:

Я не обратил внимание на эти термины. НЕТ у сигнала АЧХ и ФЧХ, есть амплитудный и фазовый спектры! Курс радиотехнические цепи и сигналы в помощь, учебники есть. Заодно и книгу Сергиенко А.Б. по цифровой обработке сигналов надо почитать.
Функцию unwrap в данном случае использовать смысла нет, она как раз для построения ФЧХ и нужна.

По фазовому спектру: посмотрите данные FFT и обнаружите, что результаты представляют собой весьма малые числа, как в вещественной, так и в мнимой частях. При вычислении аргумента комплексного числа в этом случае потеря точности практически неизбежна. Откуда Ваш изрезанный график. К тому же, у Вас кратна частота сигнала частоте дискретизации. Поставьте частоту 801, 805 и графики сильно изменятся. Будут, на мой взгляд, более осмысленными для первоначального понимания. Ничего фантастического у Вас нет, все объяснимо.


--------------------
Наука изощряет ум; ученье вострит память. Козьма Прутков
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- Грендайзер   Как задать значения ф-ций   Aug 31 2016, 12:20
- - MiPe   Цитата(Грендайзер @ Aug 31 2016, 15:20) д...   Aug 31 2016, 12:27
|- - Грендайзер   Ничё се Цитата(MiPe @ Aug 31 2016, 15...   Aug 31 2016, 13:01
|- - MiPe   Цитата(Грендайзер @ Aug 31 2016, 16:01) В...   Aug 31 2016, 13:09
- - Грендайзер   Ух ты, спасибо за помощь Правда с функцией tf не...   Aug 31 2016, 13:19
|- - MiPe   Цитата(Грендайзер @ Aug 31 2016, 16:19) У...   Aug 31 2016, 13:24
- - Tpeck   Цитата(Грендайзер @ Aug 31 2016, 15:20) З...   Aug 31 2016, 14:27
|- - alex_os   Если передаточная ф-я фильтра то можно freqz: fr...   Aug 31 2016, 14:53
- - Grizzzly   Цитата(Грендайзер @ Aug 31 2016, 15:20) 2...   Aug 31 2016, 15:45
- - Грендайзер   ЦитатаПочитайте книгу Matlab для студента. Спасибо...   Sep 1 2016, 08:05
|- - Grizzzly   Цитата(Грендайзер @ Sep 1 2016, 11:05) do...   Sep 1 2016, 08:08
- - Грендайзер   ЦитатаДа, это обычный синк. Ясно, спасибо. Правда ...   Sep 1 2016, 08:21
- - Грендайзер   Снова здравствуйте. Продолжаю осваивать Matlab и т...   Sep 15 2016, 12:31
|- - mdmitry   Цитата(Грендайзер @ Sep 15 2016, 16:31) С...   Sep 15 2016, 12:39
||- - Грендайзер   Цитата(mdmitry @ Sep 15 2016, 17:54) Я не...   Sep 15 2016, 16:57
||- - Grizzzly   Цитата(Грендайзер @ Sep 15 2016, 19:57) И...   Sep 15 2016, 18:01
|- - Grizzzly   Цитата(Грендайзер @ Sep 15 2016, 16:35) Я...   Sep 15 2016, 15:17
- - Грендайзер   Сейчас ситуация поменялась. Всёравно не совсем пон...   Sep 15 2016, 18:19


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 29th July 2025 - 21:48
Рейтинг@Mail.ru


Страница сгенерированна за 0.02286 секунд с 7
ELECTRONIX ©2004-2016