Как сделать единицу с нуля?

Question

[novichok151413]

Добрый вечер, прошу помощи в написании скрипта для импульсной реакции рекурсивного ЦФ фильтра. Написал следующее:

q = 1;
a = [];
for i = 0:30
    q = 0.71 * q;
    a(end + 1) = q;
end
nexttile
stem(a,'-k', 'LineWidth', 2)
set(gca, 'FontName', 'Helvetica', 'FontSize', 12); 
grid on;
box off;
xlabel('i (samples)');
ylabel('Amplitude');

но не могу придумать, как поставить единицу в 0, в данном варианте скрипта её нет вообще, но если добавлять её как a = [1], то она встаёт в под единицу на графике, а не в ноль. Помогите исправить

Answer 1

[Александр Скуснов]

Под рукой Матлаба нет, но вы выводите график в режиме f(y), а надо в режиме f(x,y...) (где f это plot, stem и т.п.). Первый режим организует x как индексы y (целые числа, начиная с 1).
Организуйте x, начиная с 0, шаг любой (тоже может быть равен 1).

Answer 2

[Сергей Галкин]

x=0:30;
i=0:30;
a=0.71.^(i+1);
stem(i,a,'k','LineWidth',2)
set(gca, 'FontName', 'Helvetica', 'FontSize', 12); 
grid on
box off
xlabel('i (samples)')
ylabel('Amplitude')