资源描述:
《arma模型程序》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、clear,clccloseallset(0,'defaultaxeslinestyleorder',{'',':+','-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
2、);%画图并且取得二维图形中最大值和最小值的索引号gridonset(gca,'layer','top')%将栅格放在最上层,防止遮盖axisequal%使每一格的刻度相同,同时也使长度一样xlabel('-pileqthetaleqpi','fontsize',20)ylabel('tan(sin(theta))-sin(tan(theta))','fontsize',20)title('it{-pileqthetaleqpi的三角函数图}','fontsize',20)text(0,0,'bulletit原点','fontsize'
3、,25)%文本的精确定位%标注曲线的最大值和最小值x=get(h,'xdata');%获得二维曲线的数据y=get(h,'ydata');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))],...'
4、horizontalalignment','center','verticalalignment','top','fontsize',17)宋哲2011/8/1611:07:42%AR,MAI型clear,clcclosealla=load('');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(1
5、22),parcorr(b)cs1=ar(b,2)%AR模型,b必须是列向量bhat=predict(cs1,[b;20])%Discrete-timeIDPOLYmodel:%A(q)=1-qA-1-qA-2b15=bhat{1}(end)bhat2=predict(cs1,[b;b15;-1])宋哲2011/8/1611:27:00%用逻辑斯特预测clear,clca=load('');a=a([2,4],:);a=a';a=a(:);%将数据税收展成按时间序列排列的x0=a(1);tt=[1:14]';t0=1;xt=@(cs,t)cs(1)/(1+(cs(1
6、)/x0-1)exp(-cs(2)(t-t0)));cs=lsqcurvefit(xt,rand(2,1),tt(2,end),a(2:end))%以上程序为三维视图和等高线画图程序clc,cleara=load('');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,c=contour(x0,y0,
7、a,v);clabel(c)宋哲2011-08-179:32:44clc,clear%该程序计算最小面积a=load('');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;fori=1:m-1forj=1:n-1p1=[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)];p4=[x(i),y(j+1),z(i,j