matlab simscape stepper motor control in matlab

main.m
% Define motor parameters 
R = 3.3;        % Resistance [ohms]
L = 25e-3;      % Inductance [H]
J = 2.8e-6;     % Rotor moment of inertia [kg*m^2]
B = 3.43e-6;    % Viscous friction [N*m*s]
Kt = 0.0274;    % Torque constant [N*m/A]
Ke = 0.0274;    % Voltage constant [V*s/rad]
Tp = 2e-3;      % Pole pitch [m]
NP = 50;        % Number of poles
PM = 1.5;       % Permanent magnet flux [Wb]

% Define simulation parameters
t_stop = 1;     % Simulation time [s]
dt = 1e-6;      % Time step [s]

% Create Simulink model
model = 'stepper_motor';
open_system(model);

% Set up the stepper motor subsystem
subsys = [model '/Stepper Motor'];
set_param(subsys, 'Solver', 'ode14x');

% Set the motor parameters
set_param([subsys '/Motor'], 'R', num2str(R));
set_param([subsys '/Motor'], 'L', num2str(L));
set_param([subsys '/Motor'], 'J', num2str(J));
set_param([subsys '/Motor'], 'B', num2str(B));
set_param([subsys '/Motor'], 'Kt', num2str(Kt));
set_param([subsys '/Motor'], 'Ke', num2str(Ke));
set_param([subsys '/Motor'], 'P', num2str(NP));
set_param([subsys '/Motor'], 'PM', num2str(PM));
set_param([subsys '/Motor'], 'Tp', num2str(Tp));

% Set the simulation time
set_param(model, 'StopTime', num2str(t_stop));

% Run the simulation
sim(model);

% Plot the results
figure;
plot(Theta.Time, Theta.Data);
xlabel('Time [s]');
ylabel('Angle [rad]');
title('Rotor Angle vs. Time');

figure;
plot(Current.Time, Current.Data);
xlabel('Time [s]');
ylabel('Current [A]');
title('Stator Current vs. Time');

figure;
plot(Torque.Time, Torque.Data);
xlabel('Time [s]');
ylabel('Torque [N*m]');
title('Motor Torque vs. Time');
1610 chars
59 lines