Курсовая работа на тему: Исследование режима динамического равновесия спектра короткого солитона в неоднородных средах. Выполнила: Студентка 11 -ПМИ-1 Наседкина Татьяна.
Введение Цель работы: Определение множества значений параметров автономной системы дифференциальных уравнений второго порядка, отвечающих замкнутым кривым (траекториям) первого интеграла системы.
Задачи работы: 1. Исследовать уравнение 4 -ой степень. 2. Найти значения параметров уравнения, при которых траектория локализована (имеет хотя бы один положительный корень). 3. Визуализировать фрагментацию плоскости по полученным результатам. 4. Построить фазовые плоскости в различных точках плоскости (локализованных и нелокализованных). 5. Подвести итоги.
Автономная система дифференциальных уравнений второго порядка.
Первый интеграл автономной системы дифференциальных уравнений второго порядка.
Первый интеграл при k=0.
Численная схема. #include<stdio. h> #include<math. h> float func(float n, float a, float b, float c); int main(int argc, char * argv[]) { if(argc!=3) { puts("Wrong format!!!"); return -1; } FILE *fw, *wf; fw = fopen(argv[1], "w");
![if(!fw) { puts("Error create file!!!"); return -3; } wf = fopen(argv[2], "w"); if(!wf) {](https://present5.com/presentation/22155943_175834600/image-8.jpg "if(!fw) { puts(\"Error create file!!!\"); return -3; } wf = fopen(argv[2], \"w\"); if(!wf) {")
if(!fw) { puts("Error create file!!!"); return -3; } wf = fopen(argv[2], "w"); if(!wf) { puts("Error create file!!!"); return -4; } float n, a, c, f=0, x=0, y=0; int k=0; a=-5;
while(a<=5) { n=0; c=-15; while(c<=15) { n=0; while(n<=10) {f=func(n, a, 1 -a, c); if(f<0) x=f; if(f>0) y=f; if(y>0 && x<0) { fprintf(fw, "%fn", a); fprintf(wf, "%fn", c); k++; y=0; x=0; }
n+=0. 02; } if(k==0) { c+=0. 02; } a+=0. 02; } fprintf(fw, "Non"); fprintf(fw, "*****************n"); fprintf(wf, "Non"); fprintf(wf, "*****************n"); } else printf("%dn", k); printf("The endna"); fclose(fw); fclose(wf); return 0; }
![Визуализация плоскости параметров. function [] = Plane () a=load('res_a. txt'); b=load('res_z 0 p. txt');](https://present5.com/presentation/22155943_175834600/image-11.jpg "Визуализация плоскости параметров. function [] = Plane () a=load('res_a. txt'); b=load('res_z 0 p. txt');")
Визуализация плоскости параметров. function [] = Plane () a=load('res_a. txt'); b=load('res_z 0 p. txt'); scatter(a, b, 20); grid on; xlabel('z 0/p'); ylabel('a'); end
Построения фазовых плоскостей
Итоги работы. Проведена фрагментация плоскости параметров с целью выяснения границ динамического равновесия интегральных характеристик солитона.
Спасибо за внимание.
Скачать презентацию Курсовая работа на тему Исследование режима динамического равновесия
курсовая работа.pptx