gusucode.com > 双凸轮传动变速器的MATLAB模型 6速双凸轮，带经济燃油评估，控制系统matlab源码程序 > Param_Sweep/DCT_Parameter_Sweep_2D_norm.m

    %% PERFORM PARAMETER SWEEP, NORMAL MODE
% Copyright 2011 The MathWorks(TM), Inc.
mdl = 'Dual_Clutch_Trans';
setup_model_param_sweep2

%% BUILD TARGET
%rtp = Simulink.BlockDiagram.buildRapidAcceleratorTarget(mdl,'AddTunableParameterInfo','on');

%% GENERATE PARAMETER SETS

% QUICK TEST
%mingeardiff = 16:1.25:20;
%rampconst = 0:15:60;

% FOR COMPARING METHODS/SOLVERS
mingeardiff = 16:1:20;
rampconst = 0:5:60;

% DETAILED SWEEP
%mingeardiff = 16:0.2:20;
%rampconst = 0:1:60;

%[SimSettings,rc,mg] = Generate_Sim_Settings(rampconst,mingeardiff,rtp);
pSetNum = 0;
for rc_i = 1:length(rampconst)
    for mgd_i = 1:length(mingeardiff)
        pSetNum = pSetNum+1;
        rc(pSetNum) = rampconst(rc_i);
        mg(pSetNum) = mingeardiff(mgd_i);
    end
end
%numSims = length(SimSettings);
%out = cell(1, numSims);
numSims = pSetNum

% VARIABLE FOR SAVING RESULTS
clear FuelUsedSET

%% START MATLAB POOL
%matlabpool;
%Initialize_MLPool

%% SIMULATE
tic;
for run_i = 1:numSims
    % CALCULATE NEW SHIFT MAP
    [Upshift_Speeds Downshift_Speeds Pedal_Positions] = Calc_Shift_Map_RO(rc(run_i),mg(run_i));
    
    figure(2)
    Plot_Gear_Shift_Schedule(Pedal_Positions, Upshift_Speeds, Downshift_Speeds);
    set(gca,'XLim',[0 100]);
    
    %disp(['Simulating with rampconst = ' num2str(rampconst(run_i))]);
    %out{run_i} = sim(mdl,SimSettings{run_i});
    %disp(['FINISHED with rampconst = ' num2str(rampconst(run_i)) ' Fuel Used = ' num2str(FuelUsedLiters)]);
    sim(mdl)
    FuelUsedSET_v(run_i) = FuelUsedLiters;
end
Elapsed_Sim_Time = toc;
disp(['Elapsed Sim Time = ' num2str(Elapsed_Sim_Time)]);

%% PLOT RESULTS
out_i = 0;
for rc_i=1:length(rampconst)
    for mg_i=1:length(mingeardiff)
        out_i = out_i+1;
        %FuelUsedSET(mg_i,rc_i) = out{out_i}.find('FuelUsedLiters');
        FuelUsedSET(mg_i,rc_i) = FuelUsedSET_v(out_i);
    end
end

figure;
set(gcf,'Position',[104   78   560   420]);
[rc_array mgd_array] = meshgrid(rampconst,mingeardiff);
mesh(mgd_array,rc_array,FuelUsedSET);
%surf(mgd_array,rc_array,FuelUsedSET);
title('Fuel Use With Varying Shift Schedule','FontWeight','Bold','FontSize',14);
xlabel('Offset Parameter','FontWeight','Bold','FontSize',12);
ylabel('Ramp Parameter','FontWeight','Bold','FontSize',12);
zlabel('Fuel Used (L)','FontWeight','Bold','FontSize',12);
legend('Parameter Sweep','Location','NorthWest');

%save FuelUsedSET_Sweep2D_AbsDes FuelUsedSET rampconst mingeardiff
%clear rampconst mingeardiff FuelUsedSET


%% SAVE RESULTS
trans_var = get_param([mdl '/Dual Clutch Transmission'],'BlockChoice');
veh_var = get_param([mdl '/Vehicle Dynamics'],'BlockChoice');
solver_setting = get_param(mdl,'Solver');
if (strcmp(solver_setting,'ode23t'))
    solver_var = 'Desktop';
else
    solver_var = 'Real Time';
end

SaveVarName = ['FuelUsedSweep2D_' trans_var(1:3) solver_var(1:3) '_' num2str(numel(FuelUsedSET))];
eval([SaveVarName '.FuelUsedSET = FuelUsedSET;']);
eval([SaveVarName '.rampconst = rampconst;']);
eval([SaveVarName '.mingeardiff = mingeardiff;']);
eval([SaveVarName '.Elapsed_Sim_Time = Elapsed_Sim_Time;']);
eval([SaveVarName '.Variants.Transmission = ''' trans_var ''';']);
eval([SaveVarName '.Variants.Vehicle = ''' veh_var ''';']);
eval([SaveVarName '.Variants.SolverSettings = ''' solver_var ''';']);

save(SaveVarName, SaveVarName);

clear rampconst mingeardiff FuelUsedSET

%% RESET MODEL AND CLOSE POOL
%matlabpool close
reset_model_param_sweep