%突变量算法程序
Ts=0.0;Te=0.2;N=200; %仿真起始(start)、结束(end)时间,一周期采样点数
f=50; %仿真频率
T=1/f; %一个周期时长
ts=T/N; %一个周期内采样间隔
ia=out.Ia{1}.Values; %A相电流采样值
ib=out.Ib{1}.Values; %B相电流采样值
iab=ia-ib; %计算得iab
n=1:1:2001; %n可从1连续取至2001
t=(n-1)*ts; %使计算值与时间相对应
Delta_ia=[]; %定义突变量ia为一个数组
Delta_iab_n=[]; %定义突变量iab_n为一个数组
Delta_Iab=[]; %定义突变量Iab为一个数组
for i=401:1:2001 %n从401连续取至2001
Delta_ia(n)=abs(ia(i)-ia(i-N)); %计算得数组ia内的值
Delta_iab_n(n)=abs(iab(i)-iab(i-N)); %计算得数组iab_n内的值
Delta_Iab(n)=abs(Delta_iab_n(i)- Delta_iab_n(i-N)); %计算得数组Iab内的值
end
plot(t, Delta_ia,'b-',t, Delta_iab_n,'g--',t, Delta_Iab,'r-.') ; legend('Delta_ia', 'Delta_iab_n',
'Delta_Iab') %画图
报错:Index exceeds the number of array elements. Index must not exceed 1.