如何在m文件中编程调用simulink输出的数组变量

fmx2000

我做了一个小算法，需要用matlab编程验证。学了一段时间编了如下这个程序，有个问题，自己一直解决不了。继续来求救，希望大侠们的指导和帮助！
问题关键是如何在m文件中编程调用pe这个数组变量，其为simulink用ode23tb(变步长)仿真输出存储在workplace 中的double array
function dy=generator(y,pe,Time);
%定义常量
Ws=377;
Tj=6.4;
D=2;
Pm=0.7516;
dy==zeros(2,1); %初始化
dy(1)=[Pm-Pe-D*y(1)]/Tj;
dy(2)=Ws*y(1);

function [y1,y2]=generatorsimu;
y0 = [0; 19.25];           % 取初值
options = optimset('Display','off');
i=1
while i<=size(pe,1)          %仿真输出的pe是大小为723的Array存储在workplace里                        
x = fsolve(@generator, y0, optionspe,Time(:,1),pe(:,1));%如何调用数组pe 及其对应的时间
y1(i)=y(1);
y2(i)=y(2);
i=i+1;
end
plot(k,y1,'-b',k,y2,'-r');
xlabel('peTime')
legend('y1','y2')
simulink仿真时采取的是变步长ode23tb法，仿真得到的pe是大小为723的Array。在用fsolve求解程序中应该如何调用Pe的值，以及如何控制t的步长，与pe对应的时间Time变量一致？