50628

Метод покоординатного спуска

Лабораторная работа

Математика и математический анализ

Напишем программу минимизации функции методом покоординатного спуска: Минимизируем исходную функцию стандартными средствами MatLab

Русский

2014-01-27

56.5 KB

77 чел.

Лабораторные работы

№8

Тема

Метод покоординатного спуска

Ф.И.О.

Пастухова Светлана Владимировна

Группа

403

Вариант

15

Минимизировать функцию f(x,y)=ax + by + exp(cx2 + dy2 )

 N

a

b

c

d

9

22

-0.2

1.44

0.23

Порядок выполнения работы:

  1.  Построим график функции:

ezsurf('13*x-0.2*y+exp(1.44*x^2+0.23*y^2)')

  1.  Минимизируем исходную функцию стандартными средствами MatLab:

x = [0,-6];

min = fminsearch(@my_fun,x)

где x=[0,-6] – начальное приближение, а @my_fun:

function f = my_fun(x)

f =13*x(1)-0.2*x(2)+exp(1.44*x(1)^2+0.23*x(2)^2); 

Получим:

min =[-1.0167, 0.0974]

f(xmin)=  -8.7965

  1.  Напишем программу минимизации функции методом покоординатного спуска:
  2.  %Минимизируем исходную функцию стандартными средствами MatLab:
  3.  x = [0,-6];
  4.  min = fminsearch(@my_fun,x)
  5.  %где
  6.  x=[0,-6] %– начальное приближение, а @my_fun:
  7.  function f = my_fun(x)
  8.  f =13*x(1)-0.2*x(2)+exp(1.44*x(1)^2+0.23*x(2)^2);
  9.  %Получим:
  10.  min =[-1.0167, 0.0974]
  11.  f(xmin)=  -8.7965
  12.  %Напишем программу минимизации функции методом покоординатного спуска:
  13.  aaa=13; bbb=-0.2; ccc=1.44; ddd=0.23;
  14.  x1=0;
  15.  x2=-6;
  16.  e=0.0001;
  17.  k=2;
  18.  h=0.2;
  19.  e1=e/k;
  20.  l=0;
  21.  d=e+1;
  22.  while(d>e)
  23.      d=abs(h);
  24.      h1=h;
  25.      z=aaa*x1+bbb*x2+exp(ccc*x1^2+ddd*x2^2);
  26.      z1=z+1;
  27.      d1=abs(h1);
  28.      while (d1>e)
  29.          d1=abs(h1);
  30.          x1=x1+h1;
  31.          z1=aaa*x1+bbb*x2+exp(ccc*x1^2+ddd*x2^2);
  32.          if (z1>=z) h1=-h1/k; end
  33.          z=z1;
  34.      end
  35.      h1=h;
  36.      z=aaa*x1+bbb*x2+exp(ccc*x1^2+ddd*x2^2);
  37.      d1=e1+1;
  38.      z1=z+1;
  39.      d1=abs(h1);
  40.      while (d1>e)
  41.          d1=abs(h1);
  42.          x2=x2+h1;
  43.          z1=aaa*x1+bbb*x2+exp(ccc*x1^2+ddd*x2^2);
  44.          if (z1>=z)  h1=-h1/k; end
  45.          z=z1;
  46.      end
  47.      h=h/k;
  48.      l=l+1;
  49.  end
  50.  x1
  51.  x2
  52.  aaa*x1+bbb*x2+exp(ccc*x1^2+ddd*x2^2)
  53.  e
  54.  %В результате получим:
  55.  min =[ -1.0165, 0.0981]
  56.  fmin = -8.7965
  57.  n=12
  58.   
  59.  

В результате получим:

min =[ -1.0165, 0.0981]

fmin = -8.7965

n=12


 

А также другие работы, которые могут Вас заинтересовать

15384. Телевизионная реклама батончика NUTS 167 KB
  Введение У нас в стране реклама приобретает все большее значение. Ее целью является привлечение внимания покупателей к фирме ее возможностям и выпускаемым товарам подчеркивая их высокое качество новизну надежность удобство использования и приемлемую цену. Во в
15385. Толерантность и образование Теория 196.5 KB
  Помни что ты только человек Слова эпиграфа заимствованы из триумфального ритуала римских полководцев. Их повторял триумфатору специально приставленный для этой цели глашатай. Тогда они звучали как предостережение против установления тирании.Сегодня так можно сформу...
15386. Уголовное право. Незаконное хранение, приобретение и распространение наркотических средств 195.5 KB
  Введение Повышение эффективности уголовноправовой борьбы с незаконным хранением распространением наркотиков во многом связано со своевременным совершенствованием уголовного законодательства и практики его применения. Вопрос о путях совершенствования уголовног
15387. Уголовное право. Преступления против собственности, совершенные путем обмана или злоупотребления доверием 293 KB
  Собственность - это исторически определенная общественная форма присвоения материальных благ, прежде всего, средств производства. Отношения собственности составляют основу экономики. Собственность представляет собой одну из важнейших социальных ценностей...
15388. Ценообразование и ценовая политика на предприятии Боготолмолоко 184 KB
  Теоретические аспекты ценообразования на предприятии 1.1 Сущность и функции цен в рыночной экономике 1.2 Процесс и факторы ценообразования 1.3 Цели ценовой политики и принципы ценообразования 1.4 Информация необходимая при принятии у
15389. Динамика показателей денежного обращения в РФ (на начало I квартала 2010 г) 242 KB
  Курсовая работа Экономическая статистика Содержание Введение Глава I. Основы статистики денежного обращения 1.1 Денежные агрегаты 1.2 Денежная база 1.3 Узкие деньги 1.4 Широкие деньги 1.5 Денежные агрегаты в России Глава II. Динамика показателей денежного ...
15390. Генератор пилообразного напряжения 189.5 KB
  Реферат Курсовой проект выполнен в объеме 32 страниц содержит двенадцать рисунков использовано 4 источника. Разработанный генератор пилообразного напряжения может быть использован в контрольноизмерительной аппаратуре в цифроаналоговых преобразователях и т. п.
15392. Эффективность производства кормов в ОАО Сибирская Нива 342.5 KB
  ОАО Сибирская Нива Автогараж Ц.Р.М. Склад ГСМ Отделение №1 Отделение №2 Отделение №3 Зерносклад Введение Кормопроизводство является основополагающей отраслью сельского хозяйства научнотехнический уровень развити