График в MATLAB, как нарисовать? Опции

[Stefan1]

Доброго времени суток.
У меня есть зависимость x от у, которая изменяется с течением времени. Как мне построить график зависимости х от у в последний момент времени.

Сообщение отредактировал Stefan1 - Nov 20 2012, 09:56

[V_G]

Ну, так пишите plot(x[последние моменты],y[последние моменты]),
где последние моменты - диапазон индексов массивов x,y, соответствующих по-вашему последним моментам времени.

[Stefan1]

Ну, так пишите plot(x[последние моменты],y[последние моменты]),
где последние моменты - диапазон индексов массивов x,y, соответствующих по-вашему последним моментам времени.

Спасибо, заработало!

Сообщение отредактировал Stefan1 - Nov 21 2012, 17:04

[Stefan1]

По ходу решения задачи опять возник вопрос: а можно ли из синусоидальной функции с наличием высших гармоник u в зависимости от времени выделить первую гармонику и построить ее на графике также в зависимости от времени?
Я раскладываю в ряд Фурье функцию и затем выделяю из нее первую гармонику:

uf=fft(u); u1=uf(2)/256*2; uf=[];

Просто uf(2) - не строится. Как мне построить первую гармонику напряжения u в зависимости от времени (отбросив высшие гармоники)?

[Xenia]

По ходу решения задачи опять возник вопрос: а можно ли из синусоидальной функции с наличием высших гармоник u в зависимости от времени выделить первую гармонику и построить ее на графике также в зависимости от времени?
Я раскладываю в ряд Фурье функцию и затем выделяю из нее первую гармонику:
uf=fft(u); u1=uf(2)/256*2; uf=[];
Просто uf(2) - не строится. Как мне построить первую гармонику напряжения u в зависимости от времени (отбросив высшие гармоники)?

Наверное, вычислить обратное преобразование Фурье - ifft.
Например, так:
Стираем (зануляем) в массиве uf все частоты, кроме нужной вам
u1(1)=0
u1(3:length(u))=0
или в одну строку, если она вас не устрашит
u1([1 3:length(u)])=0
написала length(u1), т.к. не знаю точно, сколько у вас частот, вы же можете подставить конкретную константу,
потом
y=ifft(u1);
а дальше строим график обычным образом.

[Stefan1]

Наверное, вычислить обратное преобразование Фурье - ifft.
Например, так:
Стираем (зануляем) в массиве uf все частоты, кроме нужной вам
u1(1)=0
u1(3:length(u))=0
или в одну строку, если она вас не устрашит
u1([1 3:length(u)])=0
написала length(u1), т.к. не знаю точно, сколько у вас частот, вы же можете подставить конкретную константу,
потом
y=ifft(u1);
а дальше строим график обычным образом.

Все работает, только вот что-то не похоже на 1-ю гармонику (т.е. синусоиду). Может быть из-за того что я не то число ставлю вместо length(u1), т.к. у меня сигнал модулирован низкой частотой и первая гармоника здесь соответствует модулирующиму сигналу, а не несущей. Как бы теперь выделить первую гармонику несущей - наверно в ручную подбирать? т.к. частот модулирующего сигнала я не знаю.

Сообщение отредактировал Stefan1 - Nov 22 2012, 10:42

[Xenia]

Все работает, только вот что-то не похоже на 1-ю гармонику (т.е. синусоиду). Может быть из-за того что я не то число ставлю вместо length(u1), т.к. у меня сигнал модулирован низкой частотой и первая гармоника здесь соответствует модулирующиму сигналу, а не несущей. Как бы теперь выделить первую гармонику несущей - наверно в ручную подбирать? т.к. частот модулирующего сигнала я не знаю.

Возьмем 8-точеный пример:

1) Загоним туда косинусоиду (удобна тем, что ее спектр действителен):

Код

x = [0:7]          // команда
x =
     0     1     2     3     4     5     6     7     // отклик МатЛаба

y = cos(2*pi/8*x)          // команда
y =
    1.0000    0.7071    0.0000   -0.7071   -1.0000   -0.7071   -0.0000    0.7071     // отклик МатЛаба

2) Делаем fft

Код

u=fft(y)          // команда
u =
  -0.0000     4.0000    0.0000   -0.0000    0.0000   -0.0000    0.0000    4.0000     // отклик МатЛаба
u'          // команда напечатать в столбик
  -0.0000  // u(1) постоянная составляющая      
   4.0000  // u(2) амплитуда первой гармоники (самой медленной)
   0.0000  // u(3) амплитуда второй гармоники          
  -0.0000  // u(4) амплитуда третей гармоники
   0.0000  // u(5) амплитуда четвертой гармоники, частота Найквиста (самая быстрая)          
  -0.0000  // u(6) зеркальное комплексносопряженное от u(4)
   0.0000  // u(7) зеркальное комплексносопряженное от u(3)          
   4.0000  // u(8) зеркальное комплексносопряженное от u(2)

Ось симметрии проходит через частоту Найквиста.
Постоянная составляющая своего отражения не имеет.
"Комплексносопряженное" означает, что действительная часть числа сохраняется, а мнимое изменяет свой знак.

3) Обнуляем ненужные гармоники
Кроме u(2) надо оставить незануленной еще и ее отражение u(8).
Когда в массиве не 8 элементов, а больше, то надо рассчитать номер зеркального отражения.
Для первой гармоники он самый последний в ряду:

Код

u([1 3:length(u)-1])=0
u =
    0     4.0000      0        0        0        0        0          4.0000

4) Делаем обратное преобразование Фурье

Код

y2=ifft(u)
y2 =
    1.0000    0.7071    0.0000   -0.7071   -1.0000   -0.7071   -0.0000       0.7071

Готово! Поскольку у нас изначально была только 1-ая гармоника, то она и осталась.

[Stefan1]

Когда в массиве не 8 элементов, а больше, то надо рассчитать номер зеркального отражения.

Спасибо за пояснения! Только вот у меня этих элементов - несколько тысяч, как это зеркальное отражение рассчитать то?
И еще один вопрос: если мой сигнал модулирован низкой частотой, то 1-я гармоника будет именно этого модулирующего сигнала, может есть способы его убрать?

Сообщение отредактировал Stefan1 - Nov 22 2012, 18:07