mirror of
https://github.com/ra3xdh/qucs_s
synced 2025-03-28 21:13:26 +00:00
83 lines
1.6 KiB
Matlab
83 lines
1.6 KiB
Matlab
% synchronous_boost_converter_example.m
|
|
|
|
cd(fileparts(which('synchronous_boost_converter_example.m')));
|
|
|
|
tstart = 0;
|
|
n = 2000;
|
|
tend = 5e-4;
|
|
|
|
%% fixed-step synchronous solver test
|
|
|
|
clear qtr_sync1
|
|
|
|
% qtr_sync1 = synctrcircuit('m-test-2.net');
|
|
qtr_sync1 = synctrcircuit('boostconverter.net');
|
|
|
|
qtr_sync1.init(tstart);
|
|
|
|
% get the number of nodes
|
|
N = qtr_sync1.getn;
|
|
|
|
% get the number of voltage sources
|
|
M = qtr_sync1.getm;
|
|
|
|
% choose some time points
|
|
T1 = linspace(tstart, tend, n)';
|
|
|
|
Y1 = zeros(numel(T1), M+N);
|
|
Y1(1, 1:(N+M)) = qtr_sync1.getsolution();
|
|
|
|
for ind = 2:numel(T1)
|
|
|
|
fprintf(1, 't: %f\n', T1(ind));
|
|
|
|
qtr_sync1.stepsolve(T1(ind));
|
|
|
|
Y1(ind, 1:(N+M)) = qtr_sync1.getsolution();
|
|
|
|
qtr_sync1.acceptstep(T1(ind));
|
|
|
|
% qtr_sync.printx()
|
|
|
|
end
|
|
|
|
figure; plotyy(T1, Y1(:,1:N), T1, Y1(:,((N+1):(N+M))));
|
|
|
|
%% ode test
|
|
|
|
clear qtr_sync2
|
|
|
|
% dbstop in synctrcircuit.m at 203
|
|
% dbstop in synctrcircuit.m at 218
|
|
% dbstop in synctrcircuit.m at 220
|
|
|
|
% qtr_sync2 = synctrcircuit('m-test-2.net');
|
|
qtr_sync2 = synctrcircuit('boostconverter.net');
|
|
|
|
qtr_sync2.init(tstart);
|
|
|
|
y0 = qtr_sync2.getsolution();
|
|
|
|
odeopts = odeset( 'OutputFcn', @synctrcircuit.odeoutputfcn, ...
|
|
'InitialStep', (tend-tstart)/n/2, ...
|
|
... 'MaxStep', (tend-tstart)/n, ...
|
|
'Refine', 1, ...
|
|
'RelTol', 1e-2);
|
|
|
|
[T2, Y2] = ode15s( @synctrcircuit.odefcn, [ tstart+(0.001)*tstart, tend ], y0, odeopts, qtr_sync2);
|
|
|
|
figure; plotyy(T2, Y2(:,1:N), T2, Y2(:,((N+1):(N+M))));
|
|
|
|
|
|
%% comparison
|
|
|
|
|
|
figure; plot(T1, Y1(:,3), 'r', T2, Y2(:,3), 'b');
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|