可靠性form 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函数得到的参数，取出来

for i=1:n;

y1(i