couple.m

来自「基于MATLAB电机磁场程序」· M 代码 · 共 22 行

tic; clear;clc; close all;
aa=22.5*pi/180; %aa=aa+pi/2;
pas=0.2
nom=3*pi/pas
nom=round(nom)+1
nnn=1:nom
ener=zeros(1,nom)
for s=1:nom
    MsP10
    aa=aa+pas
    energie1
    ener(1,s)=energ_total
end
coup(1,1)=0
for s=2:nom
coup(1,s)=ener(1,s)-ener(1,s-1)
coup(1,s)=coup(1,s)/pas
end
subplot(2,1,1),plot(nnn,ener,'r'),title('energie'),ylabel('energie magnetique stockee'),xlabel('k'),hold on
subplot(2,1,2),plot(nnn,coup,'b'),title('couple'),ylabel('couple'),xlabel('k')
toc