Variable step time with series of simulations in simulink - Plotting Difficulties
3 views (last 30 days)
Show older comments
Alfred Kulldorff
on 16 Mar 2020
Answered: Arthur Roué
on 16 Mar 2020
Hallo!
I have made an model where observ diffrent aspect of a PEMFC stack depending on load and temprature.
since each time i run the simulation the number of steps is different. Wich makes the specific Efficiency for that specific current uncompareble with the next simulation.
The only alternative i see is to run with a fixed value but that makes the simulation to take hours... All help with alternative suloutions is appriciated!
i put the code in the bottom.
Thanks in advance!
Here on left i tried to set a specific value like this:
surf(EffM)
xlabel('Current (A)','FontSize',12,'FontWeight','bold')
ylabel('Temp (C)','FontSize',12,'FontWeight','bold')
zlabel('Efficiency (%)','FontSize',12,'FontWeight','bold')
colormap(jet)
and on the right;
figure('Name','Eff LUT')
surf(CurrentM(1,:),tempA(2:Sweep+1),EffM)
xlabel('Current')
ylabel('Temp')
colormap(jet)
The hole code:
Tstart=0; % Starting time for the simulation [s]
Tstop=10;
Sweep=20;
Kelvin=273.15;
PEMTemp=301;
tempA=-50:10:150;
tempA=tempA+Kelvin;
t=length( time );
EffM=zeros(Sweep,t+40);
CurrentM=zeros(10,t+40);
VoltageM=zeros(Sweep,t+40);
H2flowM=zeros(Sweep,t+40);
ttempM=zeros(Sweep,t+40);
ttempM2=zeros(Sweep,11);
for x=1:Sweep
%for y=1:100
%Speed = SpeedA(y);
% PEMTemp44= tempA(x);
PEMTemp= tempA(x);
open_system('PEM_Dynamic_Temp_LUT')
set_param('PEM_Dynamic_Temp_LUT/Constant6', 'value', 'PEMTemp')
set_param('PEM_Dynamic_Temp_LUT', 'SimulationCommand', 'update');
sim('PEM_Dynamic_Temp_LUT',[Tstart,Tstop])
l2=length( Current );
lt=t+40-l2
tt=ones(lt,1);
CurrentM(x, :)= [Current ; tt];
VoltageM(x, :)= [Voltage ; tt];
EffM(x, :)=[Eff ; tt];
H2flowM(x, :)=[H2flow ; tt];
end
0 Comments
Accepted Answer
Arthur Roué
on 16 Mar 2020
You can interpolate your results with interp1 function or a griddedInterpolant object. This way, you'll always have the same points to compare your simulations.
0 Comments
More Answers (0)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!