Генетический алгоритм - реализовать это в матлабе

@V_i_p · Регистрация: 01.12.2012

Author24 — интернет-сервис помощи студентам

Помогите реализовать это в матлабе пожалуйста. Срочно надо решить. Там есть документ это начало. В документе есть пример, по этому примеру надо сделать 10 пунктов которые приведены ниже

Порядок решения задачи
1.Создание исходной популяции: приемлемые варианты.
2.Аргументация выбора сбособа создания исходной популяции.
3.Кодирование особи-решения.
4.Выбор операторов кроссинговера и мутации.
5.Принятиее решения о методе выбора родительских особей.
6.Выбор подхода к сокращению популяции.
7.Реализация модели и графического симулятора.
8.Проверка работоспособности выбранной модификации модели ГА.
9.Определение оптимальных параметров ГА.
10.Оценка работы модели ГА и реализации.

@VSI · 15.11.2016, 14:38

V_i_p, загляните в самый низ этой страницы в раздел "Похожие темы". Можете найти для себя много полезной информации...

@V_i_p · 15.11.2016, 22:10 **[ТС]**

Я просмотрел, но ничего не нашол

Добавлено через 21 минуту
А можете помочь решить??? или хотя бы начало сделать чтобы мне было от чего оттолкнуться?

Добавлено через 1 час 11 минут

Кликните здесь для просмотра всего текста

Matlab M

n=10; 
x=(-n):(n);
//--РОДИТЕЛИ 0-10
for m=1:n 
i=10*rand(); 
roditel(1,m)=i; 
end
//--РОДИТЕЛИ целые числа--// 
roditel=randint(1,n,[-n n]); 
for m=1:n 
for z=1:p 
i=randint(1,1,[1 n]); 
k1(z,m)=roditel(1,i); 
i=randint(1,1,[1 n]); 
k2(z,m)=roditel(1,i); 
end 
end
 
for m=1:n 
for z=1:p 
x=k1(z,m); 
f1(z,m)=(1.85-x).*(cos(3.5.*x)-0.5); 
x=k2(z,m); 
f2(z,m)=(1.85-x).*(cos(3.5.*x)-0.5); 
end 
end 
[y1,Index1]=max(f1); %максимизация функции 
[y2,Index2]=max(f2);
 
for m=1:n 
i=Index1(1,m); 
roditel(1,m)=k1(i,m); 
i=Index2(1,m); 
roditel(2,m)=k2(i,m); 
end 
///кроссинговер
for m=1:n 
i=roditel(2,m)-roditel(1,m); 
for z=1:2 
qe=rand(); %случайный множитель 
potomok(z,m)=roditel(z,m)+i*qe; 
end 
end 
 
for m=1:n 
for z=1:2 
while potomok(z,m)<-10 | potomok(z,m)>10 
qe=rand(); 
potomok(z,m)=roditel(z,m)+i*qe; 
end 
end 
end 
 
 
 
 
вот есть код, но он не работает коректно..... можете подкоректировать пожалуйста????
 
[size="1"][color="grey"][I]Добавлено через 15 секунд[/I][/color][/size]
n=10; 
x=(-n):(n);
//--РОДИТЕЛИ 0-10
for m=1:n 
i=10*rand(); 
roditel(1,m)=i; 
end
//--РОДИТЕЛИ целые числа--// 
roditel=randint(1,n,[-n n]); 
for m=1:n 
for z=1:p 
i=randint(1,1,[1 n]); 
k1(z,m)=roditel(1,i); 
i=randint(1,1,[1 n]); 
k2(z,m)=roditel(1,i); 
end 
end
 
for m=1:n 
for z=1:p 
x=k1(z,m); 
f1(z,m)=(1.85-x).*(cos(3.5.*x)-0.5); 
x=k2(z,m); 
f2(z,m)=(1.85-x).*(cos(3.5.*x)-0.5); 
end 
end 
[y1,Index1]=max(f1); %максимизация функции 
[y2,Index2]=max(f2);
 
for m=1:n 
i=Index1(1,m); 
roditel(1,m)=k1(i,m); 
i=Index2(1,m); 
roditel(2,m)=k2(i,m); 
end 
///кроссинговер
for m=1:n 
i=roditel(2,m)-roditel(1,m); 
for z=1:2 
qe=rand(); %случайный множитель 
potomok(z,m)=roditel(z,m)+i*qe; 
end 
end 
 
for m=1:n 
for z=1:2 
while potomok(z,m)<-10 | potomok(z,m)>10 
qe=rand(); 
potomok(z,m)=roditel(z,m)+i*qe; 
end 
end 
end

