Моделирование реакций

@marz9n · Регистрация: 02.07.2018

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

Имеется двумерное дискретное пространство.(я предполагаю матрица любой размерности).Имеются три вещества
A,B и C (у каждого из веществ свой уникальный цвет).Концентрация веществ в каждой из точек пространства в следующий момент времени описывается уравнениями:
a(t+Δt) = (a+a*(b-c));
b(t+Δt) = (b+b*(c-a));
c(t+Δt) = (c+c*(a-b));
Для учета диффузии в этих уравнениях в качестве величин a,b и c берется среднее арифметическое от всех точек(собственно точки и восьми соседних) дискретного пространства.
Изначально концентрация веществ берется случайной.
На пространство накладываются периодические граничные условия.

Вот такая вот непростая задачка, на данный момент имею изначальную концентрацию вещества, возникли трудности с созданием секундомера,для использования выше приведенных формул. Ищу помощи, читал самоучитель matlab, но не понял ничего, синтаксис языка очень мне мало понятен, кто может накидать примером с секундомерами или еще что либо связанное с темой буду очень рад

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
N = 100;
figure();
hold on;
for z = 1: N    
    for k = 1: N
        color = [];
        randTouch = ceil((3-0)*rand(1,1));
        switch randTouch
        case 1
            color = 'r';
        case 2 
            color = 'y';
        case 3
            color = 'b';
        end
        scatter(z, k, color,'.');
    end   
end

@nuHrBuH · 21.03.2020, 16:17

Сообщение от marz9n

a(t+Δt) = (a+a*(b-c));
b(t+Δt) = (b+b*(c-a));
c(t+Δt) = (c+c*(a-b));

Получается от времени не зависит

Сообщение от marz9n

На пространство накладываются периодические граничные условия.

Какие

@marz9n · 21.03.2020, 16:40 **[ТС]**

Это я полное задание переписал, [delete]кто знает, что значит граничные условия, но я подозреваю что это ограничение размерности матрицы, на данное время я пока написал вот такой код(я не уверен что код правильный), сейчас думаю, как эта анимация должна проигрываться.

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

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
tic
Matr = 100;
N = ceil(rand*1000);
M = zeros(Matr);
figure();
hold on;
for z = 1: N
    x = ceil(rand*100);
    y = ceil(rand*100);
    for k = 1: N
         if(x > 1)&&(x < 99) && (y > 1) && (y < 99)
        color = [];
        randTouch = ceil((3-0)*rand(1,1));
        M(x,y)=1; % метод диффузии
        switch randTouch
        case 1
            color = 'r';   
            a = ceil(M(x+(-1:1),y+(-1:1)));
           % disp(a);
            scatter(x, y, color,'.');
            break
        case 2 
            color = 'y';
            b = ceil(M(x+(-1:1),y+(-1:1)));
          %   disp(b);
            scatter(x, y, color,'.');
            break
        case 3
            color = 'b';
            c =  ceil(M(x+(-1:1),y+(-1:1)));
          %   disp(c);
            scatter(x, y, color,'.');
            break
        end  
          timerVal = tic;
          a(timerVal) = (a+a*(b-c));
          b(timerVal) = (b+b*(c-a));
          c(timerVal) = (c+c*(a-b));
          disp(timerVal);
       axis([0 100 0 100])
         end
    end     
end
toc

@nuHrBuH · 21.03.2020, 17:57

Раз веществ 3, то и матриц должно быть 3 (для каждого отдельная)
Дискретное пространство в каждой своей точке тоже имеет значение
Среднее значение совсем не используете?

@marz9n · 21.03.2020, 19:07 **[ТС]**

скажи, а отображаться эти 3 матрицы как могут?, накладываться друг на друга или как?

Добавлено через 24 минуты
Да, я проверил, реально можно 3 матрицы наложить друг на друга

@nuHrBuH · 21.03.2020, 19:09

Вещество имеет значение в каждой точке пространства (матрица)
Я бы использовал pcolor или imagesc для отображения
А отображать каждую матрицу отдельно
Или рядом

Matlab M
1
imagesc([A B C])

@marz9n · 21.03.2020, 19:10 **[ТС]**

а где тут можно использовать среднее значение?, подскажите я немного не понимаю, что нужно делать

@nuHrBuH · 21.03.2020, 19:36

