MFC学习笔记之020matlab混合编程之index out of bounds

MFC调用matlab时，我用两个按钮实现调用matlab计算的，调用的是同一个函数，只是参数不一样而已。问题就来了：第一次调用正确，第二次就不对了。错误是：index out of bounds(索引越界)。
我把m文件的代码贴出来：

C/C++ code

?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64

function [ts,sigma,Os,wc,kpp,tdd]=FindOptimal2(wn,zeta,Pm,sigma0,ts0,choose) %先搜索ts最小的PD控制参数和性能指标 %wn=6;zeta=0.3;Wc=3;Pm=55;sigma0=0.1;ts0=2; %choose是用来选择最优的判定标准：choose=1，ts最小；choose=2，sigma最小 G0=tf(wn^2,[1 2*zeta*wn wn^2 0]); kp_matrix=ones(1,5000);ti_matrix=ones(1,5000);td_matrix=ones(1,5000); Kp_matrix=ones(1,5000);Kd_matrix=ones(1,5000);wc_matrix=[]; Os_matrix=ones(1,5000);ts_matrix=[];sigma_matrix=[];  %Os是振荡次数   j=1;space=0.05;Wc_max=10; for  Wc=0.1:space:Wc_max      [Kp_matrix(j),Ki_matrix(j),Kd_matrix(j),~,~]=OnlyPD(wn,zeta,Wc,Pm,sigma0,ts0);      kp_matrix(j)=Kp_matrix(j);td_matrix(j)=Kd_matrix(j);      Gc=kp_matrix(j)*tf([td_matrix(j) 1],[0 1]);      [sigma_matrix(j),ts_matrix(j),Os_matrix(j)]=StepPara(G0,Gc);      wc_matrix(j)=Wc;      j=j+1; end   %ts最小时性能指标和控制参数 if  choose==1 for  j=1:(Wc_max/space-1)      if  ts_matrix(j)==min(ts_matrix)          k=j;          ts=ts_matrix(j);          sigma=sigma_matrix(k);Os=Os_matrix(k);          wc=wc_matrix(k);kpp=kp_matrix(k);tdd=td_matrix(k);          %ts_min=[ts,sigma,Os,wc,kpp,tdd];      end end end   %sigma最小时的性能指标和控制参数 %用这个作为性能指标时，不能使用前面计算出来的kp td,而要引入Ki，在Kd转td时用到Ki，最后再放弃算出来的ti if  choose==2      j=1;space=0.05;Wc_max=9; for  Wc=0.1:space:Wc_max      [Kp_matrix(j),Ki_matrix(j),Kd_matrix(j),~,~]=OnlyPD(wn,zeta,Wc,Pm,sigma0,ts0);      kp_matrix(j)=(Kd_matrix(j)*Ki_matrix(j)+1)*Kp_matrix(j);       ti_matrix(j)=(Kd_matrix(j)*Ki_matrix(j)+1)/Ki_matrix(j);      td_matrix(j)=Kd_matrix(j)/(Kd_matrix(j)*Ki_matrix(j)+1); %Gc=kp*(1+td*s+1/(ti*s))      Gc=kp_matrix(j)*(tf([td_matrix(j) 1],[0 1]));      [sigma_matrix(j),ts_matrix(j),Os_matrix(j)]=StepPara(G0,Gc);      wc_matrix(j)=Wc;      j=j+1; end            for  j=1:(Wc_max/space-1)          if  sigma_matrix(j)==min(sigma_matrix)              k=j;          ts=ts_matrix(j);          sigma=sigma_matrix(k);Os=Os_matrix(k);          wc=wc_matrix(k);kpp=kp_matrix(k);tdd=td_matrix(k);                        end      end end   %Wc=0.1:space:Wc_max; %plot(Wc,ts_matrix, '-' ,Wc,sigma_matrix, ':' );legend( '随着wc的增大ts变化曲线' , '随着wc的增大sigma变化曲线' );grid on; %plot(Wc,sigma_matrix, ':' );legend( '随着wc的增大sigma变化曲线' ); end

第一次调用是用[ts,sigma,Os,wc,kpp,tdd]=FindOptimal2(6,0.3,55,0.1,5,1)，第二次调用是用[ts,sigma,Os,wc,kpp,tdd]=FindOptimal2(6,0.3,55,0.1,5,2).

谁能指点一下这是为什么？

解答：是自己放了一个超级低级的错误：MFC调用MTATLAB是，输入choose是一个数值我弄成了二维的数组。真是该死。