��ǰλ�ã��Ŀ����� > ���з��� > �ߵȽ��� > ��ѧ > ����糡����ϵͳ�Ŀɿ�������Matlab����
������ش��ĵ���ȨͶ��

����糡����ϵͳ�Ŀɿ�������Matlab����

��Ҫ�������糡����ϵͳ�Ŀɿ�����������,����ʱ�����ؿ���ģ�ⷨ��

clc
clear

year=50 %ģ�������

for l=1:year;

%%%%%%%%%%%%%%%%%%%%%%%%% ��һ�������ݵ�����Ԥ���� %%%%%%%%%%%%%%%%%%%

%SW0=load('windspeed.txt'); %����ԭʼ��������
SW0=xlsread('windspeed.xls');
SW0=SW0';
SW0=SW0/10*3.6; %ԭʼ���ݵķ��ٵ�λΪ0.1m/s������ת��Ϊkm/h
N=size(SW0,2);
mu=mean(SW0);
sigma=var(SW0);
sigma=sigma^0.5; %��������ƽ��ֵ�ͱ�׼��
y=(SW0-mu)./sigma; %����Ԥ����

% figure(1);
% subplot(211);
% autocorr(y); %���������ͼ
% title('�����ͼ');
% subplot(212);
% parcorr(y); %����ƫ�����ͼ
% title('ƫ���ͼ');


%%%%%%%%%%%%%% �ڶ���������AIC׼��ȷ��ARMAģ�͵Ľ��� %%%%%%%%%%%%%%%%%%

for n=2:7;
m=armax(y',[n,n-1]);
fai=-m.a;
theta=m.c; %��armax�����õ��IJ�����ȡ����

for i=1:n;
y1(i)=y(i);
end

Noise=m.NoiseVariance^0.5;
e=normrnd(0,Noise,1,N);
for t=n+1:1:N;
y1(t)=0;
for j=2:n+1;
y1(t)=y1(t)+fai(j)*y1(t-(j-1));
end
for k=1:n;
y1(t)=y1(t)+theta(k)*e(t-(k-1));
end
end %y1(t)ΪԤ��ֵ

s(n)=0;
for i1=1:N;
residual=y1(i1)-y(i1);
s(n)=s(n)+residual^2; %��ȡ�в�ƽ����
end

AIC(n)=N*log(s(n))+2*n-1; %��AIC
end

arma=AIC(1,2:7); %n=2-7ʱ����AIC��ֵ
[AIC,n1]=min(arma);
n1=n1+1; %n1Ϊ�õ���ARMAģ�͵Ľ���

%%%%%%%%%%%%%%%%% ����������ARMAģ��Ԥ����ٲ�ȷ�����ٷֲ� %%%%%%%%%%%%%%%%%

m=armax(y',[n1,n1-1]);
fai=-m.a;
theta=m.c;

for i=1:n1;
y2(i)=y(i);
end

Noise=m.NoiseVariance^0.5;
e=normrnd(0,Noise,1,8736);
for t=n1+1:1:8736;
y2(t)=0;
for j=2:n1+1;
y2(t)=y2(t)+fai(j)*y2(t-(j-1));
end
for k=1:n1;
y2(t)=y2(t)+theta(k)*e(t-(k-1));
end
end %y1(t)ΪԤ��ֵ
y2;
SW=y2.*sigma+mu;

% figure(2)
% subplot(121);
% hist(SW0,100);
% xlabel('����');ylabel('Ƶ��');title('ԭʼ���ٵķֲ�');
% subplot(122);
% hist(SW,100);
% xlabel('����');ylabel('Ƶ��');title('Ԥ����ٵķֲ�');


%�õ�8736��Сʱ��Ԥ����٣�SW 1*8736


%%%%%%% ���IJ�����糡��ת�ƹ��̣�ȷ����һ��������״̬�ֱ���ڵ�ʱ��%%%%%%%

%���������״̬ģ��
l


ambdaRD=5.84;lambdaRF=7.96;lambdaDR=48.3;lambdaFR=58.4;lambdaDF=0;lambdaFD=0; %�����������3��״̬��ת����
T=[1 1 1;lambdaRD -lambdaDR-lambdaDF lambdaFD;

��1ҳ

�������Word�ĵ�������أ�����糡����ϵͳ�Ŀɿ�������Matlab����

������1-4ҳ����4ҳ��

����ϲ��

���ض���