zfhgsedbfvsehkbnvcuizdh

%{
% Obwód magnetyczny liniowy przykład 1
%}

% stałe
m0 = 4*pi*10^-7;

% zasilania

I1 = 2; %Amper
z1 = 30; %Ohm

I3 = 5; %Amper
z3 = 100; %Ohm

% Odcinek l1 - blacha transformatorowa

l1 = 0.03;      % długość [m]
s1 = 2e-4;      % powierzchnia przekroju [m^2]
mr1 = 600;      % przenikalność magnetyczna względna

% Odcinek l2 - miedź

l2 = 0.05;      % długość [m]
s2 = 2e-4;      % powierzchnia przekroju [m^2]
mr2 = 0.9999;   % przenikalność magnetyczna względna

% Odcinek l3 - stop Ni+Fe

l3 = 0.1;           % długość [m]
s3 = 2e-4;      % powierzchnia przekroju [m^2]
mr3 = 2700;     % przenikalność magnetyczna względna

% Reluktancje

Rm1 = l1/(m0*mr1*s1);   % reluktancja odcinka l1
Rm2 = l2/(m0*mr2*s2);   % reluktancja odcinka l2
Rm3 = l3/(m0*mr3*s3);   % reluktancja odcinka l3

% z magnetycznych praw Kirchhoffa
A = [
    1   -1      -1;
    Rm1 Rm2     0;
    0       -Rm2    Rm3
];

b = [
    0;
    I1*z1;
    -I3*z3
];

% obliczamy rozpływ strumieni indukcji magnetycznej
phi = linsolve(A, b);