想问一下,我的.mat文件在bin目录里。运行程序时还是显示未定义变量
load xishu
a=zeros(2611,1);
b=zeros(2611,1);
a=xishu(1:2601,1);
b=xishu(2602:end,1);
diffgold1=zeros(2601,1);
gold1=b;
diffgold=diff(gold1,1);
figure(1)
autocorr(diffgold)
% figure(2)
% parcorr(diffgold)
maxLags=4;
AICSet=zeros(maxLags,maxLags);
parfor i=1:maxLags
for j=1:maxLags
mdl=arima('ARLags',[1:i],'MALags',[1:j]);
[EstMdl,EstParamCov,logL,info]=estimate(mdl,diffgold,'display','off');
AICSet(i,j)=aicbic(logL,length(info.X))
end
end
figure
heatmap(AICSet/1000)
Y0=gold1(1:end);
mdl=arima(3,1,3)
N=length(gold1)
[fit,VarCov]=estimate(mdl,gold1(1:end),'Y0',Y0)
[Yf1,YMSE]=forecast(fit,10,'Y0',Y0)
figure
plot(1:N,gold1,'b','LineWidth',1)
hold on
plot(2602:2611,Yf1,'r--','LineWidth',1.5)
x1=1:5216;
x=x1';
a(2602:2611,1)=[150.4079;149.0753;145.5066;141.7708;139.8137;140.7303;144.1829;148.5747;151.8423;152.4808]
b(2602:2611,1)=[0.2854;-0.7522;0.2101;0.4024;-0.5769;-0.1477;0.5064;-0.3183;-0.407;0.4274]
C=idwt(a,b,'db4')
plot(x(1:5200,1),C(1:5200))
hold on
plot(x(5201:end,1),C(5201:end,1),'r')
运行之后就成了这样
有没有大佬帮我看一下
你的xishu来自load xishu,那应该load xishu这个还没生效
您好,我是有问必答小助手,您的问题已经有小伙伴解答了,您看下是否解决,可以追评进行沟通哦~
如果有您比较满意的答案 / 帮您提供解决思路的答案,可以点击【采纳】按钮,给回答的小伙伴一些鼓励哦~~
ps:问答VIP仅需29元,即可享受5次/月 有问必答服务,了解详情>>>https://vip.csdn.net/askvip?utm_source=1146287632
非常感谢您使用有问必答服务,为了后续更快速的帮您解决问题,现诚邀您参与有问必答体验反馈。您的建议将会运用到我们的产品优化中,希望能得到您的支持与协助!
速戳参与调研>>>https://t.csdnimg.cn/Kf0y