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

 
 
 
Reply to this topicStart new topic
> Функция неопределенности сигнала, вычисление ФН в MatLAB
petr0ff
сообщение May 22 2010, 19:13
Сообщение #1





Группа: Участник
Сообщений: 7
Регистрация: 22-05-10
Пользователь №: 57 458



возник следующий вопрос.
функция неопределенности ЛЧМ-сигнала имеет вытянутую форму по диагонали, если смотреть сверху.

однако столкнулся со следующей проблемой. скачал файл для вычисления ФН по ссылке
Matlab central
аналогично тут еще вариант

если по ней вычислять ФН ЛЧМ сигнала, то ФН получается в форме креста, если смотреть сверху.
На этой же странице приведен скрин скрипта. Там тоже ФН в форме креста. Почему так и какой должна быть ФН ЛЧМ-сигнала????
Go to the top of the page
 
+Quote Post
petr0ff
сообщение May 23 2010, 08:19
Сообщение #2





Группа: Участник
Сообщений: 7
Регистрация: 22-05-10
Пользователь №: 57 458



вот результат вычисления ФН ЛЧМ-сигнала

тут тоже приведен скрин работы функции
http://www.mathworks.com/matlabcentral/fil...y-function-plot

Сообщение отредактировал petr0ff - May 23 2010, 08:23
Go to the top of the page
 
+Quote Post
Serg76
сообщение May 23 2010, 11:52
Сообщение #3


Профессионал
*****

Группа: Участник
Сообщений: 1 050
Регистрация: 4-04-07
Пользователь №: 26 775



Цитата(petr0ff @ May 22 2010, 23:13) *
если по ней вычислять ФН ЛЧМ сигнала, то ФН получается в форме креста, если смотреть сверху.
На этой же странице приведен скрин скрипта. Там тоже ФН в форме креста. Почему так и какой должна быть ФН ЛЧМ-сигнала????

По-видимому у вас сигнал представлен суммой двух ЛЧМ импульсов - со спадающей и возрастающей частотой
Go to the top of the page
 
+Quote Post
petr0ff
сообщение May 23 2010, 20:49
Сообщение #4





Группа: Участник
Сообщений: 7
Регистрация: 22-05-10
Пользователь №: 57 458



Цитата(Serg76 @ May 23 2010, 14:52) *
По-видимому у вас сигнал представлен суммой двух ЛЧМ импульсов - со спадающей и возрастающей частотой

проверял на двух сигналах - chirp и свой ЛЧМ - тот же результат...

вот полученные квадратуры:

по ним - вроде бы действительно частота возрастает, потом спадает. но использую стандартную функцию chirp - она ведь формирует линейно возрастающую частоту...

Сообщение отредактировал petr0ff - May 23 2010, 21:28
Go to the top of the page
 
+Quote Post
petr0ff
сообщение May 23 2010, 21:57
Сообщение #5





Группа: Участник
Сообщений: 7
Регистрация: 22-05-10
Пользователь №: 57 458



на предыдущем графике перепутана одни из квадратур

вроде бы получил какие-то результаты.
но они оказались разными. Находил ФН ЛЧМ-сигнала с помощью двух функций:
ссылка 1
ссылка 2
получил следующие результаты:




какая из этих функций работает верно????
Go to the top of the page
 
+Quote Post
Serg76
сообщение May 24 2010, 07:57
Сообщение #6


Профессионал
*****

Группа: Участник
Сообщений: 1 050
Регистрация: 4-04-07
Пользователь №: 26 775



Цитата(petr0ff @ May 24 2010, 00:49) *
проверял на двух сигналах - chirp и свой ЛЧМ - тот же результат...

вот полученные квадратуры:

по ним - вроде бы действительно частота возрастает, потом спадает. но использую стандартную функцию chirp - она ведь формирует линейно возрастающую частоту...

