RBF_model.rar_Matlab两相_两相流_两相流 MATLAB_神经网络_网络流

%数据载入 cl_1_2_01=load('C:\6.1data\cl_1_2_01.TXT'); cl_1_2_02=load('C:\6.1data\cl_1_2_02.TXT'); cl_1_2_03=load('C:\6.1data\cl_1_2_11.TXT'); cl_1_2=[cl_1_2_01;cl_1_2_02;cl_1_2_03]; cl_1_3_01=load('C:\6.1data\cl_1_3_01.TXT'); cl_1_3_02=load('C:\6.1data\cl_1_3_02.TXT'); cl_1_3_03=load('C:\6.1data\cl_1_3_11.TXT'); cl_1_3=[cl_1_3_01;cl_1_3_02;cl_1_3_03]; cl_2_3_01=load('C:\6.1data\cl_2_3_01.TXT'); cl_2_3_02=load('C:\6.1data\cl_2_3_02.TXT'); cl_2_3_03=load('C:\6.1data\cl_2_3_11.TXT'); cl_2_3=[cl_2_3_01;cl_2_3_02;cl_2_3_03]; cl=[cl_1_2;cl_1_3;cl_2_3]; hxl_1_1_01=load('C:\6.1data\hxl_01.TXT'); hxl_1_1_02=load('C:\6.1data\hxl_02.TXT'); hxl_1_1_03=load('C:\6.1data\hxl_03.TXT'); hxl_1_1_04=load('C:\6.1data\hxl_04.TXT'); hxl_1_1_05=load('C:\6.1data\hxl_05.TXT'); hxl_1_1_06=load('C:\6.1data\hxl_06.TXT'); hxl=[hxl_1_1_01;hxl_1_1_02;hxl_1_1_03;hxl_1_1_04;hxl_1_1_05;hxl_1_1_06]; hzl_1_1_01=load('C:\6.1data\hzl_01.TXT'); hzl_1_1_02=load('C:\6.1data\hzl_02.TXT'); hzl_1_1_03=load('C:\6.1data\hzl_03.TXT'); hzl_1_1_04=load('C:\6.1data\hzl_04.TXT'); hzl_1_1_05=load('C:\6.1data\hzl_05.TXT'); hzl_1_1_06=load('C:\6.1data\hzl_06.TXT'); hzl_1_1_07=load('C:\6.1data\hzl_07.TXT'); hzl_1_1_08=load('C:\6.1data\hzl_08.TXT'); hzl_1_1_09=load('C:\6.1data\hzl_09.TXT'); hzl_1_1_10=load('C:\6.1data\hzl_10.TXT'); hzl=[hzl_1_1_01;hzl_1_1_02;hzl_1_1_03;hzl_1_1_04;hzl_1_1_05;hzl_1_1_06;hzl_1_1_07;hzl_1_1_08;hzl_1_1_09;hzl_1_1_10]; %以上数据为6.1号采集数据，初始化空管值为kong_csh pzl_1_1_01=load('C:\6.1data\pzl_x2_d5_16_01.TXT'); pzl_1_1_02=load('C:\6.1data\pzl_x2_d5_16_02.TXT'); pzl_1_1_03=load('C:\6.1data\pzl_x2_d5_16_03.TXT'); pzl_1_1_04=load('C:\6.1data\pzl_x3_d7_16_01.TXT'); pzl_1_1_05=load('C:\6.1data\pzl_x5_d2_16_01.TXT'); pzl_1_1_06=load('C:\6.1data\pzl_x5_d2_16_02.TXT'); pzl_1_1_07=load('C:\6.1data\pzl_x5_d2_16_03.TXT'); pzl_1_1_08=load('C:\6.1data\pzl_x7_d3_16_01.TXT'); pzl=[pzl_1_1_01;pzl_1_1_02;pzl_1_1_03;pzl_1_1_04;pzl_1_1_05;pzl_1_1_06;pzl_1_1_07;pzl_1_1_08]; kong_01=load('C:\6.1data\kong_01.TXT'); kong_02=load('C:\6.1data\kong_02.TXT'); kong_03=load('C:\6.1data\kong_03.TXT'); kong_04=load('C:\6.1data\kong_04.TXT'); kong_05=load('C:\6.1data\kong_05.TXT'); kong_06=load('C:\6.1data\kong_06.TXT'); kong=[kong_01;kong_02;kong_03;kong_04;kong_05;kong_06]; man_01=load('C:\6.1data\man_01.TXT'); man_02=load('C:\6.1data\man_02.TXT'); man_03=load('C:\6.1data\man_03.TXT'); man_04=load('C:\6.1data\man_04.TXT'); man_05=load('C:\6.1data\man_05.TXT'); %man_06=load('C:\6.1data\man_06.TXT'); man=[man_01;man_02;man_03;man_04;man_05]; %求初始化的值(5.29号采集的数据)，6.1号采集的数据初始化空管值为kong_csh kong_csh=load('C:\6.1data\kong_csh_61.txt'); kong2_csh=mean(kong); %求变化值 %diff_cl_1_2=dif(cl_1_2,kong_csh); %diff_cl_1_3=dif(cl_1_3,kong_csh); %diff_cl_2_3=dif(cl_2_3,kong_csh); diff_cl=dif(cl,kong_csh); diff_hxl=dif(hxl,kong_csh); diff_hzl=dif(hzl,kong_csh); diff_pzl=dif(pzl,kong2_csh); diff_kong=dif(kong,kong2_csh); diff_man=dif(man,kong2_csh); %求特征参数 %adm_diff_cl_1_2=adm(diff_cl_1_2); %adm_diff_cl_1_3=adm(diff_cl_1_3); %adm_diff_cl_2_3=adm(diff_cl_2_3); adm_diff_cl=adm(diff_cl); adm_diff_hxl=adm(diff_hxl); adm_diff_hzl=adm(diff_hzl); adm_diff_pzl=adm(diff_pzl); adm_diff_kong=adm(diff_kong); adm_diff_man=adm(diff_man); %ads_diff_cl_1_2=ads(diff_cl_1_2); %ads_diff_cl_1_3=ads(diff_cl_1_3); %ads_diff_cl_2_3=ads(diff_cl_2_3); ads_diff_cl=ads(diff_cl); ads_diff_hxl=ads(diff_hxl); ads_diff_hzl=ads(diff_hzl); ads_diff_pzl=ads(diff_pzl); ads_diff_kong=ads(diff_kong); ads_diff_man=ads(diff_man); %avm_diff_cl_1_2=avm(diff_cl_1_2); %avm_diff_cl_1_3=avm(diff_cl_1_3); %avm_diff_cl_2_3=avm(diff_cl_2_3); avm_diff_cl=avm(diff_cl); avm_diff_hxl=avm(diff_hxl); avm_diff_hzl=avm(diff_hzl); avm_diff_pzl=avm(diff_pzl); avm_diff_kong=avm(diff_kong); avm_diff_man=avm(diff_man); %sfm1_diff_cl_1_2=sfm1(diff_cl_1_2); %sfm1_diff_cl_1_3=sfm1(diff_cl_1_3); %sfm1_diff_cl_2_3=sfm1(diff_cl_2_3); sfm1_diff_cl=sfm1(diff_cl); sfm1_diff_hxl=sfm1(diff_hxl); sfm1_diff_hzl=sfm1(diff_hzl); sfm1_diff_pzl=sfm1(diff_pzl); sfm1_diff_kong=sfm1(diff_kong); sfm1_diff_man=sfm1(diff_man); %sfm2_diff_cl_1_2=sfm2(diff_cl_1_2); %sfm2_diff_cl_1_3=sfm2(diff_cl_1_3); %sfm2_diff_cl_2_3=sfm2(diff_cl_2_3); sfm2_diff_cl=sfm2(diff_cl); sfm2_diff_hxl=sfm2(diff_hxl); sfm2_diff_hzl=sfm2(diff_hzl); sfm2_diff_pzl=sfm2(diff_pzl); sfm2_diff_kong=sfm2(diff_kong); sfm2_diff_man=sfm2(diff_man); %tbr_diff_cl=tbr(diff_cl); %tbr_diff_hxl=tbr(diff_hxl); %tbr_diff_hzl=tbr(diff_hzl); %tbr_diff_pzl=tbr(diff_pzl); %tbr_diff_kong=tbr(diff_kong); %tbr_diff_man=tbr(diff_man); %rm_diff_cl_1_2=rm(diff_cl_1_2); %rm_diff_cl_1_3=rm(diff_cl_1_3); %rm_diff_cl_2_3=rm(diff_cl_2_3); %rm_diff_hxl=rm(diff_hxl); %rm_diff_hzl=rm(diff_hzl); %rm_diff_pzl=rm(diff_pzl); %rm_diff_kong=rm(diff_kong); %rm_diff_man=rm(diff_man); %sf_diff_cl_1_2=sf(diff_cl_1_2); %sf_diff_cl_1_3=sf(diff_cl_1_3); %sf_diff_cl_2_3=sf(diff_cl_2_3); %sf_diff_hxl_1_1=sf(diff_hxl_1_1); %sf_diff_hzl_1_1=sf(diff_hzl_1_1); %sf_diff_pzl_1_1=sf(diff_pzl_1_1); %sf_diff_kong=sf(diff_kong); %sf_diff_man=sf(diff_man); %fm_diff_cl_1_2=fm(diff_cl_1_2); %fm_diff_cl_1_3=fm(diff_cl_1_3); %fm_diff_cl_2_3=fm(diff_cl_2_3); fm_diff_cl=fm(diff_cl); fm_diff_hxl=fm(diff_hxl); fm_diff_hzl=fm(diff_hzl); fm_diff_pzl=fm(diff_pzl); fm_diff_kong=fm(diff_kong); fm_diff_man=fm(diff_man); %fs_diff_cl_1_2=fs(diff_cl_1_2); %fs_diff_cl_1_3=fs(diff_cl_1_3); %fs_diff_cl_2_3=fs(diff_cl_2_3); fs_diff_cl=fs(diff_cl); fs_diff_hxl=fs(diff_hxl); fs_diff_hzl=fs(diff_hzl); fs_diff_pzl=fs(diff_pzl); fs_diff_kong=fs(diff_kong); fs_diff_man=fs(diff_man); %sam_diff_cl_1_2=sam(diff_cl_1_2); %sam_diff_cl_1_3=sam(diff_cl_1_3); %sam_diff_cl_2_3=sam(diff_cl_2_3); sam_diff_cl=sam(diff_cl); sam_diff_hxl=sam(diff_hxl); sam_diff_hzl=sam(diff_hzl); sam_diff_pzl=sam(diff_pzl); sam_diff_kong=sam(diff_kong); sam_diff_man=sam(diff_man); %sas_diff_cl_1_2=sas(diff_cl_1_2); %sas_diff_cl_1_3=sas(diff_cl_1_3); %sas_diff_cl_2_3=sas(diff_cl_2_3); sas_diff_cl=sas(diff_cl); sas_diff_hxl=sas(diff_hxl); sas_diff_hzl=sas(diff_hzl); sas_diff_pzl=sas(diff_pzl); sas_diff_kong=sas(diff_kong); sas_diff_man=sas(diff_man); %特征参数组成的辨识矩阵 bsjz_cl=[adm_diff_cl,ads_diff_cl,avm_diff_cl,sfm1_diff_cl,sfm2_diff_cl,fm_diff_cl,fs_diff_cl,sam_diff_cl,sas_diff_cl]'; %bsjz_cl_1_2=[adm_diff_cl_1_2,ads_diff_cl_1_2,avm_diff_cl_1_2,sfm1_diff_cl_1_2,sfm2_diff_cl_1_2,fm_diff_cl_1_2,fs_diff_cl_1_2,sam_diff_cl_1_2,sas_diff_cl_1_2,rm_diff_cl_1_2]'; %bsjz_cl_1_3=[adm_diff_cl_1_3,ads_diff_cl_1_3,avm_diff_cl_1_3,sfm1_diff_cl_1_3,sfm2_diff_cl_1_3,fm_diff_cl_1_3,fs_diff_cl_1_3,sam_diff_cl_1_3,sas_diff_cl_1_3,rm_diff_cl_1_3]'; %bsjz_cl_2_3=[adm_diff_cl_2_3,ads_diff_cl_2_3,avm_diff_cl_2_3,sfm1_diff_cl_2_3,sfm2_diff_cl_2_3,fm_diff_cl_2_3,fs_diff_cl_2_3,sam_diff_cl_2_3,sas_diff_cl_2_3,rm_diff_cl_2_3]'; bsjz_hxl=[adm_diff_hxl,ads_diff_hxl,avm_diff_hxl,sfm1_diff_hxl,sfm2_diff_hxl,fm_diff_hxl,fs_diff_hxl,sam_diff_hxl,sas_diff_hxl]'; bsjz_hzl=[adm_diff_hzl,ads_diff_hzl,avm_diff_hzl,sfm1_diff_hzl,sfm2_diff_hzl,fm_diff_hzl,fs_diff_hzl,sam_diff_hzl,sas_diff_hzl]'; bsjz_pzl=[adm_diff_pzl,ads_diff_pzl,avm_diff_pzl,sfm1_diff_pzl,sfm2_diff_pzl,fm_diff_pzl,fs_diff_pzl,sam_diff_pzl,sas_diff_pzl]'; bsjz_kong=[adm_diff_kong,ads_diff_kong,avm_diff_kong,sfm1_diff_kong,sfm2_diff_kong,fm_diff_kong,fs_diff_kong,sam_diff_kong,sas_diff_kong]'; bsjz_man=[adm_diff_man,ads_diff_man,avm_diff_man,sfm1_diff_man,sfm2_diff_man,fm_diff_man,fs_diff_man,sam_diff_man,sas_diff_man]'; %随即抽取80组作为训练，20组作为预测 %调用自己编写的随机抽取函数：sjcq %[train_cl_1_2,predict_cl_1_2]=sjcq(bsjz_cl_1_2,50,12); %[train_cl_1_3,predict_cl_1_3]=sjcq(bsjz_cl_1_3,50,7); %[train_cl_2_3,predict_cl_2_3]=sjcq(bsjz_cl_2_3,50,20); [train_cl,predict_cl]=sjcq(bsjz_cl,80,20); [train_hxl,predict_hxl]=sjcq(bsjz