Updates for R2025a (#40)

Author: isaacito12

BEV/BEV_main_script.html

@@ -1,5 +1,5 @@
 function parent = BEV_ResultsCompactPlot(NameValuePair)
-%% Creates plots of simulation results
+%% Create plots of simulation results
 % To see how this function creates a plot, simply run this function without arguments.
 
 % Copyright 2022-2025 The MathWorks, Inc.
@@ -13,8 +13,8 @@
   % After simulation, create a timetable with extractTimetable.
   NameValuePair.SimData timetable = timetable(seconds([0;10]), [0;0], [0;0], [0;0], [0;0], [0;0], [0;0], [0;0], [0;0], [0;0], ...
     VariableNames = [ ...
-    "Vehicle Speed (km/hr)", ...
-    "Reference Vehicle Speed (km/hr)", ...
+    "Vehicle Speed kph", ...
+    "Reference Vehicle Speed kph", ...
     "G-Force", ...
     "Motor Torque Command", ...
     "Motor Temperature", ...
@@ -57,11 +57,11 @@
 
 ax = nexttile(tl);
 
-plot(ax, simData, "Time", "Vehicle Speed (km/hr)", LineWidth=2)
+plot(ax, simData, "Time", "Vehicle Speed kph", LineWidth=2)
 hold(ax, "on")
 grid(ax, "on")
-plot(ax, simData, "Time", "Reference Vehicle Speed (km/hr)", LineWidth=2)
-setMinimumYRange(ax, simData.("Reference Vehicle Speed (km/hr)"), dy_threshold=2);
+plot(ax, simData, "Time", "Reference Vehicle Speed kph", LineWidth=2)
+setMinimumYRange(ax, simData.("Reference Vehicle Speed kph"), dy_threshold=2);
 ylabel(ax, "")  % Hide variable name defined in the timetable.
 xlim(ax, "tight")
 xlabel(ax, "")

BEV/BEV_main_script.m

@@ -21,19 +21,19 @@
 This is a simple, fast running BEV model which can estimate the electrical efficiency of the vehicle. It is also suitable for further customizations for more focused analysis of individual components at vehicle system level.
 
 
-To open the model, navigate Toolstrip > Project Shortcut tab, and click "**BEV system level**" button.
+To open the model, navigate Toolstrip > Project Shortcut tab, and click the "**BEV model**" button.
 
 
 This script shows an example workflow to programmatically open model, run simulation, collect simulation data, visualize result, save data to text\-format file, read saved data file, and analyze. Feel free to modify this script (probably saving as another file first) and experiment to get new results.
 
 
-You can find more scripts demonstrating other simulation cases in the **SimulationCases** folder.
+You can find more scripts demonstrating other simulation cases in the BEV > Model\-\* > SimulationCases folders.
 
 <a id="H_EC484CEF"></a>
 
 # Run Simulation
 
-This section sets up the model and runs simulation. To run this script at once, navigate Toolstrip > Live Editor tab, and click "Run" button. You can also run this section only by clicking "Run Section" button.
+This section sets up the model and runs simulation. To run this script at once, navigate Toolstrip > Live Editor tab, and click the Run button. You can also run this section only by clicking the Run Section button.
 
 ```matlab
 modelName = "BEV_system_model";
@@ -82,9 +82,8 @@ simOut = sim(modelName);
 simData = extractTimetable(simOut.logsout);
 fig = BEV_ResultsCompactPlot( SimData=simData, PlotTemperature=false );
 % Save the plot to a PNG file.
-prjRoot = currentProject().RootFolder;
 imgFilename = "BEV_SimulationResultPlot.png";
-exportgraphics(fig, fullfile(prjRoot, "BEV", "simulation-results", imgFilename))
+exportgraphics(fig, fullfile(currentProject().RootFolder, "BEV", "simulation-results", imgFilename))
 ```
 
 <center><img src="media/BEV_main_script_media/figure_0.png" width="702" alt="figure_0.png"></center>
@@ -100,23 +99,21 @@ Save logged signals to a CSV file for later analysis. CSV text format is used ra
 simData = extractTimetable(simOut.logsout);
 
 % Adjust time format so as not to lose the subsecond information.
-simData.Time.Format = 'hh:mm:ss.SSSS';
+simData.Time.Format = "hh:mm:ss.SSSS";
 
 % Signal names to save in file.
 % These must be specified in the model as signal logging names.
 % For example, in the BEV system model, see the Measurement subsystem.
 dataColumns = [ ...
   "HV Battery SOC", "HV Battery Power", "HV Battery Current", ...
-  "G-Force", "Vehicle Speed (km/hr)" ];
+  "G-Force", "Vehicle Speed kph" ];
 
 % Select the logged signals to save.
 simData = simData(:, dataColumns);
 
 % Add unit information to signal names.
 varNames = string(simData.Properties.VariableNames');
-% disp(varNames)
 varUnits = string(simData.Properties.VariableUnits');
-% disp(varUnits)
 varNames2 = varNames + " (" + varUnits + ")";
 disp(varNames2)
 ```
@@ -126,15 +123,15 @@ disp(varNames2)
     "HV Battery Power (kW)"
     "HV Battery Current (A)"
     "G-Force (1)"
-    "Vehicle Speed (km/hr) (km/hr)"
+    "Vehicle Speed kph (km/hr)"
 ```
 
 ```matlab
 simData.Properties.VariableNames = varNames2;
 
 % Save data to CSV file.
 simResultFilename  = "BEV_SimulationResult_1.csv";
-simResultFile_FullPath = fullfile(prjRoot, "BEV", "simulation-results", simResultFilename);
+simResultFile_FullPath = fullfile(currentProject().RootFolder, "BEV", "simulation-results", simResultFilename);
 writetimetable(simData, simResultFile_FullPath)
 ```
 
@@ -147,41 +144,31 @@ Open the saved CSV file in text editor and check that the variable names are sav
 This section reads a simulation result CSV file which was saved in the previous section, and do some analysis on the data. This section should work independently without running the previous section as long as the result file exists. At the end of this section, you get the electric efficiency of the vehicle according to the drive cycle for which the data was collected.
 
 ```matlab
-prjRoot = currentProject().RootFolder;
 simResultFilename  = "BEV_SimulationResult_1.csv";
-simResultFile_FullPath = fullfile(prjRoot, "BEV", "simulation-results", simResultFilename);
+simResultFile_FullPath = fullfile(currentProject().RootFolder, "BEV", "simulation-results", simResultFilename);
 
 % Read a CSV file containing simulation result and store it to a timetable.
 data = readtimetable(simResultFile_FullPath, VariableNamingRule="preserve");
 
 % Adjust time format.
-data.Time.Format = 's';
+data.Time.Format = "s";
 
 % Separate variable names and unit strings.
-varNames_with_unit = data.Properties.VariableNames';
-% disp(varNames_with_unit)
+varNames_with_unit = string(data.Properties.VariableNames');
 varNames = extractBefore(varNames_with_unit, " (");
 disp(varNames)
 ```
 
 ```matlabTextOutput
-    {'HV Battery SOC'    }
-    {'HV Battery Power'  }
-    {'HV Battery Current'}
-    {'G-Force'           }
-    {'Vehicle Speed'     }
+    "HV Battery SOC"
+    "HV Battery Power"
+    "HV Battery Current"
+    "G-Force"
+    "Vehicle Speed kph"
 ```
 
 ```matlab
-varUnits = strings(5, 1);
-for idx = 1 : numel(varNames_with_unit)
-  % Return value from extractBetween is an N-by-1 cell array
-  % for each varNames_with_unit where N could be different.
-  % We keep the first element only,
-  % we cannot assign them to varUnits at once outside for-loop.
-  tmpUnits = extractBetween(varNames_with_unit{idx}, "(", ")");
-  varUnits{idx} = tmpUnits{1};
-end
+varUnits = extractBetween(varNames_with_unit, "(", ")");
 disp(varUnits)
 ```
 
@@ -196,10 +183,6 @@ disp(varUnits)
 ```matlab
 data.Properties.VariableNames = varNames;
 data.Properties.VariableUnits = varUnits;
-
-% Define a little utility to manipulate the unit representation
-% in Simscape's physical values.
-tidyUnit = @(x,u) simscape.Value( value(x,u), u );
 ```
 
 Time
@@ -212,50 +195,54 @@ dt = diff(t);
 Travelled Distance
 
 ```matlab
-dataVehSpd = data.("Vehicle Speed");
-unitStr = varUnits(varNames == "Vehicle Speed");
-vehicleSpeed = simscape.Value(dataVehSpd, unitStr{:});
-averageSpeed = sum(vehicleSpeed)/numel(vehicleSpeed)
+dataVehSpd = data.("Vehicle Speed kph");
+unitStr = varUnits(varNames == "Vehicle Speed kph");
+vehicleSpeed = simscape.Value(dataVehSpd, unitStr);
+averageSpeed = sum(vehicleSpeed)/numel(vehicleSpeed);
+disp("Average speed: " + value(averageSpeed) + " " + string(unit(averageSpeed)))
 ```
 
 ```matlabTextOutput
-averageSpeed = 
-   37.2886 (km/hr)
-
+Average speed: 37.2886 km/hr
 ```
 
 ```matlab
-maxSpeed = max(vehicleSpeed)
+maxSpeed = max(vehicleSpeed);
+disp("Maximum speed: " + value(maxSpeed) + " " + string(unit(maxSpeed)))
 ```
 
 ```matlabTextOutput
-maxSpeed = 
-   70.0078 (km/hr)
-
+Maximum speed: 70.0078 km/hr
 ```
 
 ```matlab
 tmpDistance = sum(vehicleSpeed(2:end).*dt);
-travelledDistance = tidyUnit( tmpDistance, "km" )
+travelledDistance = convert( tmpDistance, "km" );
+disp("Travelled distance: " + value(travelledDistance) + " " + string(unit(travelledDistance)))
 ```
 
 ```matlabTextOutput
-travelledDistance = 
-    0.9626 (km)
-
+Travelled distance: 0.96263 km
 ```
 
 
 G Force
 
 ```matlab
 G = data.("G-Force");
-min(G), max(G)
+disp("Minimun G: " + min(G))
+```
+
+```matlabTextOutput
+Minimun G: -0.077098
+```
+
+```matlab
+disp("Maximum G: " + max(G))
 ```
 
 ```matlabTextOutput
-ans = -0.0771
-ans = 0.1985
+Maximum G: 0.19849
 ```
 
 
@@ -264,32 +251,41 @@ Battery Power
 ```matlab
 dataBattPwr = data.("HV Battery Power");
 unitStr = varUnits(varNames == "HV Battery Power");
-batteryPower = simscape.Value(dataBattPwr, unitStr{:});  % J/s
-batteryEnergyUsed = sum(batteryPower(2:end).*dt)
+batteryPower = simscape.Value(dataBattPwr, unitStr);  % J/s
+batteryEnergyUsed = sum(batteryPower(2:end).*dt);
+batteryEnergyUsed = convert(batteryEnergyUsed, "kWh");
+disp("Battery energy used: " + value(batteryEnergyUsed) + " " + string(unit(batteryEnergyUsed)))
 ```
 
 ```matlabTextOutput
-batteryEnergyUsed = 
-  518.4123 (kW*s)
-
+Battery energy used: 0.144 kWh
 ```
 
 ```matlab
-batteryEnergyUsed = tidyUnit(batteryEnergyUsed, "kW*hr")
+energyEfficiency_kWh_per_100km = 100 * value(batteryEnergyUsed / travelledDistance, "kWh/km");
+disp("Energy efficiency: " + energyEfficiency_kWh_per_100km + " kWh per 100 km")
 ```
 
 ```matlabTextOutput
-batteryEnergyUsed = 
-    0.1440 (hr*kW)
+Energy efficiency: 14.9594 kWh per 100 km
+```
 
+```matlab
+energyEfficiency_km_per_kWh = convert(travelledDistance / batteryEnergyUsed, "km/kWh");
+disp("Energy efficiency: " + value(energyEfficiency_km_per_kWh) + " " + string(unit(energyEfficiency_km_per_kWh)))
+```
+
+```matlabTextOutput
+Energy efficiency: 6.6847 km/kWh
 ```
 
 ```matlab
-energyEfficiency_kWh_per_100km = 100 * value(batteryEnergyUsed / travelledDistance, "kW*hr/km")
+energyEfficiency_mi_per_kWh = convert(travelledDistance / batteryEnergyUsed, "mi/kWh");
+disp("Energy efficiency: " + value(energyEfficiency_mi_per_kWh) + " " + string(unit(energyEfficiency_mi_per_kWh)))
 ```
 
 ```matlabTextOutput
-energyEfficiency_kWh_per_100km = 14.9594
+Energy efficiency: 4.1537 mi/kWh
 ```
 
 

BEV/BEV_system_model.mdl

@@ -1,4 +1,4 @@
-Time,HV Battery SOC (%),HV Battery Power (kW),HV Battery Current (A),G-Force (1),Vehicle Speed (km/hr) (km/hr)
+Time,HV Battery SOC (%),HV Battery Power (kW),HV Battery Current (A),G-Force (1),Vehicle Speed kph (km/hr)
 00:00:00.0000,70,0,0,-0,0
 00:00:00.0000,70,0,0,-0,0
 00:00:00.1000,70,0,0,-0,0

BEV/Utility/BEV_ResultsCompactPlot.m

@@ -0,0 +1,22 @@
+function App = BEVProjectNavigationApp()
+%% BEV project navigation app
+% This is a project entry-point app to find some key models, scripts, etc. in the project.
+
+% Copyright 2024-2025 The MathWorks, Inc.
+
+arguments (Output)
+  App (1,:) BEVProjectNavigationAppMain
+end  % arguments
+
+disp("Starting: BEV Project Navigation App")
+
+project_navigator = BEVProjectNavigationAppMain;
+
+project_navigator.Window.HeaderUI.AppSourceName = mfilename;
+
+if nargout > 0
+  % Returning the App variable is optional because the app class object persists
+  % even if it is not returned.
+  App = project_navigator;
+end  % if
+end  % function