当前位置： 首页 > news >正文

第68步 时间序列建模实战：ARIMA建模（Matlab）

news 2026/2/9 22:01:09

基于WIN10的64位系统演示

一、写在前面

这一期，我们使用Matlab进行SARIMA模型的构建。

不同样，这里使用另一个数据：

采用《PLoS One》2015年一篇题目为《Comparison of Two Hybrid Models for Forecasting the Incidence of Hemorrhagic Fever with Renal Syndrome in Jiangsu Province, China》文章的公开数据做演示。数据为江苏省2004年1月至2012年12月肾综合症出血热月发病率。运用2004年1月至2011年12月的数据预测2012年12个月的发病率数据。

Matlab语言，估计很多人没听说过。

所以，直接上我封装好的Matlab小程序即可。

二、Matlab建立SARIMA实战

1.界面介绍

1.1.数据输入 输入用于预测的原始数据，如图一。

1.2.原始数据绘图 点击弹出原始数据曲线图，以便观察数据的消长趋势以及周期性。

1.3.实际值输入 如果用于测试模型的预测精度，此时要预测的数据已知，则输入该数据集。如果用于预测数据，则要预测的数据未知，此时全部输入0。例如要预测6个数据，则需输入6个0。

1.4.平稳检验和相关图模块

1.4.1.原始数据 点击原始数据按钮，在相关图模块显示原始数据的自相关图以及偏相关图，并在显示检验结果（是否平稳）。

1.4.2.一般拆分 填入拆分次数，点击一般拆分按钮，在相关图模块显示经过一般拆分数据的自相关图以及偏相关图，并在显示检验结果（是否平稳）；拆分次数即为d值。

1.4.2.季节性拆分 填入拆分次数，点击季节性拆分按钮，在相关图模块显示经过季节性拆分数据的自相关图以及偏相关图，并在显示检验结果（是否平稳）；拆分次数即为D值。

1.4.3.一般+季节性拆分 点击一般+季节性拆分按钮，在相关图模块显示经过一般和季节性拆分数据的自相关图以及偏相关图，并在显示检验结果（是否平稳）。

1.5.自动参数寻找模块

1.5.1.Log 勾选后，原始数据经过Log变换。此选项影响所有步骤，包括平稳性检验、模型构建、拟合和预测。

1.5.2寻找 必须手动输入参数d和D，周期s默认为12，可以根据实际情况修改。点击寻找按钮后，程序自动寻找最优p、q、P、Q参数（四个参数取值均从0至3）；并将最优参数值输出。

1.6.预测模块

1.6.1.预测数目（个） 显示需要预测数据的个数，程序自动从实际值输入文本框识别个数并输出。

1.6.2.模型参数 输入模型7个参数，其中周期s默认为12，可以根据需要修改。

1.6.3.预测 点击预测按钮后，输出AIC值、BIC值、预测值、实际值，弹出模型参数检验的窗口，并在曲线图模块显示实际值和预测值的曲线图，在误差模块显示预测误差的四个指标。

1.7.曲线图模块

1.7.1.拟合 点击拟合按钮后，显示模型拟合效果曲线图，并在误差模块输出拟合误差的四个指标，弹出显示模型具体拟合值的窗口。

1.7.2.预测 点击拟合按钮后，显示模型拟合效果曲线图，并在误差模块输出拟合误差的四个指标。

1.8.误差模块

1.8.1.拟合误差 MAE 输出模型预测的平均绝对误差；

          MAPE 输出模型预测的平均相对误差；

          MSE 输出模型预测的均方误差；

          RMSE 输出模型预测的均方根误差。

1.8.2.预测误差 MAE 输出模型预测的平均绝对误差；

          MAPE 输出模型预测的平均相对误差；

          MSE 输出模型预测的均方误差；

          RMSE 输出模型预测的均方根误差。

2.实例演示

2.1.打开软件，输入原始数据和实际值，点击原始数据绘图（这里我们用于测试模型精度，故实际值已知），如图二。

2.2.点击原始数据检验按钮，检验结果显示不平稳，如图三；

2.3.可以看到在12logs处数值很高，显示具有季节性，需进行一次季节性拆分；在季节性拆分按钮后的文本框输入1，点击季节性拆分按钮，检验结果显示平稳，如图四；

2.4.此时可在进行一次一般拆分，在一般拆分按钮后的文本框输入1，点击一般+季节性拆分按钮，检验结果显示平稳，如图五；

2.5.输入参数d、D，点击寻找按钮，程序开始寻找最优参数并输出，如图六、七；

2.6.输入参数七个模型参数，点击预测按钮，输出AIC值、BIC值、预测值、实际值，弹出模型参数检验的窗口，并在曲线图模块显示实际值和预测值的曲线图，在误差模块显示预测误差的四个指标，如图八；

2.7.在曲线图模块点击拟合按钮，显示模型拟合效果曲线图，并在误差模块输出拟合误差的四个指标，弹出显示模型具体拟合值的窗口，如图九。

2.8.拟合预测结果比较

模型

预测误差

拟合误差

MAE

MAPE

MAE

MAPE

Matlab ARIMA

0.0074

0.4148

0.0085

0.4246

SPSS ARIMA

0.0090

0.4859

0.0089

0.4219

由以上图表可知，Matlab所构建的ARIMA模型拟合以及预测效果跟SPSS基本一致。

三、小软件安装在此处

1. 双击MCRInstaller.exe安装

2. 自动解压：

3. 无脑点击下一步

4. 安装完成。

四、底层代码

可能有需要，我就贴上吧：

