Математический форум Math Help Planet
Обсуждение и решение задач по математике, физике, химии, экономике Теоретический раздел |
Часовой пояс: UTC + 3 часа [ Летнее время ] |
новый онлайн-сервис число, сумма и дата прописью |
|
Часовой пояс: UTC + 3 часа [ Летнее время ] |
Страница 1 из 3 |
[ Сообщений: 21 ] | На страницу 1, 2, 3 След. |
|
Автор | Сообщение | |
---|---|---|
BJIaquMup |
|
|
|
||
Вернуться к началу | ||
VSI |
|
|
BJIaquMup
http://www.math.rsu.ru/mexmat/kvm/MME/dsarch/logist.html Цитата: Простенькую программу... А своими словами можно? |
||
Вернуться к началу | ||
BJIaquMup |
|
|
Ну вот код с матлаба
Код: % x(i+1)=L(k)/(x(i)^2(exp(-x(i))+alpha)) r=130; L=zeros(1,2000); y=zeros(1,r); alpha=0.00729735253765; dL=.007; for k=1:length(L) for i=1:length(y) if (k==1) L(k)=1; y(k,:)=1; else if (i==1) L(k)=L(k-1)+dL; V=y(k-1,r)*y(k-1,r)*(exp(-y(k-1,r))+alpha); y(k,i)=L(k)/V; else V1=y(k,i-1)*y(k,i-1)*(exp(-y(k,i-1))+alpha); y(k,i)=L(k)/V1; end end end end for i=1:r plot(L,y(:,i),'.k') % hold on title('D. Volov'); xlabel('r');ylabel('x'); end У меня что-то ума не хватает переписать его для Математики. В смысле, график не знаю как сделать. |
||
Вернуться к началу | ||
VSI |
|
|
BJIaquMup
Сорри, в MathLab'е я, к сожалению, Может попробуете своими словами обрисовать Вашу проблему? |
||
Вернуться к началу | ||
BJIaquMup |
|
|
VSI писал(а): Сорри, в MathLab'е я, к сожалению, Может попробуете своими словами обрисовать Вашу проблему? Дык тут вроде бы ничего сложного... Вот по этой формуле: [math]x_{n+1}=a*x_n*(1-x_n)[/math] x - по оси ординат, a - по оси абцисс. Хитрость только в том, что x ищется на протяжении где-то штук 100 -- 200 итераций. Потому что вначале значения x хаотичные, а потом всё это дело конденсируется в одной или нескольких точках. Вот их-то (эти точки) и надо отобразить на графике. Блин, чё-то не знаю как на Математике сделать... Вот что должно получиться: |
||
Вернуться к началу | ||
VSI |
|
|
BJIaquMup
Вот попытка "перевести" Вашу программу из MathLAB'а в Mathcad v.14... Вот что получилось (в первом приближении): 1) Все четные (и все нечетные) столбцы массива Ау на графике совпадают друг с другом!!! Поэтому ограничился построением всего 2-х кривулек... (В Mathcad'e в одной графической области можно построить одновременно не более 16 графиков - все равно маловато будет...) 2) В программах при вычислении V вынужден был добавить весьма маленький поправочный коэффициент (т.к. V равнялось 0 и, естественно, при делении на V возникала ошибка... P.S. Эти "красные" и "синие" что-то мне напоминают... (школьные уроки по зоологии..., биологии...) О, вспомнил! Инфузории-туфельки... Точно! |
||
Вернуться к началу | ||
BJIaquMup |
|
|
Программа не моя. Она взята вот отсюда: http://www.sciteclibrary.ru/cgi-bin/yabb2/YaBB.pl?num=1327585698/510. Мне бы такой отродясь не написать.
Здесь вы не правильно закрасили области синей и красной краской. Там пустоты. Но получились две крыски. А то, что заштриховано - это глаза крыски. И каждый "глаз" очередной "крыски" так же виден, как крыска. Млин, но мне надо получить это на Математике. |
||
Вернуться к началу | ||
aspirant2007 |
|
|
вот, по-быстрому накидал
рекурсия в 1000 шагов f[a_, x_] := Nest[a*#*(1 - #) &, x, 1000]; построение графика. как я понял, чтобы получить картинку с бифуркационной диаграммой, подобной приведенной выше, необходимо в координатах (а, х) построить семейство кривых при различных начальных значениях х. я задаю диапазон начальных значений х: {0, 1, 0.1} и визуализирую это семейство Plot[ Table[ f[a, i], {i, 0., 1, 0.1} ], {a, -2, 4}, PlotRange -> All ] Вложение: bfd.jpg [ 8.23 Кб | Просмотров: 2892 ] можно еще посмотреть на это семейство в анимашке: Animate[Plot[f[a, x], {a, -2, 4}], {x, 0, 1}, AnimationRunning -> False] ПС только вот уж не понимаю, почему моя картинка отличается от других (и приведенных тут, и в инете) |
||
Вернуться к началу | ||
BJIaquMup |
|
|
А вы полностью не можете здесь программу привести? Что-то у меня ничего не получается. Вроде бы и прога невелика...
У меня Mathematica 5.0 |
||
Вернуться к началу | ||
aspirant2007 |
|
|
Так я же полностью и привел (для 8 версии)
В 5.2 будут незначительные изменения в аргументе ф-ции Plot В данном примере при определении функции [math]f(a,x)[/math] установил глубину рекурсии 10 (чтобы машина не долго думала) Ну и все Отображаю на графике одновременно десять функций [math]f(a,x)[/math] при различных начальных значениях [math]x[/math] (0, 0.1, 0.2, ... , 1) Отображаются они все одним цветом А вот я их отображаю в цвете (PlotStyle -> Table[Hue[i], {i, 0, 1, 0.1}] - от красного Hue[0] через всю цветовую гамму до красного Hue[1]) |
||
Вернуться к началу | ||
На страницу 1, 2, 3 След. | [ Сообщений: 21 ] |
Похожие темы | Автор | Ответы | Просмотры | Последнее сообщение |
---|---|---|---|---|
Отображение
в форуме Численные методы |
0 |
370 |
12 апр 2014, 06:04 |
|
Отображение
в форуме Функциональный анализ, Топология и Дифференциальная геометрия |
2 |
371 |
28 ноя 2017, 00:28 |
|
Отображение 7->3
в форуме Размышления по поводу и без |
5 |
238 |
01 июл 2019, 23:39 |
|
Отображение | 12 |
296 |
30 апр 2020, 20:54 |
|
Отображение | 1 |
156 |
26 май 2020, 12:43 |
|
Линейное отображение
в форуме Линейная и Абстрактная алгебра |
1 |
268 |
24 ноя 2015, 21:02 |
|
Липшицево отображение
в форуме Функциональный анализ, Топология и Дифференциальная геометрия |
19 |
510 |
17 ноя 2022, 17:07 |
|
Сжимающее отображение
в форуме Функциональный анализ, Топология и Дифференциальная геометрия |
12 |
832 |
23 май 2014, 22:48 |
|
Конформное отображение | 2 |
548 |
29 ноя 2014, 10:51 |
|
Конформное отображение | 2 |
195 |
25 окт 2020, 10:44 |
Часовой пояс: UTC + 3 часа [ Летнее время ] |
Кто сейчас на конференции |
Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 6 |
Вы не можете начинать темы Вы не можете отвечать на сообщения Вы не можете редактировать свои сообщения Вы не можете удалять свои сообщения Вы не можете добавлять вложения |