2022年ARMA模型程序[归 .pdf
《2022年ARMA模型程序[归 .pdf》由会员分享,可在线阅读,更多相关《2022年ARMA模型程序[归 .pdf(6页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、clear,clc close all set(0,defaultaxeslinestyleorder,-.p,:+,-s); %设置默认的线性和标识符set(0,defaultaxescolororder,0,0,1)% 设置默认的线条颜色x=-pi:pi/10:pi; y=tan(sin(x)-sin(tan(x); set(gca,xtick,min(x),max(x) %设置 x 轴的范围set(gca,ylim, min(y) ,max(y),layer,top)%设置 y 轴的范围h=plot(x,y,markersize,9); % 画图并且取得二维图形中最大值和最小值的索引号g
2、rid on set(gca,layer,top) %将栅格放在最上层,防止遮盖axis equal %使每一格的刻度相同,同时也使长度一样xlabel(-pileq theta leqpi,fontsize,20) ylabel(tan(sin(theta)-sin(tan(theta),fontsize,20) title(it-pileq thetaleqpi的三角函数图 ,fontsize,20) text(0,0,bulletit 原点 ,fontsize,25) %文本的精确定位%标注曲线的最大值和最小值x=get(h,xdata); % 获得二维曲线的数据y=get(h,ydat
3、a); imin=find(min(y)=y); imax=find(max(y)=y); text(x(imin),y(imin),最小值 =,num2str(y(imin),. horizontalalignment,center,verticalalignment,bottom,fontsize,17) text(x(imax),y(imax),最大值 =,num2str(y(imax),. horizontalalignment,center,verticalalignment,top,fontsize,17) 宋哲 2011/8/16 11:07:42 %AR,MA模型clear,c
4、lc close all a=load(22.txt); a=a(2,4,:);a=a; a=a(:);%将数据税收展成按时间序列排列的% plot(a,-p) b=diff(a); % plot(b,-p) r1=autocorr(b)% 计算自相关系数r2=parcorr(b) %计算偏自相关系数figure, subplot(121),autocorr(b) subplot(122), parcorr(b) cs1=ar(b,2) %AR 模型, b 必须是列向量bhat=predict(cs1,b;20) %Discrete-time IDPOLY model: %A(q) = 1 -
5、 0.5684 q-1 - 0.2493 q-2 b15=bhat1(end) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页,共 6 页 - - - - - - - - - bhat2=predict(cs1,b;b15;-1) 宋哲 2011/8/16 11:27:00 %用逻辑斯特预测clear,clc a=load(22.txt); a=a(2,4,:);a=a; a=a(:);%将数据税收展成按时间序列排列的x0=a(1);tt=1:14; t0=1; xt=(cs
6、,t)cs(1)/(1+(cs(1)/x0-1)*exp(-cs(2)*(t-t0); cs=lsqcurvefit(xt,rand(2,1),tt(2,end),a(2:end) %以上程序为三维视图和等高线画图程序clc, clear a=load(data5B.txt); amin=min(min(a), amax=max(max(a) x0=0:5:4000; y0=0:5:3000; xx0,yy0=meshgrid(x0,y0); mesh(xx0,yy0,a) v=-49,1,12,44,75,107,171,184:80:357,357:100:614,614; figure,
7、 c=contour(x0,y0,a,v); clabel(c) 宋哲 2011-08-17 9:32:44 clc, clear %该程序计算最小面积a=load(data5B.txt); a=a; x0=0:5:4000; y0=0:5:3000; pp=csape(x0,y0,a) x=0:4000; y=0:3000; z=fnval(pp,x,y); m,n=size(z); s=0; for i=1:m-1 for j=1:n-1 p1=x(i),y(j),z(i,j); p2=x(i+1),y(j),z(i+1,j); p3=x(i+1),y(j+1),z(i+1,j+1); p
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022年ARMA模型程序归 2022 ARMA 模型 程序
限制150内