function varargout = ARIMAb(varargin)
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...'gui_Singleton',  gui_Singleton, ...'gui_OpeningFcn', @ARIMA2_OpeningFcn, ...'gui_OutputFcn',  @ARIMA2_OutputFcn, ...'gui_LayoutFcn',  [] , ...'gui_Callback',   []);
if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});
end
if nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
elsegui_mainfcn(gui_State, varargin{:});
end
function ARIMA2_OpeningFcn(hObject, eventdata, handles, varargin)
global objs
handles.output = hObject;
movegui(hObject,'north')
objs.log=handles.logp;
objs.axes = [handles.axes1;handles.axes2;handles.axes3];
objs.disps1=[handles.jianyanjieguo;handles.yibanchaifen;handles.jijiexingchaifen];
objs.disps2=[handles.ppx;handles.d;handles.q;handles.pp;handles.dd;handles.qq];
objs.disps3 =[handles.ppp;handles.d2;handles.q2;handles.pp2;handles.dd2;handles.qq2];
objs.disps4=[handles.aic2;handles.bic2;handles.yucezhi2;handles.shijizhi2];
objs.disps5=[handles.mae;handles.mse;handles.mape;handles.rmse;handles.mae2;handles.mse2;handles.mape2;handles.rmse2];
objs.disps6=[handles.yuceshumu2;handles.s;handles.s2];
guidata(hObject, handles);
function varargout = ARIMA2_OutputFcn(hObject, eventdata, handles) 
varargout{1} = handles.output;
function shujushuru_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
guidata(hObject,handles);
function shujushuru_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function yuanshishujujianyan_Callback(hObject, eventdata, handles)
global X objs
X = str2num(get(handles.shujushuru,'String'));
if isempty(X)return
end
if eq(get(objs.log,'Value'),1)X = log(X);
end
h = adftest (X);
if h == 1U = {'平稳'}set (handles.jianyanjieguo,'String',U);
elseU = {'不平稳'}set (handles.jianyanjieguo,'String',U);
end
axes(handles.axes1)
autocorr(X);
axes(handles.axes2)
parcorr(X);
guidata(hObject,handles);
function chaifenshujujianyan_Callback(hObject, eventdata, handles)
global X T W Z objs
X = str2num(get(handles.shujushuru,'String'));
W = str2num(get(handles.jijiexingchaifen,'String'));
Z = str2num(get(handles.zhouqi,'String'));
if isempty(X) || isempty(T) || isempty(W) || isempty(Z) return
end
if eq(get(objs.log,'Value'),1)X = log(X);
end
o = length (X')/Z;
for R = 1:1:o;
Q = X (:,((R-1)*Z+1):Z*R)
data{R} = Q
end;
for S= 1:1:W;
test =[];
for R = 1:1:(o-S);
sea = data{R+1} - data{R}
test=[test,sea] 
end
end;
h = adftest(test);
if h == 1U={'ƽÎÈ'}set(handles.jianyanjieguo,'String',U)
elseU={'²»Æ½ÎÈ'}set(handles.jianyanjieguo,'String',U)
end
axes(handles.axes1)
autocorr(test);
axes(handles.axes2)
parcorr(test);
guidata(hObject,handles);
function yuceshumu1_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','0')
end
guidata(hObject,handles);
function yuceshumu1_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function ppx_Callback(hObject, eventdata, handles)
function ppx_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction d_Callback(hObject, eventdata, handles)
function d_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function q_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','0')
end
guidata(hObject,handles);
function q_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function pp_Callback(hObject, eventdata, handles)
function pp_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function dd_Callback(hObject, eventdata, handles)
function dd_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function qq_Callback(hObject, eventdata, handles)
function qq_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction s_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','0')
end
guidata(hObject,handles);
function s_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function pushbutton5_Callback(hObject, eventdata, handles)
global X objs
X = str2num(get(handles.shujushuru,'String'));
ss = str2num(get(handles.s,'String'));
value1 = str2num(get(handles.d,'String'));
value2 = str2num(get(handles.dd,'String'));
if isempty(X) || isempty(value1) || isempty(value2) || isempty(ss)return
end
X = X';
if eq(get(objs.log,'Value'),1)X = log(X);
end
test =[];
ind = [];
for p =0:1:3 
for d = value1 
for q = 0:1:3
for dap =0:1:3 
for dad =value2; 
for daq = 0:1:3 ind = [ind;p d q dap dad daq];
end
end
end
end
end
end
T = size(ind,1);
tl = length(X);
v1 = zeros(T,1);
v2 = zeros(T,1);
h = waitbar(0,'1','Name','优化ARIMA参数',...'CreateCancelBtn',...'setappdata(gcbf,''canceling'',1)');
setappdata(h,'canceling',0)
for i = 1:Tif getappdata(h,'canceling')breakendp = ind(i,1);d = ind(i,2);q = ind(i,3);dap=ind(i,4);dad=ind(i,5);daq=ind(i,6);if eq(p,0) && eq(q,0) && eq(dap,0) && eq(daq,0)amodel = arima('D',d,'Seasonality',ss*dad,'Constant',0);elseif ~eq(p,0) && eq(q,0) && eq(dap,0) && eq(daq,0)amodel = arima('ARLags',1:p,'D',d,'Seasonality',ss*dad,'Constant',0);elseif  eq(p,0) && ~eq(q,0) && eq(dap,0) && eq(daq,0)amodel = arima('MALags',1:q,'D',d,'Seasonality',ss*dad,'Constant',0);elseif  eq(p,0) && eq(q,0) && ~eq(dap,0) && eq(daq,0)amodel = arima('SARLags',1:dap,'D',d,'Seasonality',ss*dad,'Constant',0);elseif  eq(p,0) && eq(q,0) && eq(dap,0) && ~eq(daq,0)amodel = arima('SMALags',1:daq,'D',d,'Seasonality',ss*dad,'Constant',0);elseif ~eq(p,0) && ~eq(q,0) && eq(dap,0) && eq(daq,0)amodel = arima('ARLags',1:p,'MALags',1:q,'D',d,'Seasonality',ss*dad,'Constant',0);elseif ~eq(p,0) && eq(q,0) && ~eq(dap,0) && eq(daq,0)amodel = arima('ARLags',1:p,'SARLags',(1:dap)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);elseif ~eq(p,0) && eq(q,0) && eq(dap,0) && ~eq(daq,0)amodel = arima('ARLags',1:p,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);elseif  eq(p,0) && ~eq(q,0) && ~eq(dap,0) && eq(daq,0)amodel = arima('MALags',1:q,'SARLags',(1:dap)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);elseif  eq(p,0) && ~eq(q,0) && eq(dap,0) && ~eq(daq,0)amodel = arima('MALags',1:q,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);elseif  eq(p,0) && eq(q,0) && ~eq(dap,0) && ~eq(daq,0)amodel = arima('SARLags',(1:dap)*ss,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);elseif ~eq(p,0) && ~eq(q,0) && ~eq(dap,0) && eq(daq,0)amodel = arima('ARLags',1:p,'MALags',1:q,'SARLags',(1:dap)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);elseif ~eq(p,0) && ~eq(q,0) && eq(dap,0) && ~eq(daq,0)amodel = arima('ARLags',1:p,'MALags',1:q,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);elseif ~eq(p,0) && eq(q,0) && ~eq(dap,0) && ~eq(daq,0)amodel = arima('ARLags',1:p,'SARLags',(1:dap)*ss,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);elseif  eq(p,0) && ~eq(q,0) && ~eq(dap,0) && ~eq(daq,0)amodel = arima('MALags',1:q,'SARLags',(1:dap)*ss,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);elseamodel = arima('ARLags',1:p,'MALags',1:q,'SARLags',(1:dap)*ss,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);endif numel(X) >= max(amodel.P + amodel.Q + 1);[EstMdl,VarCov,logL,info]=estimate(amodel,X,'Display','off');    f = length(X)-EstMdl.P;errors = sqrt(diag(VarCov));parameters = info.X;tvalues = parameters./errors;pvalues = zeros(size(tvalues));for ii = 1:length(pvalues)pvalues(ii) = (1-tcdf(abs(tvalues(ii)),f))*2;endn = size (VarCov,1);    [aic,bic] = aicbic (logL,n,tl);v1(i) = aic;v2(i) = bic;if length(pvalues)>2pvalues = pvalues(2:end-1);if any(pvalues>0.05)v1(i) = 3*10e8;endendelsev1(i) = 4*10e8;v2(i) = 4*10e8;endwaitbar(i/T,h,['Runing ',num2str(i),' out of ',num2str(T)]);
end
delete(h);
[a,b] = min(v1);
re = ind(b,:);
set(handles.ppx,'String',num2str(re(1)));
set(handles.d,'String',num2str(re(2)));
set(handles.q,'String',num2str(re(3)));
set(handles.pp,'String',num2str(re(4)));
set(handles.dd,'String',num2str(re(5)));
set(handles.qq,'String',num2str(re(6)));
set(handles.s,'String',num2str(12));
guidata(hObject,handles);
function yuanshishujuxiangguantu_Callback(hObject, eventdata, handles)
global X objs
if eq(get(objs.log,'Value'),1)X = log(X);
end
axes(handles.axes1)
autocorr(X);
axes(handles.axes2)
parcorr(X);
function chaifenshujuxiangguantu_Callback(hObject, eventdata, handles)
global X T W Z objs
X = str2num(get(handles.shujushuru,'String'));
T = str2num(get(handles.yibanchaifend,'String'));
W = str2num(get(handles.jijiexingchaifen,'String'));
Z = str2num(get(handles.zhouqi,'String'));
if eq(get(objs.log,'Value'),1)X = log(X);
end
o = length (X')/Z;
for R = 1:1:o;
Q= X (:,((R-1)*Z+1):Z*R)
data{R}=Q
end;
for S= 1:1:W;
test =[];
for R = 1:1:(o-S);
sea = data{R+1} - data{R};
test=[test,sea] ;
end
end;
Y = diff(test,T);
axes(handles.axes1)
autocorr(Y);
axes(handles.axes2)
parcorr(Y);
function yibanchaifend_Callback(hObject, eventdata, handles)
global X T objs
X = str2num(get(handles.shujushuru,'String'));
T = str2num(get(handles.yibanchaifen,'String'));
if isempty(X) || isempty(T) return
end
if eq(get(objs.log,'Value'),1)X = log(X);
end
X =X';
Y = diff (X,T);
h = adftest (Y);
if h == 1U={'ƽÎÈ'};set(handles.jianyanjieguo,'String',U)
elseU={'²»Æ½ÎÈ'};set(handles.jianyanjieguo,'String',U)
end
axes(handles.axes1)
autocorr(Y);
axes(handles.axes2)
parcorr(Y);
guidata(hObject,handles);
function yibanchaifend_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function jijiexingchaifen_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','0')
end
guidata(hObject,handles);
function jijiexingchaifen_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction zhouqi_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','0')
end
guidata(hObject,handles);
function zhouqi_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction jianyanjieguo_Callback(hObject, eventdata, handles)
function jianyanjieguo_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function shijizhi_Callback(hObject, eventdata, handles)
function shijizhi_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function aic_Callback(hObject, eventdata, handles)
function aic_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction yucezhi_Callback(hObject, eventdata, handles)
function yucezhi_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function yuceshumu2_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','0')
end
guidata(hObject,handles);
function yuceshumu2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction canshup_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','0')
end
guidata(hObject,handles);
function canshup_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function d2_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','NaN')
end
guidata(hObject,handles);
function d2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function q2_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','NaN')
end
guidata(hObject,handles);
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction pp2_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','NaN')
end
guidata(hObject,handles);
function pp2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function dd2_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','NaN')
end
guidata(hObject,handles);
function dd2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function qq2_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','NaN')
end
guidata(hObject,handles);
function qq2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function s2_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','NaN')
end
guidata(hObject,handles);
function s2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function yuce2_Callback(hObject, eventdata, handles)
global X K B C A E F G H I L M T objs model
X = str2num(get(handles.shujushuru,'String'));
K = str2num(get(handles.ppp,'String'));
B = str2num(get(handles.d2,'String'));
C = str2num(get(handles.q2,'String'));
A = str2num(get(handles.pp2,'String'));
E = str2num(get(handles.dd2,'String'));
F = str2num(get(handles.qq2,'String'));
G = str2num(get(handles.s2,'String'));
H = str2num(get(handles.yuceshumu2,'String'));
M = str2num(get(handles.shijizhishuru,'String'));
if isempty(X) || isempty(K) || isempty(B) || isempty(C) || isempty(A) || isempty(E) || isempty(F) || isempty(G) || isempty(H) ||isempty(M)return
end
if eq(get(objs.log,'Value'),1)I = log(X');
elseI = X';
end
p = K;d=B;q = C;dap=A;daq=F;dad=E;ss=G;
if eq(p,0) && eq(q,0) && eq(dap,0) && eq(daq,0)amodel = arima('D',d,'Seasonality',ss*dad,'Constant',0);
elseif ~eq(p,0) && eq(q,0) && eq(dap,0) && eq(daq,0)amodel = arima('ARLags',1:p,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif  eq(p,0) && ~eq(q,0) && eq(dap,0) && eq(daq,0)amodel = arima('MALags',1:q,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif  eq(p,0) && eq(q,0) && ~eq(dap,0) && eq(daq,0)amodel = arima('SARLags',1:dap,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif  eq(p,0) && eq(q,0) && eq(dap,0) && ~eq(daq,0)amodel = arima('SMALags',1:daq,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif ~eq(p,0) && ~eq(q,0) && eq(dap,0) && eq(daq,0)amodel = arima('ARLags',1:p,'MALags',1:q,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif ~eq(p,0) && eq(q,0) && ~eq(dap,0) && eq(daq,0)amodel = arima('ARLags',1:p,'SARLags',(1:dap)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif ~eq(p,0) && eq(q,0) && eq(dap,0) && ~eq(daq,0)amodel = arima('ARLags',1:p,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif  eq(p,0) && ~eq(q,0) && ~eq(dap,0) && eq(daq,0)amodel = arima('MALags',1:q,'SARLags',(1:dap)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif  eq(p,0) && ~eq(q,0) && eq(dap,0) && ~eq(daq,0)amodel = arima('MALags',1:q,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif  eq(p,0) && eq(q,0) && ~eq(dap,0) && ~eq(daq,0)amodel = arima('SARLags',(1:dap)*ss,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif ~eq(p,0) && ~eq(q,0) && ~eq(dap,0) && eq(daq,0)amodel = arima('ARLags',1:p,'MALags',1:q,'SARLags',(1:dap)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif ~eq(p,0) && ~eq(q,0) && eq(dap,0) && ~eq(daq,0)amodel = arima('ARLags',1:p,'MALags',1:q,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif ~eq(p,0) && eq(q,0) && ~eq(dap,0) && ~eq(daq,0)amodel = arima('ARLags',1:p,'SARLags',(1:dap)*ss,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);
elseif  eq(p,0) && ~eq(q,0) && ~eq(dap,0) && ~eq(daq,0)amodel = arima('MALags',1:q,'SARLags',(1:dap)*ss,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);
elseamodel = arima('ARLags',1:p,'MALags',1:q,'SARLags',(1:dap)*ss,'SMALags',(1:daq)*ss,'D',d,'Seasonality',ss*dad,'Constant',0);
endmodel3=amodel;
if numel(X) <= max(amodel.P + amodel.Q + 1);msgbox('待估计变量数目大于样本数目，无法计算！');return
end
[fit3,VarCov,logL,info] = estimate(model3,I);
f = length(I)-fit3.P;
errors = sqrt(diag(VarCov));
parameters = info.X;
tvalues = parameters./errors;
pvalues = zeros(size(tvalues));
for i = 1:length(pvalues)pvalues(i) = (1-tcdf(abs(tvalues(i)),f))*2;
end
strs = [];
if ~eq(p,0)for i = 1:pstrs = [strs,{['AR(',num2str(i),')']}];%%%%end
end
if ~eq(dap,0)for i = 1:dapstrs = [strs,{['SAR(',num2str(i*ss),')']}];%%%%end
end
if ~eq(q,0)for i = 1:qstrs = [strs,{['MA(',num2str(i),')']}];%%%%end
end
if ~eq(daq,0)for i=1:daqstrs = [strs,{['SMA(',num2str(i*ss),')']}];end
end
info1 = {'Parameter','Value','Standard Error','t Statistic','p Statistic'};
info2 = [{'Constant'},num2cell(0),{'Fixed'},{'Fixed'},{'Fixed'}];
info3 = [strs,{'Variance'}]';
info4 = [info.X(2:end),errors(2:end),tvalues(2:end),pvalues(2:end)];
info5 = [info1;info2;info3 num2cell(info4)];
f = figure;
setpixelposition(f,[680,558,560,220]);
uitable(f,'Data',info5,'unit','normalized','position',[0 0 1 1],'ColumnWidth',{100});
[yF,yMSE] = forecast(fit3,H,'Y0',I);
if eq(get(objs.log,'Value'),1)L = exp(yF);
end
M = M';
MAE = mean(abs(L - M));
MAPE = mean(abs(L - M)./M);
MSE = mse (L - M);
RMSE = sqrt(mse(L - M));
[hang,lie] = size (VarCov);
T = length(X);
[aic,bic] = aicbic (logL,hang,T);
set (handles.aic2,'String',num2str(aic));
set (handles.bic2,'String',num2str(bic));
set (handles.mae2,'String',num2str(MAE));
set (handles.mape2,'String',num2str(MAPE));
set (handles.mse2,'String',num2str(MSE));
set (handles.rmse2,'String',num2str(RMSE));
set (handles.yucezhi2,'String',num2str(L),'Max',2);
set (handles.shijizhi2,'String',num2str(M),'Max',2);
axes(handles.axes3)
plot(M,'r','LineWidth',2);
hold on
plot(L,'k--','LineWidth',1.5);
xlim([0,H])
title('预测效果拟合曲线')
legend('实际值','预测值','Location','NorthEast')
hold off
model = fit3;
guidata(hObject,handles);
function shijizhi2_Callback(hObject, eventdata, handles)
function shijizhi2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function aic2_Callback(hObject, eventdata, handles)
function aic2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction yucezhi2_Callback(hObject, eventdata, handles)
function yucezhi2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function nihexiaoguo_Callback(hObject, eventdata, handles)
global model X objs
if isempty(X) || isempty(model) return
end
if eq(get(objs.log,'Value'),1)y = log(X');
elsey = X';
end
yv = zeros(size(y));
h = waitbar(0,'ÇëµÈ´ý');
for i = model.P:length(yv)-1vF1 = forecast(model,1,'Y0',y(1:i));yv(i+1) = vF1;waitbar(i/(length(yv)-1),h);
end
delete(h);
yv(1:model.P) =[];
y(1:model.P)=[];
if eq(get(objs.log,'Value'),1)y = exp(y);yv=exp(yv);
end
f = figure;
dispdata = [{'y','yp'};num2cell([y yv])];
uitable(f,'Data',dispdata,'unit','normalized','position',[0 0 1 1],'ColumnWidth',{100});
MAE = mean(abs(y - yv));
MAPE = mean(abs(y - yv)./y);
MSE = mse(y - yv);
RMSE = sqrt(mse(y - yv));
set(objs.disps5(1),'String',MAE);
set(objs.disps5(2),'String',MSE);
set(objs.disps5(3),'String',MAPE);
set(objs.disps5(4),'String',RMSE);
axes(handles.axes3)
plot(y,'r','LineWidth',2);
hold on
plot(yv,'k--','LineWidth',1.5);
xlim([0,length(y)])
title('拟合曲线')
legend('实际值','拟合值','Location','NorthEast')
hold off
guidata(hObject,handles);
function yucexiaoguo_Callback(hObject, eventdata, handles)
global X K B C A E F G H I L M T objs model
X = str2num(get(handles.shujushuru,'String'));
K = str2num(get(handles.ppp,'String'));
B = str2num(get(handles.d2,'String'));
C = str2num(get(handles.q2,'String'));
A = str2num(get(handles.pp2,'String'));
E = str2num(get(handles.dd2,'String'));
F = str2num(get(handles.qq2,'String'));
G = str2num(get(handles.s2,'String'));
H = str2num(get(handles.yuceshumu2,'String'));
M = str2num(get(handles.shijizhishuru,'String'));
if isempty(X) || isempty(K) || isempty(B) || isempty(C) || isempty(A) || isempty(E) || isempty(F) || isempty(G) || isempty(H) ||isempty(M) ||isempty(model)return
end
if eq(get(objs.log,'Value'),1)I = log(X');
elseI = X';
end
fit3=model;
[yF,yMSE] = forecast(fit3,H,'Y0',I);
if eq(get(objs.log,'Value'),1)L = exp(yF);
endaxes(handles.axes3)
plot(M,'r','LineWidth',2);
hold on
plot(L,'k--','LineWidth',1.5);
xlim([0,H])
title('预测效果拟合曲线')
legend('实际值','预测值','Location','NorthEast')
hold off
model = fit3;
guidata(hObject,handles);
function mae_Callback(hObject, eventdata, handles)
function mae_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function mape_Callback(hObject, eventdata, handles)
function mape_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction mse_Callback(hObject, eventdata, handles)
function mse_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction rmse_Callback(hObject, eventdata, handles)
function rmse_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function mae2_Callback(hObject, eventdata, handles)
function mae2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function mape2_Callback(hObject, eventdata, handles)
function mape2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction mse2_Callback(hObject, eventdata, handles)
function mse2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction rmse2_Callback(hObject, eventdata, handles)
function rmse2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction ppp_Callback(hObject, eventdata, handles)
input = str2num(get(hObject,'String'));
if (isempty(input))set(hObject,'String','NaN')
end
guidata(hObject,handles);
function ppp_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function shijizhishuru_Callback(hObject, eventdata, handles)
global objs
num = str2num(get(hObject,'String'));
if length(num)>1set(objs.disps6(1),'String',num2str(length(num)));
end
guidata(hObject,handles);
function shijizhishuru_CreateFcn(hObject, eventdata, handles)if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function yuanshishujuhuitu_Callback(hObject, eventdata, handles)
global X T
X = str2num(get(handles.shujushuru,'String'));
if isempty(X)return
end
T = length (X);
figure
hold on
plot(X,'k-','LineWidth',1.5)
xlim([0,T])
title('原始数据曲线图')
hold off
guidata(hObject,handles);
function yibanjijiexingchaifen_Callback(hObject, eventdata, handles)
global X T W Z objs
X = str2num(get(handles.shujushuru,'String'));
T = str2num(get(handles.yibanchaifen,'String'));
W = str2num(get(handles.jijiexingchaifen,'String'));
Z = str2num(get(handles.zhouqi,'String'));
if isempty(X) || isempty(T) || isempty(W) || isempty(Z) return
end
if eq(get(objs.log,'Value'),1)X = log(X);
end
o = length (X')/Z;
for R = 1:1:o;
Q = X (:,((R-1)*Z+1):Z*R)
data{R} = Q
end;
for S= 1:1:W;
test =[];
for R = 1:1:(o-S);
sea = data{R+1} - data{R}
test=[test,sea] 
end
end;
Y = diff(test,T);
h = adftest(Y);
if h == 1U={'平稳'}set(handles.jianyanjieguo,'String',U)
elseU={'不平稳'}set(handles.jianyanjieguo,'String',U)
end
axes(handles.axes1)
autocorr(Y);
axes(handles.axes2)
parcorr(Y);
guidata(hObject,handles);
function bic2_Callback(hObject, eventdata, handles)
function bic2_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function bic_Callback(hObject, eventdata, handles)
function bic_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end
function pushbutton12_Callback(hObject, eventdata, handles)
global X K B C A E F G H I L M T
X = str2num(get(handles.shujushuru,'String'));
K = str2num(get(handles.ppx,'String'));
B = str2num(get(handles.d,'String'));
C = str2num(get(handles.q,'String'));
A = str2num(get(handles.pp,'String'));
E = str2num(get(handles.dd,'String'));
F = str2num(get(handles.qq,'String'));
G = str2num(get(handles.s,'String'));
H = str2num(get(handles.yuceshumu,'String'));
M = str2num(get(handles.shijizhishuru,'String'));
I = log(X');
model3 = arima('ARLags',K,'MALags',C,'D',B,'SARLags',A*G,'SMALags',F*G,'Seasonality',G,'Constant',0);
[fit3,VarCov,logL,info] = estimate(model3,I);
[yF,yMSE] = forecast(fit3,H,'Y0',I);
L = exp (yF);
M = M';
MAE = mean(abs(L - M));
MAPE = mean(abs(L - M)./M);
MSE = mse (L - M);
RMSE = sqrt(mse(L - M));
[hang,lie] = size (VarCov);
T = length(X);
[aic,bic] = aicbic (logL,hang,T);
set (handles.aic,'String',num2str(aic));
set (handles.bic,'String',num2str(bic));
set (handles.mae2,'String',num2str(MAE));
set (handles.mape2,'String',num2str(MAPE));
set (handles.mse2,'String',num2str(MSE));
set (handles.rmse2,'String',num2str(RMSE));
set (handles.yucezhi,'String',num2str(L),'Max',2);
set (handles.shijizhi,'String',num2str(M),'Max',2);
axes(handles.axes3)
plot(M,'r','LineWidth',2);
hold on
plot(L,'k--','LineWidth',1.5);
xlim([0,H])
title('Ô¤²âЧ¹ûÄâºÏÇúÏß')
legend('ʵ¼ÊÖµ','Ô¤²âÖµ','Location','NorthEast')
hold off
guidata(hObject,handles);
function yibanchaifen_Callback(hObject, eventdata, handles)
function yibanchaifen_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction logp_Callback(hObject, eventdata, handles)
global objs
for i = 1:length(objs.axes)cla(objs.axes(i))
end
set(objs.disps1,'String','');
set(objs.disps2,'String','');
set(objs.disps3,'String','');
set(objs.disps4,'String','');
set(objs.disps5,'String','');
function tStatistic= print1(Mdl, covariance)
if size(covariance,1) ~= size(covariance,2)error(message('econ:arima:print:NonSquareCovarianceMatrix'))
end
constant = Mdl.Constant;            
beta     = Mdl.Beta;                
beta     = beta(:)';                
AR  = reflect(getLagOp(Mdl, 'AR'));
SAR = reflect(getLagOp(Mdl, 'SAR'));
MA  = getLagOp(Mdl, 'MA');
SMA = getLagOp(Mdl, 'SMA');LagsAR  = AR.Lags;                  
LagsSAR = SAR.Lags;                 
LagsMA  = MA.Lags;                  
LagsSMA = SMA.Lags;                 
LagsAR  = LagsAR(LagsAR > 0);       
LagsSAR = LagsSAR(LagsSAR > 0);
LagsMA  = LagsMA(LagsMA > 0);
LagsSMA = LagsSMA(LagsSMA > 0);if isempty(LagsAR)AR = [];
elseAR = AR.Coefficients;AR = [AR{LagsAR}];               
endif isempty(LagsSAR)SAR = []; 
elseSAR = SAR.Coefficients;SAR = [SAR{LagsSAR}];          
endif isempty(LagsMA)MA = [];
elseMA = MA.Coefficients;MA = [MA{LagsMA}];               
end
if isempty(LagsSMA)SMA = []; 
elseSMA = SMA.Coefficients;SMA = [SMA{LagsSMA}];            
endisDistributionT    =  strcmpi(Mdl.Distribution.Name, 'T');
isVarianceConstant = ~isa(Mdl.Variance, 'internal.econ.LagIndexableTimeSeries');  nARIMA = 1 + numel(LagsAR) + numel(LagsSAR) + ...numel(LagsMA) + numel(LagsSMA) + numel(beta);
if isVarianceConstantif isDistributionTparameters          = zeros(nARIMA + 2,1);parameters(end - 1) = Mdl.Variance;parameters(end)     = Mdl.Distribution.DoF;elseparameters          = zeros(nARIMA + 1,1);parameters(end)     = Mdl.Variance;endparameters(1:nARIMA) = [constant AR  SAR  MA  SMA  beta];if numel(parameters) ~= numel(solve)error(message('econ:arima:print:ModelCovarianceInconsistency'))endelse  if isDistributionTparameters(nARIMA + 1) = Mdl.Distribution.DoF;parameters(1:nARIMA)   = [constant AR  SAR  MA  SMA  beta];elseparameters = [constant AR  SAR  MA  SMA  beta];end
endFix    = ~solve;
errors = sqrt(diag(covariance));
tStatistic = parameters ./ errors;end
function pushbutton13_Callback(hObject, eventdata, handles)
close ARIMAb
main
function pushbutton14_Callback(hObject, eventdata, handles)
set(handles.shujushuru,'String',[]);
set(handles.shijizhishuru,'String',[]);
set(handles.jianyanjieguo,'String',[]);
set(handles.yibanchaifen,'String',[]);
set(handles.jijiexingchaifen,'String',[]);
set(handles.zhouqi,'String',[]);
set(handles.ppx,'String',[]);
set(handles.d,'String',[]);
set(handles.p,'String',[]);
set(handles.pp,'String',[]);
set(handles.dd,'String',[]);
set(handles.qq,'String',[]);
set(handles.s,'String',[]);
set(handles.yuceshumu2,'String',[]);
set(handles.ppp,'String',[]);
set(handles.d2,'String',[]);
set(handles.p2,'String',[]);
set(handles.pp2,'String',[]);
set(handles.dd2,'String',[]);
set(handles.qq2,'String',[]);
set(handles.s2,'String',[]);
cla(handles.axes2);
cla(handles.axes3);
set(handles.aic2,'String',[]);
set(handles.bic2,'String',[]);
set(handles.yucezhi2,'String',[]);
set(handles.shijizhi2,'String',[]);
set(handles.mae,'String',[]);
set(handles.mape,'String',[]);
set(handles.mse,'String',[]);
set(handles.rmse,'String',[]);
cla(handles.axes1);
set(handles.mae2,'String',[]);
set(handles.mape2,'String',[]);
set(handles.mse2,'String',[]);
set(handles.rmse2,'String',[]);

五、软件和数据

链接：https://pan.baidu.com/s/1ZZtT6VTUnqGtQaAZ9NlyCg?pwd=hk4p

提取码：hk4p

相关文章：

第68步 时间序列建模实战：ARIMA建模（Matlab）

基于WIN10的64位系统演示 一、写在前面 这一期&#xff0c;我们使用Matlab进行SARIMA模型的构建。 不同样&#xff0c;这里使用另一个数据&#xff1a; 采用《PLoS One》2015年一篇题目为《Comparison of Two Hybrid Models for Forecasting the Incidence of Hemorrhagic …...

编程日记 2023/9/11 22:52:24

Gin学习记录3——模版与渲染

模版与渲染 一. 返回二. 模版2.1 基础模版2.2 同名模版2.3 模版继承2.4 模版语法 一. 返回 如果只是想返回数据&#xff0c;可以使用以下函数&#xff1a; func (c *Context) JSON(code int, obj any) func (c *Context) JSONP(code int, obj any) func (c *Context) String(…...

编程日记 2023/9/11 22:51:23

Python算法练习 9.11

leetcode 392 判断子序列 给定字符串 s 和 t &#xff0c;判断 s 是否为 t 的子序列。 字符串的一个子序列是原始字符串删除一些&#xff08;也可以不删除&#xff09;字符而不改变剩余字符相对位置形成的新字符串。&#xff08;例如&#xff0c;"ace"是"abcd…...

编程日记 2023/9/11 22:50:22

2023年中秋节和国庆节放假几天？用待办软件记录放假安排并提醒

进入公历9月&#xff0c;我们都期待着下个长假的到来。那么2023年中秋节和国庆节放假几天呢&#xff1f;因为今年的中秋节是公历的9月29日&#xff0c;所以今年的中秋节和国庆节是连在一起放假的。放假时间安排是9月29日至10月6日&#xff0c;一共放假8天。而10月7日和8日则是调…...

编程日记 2023/9/11 22:49:20

使用Python实现一个完整的声音采样和模拟，使用采样声音播放输入的文字，实现代码进行详细注释，并进行测试

目录 1.功能概述 2.原理介绍 2.1.声音采样原理 2.2.PCM系统原理 2.3.声音学习与训练...

编程日记 2023/9/11 22:48:19

测试----计算机网络

文章目录 计算机网络的历史OSI/RM 协议TCP/IP协议IP地址 计算机网络的历史 50-60年代 内部通讯功能&#xff08;连接的是同一台主机&#xff0c;只能主机和终端之间通信&#xff0c;终端和终端之间的通讯只能依靠主机来传输&#xff09;60-70年代 主机和主机之间能通讯70年代-…...

编程日记 2023/9/11 22:47:18

SVN 索引版本与打包版本号不匹配

今天突然遇到了一个问题&#xff0c;SVN上传不了&#xff0c;错误提示如下&#xff1a; 解决方法&#xff1a; 1.其实&#xff0c;这是SVN库不小心搞坏了&#xff0c;只能重新再创建一个SVN仓库了。...

编程日记 2023/9/11 22:46:17

HummerRisk V1.4.1 发布

HummerRisk V1.4.1发布&#xff1a; 大家好&#xff0c;增加检测整合报告下载&#xff0c;定制多云整合报告并下载PDF&#xff0c;增加K8s 检测规则组&#xff0c;Kubernetes、Rancher、KubeSphere 检测规则组以及规则。新增云账号管理页面关联菜单&#xff0c;新增资源同步日…...

编程日记 2023/9/11 22:45:16

php的html实体和字符之间的转换

html_entity_decode() 函数是 htmlentities() 函数的反函数。用于把HTML实体转换为字符。 html_entity_decode() 函数把 HTML 实体转换为字符。 $str "<© W3CS&ccedil;h&deg;&deg;&brvbar;&sect;>"; echo html_entity_decode($str…...

编程日记 2023/9/11 22:44:14

docker-compose deploy 高可用 elasticsearch TLS

文章目录 1.sysctl2. swap3. hosts4. 配置 instances.yaml5. 创建证书6. 部署7. 修改 kibanna 密码8. 清理 1.sysctl [rootgithub es_tls]# cat /etc/sysctl.conf # sysctl settings are defined through files in # /usr/lib/sysctl.d/, /run/sysctl.d/, and /etc/sysctl.d/…...

编程日记 2023/9/11 22:43:13

让GPT成为您的科研加速器丨GPT引领前沿与应用突破之GPT4科研实践技术与AI绘图

GPT对于每个科研人员已经成为不可或缺的辅助工具&#xff0c;不同的研究领域和项目具有不同的需求。如在科研编程、绘图领域&#xff1a;1、编程建议和示例代码:无论你使用的编程语言是Python、R、MATLAB还是其他语言&#xff0c;都可以为你提供相关的代码示例。​2、数据可视化…...

编程日记 2023/9/11 22:42:12

工业互联网的破局密钥——低代码开发

纵观历史上三次工业革命的演进过程&#xff0c;以蒸汽机、电力、计算机为相应时代的“新型”基础设施建设&#xff0c;依托“通用技术”的广泛应用&#xff0c;带动了整个行业与市场的发展。 在发展数字经济成为必选题的今天&#xff0c;同样需要一个新型的、基础性的、通用型…...

编程日记 2023/9/11 22:41:11

JavaScript知识系列（2）每天10个小知识点

目录 系列文章目录JavaScript知识系列&#xff08;1&#xff09;每天10个小知识点 知识点**11. 如果 new 一个箭头函数的会怎么样****12. 箭头函数的 this 指向哪⾥&#xff1f;****13. 扩展运算符**的概念、作用、原理、特性、优点、缺点、区别、使用场景**14. Proxy**的概念、…...

编程日记 2023/9/11 22:40:10

Kotlin面向对象基础使用方法（继承、接口、Lambda、空指针检查机制等）

三、面向对象 1、继承 1.1 open改变类的继承属性 在kotlin设计时默认所有的非抽象类是无法被继承的&#xff0c;如果想要使得一个非抽象类可以被继承&#xff0c;我们需要使用open关键字。 open class Person {var name "";var age 0;fun eat() {println(name …...

编程日记 2023/9/11 22:39:09

Android USB电源管理

The USB peripheral detects the lack of 3 consecutive SOF packets as a suspend request from the USB host. 1 驱动shutdown顺序 系统关机或重启的过程中&#xff0c;会调用设备驱动的shutdown函数来完成设备的关闭操作&#xff0c;有需要的设备可以在驱动中定义该函数。其…...

编程日记 2023/9/11 22:38:08

YOLO目标检测——路标数据集+已标注voc和yolo格式标签下载分享

实际项目应用&#xff1a;自动驾驶、视频监控和安防、物体识别和分类、城市规划和地理信息系统等等数据集说明&#xff1a;YOLO路标目标检测数据集&#xff0c;真实场景的高质量图片数据&#xff0c;数据场景丰富&#xff0c;图片格式为jpg&#xff0c;共900张图片&#xff0c;…...

编程日记 2023/9/11 22:37:07

Item-Based Recommendations with Hadoop

Mahout在MapReduce上实现了Item-Based Collaborative Filtering&#xff0c;这里我尝试运行一下。 安装Hadoop 从下载Mahout并解压 准备数据 下载1 Million MovieLens Dataset&#xff0c;解压得到ratings.dat&#xff0c;用 sed ‘s/:&#x1f626;[0-9]{1,}):&#x1f626;…...

编程日记 2023/9/11 22:36:05

基于物理层网络编码的相位同步算法matlab仿真

目录 1.算法运行效果图预览 2.算法运行软件版本 3.部分核心程序 4.算法理论概述 5.算法完整程序工程 1.算法运行效果图预览 2.算法运行软件版本 matlab2022a 3.部分核心程序 ..........................................................................%数据长度 Len…...

编程日记 2023/9/11 22:35:03

数据结构——七大排序[源码+动图+性能测试]

本章代码gitee仓库&#xff1a;排序 文章目录 &#x1f383;0. 思维导图&#x1f9e8;1. 插入排序✨1.1 直接插入排序✨1.2 希尔排序 &#x1f38a;2. 选择排序&#x1f38b;2.1 直接选择排序&#x1f38b;2.2 堆排序 &#x1f38f;3. 交换排序&#x1f390;3.1 冒泡排序&#…...

编程日记 2023/9/11 22:34:02

G. The Morning Star

Problem - G - Codeforces 思路&#xff1a;想了挺长时间的&#xff0c;一直没想到一个简便的方法在瞎搞。我们发现对于某个点来说&#xff0c;其他的点如果能够跟他匹配&#xff0c;那么一定在这8个方向上&#xff0c;而同时这8个方向其实对应这4条直线&#xff0c;假设点为(x…...

编程日记 2023/9/11 22:32:57

深度学习在微纳光子学中的应用

深度学习在微纳光子学中的主要应用方向 深度学习与微纳光子学的结合主要集中在以下几个方向&#xff1a; 逆向设计 通过神经网络快速预测微纳结构的光学响应&#xff0c;替代传统耗时的数值模拟方法。例如设计超表面、光子晶体等结构。 特征提取与优化 从复杂的光学数据中自…...

编程新知 2026/2/8 15:03:04

【Python】 -- 趣味代码 - 小恐龙游戏

文章目录 文章目录 00 小恐龙游戏程序设计框架代码结构和功能游戏流程总结01 小恐龙游戏程序设计02 百度网盘地址00 小恐龙游戏程序设计框架 这段代码是一个基于 Pygame 的简易跑酷游戏的完整实现，玩家控制一个角色（龙）躲避障碍物（仙人掌和乌鸦）。以下是代码的详细介绍：…...

编程新知 2026/2/8 20:43:07

【网络安全产品大调研系列】2. 体验漏洞扫描

前言 2023 年漏洞扫描服务市场规模预计为 3.06&#xff08;十亿美元&#xff09;。漏洞扫描服务市场行业预计将从 2024 年的 3.48&#xff08;十亿美元&#xff09;增长到 2032 年的 9.54&#xff08;十亿美元&#xff09;。预测期内漏洞扫描服务市场 CAGR&#xff08;增长率&…...

编程新知 2026/2/4 12:43:08

跨链模式：多链互操作架构与性能扩展方案

跨链模式&#xff1a;多链互操作架构与性能扩展方案 ——构建下一代区块链互联网的技术基石 一、跨链架构的核心范式演进 1. 分层协议栈&#xff1a;模块化解耦设计 现代跨链系统采用分层协议栈实现灵活扩展&#xff08;H2Cross架构&#xff09;&#xff1a; 适配层&#xf…...

编程新知 2026/2/5 14:15:13

sqlserver 根据指定字符 解析拼接字符串

DECLARE LotNo NVARCHAR(50)A,B,C DECLARE xml XML ( SELECT <x> REPLACE(LotNo, ,, </x><x>) </x> ) DECLARE ErrorCode NVARCHAR(50) -- 提取 XML 中的值 SELECT value x.value(., VARCHAR(MAX))…...

编程新知 2025/10/29 4:33:03

新能源汽车智慧充电桩管理方案：新能源充电桩散热问题及消防安全监管方案

随着新能源汽车的快速普及&#xff0c;充电桩作为核心配套设施&#xff0c;其安全性与可靠性备受关注。然而&#xff0c;在高温、高负荷运行环境下&#xff0c;充电桩的散热问题与消防安全隐患日益凸显&#xff0c;成为制约行业发展的关键瓶颈。 如何通过智慧化管理手段优化散…...

编程新知 2026/1/29 3:00:56

QT： `long long` 类型转换为 `QString` 2025.6.5

在 Qt 中&#xff0c;将 long long 类型转换为 QString 可以通过以下两种常用方法实现&#xff1a; 方法 1&#xff1a;使用 QString::number() 直接调用 QString 的静态方法 number()&#xff0c;将数值转换为字符串&#xff1a; long long value 1234567890123456789LL; …...

编程新知 2026/2/1 6:30:04

HDFS分布式存储 zookeeper

hadoop介绍 狭义上hadoop是指apache的一款开源软件 用java语言实现开源框架&#xff0c;允许使用简单的变成模型跨计算机对大型集群进行分布式处理&#xff08;1.海量的数据存储 2.海量数据的计算&#xff09;Hadoop核心组件 hdfs&#xff08;分布式文件存储系统&#xff09;&a…...

编程新知 2026/1/24 4:10:44

HarmonyOS运动开发：如何用mpchart绘制运动配速图表

##鸿蒙核心技术##运动开发##Sensor Service Kit&#xff08;传感器服务&#xff09;# 前言 在运动类应用中&#xff0c;运动数据的可视化是提升用户体验的重要环节。通过直观的图表展示运动过程中的关键数据&#xff0c;如配速、距离、卡路里消耗等&#xff0c;用户可以更清晰…...

编程新知 2026/2/2 7:49:59

动态 Web 开发技术入门篇

一、HTTP 协议核心 1.1 HTTP 基础 协议全称 &#xff1a;HyperText Transfer Protocol&#xff08;超文本传输协议&#xff09; 默认端口 &#xff1a;HTTP 使用 80 端口&#xff0c;HTTPS 使用 443 端口。 请求方法 &#xff1a; GET &#xff1a;用于获取资源&#xff0c;…...

编程新知 2025/10/3 0:42:17