ex1_1.m

来自「电力系统书本的范例 想学好电力系统很有帮助 欢迎大家来下载」· M 代码 · 共 56 行

% M-file: ex1_10.m
% M-file to calculate and plot the velocity and 
% power of a linear motor as a function of load.
VB = 120;                   % Battery voltage (V)
r = 0.3;                    % Resistance (ohms)
l = 1;                      % Bar length (m)
B = 0.6;                    % Flux density (T)

% Select the forces to apply to the bar
F = 0:10:50;                % Force (N)

% Calculate the currents flowing in the motor.
i = F ./ (l * B);

% Calculate the induced voltages on the bar.
eind = VB - i .* r;

% Calculate the velocities of the bar.
v_bar = eind ./ (l * B);

% Calculate the power supplied by the bar.
power = eind .* i;

% Plot the velocity of the bar versus load.


l1 = 0.45;                  % Length of region 1
l2 = 1.3;                   % Length of region 2
a1 = 0.01;                  % Area of region 1
a2 = 0.015;                 % Area of region 2
ur = 2500;                  % Relative permeability
u0 = 4*pi*1E-7;             % Permeability of free space
n = 200;                    % Number of turns on core
i = 1;                      % Current in amps

% Calculate the first reluctance 
r1 = l1 / (ur * u0 * a1);
disp (['r1 = ' num2str(r1)]);

% Calculate the second reluctance 
r2 = l2 / (ur * u0 * a2);
disp (['r2 = ' num2str(r2)]);

% Calculate the total reluctance
rtot = r1 + r2;

% Calculate the mmf
mmf = n * i;

% Finally, get the flux in the core
flux = mmf / rtot;
b1=flux/a1
b2=flux/a2
% Display result
disp (['Flux = ' num2str(flux)]);