preisach算法模型

preisach算法模型，用于simulink建模仿真分析，可以用于电池模型了解

function varargout = preis_sum_anev（varargin）
% Preis_clas_AnEv Application M-file for Preis_clas_AnEv.fig
%    FIG = Preis_clas_AnEv launch preis_sum_anev GUI.
%    Preis_clas_AnEv（‘callback_name‘ ...） invoke the named callback.

% Last Modified by GUIDE v2.0 19-May-2005 16:21:33

if nargin == 0  % LAUNCH GUI
    global H M H1 M1 H2 jel N_s pathname filename Lj nLa H_comp M_comp cik n_m Ev jel_Ev
    H = []; M = []; H1 = []; M1 = []; H2 = []; jel = []; N_s = [];  
    pathname = []; filename = []; H_comp = []; M_comp = []; Ev = [];
    
    fig = openfig（mfilename‘reuse‘）;

	% Use system color scheme for figure:
	set（fig‘Color‘get（0‘defaultUicontrolBackgroundColor‘））;

	% Generate a structure of handles to pass to callbacks and store it. 
	handles = guihandles（fig）;
	guidata（fig handles）;

	if nargout > 0
		varargout{1} = fig;
	end
    
    set（handles.radiobutton_table_h1h2‘Value‘1）;
    set（handles.radiobutton_table_rev‘Value‘0）;
    
    %the parameters of mu（h1h2）=sum（fi（-h1）*fi（h2））
    set（handles.edit_n‘String‘num2str（1））;
    x_h1h2（11） = 0.7;    x_h1h2（12） = 0.3;    x_h1h2（13） = 0.4;
    set（handles.edit_ai‘String‘num2str（x_h1h2（11）））;
    set（handles.edit_bi‘String‘num2str（x_h1h2（12）））;
    set（handles.edit_ci‘String‘num2str（x_h1h2（13）））;
    set（handles.listbox_param_h1h2‘String‘num2str（x_h1h25））;
    %the parameters of the reversible component
    set（handles.edit_d‘String‘num2str（0.001））;
    x_rev（11） = 0.1;     x_rev（12） = 1;
    set（handles.listbox_param_rev‘String‘num2str（x_rev5））;

    Str = ‘f（x） = ai*exp（ -（x-bi）/ci ）/（ 1 + exp（ -（x-bi）/ci ） ）^2‘;
    set（handles.text_distrib‘String‘Str）;
    
    Hsat = 1;
    n_H = 100;
    set（handles.edit_n_H‘String‘num2str（n_H））;
    set（handles.edit_Hsat‘String‘num2str（Hsat））;

    set（handles.edit_MFE‘String‘num2str（500））;
    
    set（handles.pushbutton_optim‘Enable‘‘off‘）;
  
    set（handles.slider_H‘Min‘-Hsat）;
    set（handles.slider_H‘Max‘Hsat）;
    set（handles.edit_H‘String‘num2str（0））;
    
    n_m=10;
    Lj=zeros（n_m2）; %a lepcsosgorbe inicializalasa
    nLa=1;
    cik = 1;
    jel_Ev = 0;
elseif ischar（varargin{1}） % INVOKE NAMED SUBFUNCTION OR CALLBACK
 
	try
		[varargout{1:nargout}] = feval（varargin{:}）; % FEVAL switchyard
	catch
		disp（lasterr）;
	end

end


