Моделирование движения электрических зарядов в постоянном магнитном поле

@Mia07 · Регистрация: 14.12.2015

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

у меня есть код программы
/////////

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
function f=LorenzForce(t,z)
global H
f=zeros(6,1);
 
f(1) = z(2);
f(2) = 2*pi*(H(3)*z(4)-H(2)*z(6));
f(3) = z(4);
f(4) = -2*pi*(H(3)*z(2)-H(1)*z(6));
f(5) = z(6);
f(6) = 2*pi*(z(2)*H(2)-z(4)*H(1));
 
////////
clear all;
clc;
close all;
 
global H
H=[0 0 1];
x0=0;
y0=0;
z0=0;
vx0=0.1;
vy0=0;
vz0=0.01;
 
[t R] = ode45('LorenzForce',[0:10/1024:10],[x0 vx0 y0 vy0 z0 vz0]);
 
plot3(R(:,1),R(:,3),R(:,5),x0,y0,z0,'o','MarkerSize',6);
axis square; grid on

его нужно редактировать так чтобы вывести график зависимости координат
примерно

Моделирование движения электрических зарядов в постоянном магнитном поле

@SSC · 15.12.2015, 07:41

Уточняйте задачу. Ничего не понятно.
В программе проведено решение 3-х мерного движения 1 объекта в интервале времени 0-10.
На графиках есть переменная Nexp в интервале 0-100 - это что такое?
Есть гистограмма распределения, близкого к нормальному. Распределение чего? По оси X это очень похоже на x(t), но эта величина очень зависит от шага интегрирования. Причем шаг вывода результатов у Вас задан фиксированно, как интервал деленный на 1024, то есть Вы в массиве результатов получаете 1025 точек траектории.

@Mia07 · 15.12.2015, 15:26 **[ТС]**

Изобразить гистограмму частот для выборки координат. Выдвинуть гипотезу H0 о разделе данной выборки и, применяя критерий Пирсона, подтвердить или опровергнуть данную гипотезу.

этот график не для этой программы (это подано как пример вывода результата)

@SSC · 16.12.2015, 09:25

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
function reshenie
clear all;
clc;
close all;
 
global H
H=[0 0 1];
x0=0;
y0=0;
z0=0;
vx0=0.1;
vy0=0;
vz0=0.01;
 
[t R] = ode45(@LorenzForce,[0:10/1024:10],[x0 vx0 y0 vy0 z0 vz0]);
 
plot3(R(:,1),R(:,3),R(:,5),x0,y0,z0,'o','MarkerSize',6);
axis square; grid on
 
end
 
function f=LorenzForce(t,z)
global H
f=zeros(6,1); 
f(1) = z(2);
f(2) = 2*pi*(H(3)*z(4)-H(2)*z(6));
f(3) = z(4);
f(4) = -2*pi*(H(3)*z(2)-H(1)*z(6));
f(5) = z(6);
f(6) = 2*pi*(z(2)*H(2)-z(4)*H(1));
end

Спиральку моделирует, делайте выборку и обрабатывайте

@Mia07 · 16.12.2015, 11:29 **[ТС]**

а можете подсказать как ее делать? я понятия не имею как ее делать

@SSC · 16.12.2015, 13:14

В том то и вопрос - выборку чего?

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

гистограмму частот для выборки координат

Координат x y z, их комбинаций или чего еще? В какой интервал времени, с какой периодичностью?

@Mia07 · 16.12.2015, 18:27 **[ТС]**

задаете флуктуацию магнитного поля, от оптимального значения плюс-минус 20% из функцию
RANDOM с выбранным распределением, например Лоренца, и вычисляете радиусы траектории для каждого компьютерного эксперимента (взять 10000)

@SSC · 17.12.2015, 07:57

С Вами невозможно работать.
1. Общая задача не поставлена.
2. Каждый Ваш ответ, из-за п.1, вместо разъяснения задачи, порождает множество новых вопросов. Из которых самый главный: А является ли представленная Вами модель, адекватной для решения задачи (неизвестной нам)?

@Mia07 · 17.12.2015, 09:14 **[ТС]**

Извините за некорректно поставленные условия.
Если вам это поможет. Эту модель я брал из книги Компьютерное моделирование физических процессов в пакете MATLAB, С. В. Поршнев.

@Ризочка · 02.02.2019, 22:51

ошибку выдает,можете снова отправить

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11680&d=1772460536 Одним из. . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .
SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .

@SSC $Эксперт по математике/физике$ 3390 / 1913 / 571 Регистрация: 09.04.2015 Сообщений: 5,365
	15.12.2015, 07:41
	Уточняйте задачу. Ничего не понятно. В программе проведено решение 3-х мерного движения 1 объекта в интервале времени 0-10. На графиках есть переменная Nexp в интервале 0-100 - это что такое? Есть гистограмма распределения, близкого к нормальному. Распределение чего? По оси X это очень похоже на x(t), но эта величина очень зависит от шага интегрирования. Причем шаг вывода результатов у Вас задан фиксированно, как интервал деленный на 1024, то есть Вы в массиве результатов получаете 1025 точек траектории. 0

@Mia07 0 / 0 / 0 Регистрация: 14.12.2015 Сообщений: 10
	15.12.2015, 15:26 [ТС]
	Изобразить гистограмму частот для выборки координат. Выдвинуть гипотезу H0 о разделе данной выборки и, применяя критерий Пирсона, подтвердить или опровергнуть данную гипотезу. этот график не для этой программы (это подано как пример вывода результата) 0

@Mia07 0 / 0 / 0 Регистрация: 14.12.2015 Сообщений: 10
	16.12.2015, 11:29 [ТС]
	а можете подсказать как ее делать? я понятия не имею как ее делать 0

@Mia07 0 / 0 / 0 Регистрация: 14.12.2015 Сообщений: 10
	16.12.2015, 18:27 [ТС]
	задаете флуктуацию магнитного поля, от оптимального значения плюс-минус 20% из функцию RANDOM с выбранным распределением, например Лоренца, и вычисляете радиусы траектории для каждого компьютерного эксперимента (взять 10000) 0

@SSC $Эксперт по математике/физике$ 3390 / 1913 / 571 Регистрация: 09.04.2015 Сообщений: 5,365
	17.12.2015, 07:57
	С Вами невозможно работать. 1. Общая задача не поставлена. 2. Каждый Ваш ответ, из-за п.1, вместо разъяснения задачи, порождает множество новых вопросов. Из которых самый главный: А является ли представленная Вами модель, адекватной для решения задачи (неизвестной нам)? 0

@Mia07 0 / 0 / 0 Регистрация: 14.12.2015 Сообщений: 10
	17.12.2015, 09:14 [ТС]
	Извините за некорректно поставленные условия. Если вам это поможет. Эту модель я брал из книги Компьютерное моделирование физических процессов в пакете MATLAB, С. В. Поршнев. 0

@Ризочка 0 / 0 / 0 Регистрация: 02.02.2019 Сообщений: 1
	02.02.2019, 22:51
	ошибку выдает,можете снова отправить 0