как я вам и говорил, сигнал представлен сначала спадающей частотой, затем возрастающей. а вообще смоделируйте сигнал сами, так, по крайней мере, достоверно будете знать, что синтезировали
Go to the top of the page
 
+Quote Post
petr0ff
сообщение May 24 2010, 20:31
Сообщение #7





Группа: Участник
Сообщений: 7
Регистрация: 22-05-10
Пользователь №: 57 458



Цитата(Serg76 @ May 24 2010, 10:57) *
как я вам и говорил, сигнал представлен сначала спадающей частотой, затем возрастающей. а вообще смоделируйте сигнал сами, так, по крайней мере, достоверно будете знать, что синтезировали

хм... я вот тоже так думал, но строил график ЛЧМ сигнала и на нем видно, что частота только возрастает.
А это квадратуры. При расчете ФН я вместо (Q, I) взял (Q, Q). На графике видно, что эти квадратуры одинаковые. На самом деле квадратуры имеют вид:

после нахождения ошибки и еще небольшой доработки получилась диагональная ФН. т.е. все как надо. Крестообразная ФН получалась из-за того, что обрабатывал вещественный сигнал и ФН находил от вещественных сигналов т.е. p(S, S) а не p(S, S*).

Вопрос возник в следующем.
на официальном сайте матлаба разработчики скинули функции для расчета ФН. Они открываются по приведенным выше ссылкам. Если рассчитывать ФН по этим функциям, то почему-то получаются разные результаты. Графики приведены выше. Но какая из них правильная - не знаю. Мне кажется, что верхняя. Сам тоже написал функцию для расчета ФН и получил верхний график.

И ещё. Почему-то не вижу в чем причина. Использую функцию chirp и свою функцию для ЛЧМ. И они дают разный результат. Почему????
f0 = 25;
df = 20;
f1 = f0 - df;
f2 = f0 + df;
t = 0:0.0001:1;
s1 = chirp(t, f1, 1, f2, 'linear');

beta = (f2 - f1) / 1;
s2 = cos(2 * pi * (f1 + beta * t) .* t);

figure
subplot(2,1,1)
plot(t, s1), grid on
subplot(2,1,2)
plot(t, s2), grid on


видно, что на 2-ом графике в конце сигнала частота в 2 раза выше.
От моего сигнала квадратуры имеют такой же вид.

Сообщение отредактировал petr0ff - May 24 2010, 20:37
Go to the top of the page
 
+Quote Post
mdmitry
сообщение May 25 2010, 15:06
Сообщение #8


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

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



Цитата(petr0ff @ May 25 2010, 00:31) *
[sub]
beta = (f2 - f1) / 1;

Это правильное значение для beta? (деление на 1, а не на t_max)


--------------------
Наука изощряет ум; ученье вострит память. Козьма Прутков
Go to the top of the page
 
+Quote Post
petr0ff
сообщение May 25 2010, 18:30
Сообщение #9





Группа: Участник
Сообщений: 7
Регистрация: 22-05-10
Пользователь №: 57 458



Цитата(mdmitry @ May 25 2010, 18:06) *
Это правильное значение для beta? (деление на 1, а не на t_max)

в приведенном примере t_max = 1.
можно изменить код:
...
tmax = 0.5;
t = 0:0.0001:tmax;
s1 = chirp(t, f1, tmax, f2, 'linear');

beta = (f2 - f1) / tmax;
s2 = cos(2 * pi * (f1 + beta * t) .* t);
...

но результат от этого почему-то не изменился...
Go to the top of the page
 
+Quote Post
petr0ff
сообщение May 26 2010, 06:04
Сообщение #10





Группа: Участник
Сообщений: 7
Регистрация: 22-05-10
Пользователь №: 57 458



Цитата
Потому что в формуле для ЛЧМ beta пополам разделить забыли.
Go to the top of the page
 
+Quote Post

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

 


RSS Текстовая версия Сейчас: 16th June 2025 - 15:58
Рейтинг@Mail.ru


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