%| ABOUT CALLBACKS:
%| GUIDE automatically appends subfunction prototypes to this file and 
%| sets objects‘ callback properties to call them through the FEVAL 
%| switchyard above. This comment describes that mechanism.
%|
%| Each callback subfunction declaration has the following form:
%| （H EVENTDATA HANDLES VARARGIN）
%|
%| The subfunction name is composed using the object‘s Tag and the 
%| callback type separated by ‘_‘ e.g. ‘slider2_Callback‘
%| ‘figure1_CloseRequestFcn‘ ‘axis1_ButtondownFcn‘.
%|
%| H is the callback object‘s handle （obtained using GCBO）.
%|
%| EVENTDATA is empty but reserved for future use.
%|
%| HANDLES is a structure c

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-12-08 07:00  Preisach model\
     文件       48883  2015-12-04 16:23  Preisach model\Algorithm_Preis_class_AnEv.pdf
     文件       38478  2015-12-04 16:23  Preisach model\Algorithm_Preis_class_point.pdf
     文件      206170  2015-12-04 16:23  Preisach model\HB_2.bmp
     文件      762514  2015-12-04 16:22  Preisach model\Poster_SzaboFuzi_compumag_2015.pdf
     目录           0  2015-12-04 19:44  Preisach model\Preis_class_AnEv\
     目录           0  2015-12-04 19:44  Preisach model\Preis_class_AnEv\Preis_class_AnEv\
     文件       19800  2015-12-04 19:44  Preisach model\Preis_class_AnEv\Preis_class_AnEv\HB_0_11.dat
     文件         210  2015-12-04 19:44  Preisach model\Preis_class_AnEv\Preis_class_AnEv\par_HB_0_11_3_1_class.dat
     文件       53848  2015-12-04 19:44  Preisach model\Preis_class_AnEv\Preis_class_AnEv\Preis_sum_AnEv.fig
     文件       30532  2015-12-04 19:44  Preisach model\Preis_class_AnEv\Preis_class_AnEv\Preis_sum_AnEv.m
     目录           0  2015-12-04 19:44  Preisach model\Preis_class_AnEv_a\
     目录           0  2015-12-04 19:44  Preisach model\Preis_class_AnEv_a\Preis_class_AnEv_a\
     文件        2022  2015-12-04 19:44  Preisach model\Preis_class_AnEv_a\Preis_class_AnEv_a\kimenet_ClForm.m
     文件        4947  2015-12-04 19:44  Preisach model\Preis_class_AnEv_a\Preis_class_AnEv_a\main_An.m
     文件       42568  2015-12-04 19:44  Preisach model\Preis_class_AnEv_a\Preis_class_AnEv_a\mh.txt
     文件         701  2015-12-04 19:44  Preisach model\Preis_class_AnEv_a\Preis_class_AnEv_a\param_mh.dat
     文件        2343  2015-12-04 19:44  Preisach model\Preis_class_AnEv_a\Preis_class_AnEv_a\vectf_all_1.m
     目录           0  2015-12-04 19:43  Preisach model\Preis_class_point\
     目录           0  2015-12-04 19:43  Preisach model\Preis_class_point\Preis_class_point\
     文件       43672  2015-12-04 19:43  Preisach model\Preis_class_point\Preis_class_point\graf_main.fig
     文件       36542  2015-12-04 19:43  Preisach model\Preis_class_point\Preis_class_point\graf_main.m
     文件        6834  2015-12-04 19:43  Preisach model\Preis_class_point\Preis_class_point\HB1.dat
     目录           0  2015-12-04 19:45  Preisach model\Preis_dyn_AnEv_ORSI\
     目录           0  2015-12-04 19:45  Preisach model\Preis_dyn_AnEv_ORSI\Preis_dyn_AnEv_ORSI\
     文件       25800  2015-12-04 19:45  Preisach model\Preis_dyn_AnEv_ORSI\Preis_dyn_AnEv_ORSI\HB_frec_20TO500.txt
     文件        1759  2015-12-04 19:45  Preisach model\Preis_dyn_AnEv_ORSI\Preis_dyn_AnEv_ORSI\kimenet_ClForm_abc.m
     文件        1115  2015-12-04 19:45  Preisach model\Preis_dyn_AnEv_ORSI\Preis_dyn_AnEv_ORSI\kimenet_dMpdH_AnEv_abc.m
     文件        3984  2015-12-06 16:53  Preisach model\Preis_dyn_AnEv_ORSI\Preis_dyn_AnEv_ORSI\main_Preis_Dyn_HM.m
     文件        4616  2015-12-04 19:45  Preisach model\Preis_dyn_AnEv_ORSI\Preis_dyn_AnEv_ORSI\main_Preis_Dyn_HM_1.m
     文件          50  2015-12-04 19:45  Preisach model\Preis_dyn_AnEv_ORSI\Preis_dyn_AnEv_ORSI\Param_HB_Dyn_4.txt
............此处省略3个文件信息