|
|
 |
Ответов
|
Nov 27 2008, 21:10
|
Участник

Группа: Участник
Сообщений: 72
Регистрация: 7-06-08
Пользователь №: 38 128

|
Вот мой код, без претензий на робастность и скорость: Код clc;
N = 100; % кол-во отсчетов n = 1:N; d = .1; % постоянная составляющая p = rand()*pi; % случайная фаза f = 1/(4*N); % частота (в N отсчетах 1/4 периода) SNR = 40; % ОСШ в дБ
% собсно анализируемые отсчеты s = sin(2*pi*f*n+p) + d + randn(1, N)*(d+sqrt(2))*10^(-SNR/20);
st_ = [ .5 2*pi*f 0 0 ]; fo_ = fitoptions('method','NonlinearLeastSquares','Lower',[0 0 -1.7 -1],'Upper',[5 0.10000000000000001 1.7 0.10000000000000001]); set(fo_,'Startpoint',st_); ft_ = fittype('a*sin(b*x+c)+d',... 'dependent',{'y'},'independent',{'x'},... 'coefficients',{'a', 'b', 'c', 'd'}); cf_ = fit(n',s',ft_,fo_);
a = coeffvalues(cf_);
a(4) Да, у меня тоже была мысль о разложении синуса в ряд, или как-то перейти к полиному Чебышева, но тогда постоянная составляющая будет спрятана в свободном члене полинома, и как ее оттуда вычленить я не знаю. 2GetSmart: ну если в распоряжении есть больше (м.б. больше равно) одного периода, то все становится понятно
|
|
|
|
Сообщений в этой теме
vladimir_m Определение постоянной составляющей синусоиды... Nov 21 2008, 13:03 D.I.M.A Оценивать мат. ожидание процесса. Nov 21 2008, 13:16 rezident А постоянная составляющая по какому-то закону изме... Nov 21 2008, 13:24 vladimir_m ЦитатаА постоянная составляющая по какому-то закон... Nov 21 2008, 13:34 sysel По идее, для получения постоянной составляющей над... Nov 21 2008, 13:49 Самурай Цитата(vladimir_m @ Nov 21 2008, 16:03) З... Nov 21 2008, 14:12 vladimir_m Большое окно - это хорошо. Расчетом занимается и т... Nov 21 2008, 14:17 fontp Цитата(vladimir_m @ Nov 21 2008, 17:17) Б... Nov 21 2008, 15:07  Serg76 Цитата(fontp @ Nov 21 2008, 18:07) Фильтр... Nov 26 2008, 16:34 ivan219 Зделать FFT и 0 отщёт будет = пост состовляющей. Nov 21 2008, 19:57 GetSmart Цитата(ivan219 @ Nov 22 2008, 01:57) Здел... Nov 24 2008, 22:23 ivan219 1 Nov 21 2008, 19:57 Stanislav Цитата(ivan219 @ Nov 21 2008, 22:57) Здел... Nov 21 2008, 20:36 Михаил_K А чем не подходит это
http://electronix.ru/forum/i... Nov 22 2008, 10:01 Stanislav Цитата(Михаил_K @ Nov 22 2008, 13:01) А ч... Nov 22 2008, 12:53  Михаил_K Цитата(Stanislav @ Nov 22 2008, 15:53) А ... Nov 26 2008, 04:58   Stanislav Цитата(Михаил_K @ Nov 26 2008, 07:58) Что... Nov 26 2008, 09:02    Михаил_K Цитата(Stanislav @ Nov 26 2008, 12:02) Мд... Nov 26 2008, 09:04     Stanislav Цитата(Михаил_K @ Nov 26 2008, 12:04) Буд... Nov 26 2008, 09:08 vladimir_m Всем спасибо!
Лично я остановился на варианте ... Nov 22 2008, 12:31 Tanya Цитата(vladimir_m @ Nov 22 2008, 15:31) В... Nov 22 2008, 12:58  fontp Цитата(Tanya @ Nov 22 2008, 15:58) Нескол... Nov 24 2008, 08:17   Stanislav Цитата(fontp @ Nov 24 2008, 11:17) В прин... Nov 24 2008, 08:30    fontp Цитата(Stanislav @ Nov 24 2008, 11:30) Ст... Nov 24 2008, 08:41     Stanislav Цитата(fontp @ Nov 24 2008, 11:41) Просто... Nov 24 2008, 09:01      fontp Цитата(Stanislav @ Nov 24 2008, 12:01) Кр... Nov 24 2008, 09:51       alexkok Цитата(fontp @ Nov 24 2008, 12:51) Про ... Nov 24 2008, 21:03       Stanislav Цитата(fontp @ Nov 24 2008, 12:51) Да лад... Nov 25 2008, 00:31 GetSmart Цитата(Stanislav)Для вычисления параметров синусои... Nov 22 2008, 17:12 Stanislav Цитата(GetSmart @ Nov 22 2008, 20:12) На ... Nov 22 2008, 18:09 GetSmart Цитата(Stanislav)Послушайте, уважаемый Получатель ... Nov 22 2008, 23:38 Stanislav Цитата(GetSmart @ Nov 23 2008, 02:38) Тол... Nov 23 2008, 00:14 GetSmart Цитата(Stanislav)Возможность практической реализац... Nov 23 2008, 22:24 Stanislav Цитата(GetSmart @ Nov 24 2008, 01:24) Я у... Nov 23 2008, 23:20  GetSmart Цитата(Stanislav @ Nov 24 2008, 05:20) А ... Nov 23 2008, 23:45   Stanislav Цитата(GetSmart @ Nov 24 2008, 02:45) Ну ... Nov 23 2008, 23:49 GetSmart По приходу нового буфера берётся несколько последн... Nov 24 2008, 23:54 Stanislav ЗЫ. fontp, я готов прекратить спор, если будет при... Nov 25 2008, 07:30 fontp Цитата(alexkok @ Nov 25 2008, 00:03) Не о... Nov 25 2008, 07:54  Stanislav Цитата(fontp @ Nov 25 2008, 10:54) Fittin... Nov 25 2008, 08:23  alexkok Цитата(fontp @ Nov 25 2008, 10:54) Fittin... Nov 25 2008, 17:52   blackfin Цитата(alexkok @ Nov 25 2008, 20:52) Fitt... Nov 26 2008, 05:23 GetSmart Цитата(alexkok)Fitting, как Вам уже ответили, при ... Nov 25 2008, 23:17 alexkok Цитата(GetSmart @ Nov 26 2008, 02:17) А э... Nov 26 2008, 08:10 GetSmart Цитата(blackfin)Для выделения полезного сигнала на... Nov 26 2008, 07:03 alex_os Цитата(GetSmart @ Nov 26 2008, 10:03) При... Nov 26 2008, 08:25 RadioJunior пост не в тему...
Stanislav, не обижайтесь, но эт... Nov 26 2008, 09:11 Stanislav Цитата(RadioJunior @ Nov 26 2008, 12:11) ... Nov 26 2008, 21:29 Михаил_K Цитата(RadioJunior @ Nov 26 2008, 12:11) ... Nov 28 2008, 05:00  Stanislav Цитата(mvb @ Nov 28 2008, 00:10) Вот мой ... Nov 28 2008, 06:55 GetSmart Цитата(RadioJunior)Со мной кто-нибудь согласен по ... Nov 26 2008, 10:43 mvb Господа, сразу прошу прощения за то, что вмешиваюс... Nov 27 2008, 00:07 Stanislav Цитата(mvb @ Nov 27 2008, 03:07) Для реал... Nov 27 2008, 00:23 GetSmart Цитата(mvb)2. Как здесь можно использовать дпф или... Nov 27 2008, 04:12 mvb Если бы у меня был бы фрейм с достаточным количест... Nov 28 2008, 10:03 Stanislav Цитата(mvb @ Nov 28 2008, 13:03) Если бы ... Nov 28 2008, 11:45 GetSmart Цитата(mvb)2. Далее скомпенсировал бы ошибку в ДПФ... Nov 28 2008, 13:08 Stanislav Цитата(GetSmart @ Nov 28 2008, 16:08) Да,... Nov 28 2008, 20:59 mvb Попробовал на досуге определить период синуса в ма... Nov 30 2008, 09:21
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|