Дискуссионный математический форумМатематический форум
Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 3 часа [ Летнее время ]
новый онлайн-сервис
число, сумма и дата прописью

Часовой пояс: UTC + 3 часа [ Летнее время ]




Начать новую тему Ответить на тему  [ 1 сообщение ] 
Автор Сообщение
 Заголовок сообщения: Применение ммп для матриц
СообщениеДобавлено: 17 окт 2015, 16:01 
Не в сети
Начинающий
Зарегистрирован:
17 окт 2015, 15:50
Сообщений: 1
Cпасибо сказано: 0
Спасибо получено:
0 раз в 0 сообщении
Очков репутации: 1

Добавить очки репутацииУменьшить очки репутации
Здравствуйте, помогите, пожалуйста, разобраться с задачей:
На вход подается бинарная матрица, каждый элемент которой Хij, зависит от элемента θi вектора θ и элемента βj вектора β. Требуется по выборке {Хij} определить вектора θ и β. (рис.1)

Определение векторов происходит на основе метода максимального правдоподобия. (β1 полагается равным нулю, т.е остальсные β будут рассчитываться относительно β1) Функция правдоподобия (рис.2), (рис.3)
С целью максимизации используется квазиньютоновский алгоритм Бройдена-Флетчера-Голдфарба-Шенно. Стартовой точкой метода полагается нулевой вектор u0 (рис.4)
Результатом выполнения алгоритма должны являться вектор δ и вектор θ

При реализации алгоритма для примера (с известным ответом) матрицы 10*10 никак не могу получить нужные вектора.

f.m
function L = f(u)%f(u_theta, u_delta)
L = 0;
theta = 1:10;
delta = 1:10;
for i = 1:10
theta(i) = u(i);
end;

for i = 11:20
delta(i-10) = u(i);
end;
x = [1, 0, 1, 1, 0, 1, 1, 1, 0, 1 ; ...
0, 0, 0, 1, 1, 1, 0, 1, 1, 1 ; ...
1, 1, 1, 0, 1, 1, 1, 1, 1, 1 ; ...
1, 1, 0, 0, 0, 0, 1, 0, 1, 1 ; ...
0, 1, 1, 0, 0, 1, 0, 0, 1, 1 ; ...
1, 0, 1, 1, 1, 1, 0, 1, 0, 1 ; ...
0, 1, 1, 1, 1, 1, 1, 1, 1, 1 ; ...
0, 0, 1, 1, 1, 1, 0, 0, 0, 0 ; ...
1, 1, 1, 0, 1, 1, 1, 1, 1, 1 ; ...
1, 1, 1, 1, 1, 1, 0, 1, 1, 0 ; ];
for i = 1:10
for j = 2:10
L = L + (theta(i))*x(i,1) - log(1 + exp(theta(i))) + (theta(i)-delta(j))*x(i,j) - log(1 + exp(theta(i) - delta(j)));
end
end
L = -L;
end %f(u_theta(k2,k1), u_delta(k2,k1) f(u_theta(k2,k1+1),u_delta(k2,k1+1)zz


Un.m
clear all;
x0 = [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ];

fminunc(@f, x0)



>> Un
Warning: Gradient must be provided for trust-region algorithm;
using line-search algorithm instead.
> In fminunc at 383
In Un at 4

Local minimum found.

Optimization completed because the size of the gradient is less than
the default value of the function tolerance.

<stopping criteria details>

ans =
Columns 1 through 18
1.3808 -1.1211 2.6343 0.6940 -1.4191 1.3808 -0.3430 -1.7571 2.6343 1.8660 0 0.0000 -1.4667 -0.0000 -0.6934 -2.4800 0.6507 -0.6934

Columns 19 through 20
-0.6934 -1.4667
>>


Очень долго не могу получить требуемые вектора для конкретной матрицы. (рис5)
Пожалуйста, помогите разобраться. Очень надеюсь на вашу помощь!!!

Изображение
Изображение
Изображение
Изображение
Изображение

Вернуться к началу
 Профиль  
Cпасибо сказано 
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему      Страница 1 из 1 [ 1 сообщение ]

 Похожие темы   Автор   Ответы   Просмотры   Последнее сообщение 
Применение функций от матриц

в форуме Линейная и Абстрактная алгебра

Flow

0

208

27 дек 2018, 19:56

Вычеты и их применение

в форуме Объявления участников Форума

AnnKond

4

256

02 янв 2023, 12:58

Применение кватернионов

в форуме Линейная и Абстрактная алгебра

mono_libre

12

3272

26 апр 2015, 01:22

Применение математики

в форуме Размышления по поводу и без

Andrei762

4

966

28 апр 2014, 14:28

Применение логарифма

в форуме Алгебра

xxenoss

0

282

11 май 2021, 03:19

Применение производной

в форуме Алгебра

Nora

4

144

30 дек 2023, 18:22

Применение производных

в форуме Начала анализа и Другие разделы школьной математики

Grach

1

259

27 апр 2020, 01:33

Применение тройных интегралов

в форуме Интегральное исчисление

newak

4

198

26 апр 2021, 22:33

Задача на применение преобразований

в форуме Геометрия

nsara

16

567

13 апр 2018, 12:56

Применение формулы Байеса

в форуме Комбинаторика и Теория вероятностей

Maxersh

7

428

08 июн 2020, 22:33


Часовой пояс: UTC + 3 часа [ Летнее время ]



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 6


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Перейти:  

Яндекс.Метрика

Copyright © 2010-2023 MathHelpPlanet.com. All rights reserved