RLS最小二乘法线性拟合_最小二乘线性拟合_源码

%% I. 清空环境变量 clear all clc %% II. 训练集/测试集产生 %% % 1. 导入数据 load 'data83' Ts=1;%模型辨识采样时间 P06MCV03AI=data(:,2); % 7820*1 double 送风机A动叶调节阀位反馈 P06MCV03BI=data(:,3); % 7820*1 double 送风机B动叶调节阀位反馈 P07MCV29AI=data(:,8); % 7820*1 double 吸风机A动叶调节阀位反馈 P07MCV29BI=data(:,9); % 7820*1 double 吸风机B动叶调节阀位反馈 TOTALCOALFL=data(:,14); % 7820*1 double 总给煤瞬时量 SECAIRFLW=data(:,21); % 7820*1 double 总二次风量 P06MCV03I=(P06MCV03AI+P06MCV03BI)/2; figure(1) subplot(3,1,1) plot(P06MCV03I) subplot(3,1,2) plot(P07MCV29AI) subplot(3,1,3) plot(P07MCV29BI) figure(2) subplot(2,1,1) plot(TOTALCOALFL) subplot(2,1,2) plot(SECAIRFLW) global P06MCV03I_mean_1 P06MCV03AI_mean_1 P06MCV03BI_mean_1 P07MCV29AI_mean_1 P07MCV29BI_mean_1 TOTALCOALFL_mean_1 BLRTATOLLFLOW_mean_1 SECAIRFLW_mean_1 %% III.是否在线优化PID参数 prompt='选择是否在线优化PID参数：1：在线优化PID参数，2：使用离线优化得到的PID参数\n'; flag_0=input(prompt); if flag_0==2 fprintf('A导前区PID控制器参数：\n ') disp('kp1=0.5') disp('ki1=1.5') disp('kd1=0') disp('ka1=10') disp('k1=0.5') else %% 步骤1：模型判阶 %% prompt='选择是否自动判断模型阶数，1：自动判阶，2：人工判阶\n'; flag_panjie=input(prompt); Nf=SECAIRFLW; %滤波 b=0.04; a=[1 -0.96]; SECAIRFLW=filter(b,a,Nf); %% A导前区判阶 Udata1=[P06MCV03I P07MCV29AI P07MCV29BI TOTALCOALFL]; %输入 Ydata1=SECAIRFLW; %输出 [~,nu1]=size(Udata1); [~,ny1]=size(Ydata1); Nqs1 = 234001; %数据开始位置123001 NN = 25000; %数据长度 NN1 = Nqs1+NN; %数据结束位置 figure(1) n1 = Orderdeter(Udata1,Ydata1,Nqs1,NN,'Picturing',1); %n1:A导前区模型阶数 if flag_panjie==1 %自动判断模型阶数 else prompt='输入A导前区阶数\n'; n1=input(prompt); end %% 步骤2：模型辨识 %% prompt='选择辨识算法，1：递推最小二乘法，2：预报误差法\n'; flag_bianshi=input(prompt); %% 辨识方法1：RLS 递推最小二乘法，辨识A导前区 if flag_bianshi==1 figure(4) [theta1,fits1]= RLS_1(Udata1,Ydata1,Nqs1,NN,'order',n1,'theta1',1e-6*ones(n1*(nu1+ny1),1),'P1',1e6*eye(n1*(nu1+ny1)),'Picturing',1,'Averaging',1);%A导前区模型 title('A导前区') end %% 辨识方法2：PEM 预报误差法，辨识A导前区 if flag_bianshi==2 [theta1,fits1]=PEM(Udata1,Ydata1,Nqs1,NN,'order',n1,'theta1',1e-6*ones(n1*(nu1+ny1),1),'Picturing',1,'Averaging',1); end %% 均值计算，A导前区均值 if flag_bianshi==1 || flag_bianshi==2 P06MCV03AI_mean_1=mean(P06MCV03AI(Nqs1:NN1)); P06MCV03BI_mean_1=mean(P06MCV03BI(Nqs1:NN1)); P07MCV29AI_mean_1=mean(P07MCV29AI(Nqs1:NN1)); P07MCV29BI_mean_1=mean(P07MCV29BI(Nqs1:NN1)); TOTALCOALFL_mean_1=mean(TOTALCOALFL(Nqs1:NN1)); SECAIRFLW_mean_1=mean(SECAIRFLW(Nqs1:NN1)); P06MCV03I_mean_1=mean(P06MCV03I(Nqs1:NN1)); %% 仿真模型构造 %输入1对输出1传递函数：numerator（:,:,1）/d1; 第nu个输入对输出传递函数=numerator（:,:,nu）/dnu [numerator1,d1]=model_construction(n1,theta1,nu1); %A导前区均值 end end