Я думаю типа такого

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
for k = 1:10                % 10 итераций
    A = A1;
    for i = 2:N-1           % не доходя до края
        for j = 2:N-1       % не доходя до края
            a1 = A(...
            a = mean(a1(:));
            A1(i,j) = a+a...
        end
    end
    imagesc(A1)
    pause(0.1)
end

A, B и C - матрицы веществ на шаге (i)
A1, B1 и C1 - матрицы веществ на шаге (i+1)

@marz9n · 21.03.2020, 19:53 **[ТС]**

не могу понять, что к чему, можете пожалуйста закомментить каждую строку для более лучшего понимания проходящего?

@nuHrBuH · 22.03.2020, 15:09

Не по теме:

Я проверил свои исследования и приблизился еще на шаг к тайне мироздания

Вот как бог создавал нашу Вселенную из хаоса на поле из 40х40 пикселей

Ветер

И даже любовь

Код

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

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
clc,clear
opengl software
N = 40;
 
[A0,B0,C0] = deal(1/4+zeros(N));
A0(2:end-1,2:end-1) = rand(N-2)*0.9;
B0(2:end-1,2:end-1) = rand(N-2)*0.9;
C0(2:end-1,2:end-1) = rand(N-2)*0.9;
[A1,B1,C1] = deal(A0,B0,C0);
[A2,B2,C2] = deal(A0,B0,C0);
%%
[A1,B1,C1] = deal(A2,B2,C2);
% filename = 'diffusion2.gif';
% f = getframe(gcf); 
% [im,map] = rgb2ind(f.cdata,256);
% imwrite(im,map,filename,'DelayTime',0,'Loopcount',inf);
[Ea,Eb,Ec] = deal(0.4);
h = imagesc([A0 B0 C0]);
colorbar('location','southoutside')
pbaspect([3 1 1]*N)
for k = 1:75
    A0 = A1;
    B0 = B1;
    C0 = C1;
%     A0(820:825) = 0;
%     B0(860:865) = 0;
    for i = 2:N-1
        for j = 2:N-1
            ii = i-1:i+1;
            jj = j-1:j+1;
            
            a = A0(ii,jj);
            b = B0(ii,jj);
            c = C0(ii,jj);
            a1 = mean(a(:));
            b1 = mean(b(:));
            c1 = mean(c(:));
            A1(i,j) = a1+a1*(b1-c1);
            B1(i,j) = b1+b1*(c1-a1);
            C1(i,j) = c1+c1*(a1-b1);
        end
    end
 
    tol = cat(3,A1,B1,C1);
    if max(tol(:)) > 10
        [~,ix] = max(tol(:));
        [y,x] = ind2sub([N 3*N],ix);
        text(x+3,y+3,'темная энергия')
        break
    end
    k
%     Ea = sum(A1(:))/N^2;
%     Eb = sum(B1(:))/N^2;
%     Ec = sum(C1(:))/N^2;
%     subplot(211)
%     hold on
%     plot(k,Ea,'.b',k,Eb,'.r',k,Ec,'.g')
%     hold off
%     subplot(212)
    set(h,'cdata',[A1 B1 C1])
%     f = getframe(gcf); 
%     [im,map] = rgb2ind(f.cdata,256);
%     imwrite(im,map,filename,'DelayTime',0.1,'WriteMode','Append');
    pause(0.02)
end

@marz9n · 23.03.2020, 10:48 **[ТС]**

АФИГЕТЬ, это походу то самое, огромное вам спасибо за помощь, я бы сам не додумался такое реализовать ,вы очень очень сильно меня выручили

Новые блоги и статьи Все статьи Все блоги /
Автозаполнение реквизита при выборе элемента справочника Maks 27.03.2026 Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. При выборе "Спецтехники" (Тип Справочник. Спецтехника), заполняется. . .	Сумматор с применением элементов трёх состояний. Hrethgir 26.03.2026 Тут. https:/ / fips. ru/ EGD/ ab3c85c8-836d-4866-871b-c2f0c5d77fbc Первый документ красиво выглядит, но без схемы. Это конечно не даёт никаких плюсов автору, но тем не менее. . . всё может быть. . .	Автозаполнение реквизитов при создании документа Maks 26.03.2026 Программный код из решения ниже размещается в модуле объекта документа, в процедуре "ПриСозданииНаСервере". Алгоритм проверки заполнения реализован для исключения перезаписи значения реквизита,. . .	Команды формы и диалоговое окно Maks 26.03.2026 1. Команда формы "ЗаполнитьЗапчасти". Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. В качестве источника данных. . .
Кому нужен AOT? DevAlt 26.03.2026 Решил сделать простой ланчер Написал заготовку: dotnet new console --aot -o UrlHandler var items = args. Split(":"); var tag = items; var id = items; var executable = args;. . .	Отправка уведомления на почту при создании или изменении элементов справочника Maks 24.03.2026 Программная отправка письма электронной почты на примере типового справочника "Склады" в конфигурации БП3. Перед реализацией необходимо выполнить настройку системной учетной записи электронной. . .	модель ЗдравоСохранения 5. Меньше увольнений- больше дохода! anaschu 24.03.2026 Теперь система здравосохранения уменьшает количество увольнений. 9TO2GP2bpX4 a42b81fb172ffc12ca589c7898261ccb/ https:/ / rutube. ru/ video/ a42b81fb172ffc12ca589c7898261ccb/ Слева синяя линия -. . .	Midnight Chicago Blues kumehtar 24.03.2026 Такой Midnight Chicago Blues, знаешь?. . Когда вечерние улицы становятся ночными, а ты не можешь уснуть. Ты идёшь в любимый старый бар, и бармен наливает тебе виски. Ты смотришь на пролетающие. . .

@nuHrBuH 483 / 427 / 205 Регистрация: 04.03.2011 Сообщений: 1,259
	21.03.2020, 17:57
	Раз веществ 3, то и матриц должно быть 3 (для каждого отдельная) Дискретное пространство в каждой своей точке тоже имеет значение Среднее значение совсем не используете? 0

@marz9n 0 / 0 / 0 Регистрация: 02.07.2018 Сообщений: 23
	21.03.2020, 19:07 [ТС]
	скажи, а отображаться эти 3 матрицы как могут?, накладываться друг на друга или как? Добавлено через 24 минуты Да, я проверил, реально можно 3 матрицы наложить друг на друга 0

@marz9n 0 / 0 / 0 Регистрация: 02.07.2018 Сообщений: 23
	21.03.2020, 19:10 [ТС]
	а где тут можно использовать среднее значение?, подскажите я немного не понимаю, что нужно делать 0

@marz9n 0 / 0 / 0 Регистрация: 02.07.2018 Сообщений: 23
	21.03.2020, 19:53 [ТС]
	не могу понять, что к чему, можете пожалуйста закомментить каждую строку для более лучшего понимания проходящего? 0

@marz9n 0 / 0 / 0 Регистрация: 02.07.2018 Сообщений: 23
	23.03.2020, 10:48 [ТС]
	АФИГЕТЬ, это походу то самое, огромное вам спасибо за помощь, я бы сам не додумался такое реализовать ,вы очень очень сильно меня выручили 0

Моделирование реакций

Решение