Как задать значения ф-ций, Начал изучать матлаб и тут понеслось... Опции

[Грендайзер]

Здравствуйте. Начал изучать матлаб и появилось несколько вопросов на которые никак не найду ответов. Вопросы глупые конечно, но практики маловато так что...
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, но что то не выходит...

[MiPe]

для этого мне надо построить его дискретное преобразование фурье

Ага, а еще вам нужна миллиметровка и карандаш

CODE

coef = tf([1 1],[0.5 1 1]);
bode(coef)

2) И ещё, как задать скажем такую ф-цию:
h(n) = sin(n)/n, n /= 0;
h(n) = 1/pi, n = 0;
Пробовал при помощи if, но что то не выходит...

А вы покажите, как вы пробовали, а мы поможем...

[Грендайзер]

Ничё се

А вы покажите, как вы пробовали, а мы поможем...

Вот так:

if n == 0
h = w/pi
else
h = sin(w*n)/n*pi
end

От так

[MiPe]

Вот так:

if n == 0
h = w/pi
else
h = sin(w*n)/n*pi
end

От так

Такая запись должна работать для скалярных n. Если у вас n - вектор (MATLAB же), то лучше как-то так:

CODE

function h = myfun(n)
h = sin(n)./n;
h(n == 0) = 1/pi;

А потом можно делать так:

CODE

plot(myfun(-2:0.1:2))