【北太天元軟件上做RLC電路模擬】 https://www.bilibili.com/video/BV15P411M78h/?share_source=copy_web&;vd_source=2adc5aa7a702b808eb8b31dbd210f954

在北太天元數(shù)值計(jì)算通用軟件的代碼如下:?

% 北太天元模擬RLC 電路

% RLC 電路 包含 電容C、電感L, ?和 電阻 R

% 在t =0 的時(shí)候打開開關(guān)， 電源是直流供電

%

% 常微分方程:

% ?Vc '' (t) ?+ \frac{R}{L} Vc '(t) + \frac{1}{LC} Vc(t) = \frac{1}{LC} Vs(t)

% ?這里的 Vs(t) = 5?

%

% ?我們首先把二階常微分方程寫成一階常微分方程組

% ?引入 y(1) = Vc(t) ; y(2) = Vc'(t);?

% ?d ( y(1) ) / dt ?= ?y(2);

% ?d ( y(2) ) / dt = - R/L y(2) ?- 1/(L*C)*y(1) ?+ 1/(L*C) *Vs

Vs = 5 ; ?%直流電源的電壓?

R = 1; % ?電阻

C = 0.001; % 電容

L = 0.1; ?% 電感?

tspan = [0 1]; % 計(jì)算的時(shí)間區(qū)間

y0 = [0; 0]; %初值

odefun = @(t,y) [ y(2); ?-R/L*y(2)-1/(L*C)*y(1)+1/(L*C)*Vs];

options = odeset('RelTol', 1e-3, 'AbsTol', 1e-4);

[t, y] = ode45(odefun, tspan, y0, options);

plot(t, y(:,1), 'r-', 'LineWidth', 3);

xlabel("時(shí)間")

ylabel("Vc");

title("RLC 電路的響應(yīng)Vc (也就是電容上的電壓差)")

text(0.6,0.4, "北太天元數(shù)值計(jì)算通用軟件",'FontSize',24)