вот есть код, но он не работает коректно..... можете подкоректировать пожалуйста????

@V_i_p · 26.11.2016, 20:34 **[ТС]**

Помогите пожалуйста изменить код для 2-х переменных в пространстве xyz

Кликните здесь для просмотра всего текста

Matlab M

n=10; %размер популяции
xi=-10:0.1:10; %в дальнейшем интервал от -10 до 10
p=2;%количество элементов для турнирного отбора
 
%---Рандомно задаем значения---
for m=1:n
I=length(xi);
i=fix(I*rand());
roditel(1,m)=xi(1,i);
end
 
%---Заполняем массив случайными значениями из roditel---
%---k1 первый массив родителей, k2 второй массив родителей---
for m=1:n
for z=1:p
i=randint(1,1,[1 n]);
k1(z,m)=roditel(1,i); %k-турнирная группа
i=randint(1,1,[1 n]);
k2(z,m)=roditel(1,i);
end
end
 
%---Находим максимальные элементы из турнирной группы и записываем их---
for m=1:n
for z=1:p
x=k1(z,m);
f1(z,m)=(1.85-x).*(cos(3.5.*x)-0.5);
x=k2(z,m);
f2(z,m)=(1.85-x).*(cos(3.5.*x)-0.5);
end
end
[y1,Index1]=max(f1); %максимизация функции
[y2,Index2]=max(f2);
 
%---На основе y max записываем элементы в массив roditel--- 
for m=1:n
i=Index1(1,m);
roditel(1,m)=k1(i,m);
i=Index2(1,m);
roditel(2,m)=k2(i,m);
end
 
for m=1:n
i=roditel(2,m)-roditel(1,m);
for z=1:2
qe=rand(); %случайный множитель  
potomok(z,m)=roditel(z,m)+i*qe;
end
end
 
for m=1:n
for z=1:2
while potomok(z,m)<-10 | potomok(z,m)>10  
qe=rand(); %случайный множитель
potomok(z,m)=roditel(z,m)+i*qe;
end
end
end
 
koef_mut=0.1;
koef=[(xi(1,1))+(xi(1,I))-(xi(1,1))]/1000;
i=1;
for m=1:n
for z=1:2
mut=rand();
if mut<koef_mut
potomok_(1,i)=potomok(z,m)+koef;          %---мутация
else
potomok_(1,i)=potomok(z,m);
end
roditel_(1,i)=roditel(z,m);
i=i+1;
end
end
 
red1=n/2; %произвольно выбираем кол-во родительских особей
red2=(n*2)-red1;
 
%---Находим максимальные элементы из группы родителей и потомков и записываем их---
for m=1:n*2
x=roditel_(1,m);
y1(1,m)=(1.85-x).*(cos(3.5.*x)-0.5);
x=potomok_(1,m);
y2(1,m)=(1.85-x).*(cos(3.5.*x)-0.5); %---редукция
end
 
SredZna4_rod=mean(y1);
SredZna4_pot=mean(y2);
 
%---Записываем в новый массив red1 элементов из y1
m_=1;
m=1;
while m_<=red1 & m<=n*2
if y1(1,m)>=SredZna4_rod
i=y1(1,m);
f(1,m_)=i;
kone4_mass(1,m_)=roditel_(1,m);        %---редукция
m_=m_+1;
end
m=m+1;
end
 
%---Записываем в новый массив red2 элементов из y2
m_=red1+1;
m=1;
while m_<=n*2 & m<=n*2
if y2(1,m)>=SredZna4_pot
i=y2(1,m);
f(1,m_)=i;
kone4_mass(1,m_)=potomok_(1,m);             %---редукция
m_=m_+1;
end
m=m+1;
end
 
fplot('(1.85-x).*(cos(3.5.*x)-0.5)',[-10 10])
hold on
plot(kone4_mass,f,'r*')

@Ирина6 · 02.12.2017, 06:22

Ребят, кто может изменить этот код для двух переменных? Очень нужно

@VSI Модератор $Эксперт по математике/физике$ 5277 / 4059 / 1389 Регистрация: 30.07.2012 Сообщений: 12,429
	15.11.2016, 14:38	2
	V_i_p, загляните в самый низ этой страницы в раздел "Похожие темы". Можете найти для себя много полезной информации... 0

@Ирина6 1 / 1 / 0 Регистрация: 19.11.2012 Сообщений: 59
	02.12.2017, 06:22	5
	Ребят, кто может изменить этот код для двух переменных? Очень нужно 0