- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
您好,我正在编写一个小程序,它可以在 x 平面中使用 for 循环创建许多不同的粒子,这些粒子在各种力的作用下沿 y 方向向上移动。粒子以排斥力相互作用,因此在 ode45 计算中需要每个粒子相对于彼此的位置。我已经多次查看代码,但我不确定问题出在哪里,因为解决方案显示为“NaN”。时间间隔是为了及时打散粒子的运动,以便在ode45计算中更新所有粒子的位置。我用事件来打发时间。如果您需要更多信息,因为代码不是很清楚,请询问!我没有在代码中显示常量:
time_interval = 1000; %Time in microseconds
%The number of time intervals the particle movement is divided by
number_of_intervals = 10000/time_interval;
tstart=0;
tstop= time_interval*1e-6;
num_of_particles = 3;
w0{1} = [0;0;0;0];
sols{1} = transpose(w0{1});
x_adj_stable(1) = 0;
y_adj_stable(1) = 0;
break_vect(num_of_particles) = 0;
times{1} = 0;
%Initiate particle characteristics
for particle_num = 2:num_of_particles
w0{particle_num} = [0;0;w0{particle_num - 1}(3) + 4*1e-6;0];
x_adj_stable(particle_num) = w0{particle_num}(3);
y_adj_stable(particle_num) = 0;
sols{particle_num} = transpose(w0{particle_num});
times{particle_num} = 0;
end
%Event options, Identifying the time intervals within the solution
options = odeset('Events', @events);
%The loop which connects the intervals of time
for n = 1 : number_of_intervals
for particle_num = 1:num_of_particles
%===========================================
[t_sol, y_sol] = ode45(@pm1dwoc, [tstart, tstop], w0{particle_num}, options);
%===========================================
%The number of steps within each interval
n_steps = length(t_sol);
for interval = 2: n_steps
times{particle_num} = [times{particle_num}; t_sol(interval)];
sols{particle_num} = [sols{particle_num}; y_sol(interval, :)];
end
%Setting the initial conditions of the next interval
w0{particle_num} = y_sol(n_steps,:);
y_adj(particle_num) = y_sol(n_steps, 1);
x_adj(particle_num) = y_sol(n_steps, 3);
end
tstart = tstart + time_interval*1e-6;
tstop = tstart + time_interval*1e-6;
for n = 1:num_of_particles
y_adj_stable(n) = y_adj(n);
x_adj_stable(n) = x_adj(n);
end
end
function dwdt = pm1dwoc (t,w)
y=w(1);
vy=w(2);
x=w(3);
vx=w(4);
for particle_number = 1:num_of_particles
Dist(particle_number) = sqrt((x - x_adj_stable(particle_number))^2 + (y - y_adj_stable(particle_number))^2);
if (Dist(particle_number) > 0)
x_vect(particle_number) = (x - x_adj_stable(particle_number))/Dist(particle_number);
y_vect(particle_number) = (y - y_adj_stable(particle_number))/Dist(particle_number);
end
end
dy_component = -qw*Vd/(m*G) + qw*(-qw)/(m*4*pi*epsilon0*(2*R+2*y)^2));
dx_component = 0;
for particle_number = 1:num_of_particles
if (Dist(particle_number) > 0)
dx_component = dx_component + ((qw^2)/(0.0001*m*4*pi*epsilon0*Dist(particle_number)^2))*x_vect(particle_number);
dy_component = dy_component + ((qw^2)/(0.0001*m*4*pi*epsilon0*Dist(particle_number)^2))*y_vect(particle_number);
end
end
dwdt = [vy; dy_component; vx; dx_component];
end
function [eventvalue,stopthecalc,eventdirection] = events (t,w)
t= t*1e6;
end_time = 10000;
eventvalue = t - end_time/number_of_intervals;
stopthecalc = 1;
eventdirection = 1;
for num = 2:number_of_intervals
eventvalue = [eventvalue; t - num*(end_time/number_of_intervals)];
stopthecalc = [stopthecalc; 1];
eventdirection = [eventdirection; 1];
end
end
end
非常感谢任何帮助!
谢谢
山姆
编辑
我将发布一个简化版本,因为我删除了我知道没有任何错误的部分:
num_of_particles = 3;
w0{1} = [0;0;0;0];
sols{1} = transpose(w0{1});
x_adj_stable(1) = 0;
y_adj_stable(1) = 0;
times{1} = 0;
%Initiate particle characteristics
for particle_num = 2:num_of_particles
w0{particle_num} = [0;0;w0{particle_num - 1}(3) + 4*1e-6;0];
x_adj_stable(particle_num) = w0{particle_num}(3);
y_adj_stable(particle_num) = 0;
sols{particle_num} = transpose(w0{particle_num});
times{particle_num} = 0;
end
for particle_num = 1:num_of_particles
[t_sol, y_sol] = ode45(@pm1dwoc, [tstart, tstop], w0{particle_num});
%The number of steps within each interval
n_steps = length(t_sol);
for interval = 2: n_steps
times{particle_num} = [times{particle_num}; t_sol(interval)];
sols{particle_num} = [sols{particle_num}; y_sol(interval, :)];
end
%Setting the initial conditions of the next interval
w0{particle_num} = y_sol(n_steps,:);
y_adj(particle_num) = y_sol(n_steps, 1);
x_adj(particle_num) = y_sol(n_steps, 3);
end
%//This part is written since the y_adj and x_adj variables are constantly changing so %//stable variable version is required when used in the ode45 calculations
for n = 1:num_of_particles
y_adj_stable(n) = y_adj(n);
x_adj_stable(n) = x_adj(n);
end
end
function dwdt = pm1dwoc (t,w)
y=w(1);
vy=w(2);
x=w(3);
vx=w(4);
for particle_number = 1:num_of_particles
Dist(particle_number) = sqrt((x - x_adj_stable(particle_number))^2 + (y - y_adj_stable(particle_number))^2);
if (Dist(particle_number) > 0)
x_vect(particle_number) = (x - x_adj_stable(particle_number))/Dist(particle_number);
y_vect(particle_number) = (y - y_adj_stable(particle_number))/Dist(particle_number);
end
end
dy_component = 0;
dx_component = 0;
for particle_number = 1:num_of_particles
if (Dist(particle_number) > 0)
dx_component = dx_component + (1/(Dist(particle_number))*x_vect(particle_number);
dy_component = dy_component + (1/(Dist(particle_number))*y_vect(particle_number);
end
end
dwdt = [vy; dy_component; vx; dx_component];
end
最佳答案
一般来说,您不想分别对许多方程中的每一个进行积分,尤其是如果方程是耦合的。 ode45
及其亲戚完全能够对耦合方程组进行积分。
作为一个简单的示例,您可能有以下函数,它定义了三个变量的导数。衍生品相互依赖。
function xprime = my_ode(t, x)
xprime = [x(1) + 2*x(3);
x(2) - x(1);
3*x(3) + x(1) + x(2)];
这不会生成一个非常有趣的系统,但它展示了如何一次性集成三个不同的变量。使用 ode45
,您可以键入:
x0 = [0;0;1];
tspan = 0:0.1:1;
[T,X] = ode45(@my_ode, tspan, x0)
对于您的问题,此结构允许您计算每一步中粒子之间的分离并确定 react 。然后,您可以将这些效果作为运动方程的一部分,而不是对每个效果进行积分和调整。
关于MATLAB:在 for 循环中使用 ode45,使用动态变量进行粒子运动和交互,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/13717385/
我是 PHP 新手。我一直在脚本中使用 for 循环、while 循环、foreach 循环。我想知道 哪个性能更好? 选择循环的标准是什么? 当我们在另一个循环中循环时应该使用哪个? 我一直想知道要
我在高中的编程课上,我的作业是制作一个基本的小计和顶级计算器,但我在一家餐馆工作,所以制作一个只能让你在一种食物中读到。因此,我尝试让它能够接收多种食品并将它们添加到一个价格变量中。抱歉,如果某些代码
这是我正在学习的一本教科书。 var ingredients = ["eggs", "milk", "flour", "sugar", "baking soda", "baking powder",
我正在从字符串中提取数字并将其传递给函数。我想给它加 1,然后返回字符串,同时保留前导零。我可以使用 while 循环来完成此操作,但不能使用 for 循环。 for 循环只是跳过零。 var add
编辑:我已经在程序的输出中进行了编辑。 该程序要求估计给定值 mu。用户给出一个值 mu,同时还提供了四个不等于 1 的不同数字(称为 w、x、y、z)。然后,程序尝试使用 de Jaeger 公式找
我正在编写一个算法,该算法对一个整数数组从末尾到开头执行一个大循环,其中包含一个 if 条件。第一次条件为假时,循环可以终止。 因此,对于 for 循环,如果条件为假,它会继续迭代并进行简单的变量更改
现在我已经习惯了在内存非常有限的情况下进行编程,但我没有答案的一个问题是:哪个内存效率更高;- for(;;) 或 while() ?还是它们可以平等互换?如果有的话,还要对效率问题发表评论! 最佳答
这个问题已经有答案了: How do I compare strings in Java? (23 个回答) 已关闭 8 年前。 我正在尝试创建一个小程序,我可以在其中读取该程序的单词。如果单词有 6
这个问题在这里已经有了答案: python : list index out of range error while iteratively popping elements (12 个答案) 关
我正在尝试向用户请求 4 到 10 之间的整数。如果他们回答超出该范围,它将进入循环。当用户第一次正确输入数字时,它不会中断并继续执行 else 语句。如果用户在 else 语句中正确输入数字,它将正
我尝试创建一个带有嵌套 foreach 循环的列表。第一个循环是循环一些数字,第二个循环是循环日期。我想给一个日期写一个数字。所以还有另一个功能来检查它。但结果是数字多次写入日期。 Out 是这样的:
我想要做的事情是使用循环创建一个数组,然后在另一个类中调用该数组,这不会做,也可能永远不会做。解决这个问题最好的方法是什么?我已经寻找了所有解决方案,但它们无法编译。感谢您的帮助。 import ja
我尝试创建一个带有嵌套 foreach 循环的列表。第一个循环是循环一些数字,第二个循环是循环日期。我想给一个日期写一个数字。所以还有另一个功能来检查它。但结果是数字多次写入日期。 Out 是这样的:
我正在模拟一家快餐店三个多小时。这三个小时分为 18 个间隔,每个间隔 600 秒。每个间隔都会输出有关这 600 秒内发生的情况的统计信息。 我原来的结构是这样的: int i; for (i=0;
这个问题已经有答案了: IE8 for...in enumerator (3 个回答) How do I check if an object has a specific property in J
哪个对性能更好?这可能与其他编程语言不一致,所以如果它们不同,或者如果你能用你对特定语言的知识回答我的问题,请解释。 我将使用 c++ 作为示例,但我想知道它在 java、c 或任何其他主流语言中的工
这个问题不太可能帮助任何 future 的访问者;它只与一个小的地理区域、一个特定的时间点或一个非常狭窄的情况有关,这些情况并不普遍适用于互联网的全局受众。为了帮助使这个问题更广泛地适用,visit
我是 C 编程和编写代码的新手,以确定 M 测试用例的质因数分解。如果我一次只扫描一次,该功能本身就可以工作,但是当我尝试执行 M 次时却惨遭失败。 我不知道为什么 scanf() 循环有问题。 in
这个问题已经有答案了: JavaScript by reference vs. by value [duplicate] (4 个回答) 已关闭 3 年前。 我在使用 TSlint 时遇到问题,并且理
我尝试在下面的代码中添加 foreach 或 for 循环,以便为 Charts.js 创建多个数据集。这将允许我在此折线图上创建多条线。 我有一个 PHP 对象,我可以对其进行编码以稍后填充变量,但
我是一名优秀的程序员,十分优秀!