从数据中提取电池参数

%的电池数据加载负载ee_battery_data.matassignin (“基地”,“T1”battery_data(找到([battery_data (1:2) .T] = = 25)) .T);assignin (“基地”,《终结者2》battery_data(找到([battery_data (1:2) .T] ~ = 25)) .T);%显示电池模型模型=“ee_battery”;open_system(模型)

close_system(模型中,0);

ParsListMain = {“Vnom”,R1的,“啊”,“V1”,“AH1”};InitGuessMain = (3.6, 0.045, 2.7, 3.4, 1.4);ParsListDyn = {“一国”,“tau1”};InitGuessDyn = (0.006, 200);ParsListTemp = {“Vnom_T2”,“R1_T2”,“V1_T2”,“Rp1_T2”,“tau1_T2”};InitGuessTemp = [3.8, 0.055, 3.6, 0.006, 200);Pars0 =重塑([[ParsListMain ParsListDyn ParsListTemp);cellstr (num2str ([InitGuessMain InitGuessDyn InitGuessTemp)))), 1, []);流(' \ t % 5 s = % s \ n ',Pars0 {:});清晰的Pars0

load_system(模型);%启用快速启动加速仿真set_param(模型,“FastRestart”,“上”)Pars =重塑([ParsListMain;cellstr (num2str (InitGuessMain ') '], 1, []);为k = 1:2:长度(Pars) evalin (“基地”,(Pars {k}“=”帕尔斯{k + 1}“;”])结束帕尔斯=重塑([ParsListDyn;cellstr (num2str (InitGuessDyn ') '], 1, []);为k = 1:2:长度(Pars) evalin (“基地”,(Pars {k}“=”帕尔斯{k + 1}“;”])结束帕尔斯=重塑([ParsListTemp;cellstr (num2str (InitGuessTemp ') '], 1, []);为k = 1:2:长度(Pars) evalin (“基地”,(Pars {k}“=”帕尔斯{k + 1}“;”])结束%生成初步模型曲线和暗算数据num_lines =长度(battery_data) 1;v_model =细胞(1、num_lines);t_model =细胞(1、num_lines);legend_info_data =细胞(1、num_lines);legend_info_model =细胞(1、num_lines);为idx_data = 1: num_lines assignin (“基地”,“t_data”battery_data (idx_data) .t);assignin (“基地”,“i_data”battery_data (idx_data)。我);assignin (“基地”,“T_data”battery_data (idx_data) .T *的(长度(t_data), 1));assignin (“基地”,“T0”battery_data (idx_data) .T);assignin (“基地”,“t”t_data (2) -t_data (1));assignin (“基地”,“AH0”,啊* battery_data (idx_data) .SOC0);= sim(模型);v_model {idx_data} = out.Vo.signals.values;t_model {idx_data} = out.Vo.time;legend_info_data {idx_data} = [“临时= '…num2str (battery_data (idx_data) .T)“\ circC、数据”];legend_info_model {idx_data} = [“临时= '…num2str (battery_data (idx_data) .T)“\ circC模式”];结束情节([battery_data (1: num_lines)。t]/3600, [battery_data(1:num_lines).v],“o”,(t_model {}): / 3600, [v_model{}):)包含(的时间(小时));ylabel (“电池电压(V)”);传奇([legend_info_data legend_info_model),“位置”,“最佳”);标题(“模型与初始参数值”);

%找到环境温度数据索引idx_data =找到([battery_data (1: num_lines) .T] = = 25);assignin (“基地”,“t_data”battery_data (idx_data) .t);assignin (“基地”,“i_data”battery_data (idx_data)。我);assignin (“基地”,“T_data”battery_data (idx_data) .T *的(长度(t_data), 1));assignin (“基地”,“T0”battery_data (idx_data) .T);assignin (“基地”,“t”t_data (2) -t_data (1));assignin (“基地”,“v_data”battery_data (idx_data) .v);%优化电池的主要参数对话框选项卡assignin (“基地”,“ParsList”ParsListMain (1:4));InitGuess = InitGuessMain (1:4);OptPars = fminsearch (@ee_battery_lse InitGuess,…optimset (“TolX”1 e - 3));OptParsMain = [OptPars (1:4) InitGuessMain (5)];%更新电池块与优化参数帕尔斯=重塑([ParsListMain;cellstr (num2str (OptParsMain ') '], 1, []);为k = 1:2:长度(Pars) evalin (“基地”,(Pars {k}“=”帕尔斯{k + 1}“;”])结束%显示优化参数流([电池主要优化参数的…“对话框选项卡:\ n”]);流(' \ t % 5 s = % s \ n 'Pars {:});清晰的i_datav_datat_dataT_dataTs清晰的kInitGuess

%只使用一个电流脉冲优化电荷动力学i_pos = battery_data(1)我。* (battery_data (1), > = 0);一个=找到(diff (i_pos) > 0, 2);b =找到(diff (battery_data(1)我));c =修复((b(找到(b < (1), 1,“最后一次”))+ (1))/ 2);assignin (“基地”,“i_data”battery_data (idx_data),我(c + 1: (2)));assignin (“基地”,“v_data”,battery_data (1) .v (c + 1: (2)));assignin (“基地”,“t_data”battery_data (idx_data) .t(1:长度(i_data)));assignin (“基地”,“T_data”battery_data (idx_data) .T *的(长度(t_data), 1));assignin (“基地”,“T0”battery_data (idx_data) .T);assignin (“基地”,“t”t_data (2) -t_data (1));%之前找到电池初始充电优化电荷动力学参数assignin (“基地”,“ParsList”,{“充电”});InitGuessCharge = OptParsMain (3);OptCharge = fminsearch (@ee_battery_lse InitGuessCharge,…optimset (“TolX”1 e - 3));assignin (“基地”,“AH0”,OptCharge);%优化电池充电动力学参数assignin (“基地”,“ParsList”,(ParsListMain (2) ParsListDyn]);InitGuessDyn = [OptPars (2) InitGuessDyn];OptParsDyn = fminsearch (@ee_battery_lse InitGuessDyn,…optimset (“TolX”1 e - 3));%更新电池块优化电荷动力学参数ParsListMainDyn = [ParsListMain ParsListDyn];OptParsMainDyn = [OptPars (1) OptParsDyn (1) OptPars (3:4) InitGuessMain (5) OptParsDyn (2:3)];帕尔斯=重塑([ParsListMainDyn;cellstr (num2str (OptParsMainDyn ') '], 1, []);为k = 1:2:长度(Pars) evalin (“基地”,(Pars {k}“=”帕尔斯{k + 1}“;”])结束assignin (“基地”,“AH0”,啊* battery_data (idx_data) .SOC0);%显示优化参数流([电池的优化参数,…“包括电荷动力学:\ n”]);流(' \ t % 5 s = % s \ n 'Pars {:});清晰的i_datav_datat_dataT_dataTs清晰的i_pos一个bc清晰的k清晰的OptParsOptParsDynParsListMainDynInitGuessMainInitGuessDynParsListDynParsListMain

idx_data =找到([battery_data (1: num_lines) .T] ~ = 25);assignin (“基地”,“t_data”battery_data (idx_data) .t);assignin (“基地”,“i_data”battery_data (idx_data)。我);assignin (“基地”,“T_data”battery_data (idx_data) .T *的(长度(t_data), 1));assignin (“基地”,“T0”battery_data (idx_data) .T);assignin (“基地”,“t”t_data (2) -t_data (1)) assignin (“基地”,“v_data”battery_data (2) .v);%使用初始参数为T1, T2参数InitGuessTemp = [OptParsMainDyn (1:2) OptParsMainDyn (4) OptParsMainDyn(者);帕尔斯=重塑([ParsListTemp;cellstr (num2str (InitGuessTemp ') '], 1, []);为k = 1:2:长度(Pars) evalin (“基地”,(Pars {k}“=”帕尔斯{k + 1}“;”])结束%优化电池温度相关参数assignin (“基地”,“ParsList”ParsListTemp (1:3));OptParsTemp = fminsearch (@ee_battery_lse InitGuessTemp (1:3),…optimset (“TolX”1 e - 3));%更新电池块与优化参数帕尔斯=重塑([ParsListTemp (1:3);cellstr (num2str (OptParsTemp ') '], 1, []);为k = 1:2:长度(Pars) evalin (“基地”,(Pars {k}“=”帕尔斯{k + 1}“;”])结束assignin (“基地”,“AH0”,啊* battery_data (idx_data) .SOC0);%显示优化参数流([优化电池的温度依赖参数的…“是:\ n”]);流(' \ t % 5 s = % s \ n 'Pars {:});清晰的i_datav_datat_dataT_dataTs清晰的k清晰的OptParsMainDyn

%寻找non-room温度指数数据%只使用一个电流脉冲优化电荷动力学i_pos = battery_data (idx_data)。我。* (battery_data (idx_data), > = 0);一个=找到(diff (i_pos) > 0, 2);b =找到(diff (battery_data (idx_data)。我));c =修复((b(找到(b < (1), 1,“最后一次”))+ (1))/ 2);assignin (“基地”,“i_data”battery_data (idx_data),我(c + 1: (2)));assignin (“基地”,“v_data”battery_data (idx_data) .v (c + 1: (2)));assignin (“基地”,“t_data”battery_data (idx_data) .t(1:长度(i_data)));assignin (“基地”,“T_data”battery_data (idx_data) .T *的(长度(t_data), 1));assignin (“基地”,“T0”battery_data (idx_data) .T);assignin (“基地”,“t”t_data (2) -t_data (1))%之前找到电池初始充电优化电荷动力学参数assignin (“基地”,“ParsList”,{“充电”});InitGuessCharge = OptParsMain (3);OptCharge = fminsearch (@ee_battery_lse InitGuessCharge,…optimset (“TolX”1 e - 3));assignin (“基地”,“AH0”,OptCharge);%优化电池充电动力学参数assignin (“基地”,“ParsList”,(ParsListTemp (2) ParsListTemp (4:5)]);InitGuessTempDyn = [OptParsTemp (2) InitGuessTemp (4:5)];OptParsTempDyn = fminsearch (@ee_battery_lse InitGuessTempDyn,…optimset (“TolX”1 e - 3));%更新电池块与优化参数OptParsTempDyn = [OptParsTemp (1) OptParsTempDyn (1) OptParsTemp (3) OptParsTempDyn (2:3)];帕尔斯=重塑([ParsListTemp;cellstr (num2str (OptParsTempDyn ') '], 1, []);为k = 1:2:长度(Pars) evalin (“基地”,(Pars {k}“=”帕尔斯{k + 1}“;”])结束assignin (“基地”,“AH0”,啊* battery_data (idx_data) .SOC0);%显示优化参数流([优化电池的温度依赖参数,…“包括电荷动力学:\ n”]);流(' \ t % 5 s = % s \ n 'Pars {:});清晰的i_datav_datat_dataT_dataTs清晰的i_pos一个bc清晰的k清晰的OptChargeOptParsMainOptParsTempOptParsTempDyn清晰的帕尔斯ParsListParsListTempInitGuessChargeInitGuessTempInitGuessTempDyn

为idx_data = 1: num_lines assignin (“基地”,“t_data”battery_data (idx_data) .t);assignin (“基地”,“i_data”battery_data (idx_data)。我);assignin (“基地”,“T_data”battery_data (idx_data) .T *的(长度(t_data), 1));assignin (“基地”,“T0”battery_data (idx_data) .T);assignin (“基地”,“t”t_data (2) -t_data (1));= sim(模型);v_model {idx_data} = out.Vo.signals.values;t_model {idx_data} = out.Vo.time;结束情节([battery_data (1: num_lines)。t]/3600, [battery_data(1:num_lines).v],“o”,(t_model {}): / 3600, [v_model{}):)包含(的时间(小时));ylabel (“电池电压(V)”);传奇([legend_info_data legend_info_model),“位置”,“最佳”);标题(“模型与优化参数值”);

idx_data = 3;assignin (“基地”,“t_data”battery_data (idx_data) .t);assignin (“基地”,“i_data”battery_data (idx_data)。我);assignin (“基地”,“T_data”battery_data (idx_data) .T);assignin (“基地”,“T0”battery_data (idx_data) .T (1));assignin (“基地”,“t”t_data (2) -t_data (1));assignin (“基地”,“AH0”,啊* battery_data (idx_data) .SOC0) = sim(模型);次要情节(1,1)情节(t_data, battery_data .v (3),“o”,out.Vo。t我me, out.Vo.signals.values) xlabel(“时间(s)”);ylabel (“电池电压(V)”);传奇(“数据”,“模型”,“位置”,“最佳”);次要情节(1、2)情节(t_data i_data)包含(“时间(s)”);ylabel (“当前的需求(A)”);次要情节(3、1,3)情节(t_data, t_data)包含(“时间(s)”);ylabel (“温度(^ oC)”);标题(模型验证的);

bdclose(模型)清晰num_lineslegend_info_datalegend_info_model出v_modelt_model清晰的battery_dataidx_dataTs模型