1樓:
大家都在用ode,我就用RK4吧
%%%%%%% MATLAB %%%%%%%%%%%%%% main.m %%%%%%%dt=0.001;t
=0:dt
:100;x
=zeros(3
,length(t
));p=10
;r=28
;b=8
/3;x
(:,1)=
[-8;
8;r-
1];fortn=1
:1:length(t
)-1k1
=Lorenz(x
(:,tn),
p,r,
b);k2=
Lorenz(x
(:,tn)+
dt*k1/
2,p,
r,b);
k3=Lorenz(x
(:,tn)+
dt*k2/
2,p,
r,b);
k4=Lorenz(x
(:,tn)+
dt*k3,
p,r,
b);x(:,tn+
1)=x
(:,tn)+
dt/6*
(k1+2
*k2+2
*k3+k4
);end
plot3(x
(1,:),x(
2,:),x(
3,:),
'-');
boxon
;gridon;
drawnow
;%%%%%%% Lorenz.m %%%%%%%function
xd =
Lorenz
( x, p, r, b )xd=
zeros(3
,1);xd(1
)=-p
*x(1
)+p*
x(2);
xd(2)
=-x(
1)*x
(3)+
r*x(
1)-x
(2);xd(3
)=x(
1)*x
(2)-
b*x(
3);end
2樓:喵小嘰
在matlab裡用ode45解出x,y,z;然後直接用plot畫圖,只要係數是標準洛倫茲係數就沒有問題;
function dy=lorenzf(t,y)%%parameter definition%%p=10;
r=28;
b=(8/3);
dy=zeros(3,1);
dy(1,:)=p*(-y(1)+y(2));
dy(2,:)=(r-y(3))*y(1)-y(2);
dy(3,:)=y(1)*y(2)-b*y(3);
[t,yy]=ode45('lorenzf',[0:0.01:300],[1;1;1]);
x=yy(:,1);
y=yy(:,2);
z=yy(:,3);
figure(1);
plot3(x,y,z);
xlabel('x(t)')
ylabel('y(t)')
zlabel('z(t)')
title('Lorenz')
洛倫茲力永遠不做功嗎,那如何解釋這個問題
冬天真冷 稍偏個題,寫個稍微數學一點的回答吧 電磁學中洛倫茲力的形式是 考慮在元功 有處理一下最後一項,使用Einstein求和約定注意到,由鏈式法則,所以在我們推導出的元功表示式種,倘若固定k而交換i,j下標的取值,levi cevita符號取反,我們立刻得到求和後結果為零。因此 綜上,洛倫茲力做...
動生電動勢本質是磁通量變化導致還是洛倫茲力導致?
豆沙麵包 可以從法拉第定律考慮,磁通為B點乘面積向量S的積分,其導數有兩項,一項為動生,另一項為感生。因此在數學上表示為兩項之和,和為0可以理解為相互抵消。也就是渦旋電場對電荷運動的改變和Lorentz力對電荷運動改變相互抵消,電荷仍保持原有運動。 雪野 磁通量不變只是表面的現象,真正的原因是感生電...
磁流體發電機的非靜電力為什麼是洛倫茲力?
馬小松 假設高溫氣體帶正電荷,很明顯,正電荷是往右側金屬板方向移動,在正電荷移動的過程中,作用於該正電荷的洛倫茲力 磁場力 的方向一直與正電荷速度的方向垂直,不做功。但是偏移引起的電荷分布不均導致燈泡兩端有電勢差,從而使得燈泡有電流通過,燈泡亮了起來。但這個燈泡的能量是誰給的呢?我們假設那個燃燒器瞬...