differential-equations

1熱度

2回答

我試圖讓一個參數根據當前時間根據模型而改變。我設置了一個功能： def func(state,t): x = state[0] y = state[1] z = state[2] Iapp = 0.7 xd = ... + Iapp yd = ... zd = ... 我想有IAPP隨時間變化。也就是說，我需要它是從時

-1熱度

1回答

如何assing係數肯定和實際使用dsolve

所以我要 FullSimplify[DSolve[{iL''[t] + iL[t]/(C*L) == Vc[t]/(C*L*R), C*Vc''[t] + Vc'[t]/R + Vc[t]/L == Vin/L}, {iL[t], Vc[t]},t] {{C, L, R} \[Element] Reals && C >= 0 && L >= 0 && R >= 0}] 符號解決這些方程，我

2熱度

1回答

在Python中模擬具有自適應步長的電子運動 - 微分方程

我想模擬強激光器的振盪電場將如何推動靠近庫侖電位的+1離子電子。在y方向上，激光場是 E = Eo sin（wt）。和庫侖勢是 F =科Q1 * Q2/R^2，在R方向。強電場導致電子到tunnel ionize，所以電子的初始條件是在y方向上從原子上移動。然後，電子被激光場來回推動，並有機會與庫侖電位相互作用。我想模擬庫侖勢能如何影響電子的飛行。模擬需要三維空間，因爲我最終希望包含更復雜的激

0熱度

1回答

集成微分方程在Matlab

我面臨中，我必須解決方程系統的一個問題： \部分A \增量X = F（A）克（B） \局部乙其中f，g，h和k是A或B的一些函數，並且我有一些初始條件A（x = 0）= a，B（（）） y = 0）= b。所有以前的變量都是標量。我已經解決了這個問題，通過在一個變量中做一個小步驟，然後爲另一個變量的整個範圍（如一個ode例程）整合另一個方程，然後在前一個變量的另一個步驟等等。你有另一種想法更有

0熱度

1回答

矩陣指數

所以，我基本上不知道是否對於所述類型的矩陣A （0）（B 0）有一種簡單的方法來計算E 1在，因爲如果它是（0二）（-b 0）一會纔有rotatio n矩陣，其中項的角頻率爲b，正弦和餘弦。因爲獲得Jordan Block和Basis Change的麻煩似乎對於如此簡單的工作來說過於繁瑣。再說，我可能是錯的作爲爲例，假設A = （0 9 （-1 0）如果特徵值+ -3i。我將如何進行

0熱度

1回答

在MATLAB中模擬非線性微分問題？

我被分配了一個程序，我必須創建一個等式和數字的MATLAB模型，如下所示 - http://i.stack.imgur.com/wV0ro.png。不幸的是，我陷入了一段時間。或DH/DT =（-r^2 * SQRT（2 * G * H））/（0.5R + htan（PHI））^ 2 其中 g=386.4 in/s^2 h = 2+(34/64) in r = 1/10 in

0熱度

1回答

表格轉換成矩陣方程

我的表看起來像這樣： {975/2048 - 512 y[1] + 256 y[2]}, {175/128 + 256 y[1] - 512 y[2] + 256 y[3]}, {4095/2048 + 256 y[2] - 512 y[3] + 256 y[4]}, 我想將它轉化成一個矩陣方程： A*y=b 你有一些建議辦呢？謝謝！

0熱度

1回答

實施歐拉解決微分法方程

這個代碼是找到歐拉法MATLAB的 function [x,y]=euler_forward(f,xinit,yinit,xfinal,n) h=(xfinal-xinit)/n; % Initialization of x and y as column vectors x=[xinit zeros(1,n)]; y=[yinit zeros(1,n)]; % Calculation

1熱度

1回答

解微分方程單時間循環利用MATLAB

我有以下公式的機械系統： xdot = Ax+ Bu 我要解決的一個循環這個等式，因爲在我的每一步需要更新û但求解器如ode45或lsim解決時間間隔的微分方程。 for i = 1:10001 if x(i,:)>= Sin1 & x(i,:)<=Sout2 U(i,:) = Ueq - (K*(S/Alpha)) else U(i,:) = Ueq

0熱度

2回答

耦合difeq的數值解

所以我有一套我想解決的3個微分方程。他們可以在我的代碼中看到。我的問題是，我想這些代碼結合起來，這樣我可以有一個用於相對於迴路R（我們將會看到我有什麼： T2 = 1; [T,Y] = ode45(@ball, [0 5*T2] ,[0 0 -10]); figure plot(T,Y(:,1),'-r',T,Y(:,2),'-g',T,Y(:,3),'-b') legend('x(t)