close all, clear, clc % loading Data for Test A B1a=load('Batch1a.mat') B1b=load('Batch1b.mat') Ba=load('Batch_SCR_5_7') Bb=load('Batch_SCR_3_1p5') % Data manipulation for L error half_per=1.6; t_3p2=13:half_per:189; Ln_act_Ba=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_3p2); Ln_est_B1a=interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_3p2); Error_Ln_B1a=(Ln_act_Ba-Ln_est_B1a)./Ln_act_Ba.*100; t_int_3p2_SEG1=13:1.6:21; Ln_B1a_SEG1 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG1); Ln_act_B1a_SEG1=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG1); Err_Ln_B1a_SEG1=mean((Ln_act_B1a_SEG1-Ln_B1a_SEG1)./Ln_act_B1a_SEG1.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Ln_B1a_SEG2 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG2); Ln_act_B1a_SEG2=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG2); Err_Ln_B1a_SEG2=mean((Ln_act_B1a_SEG2-Ln_B1a_SEG2)./Ln_act_B1a_SEG2.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Ln_B1a_SEG3 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG3); Ln_act_B1a_SEG3=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG3); Err_Ln_B1a_SEG3=mean((Ln_act_B1a_SEG2-Ln_B1a_SEG3)./Ln_act_B1a_SEG3.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Ln_B1a_SEG4 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG4); Ln_act_B1a_SEG4=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG4); Err_Ln_B1a_SEG4=mean((Ln_act_B1a_SEG4-Ln_B1a_SEG4)./Ln_act_B1a_SEG4.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Ln_B1a_SEG5 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG5); Ln_act_B1a_SEG5=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG5); Err_Ln_B1a_SEG5=mean((Ln_act_B1a_SEG5-Ln_B1a_SEG5)./Ln_act_B1a_SEG5.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Ln_B1a_SEG6 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG6); Ln_act_B1a_SEG6=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG6); Err_Ln_B1a_SEG6=mean((Ln_act_B1a_SEG6-Ln_B1a_SEG6)./Ln_act_B1a_SEG6.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Ln_B1a_SEG7 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG7); Ln_act_B1a_SEG7=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG7); Err_Ln_B1a_SEG7=mean((Ln_act_B1a_SEG7-Ln_B1a_SEG7)./Ln_act_B1a_SEG7.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Ln_B1a_SEG8 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG8); Ln_act_B1a_SEG8=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG8); Err_Ln_B1a_SEG8=mean((Ln_act_B1a_SEG8-Ln_B1a_SEG8)./Ln_act_B1a_SEG8.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Ln_B1a_SEG9 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG9); Ln_act_B1a_SEG9=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG9); Err_Ln_B1a_SEG9=mean((Ln_act_B1a_SEG9-Ln_B1a_SEG9)./Ln_act_B1a_SEG9.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Ln_B1a_SEG10 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG10); Ln_act_B1a_SEG10=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG10); Err_Ln_B1a_SEG10=mean((Ln_act_B1a_SEG10-Ln_B1a_SEG10)./Ln_act_B1a_SEG10.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Ln_B1a_SEG11 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG11); Ln_act_B1a_SEG11=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG11); Err_Ln_B1a_SEG11=mean((Ln_act_B1a_SEG11-Ln_B1a_SEG11)./Ln_act_B1a_SEG11.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Ln_B1a_SEG12 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG12); Ln_act_B1a_SEG12=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG12); Err_Ln_B1a_SEG12=mean((Ln_act_B1a_SEG12-Ln_B1a_SEG12)./Ln_act_B1a_SEG12.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Ln_B1a_SEG13 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG13); Ln_act_B1a_SEG13=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG13); Err_Ln_B1a_SEG13=mean((Ln_act_B1a_SEG13-Ln_B1a_SEG13)./Ln_act_B1a_SEG13.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Ln_B1a_SEG14 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG14); Ln_act_B1a_SEG14=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG14); Err_Ln_B1a_SEG14=mean((Ln_act_B1a_SEG14-Ln_B1a_SEG14)./Ln_act_B1a_SEG14.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Ln_B1a_SEG15 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG15); Ln_act_B1a_SEG15=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG15); Err_Ln_B1a_SEG15=mean((Ln_act_B1a_SEG15-Ln_B1a_SEG15)./Ln_act_B1a_SEG15.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Ln_B1a_SEG16 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG16); Ln_act_B1a_SEG16=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG16); Err_Ln_B1a_SEG16=mean((Ln_act_B1a_SEG16-Ln_B1a_SEG16)./Ln_act_B1a_SEG16.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Ln_B1a_SEG17 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG17); Ln_act_B1a_SEG17=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG17); Err_Ln_B1a_SEG17=mean((Ln_act_B1a_SEG17-Ln_B1a_SEG17)./Ln_act_B1a_SEG17.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Ln_B1a_SEG18 = interp1(B1a.Ln_est.Time,B1a.Ln_est.Data,t_int_3p2_SEG18); Ln_act_B1a_SEG18=interp1(Ba.Ln_act.Time,Ba.Ln_act.Data,t_int_3p2_SEG18); Err_Ln_B1a_SEG18=mean((Ln_act_B1a_SEG18-Ln_B1a_SEG18)./Ln_act_B1a_SEG18.*100); Error_Ln_B1a_ALL=[Err_Ln_B1a_SEG1 Err_Ln_B1a_SEG2 Err_Ln_B1a_SEG3 Err_Ln_B1a_SEG4 Err_Ln_B1a_SEG5 Err_Ln_B1a_SEG6 Err_Ln_B1a_SEG7 Err_Ln_B1a_SEG8 Err_Ln_B1a_SEG9 Err_Ln_B1a_SEG10 Err_Ln_B1a_SEG11 Err_Ln_B1a_SEG12 Err_Ln_B1a_SEG13 Err_Ln_B1a_SEG14 Err_Ln_B1a_SEG15 Err_Ln_B1a_SEG16 Err_Ln_B1a_SEG17 Err_Ln_B1a_SEG18]; Test_A_SCR_7_Lerr=Error_Ln_B1a_ALL(1:9) Test_A_SCR_5_Lerr=Error_Ln_B1a_ALL(10:end) % Now R error Rn_act_Ba=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_3p2); Rn_est_B1a=interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_3p2); Error_Rn_B1a=(Rn_act_Ba-Rn_est_B1a)./Rn_act_Ba.*100; t_int_3p2_SEG1=13:1.6:21; Rn_B1a_SEG1 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG1); Rn_act_B1a_SEG1=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG1); Err_Rn_B1a_SEG1=mean((Rn_act_B1a_SEG1-Rn_B1a_SEG1)./Rn_act_B1a_SEG1.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Rn_B1a_SEG2 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG2); Rn_act_B1a_SEG2=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG2); Err_Rn_B1a_SEG2=mean((Rn_act_B1a_SEG2-Rn_B1a_SEG2)./Rn_act_B1a_SEG2.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Rn_B1a_SEG3 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG3); Rn_act_B1a_SEG3=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG3); Err_Rn_B1a_SEG3=mean((Rn_act_B1a_SEG2-Rn_B1a_SEG3)./Rn_act_B1a_SEG3.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Rn_B1a_SEG4 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG4); Rn_act_B1a_SEG4=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG4); Err_Rn_B1a_SEG4=mean((Rn_act_B1a_SEG4-Rn_B1a_SEG4)./Rn_act_B1a_SEG4.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Rn_B1a_SEG5 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG5); Rn_act_B1a_SEG5=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG5); Err_Rn_B1a_SEG5=mean((Rn_act_B1a_SEG5-Rn_B1a_SEG5)./Rn_act_B1a_SEG5.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Rn_B1a_SEG6 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG6); Rn_act_B1a_SEG6=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG6); Err_Rn_B1a_SEG6=mean((Rn_act_B1a_SEG6-Rn_B1a_SEG6)./Rn_act_B1a_SEG6.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Rn_B1a_SEG7 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG7); Rn_act_B1a_SEG7=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG7); Err_Rn_B1a_SEG7=mean((Rn_act_B1a_SEG7-Rn_B1a_SEG7)./Rn_act_B1a_SEG7.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Rn_B1a_SEG8 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG8); Rn_act_B1a_SEG8=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG8); Err_Rn_B1a_SEG8=mean((Rn_act_B1a_SEG8-Rn_B1a_SEG8)./Rn_act_B1a_SEG8.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Rn_B1a_SEG9 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG9); Rn_act_B1a_SEG9=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG9); Err_Rn_B1a_SEG9=mean((Rn_act_B1a_SEG9-Rn_B1a_SEG9)./Rn_act_B1a_SEG9.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Rn_B1a_SEG10 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG10); Rn_act_B1a_SEG10=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG10); Err_Rn_B1a_SEG10=mean((Rn_act_B1a_SEG10-Rn_B1a_SEG10)./Rn_act_B1a_SEG10.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Rn_B1a_SEG11 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG11); Rn_act_B1a_SEG11=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG11); Err_Rn_B1a_SEG11=mean((Rn_act_B1a_SEG11-Rn_B1a_SEG11)./Rn_act_B1a_SEG11.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Rn_B1a_SEG12 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG12); Rn_act_B1a_SEG12=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG12); Err_Rn_B1a_SEG12=mean((Rn_act_B1a_SEG12-Rn_B1a_SEG12)./Rn_act_B1a_SEG12.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Rn_B1a_SEG13 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG13); Rn_act_B1a_SEG13=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG13); Err_Rn_B1a_SEG13=mean((Rn_act_B1a_SEG13-Rn_B1a_SEG13)./Rn_act_B1a_SEG13.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Rn_B1a_SEG14 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG14); Rn_act_B1a_SEG14=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG14); Err_Rn_B1a_SEG14=mean((Rn_act_B1a_SEG14-Rn_B1a_SEG14)./Rn_act_B1a_SEG14.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Rn_B1a_SEG15 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG15); Rn_act_B1a_SEG15=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG15); Err_Rn_B1a_SEG15=mean((Rn_act_B1a_SEG15-Rn_B1a_SEG15)./Rn_act_B1a_SEG15.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Rn_B1a_SEG16 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG16); Rn_act_B1a_SEG16=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG16); Err_Rn_B1a_SEG16=mean((Rn_act_B1a_SEG16-Rn_B1a_SEG16)./Rn_act_B1a_SEG16.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Rn_B1a_SEG17 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG17); Rn_act_B1a_SEG17=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG17); Err_Rn_B1a_SEG17=mean((Rn_act_B1a_SEG17-Rn_B1a_SEG17)./Rn_act_B1a_SEG17.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Rn_B1a_SEG18 = interp1(B1a.Rn_est.Time,B1a.Rn_est.Data,t_int_3p2_SEG18); Rn_act_B1a_SEG18=interp1(Ba.Rn_act.Time,Ba.Rn_act.Data,t_int_3p2_SEG18); Err_Rn_B1a_SEG18=mean((Rn_act_B1a_SEG18-Rn_B1a_SEG18)./Rn_act_B1a_SEG18.*100); Error_Rn_B1a_ALL=[Err_Rn_B1a_SEG1 Err_Rn_B1a_SEG2 Err_Rn_B1a_SEG3 Err_Rn_B1a_SEG4 Err_Rn_B1a_SEG5 Err_Rn_B1a_SEG6 Err_Rn_B1a_SEG7 Err_Rn_B1a_SEG8 Err_Rn_B1a_SEG9 Err_Rn_B1a_SEG10 Err_Rn_B1a_SEG11 Err_Rn_B1a_SEG12 Err_Rn_B1a_SEG13 Err_Rn_B1a_SEG14 Err_Rn_B1a_SEG15 Err_Rn_B1a_SEG16 Err_Rn_B1a_SEG17 Err_Rn_B1a_SEG18]; Test_A_SCR_7_Rerr=Error_Rn_B1a_ALL(1:9) Test_A_SCR_5_Rerr=Error_Rn_B1a_ALL(10:end) % Data manipulation for SCRs 3 and 1.5 % L error Ln_act_Bb=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_3p2); Ln_est_B1b=interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_3p2); Error_Ln_B1b=(Ln_act_Bb-Ln_est_B1b)./Ln_act_Bb.*100; t_int_3p2_SEG1=13:1.6:21; Ln_B1b_SEG1 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG1); Ln_act_B1b_SEG1=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG1); Err_Ln_B1b_SEG1=mean((Ln_act_B1b_SEG1-Ln_B1b_SEG1)./Ln_act_B1b_SEG1.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Ln_B1b_SEG2 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG2); Ln_act_B1b_SEG2=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG2); Err_Ln_B1b_SEG2=mean((Ln_act_B1b_SEG2-Ln_B1b_SEG2)./Ln_act_B1b_SEG2.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Ln_B1b_SEG3 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG3); Ln_act_B1b_SEG3=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG3); Err_Ln_B1b_SEG3=mean((Ln_act_B1b_SEG2-Ln_B1b_SEG3)./Ln_act_B1b_SEG3.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Ln_B1b_SEG4 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG4); Ln_act_B1b_SEG4=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG4); Err_Ln_B1b_SEG4=mean((Ln_act_B1b_SEG4-Ln_B1b_SEG4)./Ln_act_B1b_SEG4.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Ln_B1b_SEG5 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG5); Ln_act_B1b_SEG5=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG5); Err_Ln_B1b_SEG5=mean((Ln_act_B1b_SEG5-Ln_B1b_SEG5)./Ln_act_B1b_SEG5.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Ln_B1b_SEG6 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG6); Ln_act_B1b_SEG6=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG6); Err_Ln_B1b_SEG6=mean((Ln_act_B1b_SEG6-Ln_B1b_SEG6)./Ln_act_B1b_SEG6.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Ln_B1b_SEG7 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG7); Ln_act_B1b_SEG7=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG7); Err_Ln_B1b_SEG7=mean((Ln_act_B1b_SEG7-Ln_B1b_SEG7)./Ln_act_B1b_SEG7.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Ln_B1b_SEG8 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG8); Ln_act_B1b_SEG8=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG8); Err_Ln_B1b_SEG8=mean((Ln_act_B1b_SEG8-Ln_B1b_SEG8)./Ln_act_B1b_SEG8.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Ln_B1b_SEG9 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG9); Ln_act_B1b_SEG9=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG9); Err_Ln_B1b_SEG9=mean((Ln_act_B1b_SEG9-Ln_B1b_SEG9)./Ln_act_B1b_SEG9.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Ln_B1b_SEG10 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG10); Ln_act_B1b_SEG10=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG10); Err_Ln_B1b_SEG10=mean((Ln_act_B1b_SEG10-Ln_B1b_SEG10)./Ln_act_B1b_SEG10.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Ln_B1b_SEG11 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG11); Ln_act_B1b_SEG11=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG11); Err_Ln_B1b_SEG11=mean((Ln_act_B1b_SEG11-Ln_B1b_SEG11)./Ln_act_B1b_SEG11.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Ln_B1b_SEG12 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG12); Ln_act_B1b_SEG12=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG12); Err_Ln_B1b_SEG12=mean((Ln_act_B1b_SEG12-Ln_B1b_SEG12)./Ln_act_B1b_SEG12.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Ln_B1b_SEG13 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG13); Ln_act_B1b_SEG13=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG13); Err_Ln_B1b_SEG13=mean((Ln_act_B1b_SEG13-Ln_B1b_SEG13)./Ln_act_B1b_SEG13.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Ln_B1b_SEG14 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG14); Ln_act_B1b_SEG14=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG14); Err_Ln_B1b_SEG14=mean((Ln_act_B1b_SEG14-Ln_B1b_SEG14)./Ln_act_B1b_SEG14.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Ln_B1b_SEG15 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG15); Ln_act_B1b_SEG15=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG15); Err_Ln_B1b_SEG15=mean((Ln_act_B1b_SEG15-Ln_B1b_SEG15)./Ln_act_B1b_SEG15.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Ln_B1b_SEG16 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG16); Ln_act_B1b_SEG16=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG16); Err_Ln_B1b_SEG16=mean((Ln_act_B1b_SEG16-Ln_B1b_SEG16)./Ln_act_B1b_SEG16.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Ln_B1b_SEG17 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG17); Ln_act_B1b_SEG17=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG17); Err_Ln_B1b_SEG17=mean((Ln_act_B1b_SEG17-Ln_B1b_SEG17)./Ln_act_B1b_SEG17.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Ln_B1b_SEG18 = interp1(B1b.Ln_est.Time,B1b.Ln_est.Data,t_int_3p2_SEG18); Ln_act_B1b_SEG18=interp1(Bb.Ln_act.Time,Bb.Ln_act.Data,t_int_3p2_SEG18); Err_Ln_B1b_SEG18=mean((Ln_act_B1b_SEG18-Ln_B1b_SEG18)./Ln_act_B1b_SEG18.*100); Error_Ln_B1b_ALL=[Err_Ln_B1b_SEG1 Err_Ln_B1b_SEG2 Err_Ln_B1b_SEG3 Err_Ln_B1b_SEG4 Err_Ln_B1b_SEG5 Err_Ln_B1b_SEG6 Err_Ln_B1b_SEG7 Err_Ln_B1b_SEG8 Err_Ln_B1b_SEG9 Err_Ln_B1b_SEG10 Err_Ln_B1b_SEG11 Err_Ln_B1b_SEG12 Err_Ln_B1b_SEG13 Err_Ln_B1b_SEG14 Err_Ln_B1b_SEG15 Err_Ln_B1b_SEG16 Err_Ln_B1b_SEG17 Err_Ln_B1b_SEG18]; Test_A_SCR_3_Lerr=Error_Ln_B1b_ALL(1:9) Test_A_SCR_1p5_Lerr=Error_Ln_B1b_ALL(10:end) % Now R error Rn_act_Bb=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_3p2); Rn_est_B1b=interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_3p2); Error_Rn_B1b=(Rn_act_Bb-Rn_est_B1b)./Rn_act_Bb.*100; t_int_3p2_SEG1=13:1.6:21; Rn_B1b_SEG1 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG1); Rn_act_B1b_SEG1=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG1); Err_Rn_B1b_SEG1=mean((Rn_act_B1b_SEG1-Rn_B1b_SEG1)./Rn_act_B1b_SEG1.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Rn_B1b_SEG2 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG2); Rn_act_B1b_SEG2=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG2); Err_Rn_B1b_SEG2=mean((Rn_act_B1b_SEG2-Rn_B1b_SEG2)./Rn_act_B1b_SEG2.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Rn_B1b_SEG3 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG3); Rn_act_B1b_SEG3=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG3); Err_Rn_B1b_SEG3=mean((Rn_act_B1b_SEG2-Rn_B1b_SEG3)./Rn_act_B1b_SEG3.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Rn_B1b_SEG4 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG4); Rn_act_B1b_SEG4=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG4); Err_Rn_B1b_SEG4=mean((Rn_act_B1b_SEG4-Rn_B1b_SEG4)./Rn_act_B1b_SEG4.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Rn_B1b_SEG5 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG5); Rn_act_B1b_SEG5=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG5); Err_Rn_B1b_SEG5=mean((Rn_act_B1b_SEG5-Rn_B1b_SEG5)./Rn_act_B1b_SEG5.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Rn_B1b_SEG6 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG6); Rn_act_B1b_SEG6=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG6); Err_Rn_B1b_SEG6=mean((Rn_act_B1b_SEG6-Rn_B1b_SEG6)./Rn_act_B1b_SEG6.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Rn_B1b_SEG7 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG7); Rn_act_B1b_SEG7=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG7); Err_Rn_B1b_SEG7=mean((Rn_act_B1b_SEG7-Rn_B1b_SEG7)./Rn_act_B1b_SEG7.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Rn_B1b_SEG8 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG8); Rn_act_B1b_SEG8=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG8); Err_Rn_B1b_SEG8=mean((Rn_act_B1b_SEG8-Rn_B1b_SEG8)./Rn_act_B1b_SEG8.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Rn_B1b_SEG9 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG9); Rn_act_B1b_SEG9=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG9); Err_Rn_B1b_SEG9=mean((Rn_act_B1b_SEG9-Rn_B1b_SEG9)./Rn_act_B1b_SEG9.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Rn_B1b_SEG10 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG10); Rn_act_B1b_SEG10=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG10); Err_Rn_B1b_SEG10=mean((Rn_act_B1b_SEG10-Rn_B1b_SEG10)./Rn_act_B1b_SEG10.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Rn_B1b_SEG11 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG11); Rn_act_B1b_SEG11=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG11); Err_Rn_B1b_SEG11=mean((Rn_act_B1b_SEG11-Rn_B1b_SEG11)./Rn_act_B1b_SEG11.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Rn_B1b_SEG12 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG12); Rn_act_B1b_SEG12=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG12); Err_Rn_B1b_SEG12=mean((Rn_act_B1b_SEG12-Rn_B1b_SEG12)./Rn_act_B1b_SEG12.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Rn_B1b_SEG13 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG13); Rn_act_B1b_SEG13=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG13); Err_Rn_B1b_SEG13=mean((Rn_act_B1b_SEG13-Rn_B1b_SEG13)./Rn_act_B1b_SEG13.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Rn_B1b_SEG14 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG14); Rn_act_B1b_SEG14=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG14); Err_Rn_B1b_SEG14=mean((Rn_act_B1b_SEG14-Rn_B1b_SEG14)./Rn_act_B1b_SEG14.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Rn_B1b_SEG15 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG15); Rn_act_B1b_SEG15=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG15); Err_Rn_B1b_SEG15=mean((Rn_act_B1b_SEG15-Rn_B1b_SEG15)./Rn_act_B1b_SEG15.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Rn_B1b_SEG16 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG16); Rn_act_B1b_SEG16=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG16); Err_Rn_B1b_SEG16=mean((Rn_act_B1b_SEG16-Rn_B1b_SEG16)./Rn_act_B1b_SEG16.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Rn_B1b_SEG17 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG17); Rn_act_B1b_SEG17=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG17); Err_Rn_B1b_SEG17=mean((Rn_act_B1b_SEG17-Rn_B1b_SEG17)./Rn_act_B1b_SEG17.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Rn_B1b_SEG18 = interp1(B1b.Rn_est.Time,B1b.Rn_est.Data,t_int_3p2_SEG18); Rn_act_B1b_SEG18=interp1(Bb.Rn_act.Time,Bb.Rn_act.Data,t_int_3p2_SEG18); Err_Rn_B1b_SEG18=mean((Rn_act_B1b_SEG18-Rn_B1b_SEG18)./Rn_act_B1b_SEG18.*100); Error_Rn_B1b_ALL=[Err_Rn_B1b_SEG1 Err_Rn_B1b_SEG2 Err_Rn_B1b_SEG3 Err_Rn_B1b_SEG4 Err_Rn_B1b_SEG5 Err_Rn_B1b_SEG6 Err_Rn_B1b_SEG7 Err_Rn_B1b_SEG8 Err_Rn_B1b_SEG9 Err_Rn_B1b_SEG10 Err_Rn_B1b_SEG11 Err_Rn_B1b_SEG12 Err_Rn_B1b_SEG13 Err_Rn_B1b_SEG14 Err_Rn_B1b_SEG15 Err_Rn_B1b_SEG16 Err_Rn_B1b_SEG17 Err_Rn_B1b_SEG18]; Test_A_SCR_3_Rerr=Error_Rn_B1b_ALL(1:9) Test_A_SCR_1p5_Rerr=Error_Rn_B1b_ALL(10:end) %% Test B clear % Load data B1a=load('BatchRLCRL1a.mat'); B1b=load('BatchRLCRL1b.mat'); B2a=load('BatchRLCRL2a.mat'); B2b=load('BatchRLCRL2b.mat'); B3a=load('BatchRLCRL3a.mat'); B3b=load('BatchRLCRL3b.mat'); % Load colours red = [244, 66, 66]./255; blue = [59, 104, 229]./255; green = [63, 150, 49]./255; purple = [255, 0, 233]./255; sky = [0, 212, 255]./255; orange = [255, 140, 0]./255; black = [0, 0, 0]./255; run('M12b_Resonances_in_RLCRL') close all half_per=1.6; t_3p2=13:half_per:189; t_checks = 0:0.1:190; Ln_calc_Ba_600Hz_steps=[Lg_2ap1_50Hz Lg_2ap1_50Hz Lg_2ap2_50Hz Lg_2ap2_50Hz]; % calculated in M12b_Resonances Timing_imp_step_change=[0 101 101.0001 189]; Ln_calc_Ba_600Hz=interp1(Timing_imp_step_change,Ln_calc_Ba_600Hz_steps,t_checks); Rn_calc_Ba_600Hz_steps=[Rg_2ap1_50Hz Rg_2ap1_50Hz Rg_2ap2_50Hz Rg_2ap2_50Hz]; Rn_calc_Ba_600Hz=interp1(Timing_imp_step_change,Rn_calc_Ba_600Hz_steps,t_checks); Ln_est_B2a=interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_3p2); Ln_calc_Ba_600Hz_plot=interp1(t_checks,Ln_calc_Ba_600Hz,t_3p2); Rn_est_B2a=interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_3p2); Rn_calc_Ba_600Hz_plot=interp1(t_checks,Rn_calc_Ba_600Hz,t_3p2); Eq_actual_B2a_600Hz = interp1(B2a.Eq_act.Time,B2a.Eq_act.Data,t_checks); Ed_actual_B2a_600Hz = interp1(B2a.Ed_act.Time,B2a.Ed_act.Data,t_checks); TRUE_Ln_calc_Ba600_SEG_A=Lg_2ap1_50Hz; TRUE_Ln_calc_Ba600_SEG_B=Lg_2ap2_50Hz; t_int_3p2_SEG1=13:1.6:21; Ln_B2a_SEG1 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG1); Err_Ln_B2a_SEG1=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B2a_SEG1))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Ln_B2a_SEG2 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG2); Err_Ln_B2a_SEG2=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B2a_SEG2))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Ln_B2a_SEG3 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG3); Err_Ln_B2a_SEG3=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B2a_SEG3))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Ln_B2a_SEG4 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG4); Err_Ln_B2a_SEG4=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B2a_SEG4))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Ln_B2a_SEG5 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG5); Err_Ln_B2a_SEG5=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B2a_SEG5))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Ln_B2a_SEG6 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG6); Err_Ln_B2a_SEG6=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B2a_SEG6))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Ln_B2a_SEG7 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG7); Err_Ln_B2a_SEG7=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B2a_SEG7))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Ln_B2a_SEG8 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG8); Err_Ln_B2a_SEG8=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B2a_SEG8))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Ln_B2a_SEG9 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG9); Err_Ln_B2a_SEG9=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B2a_SEG9))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Ln_B2a_SEG10 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG10); Err_Ln_B2a_SEG10=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B2a_SEG10))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Ln_B2a_SEG11 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG11); Err_Ln_B2a_SEG11=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B2a_SEG11))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Ln_B2a_SEG12 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG12); Err_Ln_B2a_SEG12=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B2a_SEG12))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Ln_B2a_SEG13 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG13); Err_Ln_B2a_SEG13=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B2a_SEG13))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Ln_B2a_SEG14 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG14); Err_Ln_B2a_SEG14=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B2a_SEG14))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Ln_B2a_SEG15 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG15); Err_Ln_B2a_SEG15=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B2a_SEG15))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Ln_B2a_SEG16 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG16); Err_Ln_B2a_SEG16=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B2a_SEG16))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Ln_B2a_SEG17 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG17); Err_Ln_B2a_SEG17=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B2a_SEG17))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Ln_B2a_SEG18 = interp1(B2a.Ln_est.Time,B2a.Ln_est.Data,t_int_3p2_SEG18); Err_Ln_B2a_SEG18=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B2a_SEG18))./TRUE_Ln_calc_Ba600_SEG_B.*100); Error_Ln_B2a_ALL=[Err_Ln_B2a_SEG1 Err_Ln_B2a_SEG2 Err_Ln_B2a_SEG3 Err_Ln_B2a_SEG4 Err_Ln_B2a_SEG5 Err_Ln_B2a_SEG6 Err_Ln_B2a_SEG7 Err_Ln_B2a_SEG8 Err_Ln_B2a_SEG9 Err_Ln_B2a_SEG10 Err_Ln_B2a_SEG11 Err_Ln_B2a_SEG12 Err_Ln_B2a_SEG13 Err_Ln_B2a_SEG14 Err_Ln_B2a_SEG15 Err_Ln_B2a_SEG16 Err_Ln_B2a_SEG17 Err_Ln_B2a_SEG18]; Test_B_SCR_5p8_Lerr=Error_Ln_B2a_ALL(1:9) Test_B_SCR_4p4_Lerr=Error_Ln_B2a_ALL(10:end) % Now R err TRUE_Rn_calc_Ba_SEG_A=Rg_2ap1_50Hz; TRUE_Rn_calc_Ba_SEG_B=Rg_2ap2_50Hz; t_int_3p2_SEG1=13:1.6:21; Rn_B2a_SEG1 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG1); Err_Rn_B2a_SEG1=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B2a_SEG1))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Rn_B2a_SEG2 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG2); Err_Rn_B2a_SEG2=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B2a_SEG2))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Rn_B2a_SEG3 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG3); Err_Rn_B2a_SEG3=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B2a_SEG3))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Rn_B2a_SEG4 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG4); Err_Rn_B2a_SEG4=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B2a_SEG4))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Rn_B2a_SEG5 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG5); Err_Rn_B2a_SEG5=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B2a_SEG5))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Rn_B2a_SEG6 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG6); Err_Rn_B2a_SEG6=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B2a_SEG6))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Rn_B2a_SEG7 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG7); Err_Rn_B2a_SEG7=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B2a_SEG7))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Rn_B2a_SEG8 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG8); Err_Rn_B2a_SEG8=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B2a_SEG8))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Rn_B2a_SEG9 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG9); Err_Rn_B2a_SEG9=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B2a_SEG9))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Rn_B2a_SEG10 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG10); Err_Rn_B2a_SEG10=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B2a_SEG10))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Rn_B2a_SEG11 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG11); Err_Rn_B2a_SEG11=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B2a_SEG11))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Rn_B2a_SEG12 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG12); Err_Rn_B2a_SEG12=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B2a_SEG12))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Rn_B2a_SEG13 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG13); Err_Rn_B2a_SEG13=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B2a_SEG13))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Rn_B2a_SEG14 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG14); Err_Rn_B2a_SEG14=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B2a_SEG14))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Rn_B2a_SEG15 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG15); Err_Rn_B2a_SEG15=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B2a_SEG15))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Rn_B2a_SEG16 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG16); Err_Rn_B2a_SEG16=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B2a_SEG16))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Rn_B2a_SEG17 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG17); Err_Rn_B2a_SEG17=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B2a_SEG17))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Rn_B2a_SEG18 = interp1(B2a.Rn_est.Time,B2a.Rn_est.Data,t_int_3p2_SEG18); Err_Rn_B2a_SEG18=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B2a_SEG18))./TRUE_Rn_calc_Ba_SEG_B.*100); Error_Rn_B2a_ALL=[Err_Rn_B2a_SEG1 Err_Rn_B2a_SEG2 Err_Rn_B2a_SEG3 Err_Rn_B2a_SEG4 Err_Rn_B2a_SEG5 Err_Rn_B2a_SEG6 Err_Rn_B2a_SEG7 Err_Rn_B2a_SEG8 Err_Rn_B2a_SEG9 Err_Rn_B2a_SEG10 Err_Rn_B2a_SEG11 Err_Rn_B2a_SEG12 Err_Rn_B2a_SEG13 Err_Rn_B2a_SEG14 Err_Rn_B2a_SEG15 Err_Rn_B2a_SEG16 Err_Rn_B2a_SEG17 Err_Rn_B2a_SEG18]; Test_B_SCR_5p8_Rerr=Error_Rn_B2a_ALL(1:9) Test_B_SCR_4p4_Rerr=Error_Rn_B2a_ALL(10:end) % Now SCR 2.3 Ln_calc_Bb_600Hz_steps=[Lg_2bp1_50Hz Lg_2bp1_50Hz Lg_2bp2_50Hz Lg_2bp2_50Hz]; Ln_calc_Bb_600Hz=interp1(Timing_imp_step_change,Ln_calc_Bb_600Hz_steps,t_checks); Rn_calc_Bb_600Hz_steps=[Rg_2bp1_50Hz Rg_2bp1_50Hz Rg_2bp2_50Hz Rg_2bp2_50Hz]; Rn_calc_Bb_600Hz=interp1(Timing_imp_step_change,Rn_calc_Bb_600Hz_steps,t_checks); Ln_est_B2b=interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_3p2); Ln_calc_Bb_600Hz_plot=interp1(t_checks,Ln_calc_Bb_600Hz,t_3p2); Rn_est_B2b=interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_3p2); Rn_calc_Bb_600Hz_plot=interp1(t_checks,Rn_calc_Bb_600Hz,t_3p2); TRUE_Ln_calc_Bb600_SEG_A=Lg_2bp1_50Hz; TRUE_Ln_calc_Bb600_SEG_B=Lg_2bp2_50Hz; t_int_3p2_SEG1=13:1.6:21; Ln_B2b_SEG1 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG1); Err_Ln_B2b_SEG1=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B2b_SEG1))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Ln_B2b_SEG2 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG2); Err_Ln_B2b_SEG2=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B2b_SEG2))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Ln_B2b_SEG3 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG3); Err_Ln_B2b_SEG3=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B2b_SEG3))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Ln_B2b_SEG4 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG4); Err_Ln_B2b_SEG4=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B2b_SEG4))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Ln_B2b_SEG5 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG5); Err_Ln_B2b_SEG5=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B2b_SEG5))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Ln_B2b_SEG6 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG6); Err_Ln_B2b_SEG6=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B2b_SEG6))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Ln_B2b_SEG7 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG7); Err_Ln_B2b_SEG7=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B2b_SEG7))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Ln_B2b_SEG8 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG8); Err_Ln_B2b_SEG8=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B2b_SEG8))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Ln_B2b_SEG9 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG9); Err_Ln_B2b_SEG9=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B2b_SEG9))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Ln_B2b_SEG10 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG10); Err_Ln_B2b_SEG10=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B2b_SEG10))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Ln_B2b_SEG11 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG11); Err_Ln_B2b_SEG11=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B2b_SEG11))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Ln_B2b_SEG12 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG12); Err_Ln_B2b_SEG12=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B2b_SEG12))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Ln_B2b_SEG13 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG13); Err_Ln_B2b_SEG13=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B2b_SEG13))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Ln_B2b_SEG14 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG14); Err_Ln_B2b_SEG14=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B2b_SEG14))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Ln_B2b_SEG15 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG15); Err_Ln_B2b_SEG15=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B2b_SEG15))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Ln_B2b_SEG16 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG16); Err_Ln_B2b_SEG16=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B2b_SEG16))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Ln_B2b_SEG17 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG17); Err_Ln_B2b_SEG17=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B2b_SEG17))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Ln_B2b_SEG18 = interp1(B2b.Ln_est.Time,B2b.Ln_est.Data,t_int_3p2_SEG18); Err_Ln_B2b_SEG18=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B2b_SEG18))./TRUE_Ln_calc_Bb600_SEG_B.*100); Error_Ln_B2b_ALL=[Err_Ln_B2b_SEG1 Err_Ln_B2b_SEG2 Err_Ln_B2b_SEG3 Err_Ln_B2b_SEG4 Err_Ln_B2b_SEG5 Err_Ln_B2b_SEG6 Err_Ln_B2b_SEG7 Err_Ln_B2b_SEG8 Err_Ln_B2b_SEG9 Err_Ln_B2b_SEG10 Err_Ln_B2b_SEG11 Err_Ln_B2b_SEG12 Err_Ln_B2b_SEG13 Err_Ln_B2b_SEG14 Err_Ln_B2b_SEG15 Err_Ln_B2b_SEG16 Err_Ln_B2b_SEG17 Err_Ln_B2b_SEG18]; Test_B_SCR_2p3_Lerr=Error_Ln_B2b_ALL(1:9) % Now R error TRUE_Rn_calc_Bb_SEG_A=Rg_2bp1_50Hz; TRUE_Rn_calc_Bb_SEG_B=Rg_2bp2_50Hz; t_int_int_3p2_SEG1=13:1.6:21; Rn_B2b_SEG1 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG1); Err_Rn_B2b_SEG1=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B2b_SEG1))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Rn_B2b_SEG2 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG2); Err_Rn_B2b_SEG2=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B2b_SEG2))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Rn_B2b_SEG3 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG3); Err_Rn_B2b_SEG3=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B2b_SEG3))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Rn_B2b_SEG4 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG4); Err_Rn_B2b_SEG4=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B2b_SEG4))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Rn_B2b_SEG5 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG5); Err_Rn_B2b_SEG5=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B2b_SEG5))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Rn_B2b_SEG6 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG6); Err_Rn_B2b_SEG6=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B2b_SEG6))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Rn_B2b_SEG7 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG7); Err_Rn_B2b_SEG7=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B2b_SEG7))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Rn_B2b_SEG8 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG8); Err_Rn_B2b_SEG8=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B2b_SEG8))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Rn_B2b_SEG9 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG9); Err_Rn_B2b_SEG9=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B2b_SEG9))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Rn_B2b_SEG10 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG10); Err_Rn_B2b_SEG10=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B2b_SEG10))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Rn_B2b_SEG11 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG11); Err_Rn_B2b_SEG11=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B2b_SEG11))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Rn_B2b_SEG12 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG12); Err_Rn_B2b_SEG12=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B2b_SEG12))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Rn_B2b_SEG13 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG13); Err_Rn_B2b_SEG13=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B2b_SEG13))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Rn_B2b_SEG14 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG14); Err_Rn_B2b_SEG14=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B2b_SEG14))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Rn_B2b_SEG15 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG15); Err_Rn_B2b_SEG15=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B2b_SEG15))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Rn_B2b_SEG16 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG16); Err_Rn_B2b_SEG16=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B2b_SEG16))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Rn_B2b_SEG17 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG17); Err_Rn_B2b_SEG17=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B2b_SEG17))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Rn_B2b_SEG18 = interp1(B2b.Rn_est.Time,B2b.Rn_est.Data,t_int_3p2_SEG18); Err_Rn_B2b_SEG18=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B2b_SEG18))./TRUE_Rn_calc_Bb_SEG_B.*100); Error_Rn_B2b_ALL=[Err_Rn_B2b_SEG1 Err_Rn_B2b_SEG2 Err_Rn_B2b_SEG3 Err_Rn_B2b_SEG4 Err_Rn_B2b_SEG5 Err_Rn_B2b_SEG6 Err_Rn_B2b_SEG7 Err_Rn_B2b_SEG8 Err_Rn_B2b_SEG9 Err_Rn_B2b_SEG10 Err_Rn_B2b_SEG11 Err_Rn_B2b_SEG12 Err_Rn_B2b_SEG13 Err_Rn_B2b_SEG14 Err_Rn_B2b_SEG15 Err_Rn_B2b_SEG16 Err_Rn_B2b_SEG17 Err_Rn_B2b_SEG18]; Test_B_SCR_2p3_Rerr=Error_Rn_B2b_ALL(1:9) %% Test C clear % Load data % This batch is for AE+RLCRL+noise+SCR7and5 B1Aa=load('BatchRL1aNOISY.mat'); % 'P_dem','P_act','Q_act','U_dem','U_act','Ln_est','Rn_est','Eq_est','Ed_est','Eq_act','Ed_act' B1Ab=load('BatchRL1bNOISY.mat'); % 'Ln_calc','Rn_calc','Ln_act','Rn_act' % This batch is for AE+RLCRL+noise+SCR3and1.5 ==== index 2 B1Ba=load('BatchRL2aNOISY.mat');% 'P_dem','P_act','Q_act','U_dem','U_act','Ln_est','Rn_est','Eq_est','Ed_est','Eq_act','Ed_act' B1Bb=load('BatchRL2bNOISY.mat'); % 'Ln_calc','Rn_calc','Ln_act','Rn_act' B1Ba_Uact=B1Ba.U_act; half_per=1.6; t_3p2=13:half_per:189; t_checks = 0:0.1:190; Ln_calc_Ba_600Hz_steps=[0.0165 0.0165 0.0231 0.0231]; % calculated in M12b_Resonances Timing_imp_step_change=[0 101 101.0001 189]; Ln_calc_Ba_600Hz=interp1(Timing_imp_step_change,Ln_calc_Ba_600Hz_steps,t_checks); Rn_calc_Ba_600Hz_steps=[0.5173 0.5173 0.7243 0.7243]; Rn_calc_Ba_600Hz=interp1(Timing_imp_step_change,Rn_calc_Ba_600Hz_steps,t_checks); Ln_est_B1a=interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_3p2); Ln_calc_Ba_600Hz_plot=interp1(t_checks,Ln_calc_Ba_600Hz,t_3p2); Rn_est_B1a=interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_3p2); Rn_calc_Ba_600Hz_plot=interp1(t_checks,Rn_calc_Ba_600Hz,t_3p2); % L error SCR 7 and 5 TRUE_Ln_calc_Ba600_SEG_A=0.0165; TRUE_Ln_calc_Ba600_SEG_B=0.0231; t_int_3p2_SEG1=13:1.6:21; Ln_B1Aa_SEG1 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG1); Err_Ln_B1Aa_SEG1=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B1Aa_SEG1))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Ln_B1Aa_SEG2 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG2); Err_Ln_B1Aa_SEG2=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B1Aa_SEG2))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Ln_B1Aa_SEG3 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG3); Err_Ln_B1Aa_SEG3=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B1Aa_SEG3))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Ln_B1Aa_SEG4 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG4); Err_Ln_B1Aa_SEG4=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B1Aa_SEG4))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Ln_B1Aa_SEG5 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG5); Err_Ln_B1Aa_SEG5=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B1Aa_SEG5))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Ln_B1Aa_SEG6 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG6); Err_Ln_B1Aa_SEG6=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B1Aa_SEG6))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Ln_B1Aa_SEG7 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG7); Err_Ln_B1Aa_SEG7=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B1Aa_SEG7))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Ln_B1Aa_SEG8 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG8); Err_Ln_B1Aa_SEG8=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B1Aa_SEG8))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Ln_B1Aa_SEG9 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG9); Err_Ln_B1Aa_SEG9=mean((TRUE_Ln_calc_Ba600_SEG_A-mean(Ln_B1Aa_SEG9))./TRUE_Ln_calc_Ba600_SEG_A.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Ln_B1Aa_SEG10 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG10); Err_Ln_B1Aa_SEG10=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B1Aa_SEG10))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Ln_B1Aa_SEG11 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG11); Err_Ln_B1Aa_SEG11=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B1Aa_SEG11))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Ln_B1Aa_SEG12 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG12); Err_Ln_B1Aa_SEG12=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B1Aa_SEG12))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Ln_B1Aa_SEG13 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG13); Err_Ln_B1Aa_SEG13=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B1Aa_SEG13))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Ln_B1Aa_SEG14 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG14); Err_Ln_B1Aa_SEG14=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B1Aa_SEG14))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Ln_B1Aa_SEG15 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG15); Err_Ln_B1Aa_SEG15=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B1Aa_SEG15))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Ln_B1Aa_SEG16 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG16); Err_Ln_B1Aa_SEG16=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B1Aa_SEG16))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Ln_B1Aa_SEG17 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG17); Err_Ln_B1Aa_SEG17=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B1Aa_SEG17))./TRUE_Ln_calc_Ba600_SEG_B.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Ln_B1Aa_SEG18 = interp1(B1Aa.Ln_est.Time,B1Aa.Ln_est.Data,t_int_3p2_SEG18); Err_Ln_B1Aa_SEG18=mean((TRUE_Ln_calc_Ba600_SEG_B-mean(Ln_B1Aa_SEG18))./TRUE_Ln_calc_Ba600_SEG_B.*100); Error_Ln_B1Aa_ALL=[Err_Ln_B1Aa_SEG1 Err_Ln_B1Aa_SEG2 Err_Ln_B1Aa_SEG3 Err_Ln_B1Aa_SEG4 Err_Ln_B1Aa_SEG5 Err_Ln_B1Aa_SEG6 Err_Ln_B1Aa_SEG7 Err_Ln_B1Aa_SEG8 Err_Ln_B1Aa_SEG9 Err_Ln_B1Aa_SEG10 Err_Ln_B1Aa_SEG11 Err_Ln_B1Aa_SEG12 Err_Ln_B1Aa_SEG13 Err_Ln_B1Aa_SEG14 Err_Ln_B1Aa_SEG15 Err_Ln_B1Aa_SEG16 Err_Ln_B1Aa_SEG17 Err_Ln_B1Aa_SEG18]; Test_C_SCR_7_Lerr=Error_Ln_B1Aa_ALL(1:9) Test_C_SCR_5_Lerr=Error_Ln_B1Aa_ALL(10:end) % R error TRUE_Rn_calc_Ba_SEG_A=0.5173; TRUE_Rn_calc_Ba_SEG_B=0.7243; t_int_3p2_SEG1=13:1.6:21; Rn_B1Aa_SEG1 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG1); Err_Rn_B1Aa_SEG1=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B1Aa_SEG1))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Rn_B1Aa_SEG2 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG2); Err_Rn_B1Aa_SEG2=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B1Aa_SEG2))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Rn_B1Aa_SEG3 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG3); Err_Rn_B1Aa_SEG3=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B1Aa_SEG3))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Rn_B1Aa_SEG4 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG4); Err_Rn_B1Aa_SEG4=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B1Aa_SEG4))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Rn_B1Aa_SEG5 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG5); Err_Rn_B1Aa_SEG5=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B1Aa_SEG5))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Rn_B1Aa_SEG6 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG6); Err_Rn_B1Aa_SEG6=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B1Aa_SEG6))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Rn_B1Aa_SEG7 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG7); Err_Rn_B1Aa_SEG7=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B1Aa_SEG7))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Rn_B1Aa_SEG8 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG8); Err_Rn_B1Aa_SEG8=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B1Aa_SEG8))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Rn_B1Aa_SEG9 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG9); Err_Rn_B1Aa_SEG9=mean((TRUE_Rn_calc_Ba_SEG_A-mean(Rn_B1Aa_SEG9))./TRUE_Rn_calc_Ba_SEG_A.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Rn_B1Aa_SEG10 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG10); Err_Rn_B1Aa_SEG10=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B1Aa_SEG10))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Rn_B1Aa_SEG11 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG11); Err_Rn_B1Aa_SEG11=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B1Aa_SEG11))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Rn_B1Aa_SEG12 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG12); Err_Rn_B1Aa_SEG12=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B1Aa_SEG12))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Rn_B1Aa_SEG13 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG13); Err_Rn_B1Aa_SEG13=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B1Aa_SEG13))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Rn_B1Aa_SEG14 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG14); Err_Rn_B1Aa_SEG14=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B1Aa_SEG14))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Rn_B1Aa_SEG15 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG15); Err_Rn_B1Aa_SEG15=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B1Aa_SEG15))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Rn_B1Aa_SEG16 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG16); Err_Rn_B1Aa_SEG16=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B1Aa_SEG16))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Rn_B1Aa_SEG17 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG17); Err_Rn_B1Aa_SEG17=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B1Aa_SEG17))./TRUE_Rn_calc_Ba_SEG_B.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Rn_B1Aa_SEG18 = interp1(B1Aa.Rn_est.Time,B1Aa.Rn_est.Data,t_int_3p2_SEG18); Err_Rn_B1Aa_SEG18=mean((TRUE_Rn_calc_Ba_SEG_B-mean(Rn_B1Aa_SEG18))./TRUE_Rn_calc_Ba_SEG_B.*100); Error_Rn_B1Aa_ALL=[Err_Rn_B1Aa_SEG1 Err_Rn_B1Aa_SEG2 Err_Rn_B1Aa_SEG3 Err_Rn_B1Aa_SEG4 Err_Rn_B1Aa_SEG5 Err_Rn_B1Aa_SEG6 Err_Rn_B1Aa_SEG7 Err_Rn_B1Aa_SEG8 Err_Rn_B1Aa_SEG9 Err_Rn_B1Aa_SEG10 Err_Rn_B1Aa_SEG11 Err_Rn_B1Aa_SEG12 Err_Rn_B1Aa_SEG13 Err_Rn_B1Aa_SEG14 Err_Rn_B1Aa_SEG15 Err_Rn_B1Aa_SEG16 Err_Rn_B1Aa_SEG17 Err_Rn_B1Aa_SEG18]; Test_C_SCR_7_Rerr=Error_Rn_B1Aa_ALL(1:9) Test_C_SCR_5_Rerr=Error_Rn_B1Aa_ALL(10:end) % SCR 3 and 1.5 Ln_calc_Bb_600Hz_steps=[0.0384 0.0384 0.0768 0.0768]; Ln_calc_Bb_600Hz=interp1(Timing_imp_step_change,Ln_calc_Bb_600Hz_steps,t_checks); Rn_calc_Bb_600Hz_steps=[1.2071 1.2071 2.4143 2.4143]; Rn_calc_Bb_600Hz=interp1(Timing_imp_step_change,Rn_calc_Bb_600Hz_steps,t_checks); Ln_est_B1Ba=interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_3p2); Ln_calc_Bb_600Hz_plot=interp1(t_checks,Ln_calc_Bb_600Hz,t_3p2); Rn_est_B1Ba=interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_3p2); Rn_calc_Bb_600Hz_plot=interp1(t_checks,Rn_calc_Bb_600Hz,t_3p2); % L error TRUE_Ln_calc_Bb600_SEG_A=0.0384; TRUE_Ln_calc_Bb600_SEG_B=0.0768; t_int_3p2_SEG1=13:1.6:21; Ln_B1Ba_SEG1 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG1); Err_Ln_B1Ba_SEG1=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B1Ba_SEG1))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Ln_B1Ba_SEG2 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG2); Err_Ln_B1Ba_SEG2=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B1Ba_SEG2))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Ln_B1Ba_SEG3 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG3); Err_Ln_B1Ba_SEG3=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B1Ba_SEG3))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Ln_B1Ba_SEG4 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG4); Err_Ln_B1Ba_SEG4=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B1Ba_SEG4))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Ln_B1Ba_SEG5 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG5); Err_Ln_B1Ba_SEG5=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B1Ba_SEG5))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Ln_B1Ba_SEG6 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG6); Err_Ln_B1Ba_SEG6=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B1Ba_SEG6))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Ln_B1Ba_SEG7 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG7); Err_Ln_B1Ba_SEG7=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B1Ba_SEG7))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Ln_B1Ba_SEG8 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG8); Err_Ln_B1Ba_SEG8=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B1Ba_SEG8))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Ln_B1Ba_SEG9 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG9); Err_Ln_B1Ba_SEG9=mean((TRUE_Ln_calc_Bb600_SEG_A-mean(Ln_B1Ba_SEG9))./TRUE_Ln_calc_Bb600_SEG_A.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Ln_B1Ba_SEG10 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG10); Err_Ln_B1Ba_SEG10=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B1Ba_SEG10))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Ln_B1Ba_SEG11 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG11); Err_Ln_B1Ba_SEG11=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B1Ba_SEG11))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Ln_B1Ba_SEG12 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG12); Err_Ln_B1Ba_SEG12=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B1Ba_SEG12))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Ln_B1Ba_SEG13 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG13); Err_Ln_B1Ba_SEG13=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B1Ba_SEG13))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Ln_B1Ba_SEG14 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG14); Err_Ln_B1Ba_SEG14=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B1Ba_SEG14))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Ln_B1Ba_SEG15 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG15); Err_Ln_B1Ba_SEG15=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B1Ba_SEG15))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Ln_B1Ba_SEG16 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG16); Err_Ln_B1Ba_SEG16=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B1Ba_SEG16))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Ln_B1Ba_SEG17 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG17); Err_Ln_B1Ba_SEG17=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B1Ba_SEG17))./TRUE_Ln_calc_Bb600_SEG_B.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Ln_B1Ba_SEG18 = interp1(B1Ba.Ln_est.Time,B1Ba.Ln_est.Data,t_int_3p2_SEG18); Err_Ln_B1Ba_SEG18=mean((TRUE_Ln_calc_Bb600_SEG_B-mean(Ln_B1Ba_SEG18))./TRUE_Ln_calc_Bb600_SEG_B.*100); Error_Ln_B1Ba_ALL=[Err_Ln_B1Ba_SEG1 Err_Ln_B1Ba_SEG2 Err_Ln_B1Ba_SEG3 Err_Ln_B1Ba_SEG4 Err_Ln_B1Ba_SEG5 Err_Ln_B1Ba_SEG6 Err_Ln_B1Ba_SEG7 Err_Ln_B1Ba_SEG8 Err_Ln_B1Ba_SEG9 Err_Ln_B1Ba_SEG10 Err_Ln_B1Ba_SEG11 Err_Ln_B1Ba_SEG12 Err_Ln_B1Ba_SEG13 Err_Ln_B1Ba_SEG14 Err_Ln_B1Ba_SEG15 Err_Ln_B1Ba_SEG16 Err_Ln_B1Ba_SEG17 Err_Ln_B1Ba_SEG18]; Test_C_SCR_3_Lerr=Error_Ln_B1Ba_ALL(1:9) Test_C_SCR_1p5_Lerr=Error_Ln_B1Ba_ALL(10:end) % R error TRUE_Rn_calc_Bb_SEG_A=1.206; TRUE_Rn_calc_Bb_SEG_B=2.413; t_int_int_3p2_SEG1=13:1.6:21; Rn_B1Ba_SEG1 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG1); Err_Rn_B1Ba_SEG1=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B1Ba_SEG1))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG2=(13:1.6:21)+10; Rn_B1Ba_SEG2 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG2); Err_Rn_B1Ba_SEG2=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B1Ba_SEG2))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG3=(13:1.6:21)+20; Rn_B1Ba_SEG3 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG3); Err_Rn_B1Ba_SEG3=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B1Ba_SEG3))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG4=(13:1.6:21)+30; Rn_B1Ba_SEG4 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG4); Err_Rn_B1Ba_SEG4=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B1Ba_SEG4))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG5=(13:1.6:21)+40; Rn_B1Ba_SEG5 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG5); Err_Rn_B1Ba_SEG5=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B1Ba_SEG5))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG6=(13:1.6:21)+50; Rn_B1Ba_SEG6 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG6); Err_Rn_B1Ba_SEG6=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B1Ba_SEG6))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG7=(13:1.6:21)+60; Rn_B1Ba_SEG7 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG7); Err_Rn_B1Ba_SEG7=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B1Ba_SEG7))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG8=(13:1.6:21)+70; Rn_B1Ba_SEG8 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG8); Err_Rn_B1Ba_SEG8=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B1Ba_SEG8))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG9=(13:1.6:21)+80; Rn_B1Ba_SEG9 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG9); Err_Rn_B1Ba_SEG9=mean((TRUE_Rn_calc_Bb_SEG_A-mean(Rn_B1Ba_SEG9))./TRUE_Rn_calc_Bb_SEG_A.*100); t_int_3p2_SEG10=(13:1.6:21)+90; Rn_B1Ba_SEG10 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG10); Err_Rn_B1Ba_SEG10=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B1Ba_SEG10))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG11=(13:1.6:21)+100; Rn_B1Ba_SEG11 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG11); Err_Rn_B1Ba_SEG11=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B1Ba_SEG11))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG12=(13:1.6:21)+110; Rn_B1Ba_SEG12 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG12); Err_Rn_B1Ba_SEG12=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B1Ba_SEG12))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG13=(13:1.6:21)+120; Rn_B1Ba_SEG13 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG13); Err_Rn_B1Ba_SEG13=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B1Ba_SEG13))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG14=(13:1.6:21)+130; Rn_B1Ba_SEG14 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG14); Err_Rn_B1Ba_SEG14=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B1Ba_SEG14))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG15=(13:1.6:21)+140; Rn_B1Ba_SEG15 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG15); Err_Rn_B1Ba_SEG15=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B1Ba_SEG15))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG16=(13:1.6:21)+150; Rn_B1Ba_SEG16 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG16); Err_Rn_B1Ba_SEG16=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B1Ba_SEG16))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG17=(13:1.6:21)+160; Rn_B1Ba_SEG17 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG17); Err_Rn_B1Ba_SEG17=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B1Ba_SEG17))./TRUE_Rn_calc_Bb_SEG_B.*100); t_int_3p2_SEG18=(13:1.6:20)+170; Rn_B1Ba_SEG18 = interp1(B1Ba.Rn_est.Time,B1Ba.Rn_est.Data,t_int_3p2_SEG18); Err_Rn_B1Ba_SEG18=mean((TRUE_Rn_calc_Bb_SEG_B-mean(Rn_B1Ba_SEG18))./TRUE_Rn_calc_Bb_SEG_B.*100); Error_Rn_B1Ba_ALL=[Err_Rn_B1Ba_SEG1 Err_Rn_B1Ba_SEG2 Err_Rn_B1Ba_SEG3 Err_Rn_B1Ba_SEG4 Err_Rn_B1Ba_SEG5 Err_Rn_B1Ba_SEG6 Err_Rn_B1Ba_SEG7 Err_Rn_B1Ba_SEG8 Err_Rn_B1Ba_SEG9 Err_Rn_B1Ba_SEG10 Err_Rn_B1Ba_SEG11 Err_Rn_B1Ba_SEG12 Err_Rn_B1Ba_SEG13 Err_Rn_B1Ba_SEG14 Err_Rn_B1Ba_SEG15 Err_Rn_B1Ba_SEG16 Err_Rn_B1Ba_SEG17 Err_Rn_B1Ba_SEG18]; Test_C_SCR_3_Rerr=Error_Rn_B1Ba_ALL(1:9) Test_C_SCR_1p5_Rerr=Error_Rn_B1Ba_ALL(10:end) %% Recursive Least Squares - Test D,E,F % Test D clear % Load % This batch is for RLS+RL+no noise+SCR7and5 B1Aaa=load('Batch1a_Rec_noisy_a_v2.mat'); % 'Batch1a_Rec_noisy_a','P_dem','P_act','Q_act','U_dem','U_act' B1Aab=load('Batch1a_Rec_noisy_b_v2.mat'); % 'Ln_act','Ln_est','Rn_act' B1Aac=load('Batch1a_Rec_noisy_c_v2.mat'); % 'Rn_est' B1Aad=load('Batch1a_Rec_noisy_d_v2.mat'); % 'Eq_est','Ed_est','Eq_act','Ed_act' U_dist2=((B1Aaa.U_dem-1.592168332809066e+05)/1.592168332809066e+05); period=0.9; simulation=9:period:128.7; t_checks = 0:0.1:130; Ln_REAL_steps=[0.0165 0.0165 0.0231 0.0231]; % insert real values Timing_imp_step_change=[0 70 70.0001 130]; Impedance_calc=interp1(Timing_imp_step_change,Ln_REAL_steps,t_checks); RN_REAL_steps=[0.5173 0.5173 0.7243 0.7243]; % insert real values Resistance_calc=interp1(Timing_imp_step_change,RN_REAL_steps,t_checks); Ln_est_B1a=interp1(B1Aab.Ln_est.Time,B1Aab.Ln_est.Data(:,(1:1:13000)),simulation);% insert real values Ln_REAL_plot=interp1(t_checks,Impedance_calc,simulation); Rn_est_B1a=interp1(B1Aac.Rn_est.Time,B1Aac.Rn_est.Data(:,(1:1:13000)),simulation);% insert real values Rn_REAL_plot=interp1(t_checks,Resistance_calc,simulation); Ln_REAL_SEG_A=0.0165; Ln_REAL_SEG_B=0.0231; period=0.9; simulation=9:period:128.7; % L error period_SEG1=18:0.9:30; Ln_B1Aa_SEG1 = interp1(B1Aab.Ln_est.Time,B1Aab.Ln_est.Data(:,(1:1:13000)),period_SEG1); Err_Ln_B1Aa_SEG1=mean((Ln_REAL_SEG_A-mean(Ln_B1Aa_SEG1))./Ln_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Ln_B1Aa_SEG2 = interp1(B1Aab.Ln_est.Time,B1Aab.Ln_est.Data(:,(1:1:13000)),period_SEG2); Err_Ln_B1Aa_SEG2=mean((Ln_REAL_SEG_A-mean(Ln_B1Aa_SEG2))./Ln_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Ln_B1Aa_SEG3 = interp1(B1Aab.Ln_est.Time,B1Aab.Ln_est.Data(:,(1:1:13000)),period_SEG3); Err_Ln_B1Aa_SEG3=mean((Ln_REAL_SEG_A-mean(Ln_B1Aa_SEG3))./Ln_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Ln_B1Aa_SEG4 = interp1(B1Aab.Ln_est.Time,B1Aab.Ln_est.Data(:,(1:1:13000)),period_SEG4); Err_Ln_B1Aa_SEG4=mean((Ln_REAL_SEG_B-mean(Ln_B1Aa_SEG4))./Ln_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Ln_B1Aa_SEG5 = interp1(B1Aab.Ln_est.Time,B1Aab.Ln_est.Data(:,(1:1:13000)),period_SEG5); Err_Ln_B1Aa_SEG5=mean((Ln_REAL_SEG_B-mean(Ln_B1Aa_SEG5))./Ln_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Ln_B1Aa_SEG6 = interp1(B1Aab.Ln_est.Time,B1Aab.Ln_est.Data(:,(1:1:13000)),period_SEG6); Err_Ln_B1Aa_SEG6=mean((Ln_REAL_SEG_B-mean(Ln_B1Aa_SEG6))./Ln_REAL_SEG_B.*100); Error_Ln_B1Aa_ALL=[Err_Ln_B1Aa_SEG1 Err_Ln_B1Aa_SEG2 Err_Ln_B1Aa_SEG3 Err_Ln_B1Aa_SEG4 Err_Ln_B1Aa_SEG5 Err_Ln_B1Aa_SEG6]; Test_D_SCR_7_Lerr=Error_Ln_B1Aa_ALL(1:3) Test_D_SCR_5_Lerr=Error_Ln_B1Aa_ALL(4:end) % R error Rn_REAL_SEG_A=0.5173; Rn_REAL_SEG_B=0.7243; period=0.9; simulation=9:period:128.7; period_SEG1=18:0.9:30; Rn_B1Aa_SEG1 = interp1(B1Aac.Rn_est.Time,B1Aac.Rn_est.Data(:,(1:1:13000)),period_SEG1); Err_Rn_B1Aa_SEG1=mean((Rn_REAL_SEG_A-mean(Rn_B1Aa_SEG1))./Rn_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Rn_B1Aa_SEG2 = interp1(B1Aac.Rn_est.Time,B1Aac.Rn_est.Data(:,(1:1:13000)),period_SEG2); Err_Rn_B1Aa_SEG2=mean((Rn_REAL_SEG_A-mean(Rn_B1Aa_SEG2))./Rn_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Rn_B1Aa_SEG3 = interp1(B1Aac.Rn_est.Time,B1Aac.Rn_est.Data(:,(1:1:13000)),period_SEG3); Err_Rn_B1Aa_SEG3=mean((Rn_REAL_SEG_A-mean(Rn_B1Aa_SEG3))./Rn_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Rn_B1Aa_SEG4 = interp1(B1Aac.Rn_est.Time,B1Aac.Rn_est.Data(:,(1:1:13000)),period_SEG4); Err_Rn_B1Aa_SEG4=mean((Rn_REAL_SEG_B-mean(Rn_B1Aa_SEG4))./Rn_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Rn_B1Aa_SEG5 = interp1(B1Aac.Rn_est.Time,B1Aac.Rn_est.Data(:,(1:1:13000)),period_SEG5); Err_Rn_B1Aa_SEG5=mean((Rn_REAL_SEG_B-mean(Rn_B1Aa_SEG5))./Rn_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Rn_B1Aa_SEG6 = interp1(B1Aac.Rn_est.Time,B1Aac.Rn_est.Data(:,(1:1:13000)),period_SEG6); Err_Rn_B1Aa_SEG6=mean((Rn_REAL_SEG_B-mean(Rn_B1Aa_SEG6))./Rn_REAL_SEG_B.*100); Error_Rn_B1Aa_ALL=[Err_Rn_B1Aa_SEG1 Err_Rn_B1Aa_SEG2 Err_Rn_B1Aa_SEG3 Err_Rn_B1Aa_SEG4 Err_Rn_B1Aa_SEG5 Err_Rn_B1Aa_SEG6]; Test_D_SCR_7_Rerr=Error_Rn_B1Aa_ALL(1:3) Test_D_SCR_5_Rerr=Error_Rn_B1Aa_ALL(4:end) % now SCR 3 and 1.5 B1Baa=load('Batch1b_Rec_noisy_a_v2.mat'); % 'P_dem','P_act','Q_act','U_dem','U_act' B1Bab=load('Batch1b_Rec_noisy_b_v2.mat'); % 'Ln_act','Ln_est','Rn_act' B1Bac=load('Batch1b_Rec_noisy_c_v2.mat'); % 'Rn_est' B1Bad=load('Batch1b_Rec_noisy_d_v2.mat'); % 'Eq_est','Ed_est','Eq_act','Ed_act' Ln_REAL_steps=[0.0384 0.0384 0.0768 0.0768]; % insert real values Timing_imp_step_change=[0 70 70.0001 130]; Impedance_calc=interp1(Timing_imp_step_change,Ln_REAL_steps,t_checks); RN_REAL_steps=[1.2071 1.2071 2.4143 2.4143]; % insert real values Resistance_calc=interp1(Timing_imp_step_change,RN_REAL_steps,t_checks); Ln_est_B1B=interp1(B1Bab.Ln_est.Time,B1Bab.Ln_est.Data(:,(1:1:13000)),simulation);% insert real values Ln_REAL_plot=interp1(t_checks,Impedance_calc,simulation); Rn_est_B1B=interp1(B1Bac.Rn_est.Time,B1Bac.Rn_est.Data(:,(1:1:13000)),simulation);% insert real values Rn_REAL_plot=interp1(t_checks,Resistance_calc,simulation); Ln_REAL_SEG_A=0.0384; Ln_REAL_SEG_B=0.0768; period=0.9; simulation=9:period:128.7; % L error period_SEG1=18:0.9:30; Ln_B1Ba_SEG1 = interp1(B1Bab.Ln_est.Time,B1Bab.Ln_est.Data(:,(1:1:13000)),period_SEG1); Err_Ln_B1Ba_SEG1=mean((Ln_REAL_SEG_A-mean(Ln_B1Ba_SEG1))./Ln_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Ln_B1Ba_SEG2 = interp1(B1Bab.Ln_est.Time,B1Bab.Ln_est.Data(:,(1:1:13000)),period_SEG2); Err_Ln_B1Ba_SEG2=mean((Ln_REAL_SEG_A-mean(Ln_B1Ba_SEG2))./Ln_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Ln_B1Ba_SEG3 = interp1(B1Bab.Ln_est.Time,B1Bab.Ln_est.Data(:,(1:1:13000)),period_SEG3); Err_Ln_B1Ba_SEG3=mean((Ln_REAL_SEG_A-mean(Ln_B1Ba_SEG3))./Ln_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Ln_B1Ba_SEG4 = interp1(B1Bab.Ln_est.Time,B1Bab.Ln_est.Data(:,(1:1:13000)),period_SEG4); Err_Ln_B1Ba_SEG4=mean((Ln_REAL_SEG_B-mean(Ln_B1Ba_SEG4))./Ln_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Ln_B1Ba_SEG5 = interp1(B1Bab.Ln_est.Time,B1Bab.Ln_est.Data(:,(1:1:13000)),period_SEG5); Err_Ln_B1Ba_SEG5=mean((Ln_REAL_SEG_B-mean(Ln_B1Ba_SEG5))./Ln_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Ln_B1Ba_SEG6 = interp1(B1Bab.Ln_est.Time,B1Bab.Ln_est.Data(:,(1:1:13000)),period_SEG6); Err_Ln_B1Ba_SEG6=mean((Ln_REAL_SEG_B-mean(Ln_B1Ba_SEG6))./Ln_REAL_SEG_B.*100); Error_Ln_B1Ba_ALL=[Err_Ln_B1Ba_SEG1 Err_Ln_B1Ba_SEG2 Err_Ln_B1Ba_SEG3 Err_Ln_B1Ba_SEG4 Err_Ln_B1Ba_SEG5 Err_Ln_B1Ba_SEG6]; Test_D_SCR_3_Lerr=Error_Ln_B1Ba_ALL(1:3) Test_D_SCR_1p5_Lerr=Error_Ln_B1Ba_ALL(4:end) % Now R error Rn_REAL_SEG_A=1.2071; Rn_REAL_SEG_B=2.4143; period=0.9; simulation=9:period:128.7; period_SEG1=18:0.9:30; Rn_B1Ba_SEG1 = interp1(B1Bac.Rn_est.Time,B1Bac.Rn_est.Data(:,(1:1:13000)),period_SEG1); Err_Rn_B1Ba_SEG1=mean((Rn_REAL_SEG_A-mean(Rn_B1Ba_SEG1))./Rn_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Rn_B1Ba_SEG2 = interp1(B1Bac.Rn_est.Time,B1Bac.Rn_est.Data(:,(1:1:13000)),period_SEG2); Err_Rn_B1Ba_SEG2=mean((Rn_REAL_SEG_A-mean(Rn_B1Ba_SEG2))./Rn_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Rn_B1Ba_SEG3 = interp1(B1Bac.Rn_est.Time,B1Bac.Rn_est.Data(:,(1:1:13000)),period_SEG3); Err_Rn_B1Ba_SEG3=mean((Rn_REAL_SEG_A-mean(Rn_B1Ba_SEG3))./Rn_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Rn_B1Ba_SEG4 = interp1(B1Bac.Rn_est.Time,B1Bac.Rn_est.Data(:,(1:1:13000)),period_SEG4); Err_Rn_B1Ba_SEG4=mean((Rn_REAL_SEG_B-mean(Rn_B1Ba_SEG4))./Rn_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Rn_B1Ba_SEG5 = interp1(B1Bac.Rn_est.Time,B1Bac.Rn_est.Data(:,(1:1:13000)),period_SEG5); Err_Rn_B1Ba_SEG5=mean((Rn_REAL_SEG_B-mean(Rn_B1Ba_SEG5))./Rn_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Rn_B1Ba_SEG6 = interp1(B1Bac.Rn_est.Time,B1Bac.Rn_est.Data(:,(1:1:13000)),period_SEG6); Err_Rn_B1Ba_SEG6=mean((Rn_REAL_SEG_B-mean(Rn_B1Ba_SEG6))./Rn_REAL_SEG_B.*100); Error_Rn_B1Ba_ALL=[Err_Rn_B1Ba_SEG1 Err_Rn_B1Ba_SEG2 Err_Rn_B1Ba_SEG3 Err_Rn_B1Ba_SEG4 Err_Rn_B1Ba_SEG5 Err_Rn_B1Ba_SEG6]; Test_D_SCR_3_Rerr=Error_Rn_B1Ba_ALL(1:3) Test_D_SCR_1p5_Rerr=Error_Rn_B1Ba_ALL(4:end) %% Test E % This batch is for RLS+RLCRL+no noise+SCR7and5 ==== index 5 B2Aaa=load('Batch2a_Rec_noisy_a_v2.mat'); B2Aab=load('Batch2a_Rec_noisy_b_v2.mat'); B2Aac=load('Batch2a_Rec_noisy_c_v2.mat'); B2Aad=load('Batch2a_Rec_noisy_d_v2.mat'); run('M12b_Resonances_in_RLCRL'); close all period=0.9; simulation=9:period:128.7; t_checks = 0:0.1:130; Ln_REAL_steps=[Lg_2ap1_50Hz Lg_2ap1_50Hz Lg_2ap2_50Hz Lg_2ap2_50Hz]; % insert real values Timing_imp_step_change=[0 70 70.0001 130]; Impedance_calc=interp1(Timing_imp_step_change,Ln_REAL_steps,t_checks); RN_REAL_steps=[Rg_2ap1_50Hz Rg_2ap1_50Hz Rg_2ap2_50Hz Rg_2ap2_50Hz]; % insert real values Resistance_calc=interp1(Timing_imp_step_change,RN_REAL_steps,t_checks); Ln_est_B2A=interp1(B2Aab.Ln_est.Time,B2Aab.Ln_est.Data(:,(1:1:13000)),simulation);% insert real values Ln_REAL_plot=interp1(t_checks,Impedance_calc,simulation); Rn_est_B2A=interp1(B2Aac.Rn_est.Time,B2Aac.Rn_est.Data(:,(1:1:13000)),simulation);% insert real values Rn_REAL_plot=interp1(t_checks,Resistance_calc,simulation); Ln_REAL_SEG_A=Lg_2ap1_50Hz; Ln_REAL_SEG_B=Lg_2ap2_50Hz; period=0.9; simulation=9:period:128.7; % L error period_SEG1=18:0.9:30; Ln_B2Aa_SEG1 = interp1(B2Aab.Ln_est.Time,B2Aab.Ln_est.Data(:,(1:1:13000)),period_SEG1); Err_Ln_B2Aa_SEG1=mean((Ln_REAL_SEG_A-mean(Ln_B2Aa_SEG1))./Ln_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Ln_B2Aa_SEG2 = interp1(B2Aab.Ln_est.Time,B2Aab.Ln_est.Data(:,(1:1:13000)),period_SEG2); Err_Ln_B2Aa_SEG2=mean((Ln_REAL_SEG_A-mean(Ln_B2Aa_SEG2))./Ln_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Ln_B2Aa_SEG3 = interp1(B2Aab.Ln_est.Time,B2Aab.Ln_est.Data(:,(1:1:13000)),period_SEG3); Err_Ln_B2Aa_SEG3=mean((Ln_REAL_SEG_A-mean(Ln_B2Aa_SEG3))./Ln_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Ln_B2Aa_SEG4 = interp1(B2Aab.Ln_est.Time,B2Aab.Ln_est.Data(:,(1:1:13000)),period_SEG4); Err_Ln_B2Aa_SEG4=mean((Ln_REAL_SEG_B-mean(Ln_B2Aa_SEG4))./Ln_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Ln_B2Aa_SEG5 = interp1(B2Aab.Ln_est.Time,B2Aab.Ln_est.Data(:,(1:1:13000)),period_SEG5); Err_Ln_B2Aa_SEG5=mean((Ln_REAL_SEG_B-mean(Ln_B2Aa_SEG5))./Ln_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Ln_B2Aa_SEG6 = interp1(B2Aab.Ln_est.Time,B2Aab.Ln_est.Data(:,(1:1:13000)),period_SEG6); Err_Ln_B2Aa_SEG6=mean((Ln_REAL_SEG_B-mean(Ln_B2Aa_SEG6))./Ln_REAL_SEG_B.*100); Error_Ln_B2Aa_ALL=[Err_Ln_B2Aa_SEG1 Err_Ln_B2Aa_SEG2 Err_Ln_B2Aa_SEG3 Err_Ln_B2Aa_SEG4 Err_Ln_B2Aa_SEG5 Err_Ln_B2Aa_SEG6]; Test_E_SCR_5p8_Lerr=Error_Ln_B2Aa_ALL(1:3) Test_E_SCR_4p4_Lerr=Error_Ln_B2Aa_ALL(4:end) % R error Rn_REAL_SEG_A=Rg_2ap1_50Hz; Rn_REAL_SEG_B=Rg_2ap2_50Hz; period=0.9; simulation=9:period:128.7; period_SEG1=18:0.9:30; Rn_B2Aa_SEG1 = interp1(B2Aac.Rn_est.Time,B2Aac.Rn_est.Data(:,(1:1:13000)),period_SEG1); Err_Rn_B2Aa_SEG1=mean((Rn_REAL_SEG_A-mean(Rn_B2Aa_SEG1))./Rn_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Rn_B2Aa_SEG2 = interp1(B2Aac.Rn_est.Time,B2Aac.Rn_est.Data(:,(1:1:13000)),period_SEG2); Err_Rn_B2Aa_SEG2=mean((Rn_REAL_SEG_A-mean(Rn_B2Aa_SEG2))./Rn_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Rn_B2Aa_SEG3 = interp1(B2Aac.Rn_est.Time,B2Aac.Rn_est.Data(:,(1:1:13000)),period_SEG3); Err_Rn_B2Aa_SEG3=mean((Rn_REAL_SEG_A-mean(Rn_B2Aa_SEG3))./Rn_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Rn_B2Aa_SEG4 = interp1(B2Aac.Rn_est.Time,B2Aac.Rn_est.Data(:,(1:1:13000)),period_SEG4); Err_Rn_B2Aa_SEG4=mean((Rn_REAL_SEG_B-mean(Rn_B2Aa_SEG4))./Rn_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Rn_B2Aa_SEG5 = interp1(B2Aac.Rn_est.Time,B2Aac.Rn_est.Data(:,(1:1:13000)),period_SEG5); Err_Rn_B2Aa_SEG5=mean((Rn_REAL_SEG_B-mean(Rn_B2Aa_SEG5))./Rn_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Rn_B2Aa_SEG6 = interp1(B2Aac.Rn_est.Time,B2Aac.Rn_est.Data(:,(1:1:13000)),period_SEG6); Err_Rn_B2Aa_SEG6=mean((Rn_REAL_SEG_B-mean(Rn_B2Aa_SEG6))./Rn_REAL_SEG_B.*100); Error_Rn_B2Aa_ALL=[Err_Rn_B2Aa_SEG1 Err_Rn_B2Aa_SEG2 Err_Rn_B2Aa_SEG3 Err_Rn_B2Aa_SEG4 Err_Rn_B2Aa_SEG5 Err_Rn_B2Aa_SEG6]; Test_E_SCR_5p8_Rerr=Error_Rn_B2Aa_ALL(1:3) Test_E_SCR_4p4_Rerr=Error_Rn_B2Aa_ALL(4:end) % This batch is for RLS+RLCRL+no noise+SCR5and1.5 ==== index 6 B2Baa=load('Batch2b_Rec_noisy_a_v2.mat'); B2Bab=load('Batch2b_Rec_noisy_b_v2.mat'); B2Bac=load('Batch2b_Rec_noisy_c_v2.mat'); B2Bad=load('Batch2b_Rec_noisy_d_v2.mat'); Ln_REAL_steps=[Lg_2bp1_50Hz Lg_2bp1_50Hz Lg_2bp1_50Hz Lg_2bp1_50Hz]; % insert real values Timing_imp_step_change=[0 70 70.0001 130]; Impedance_calc=interp1(Timing_imp_step_change,Ln_REAL_steps,t_checks); RN_REAL_steps=[Rg_2bp1_50Hz Rg_2bp1_50Hz Rg_2bp1_50Hz Rg_2bp1_50Hz ]; % insert real values Resistance_calc=interp1(Timing_imp_step_change,RN_REAL_steps,t_checks); Ln_est_B2B=interp1(B2Bab.Ln_est.Time,B2Bab.Ln_est.Data(:,(1:1:13000)),simulation);% insert real values Ln_REAL_plot=interp1(t_checks,Impedance_calc,simulation); Rn_est_B2B=interp1(B2Bac.Rn_est.Time,B2Bac.Rn_est.Data(:,(1:1:13000)),simulation);% insert real values Rn_REAL_plot=interp1(t_checks,Resistance_calc,simulation); Ln_REAL_SEG_A=Lg_2bp1_50Hz; Ln_REAL_SEG_B=Lg_2bp1_50Hz; period=0.9; simulation=9:period:128.7; % L error period_SEG1=18:0.9:30; Ln_B2Ba_SEG1 = interp1(B2Bab.Ln_est.Time,B2Bab.Ln_est.Data(:,(1:1:13000)),period_SEG1); Err_Ln_B2Ba_SEG1=mean((Ln_REAL_SEG_A-mean(Ln_B2Ba_SEG1))./Ln_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Ln_B2Ba_SEG2 = interp1(B2Bab.Ln_est.Time,B2Bab.Ln_est.Data(:,(1:1:13000)),period_SEG2); Err_Ln_B2Ba_SEG2=mean((Ln_REAL_SEG_A-mean(Ln_B2Ba_SEG2))./Ln_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Ln_B2Ba_SEG3 = interp1(B2Bab.Ln_est.Time,B2Bab.Ln_est.Data(:,(1:1:13000)),period_SEG3); Err_Ln_B2Ba_SEG3=mean((Ln_REAL_SEG_A-mean(Ln_B2Ba_SEG3))./Ln_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Ln_B2Ba_SEG4 = interp1(B2Bab.Ln_est.Time,B2Bab.Ln_est.Data(:,(1:1:13000)),period_SEG4); Err_Ln_B2Ba_SEG4=mean((Ln_REAL_SEG_B-mean(Ln_B2Ba_SEG4))./Ln_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Ln_B2Ba_SEG5 = interp1(B2Bab.Ln_est.Time,B2Bab.Ln_est.Data(:,(1:1:13000)),period_SEG5); Err_Ln_B2Ba_SEG5=mean((Ln_REAL_SEG_B-mean(Ln_B2Ba_SEG5))./Ln_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Ln_B2Ba_SEG6 = interp1(B2Bab.Ln_est.Time,B2Bab.Ln_est.Data(:,(1:1:13000)),period_SEG6); Err_Ln_B2Ba_SEG6=mean((Ln_REAL_SEG_B-mean(Ln_B2Ba_SEG6))./Ln_REAL_SEG_B.*100); Error_Ln_B2Ba_ALL=[Err_Ln_B2Ba_SEG1 Err_Ln_B2Ba_SEG2 Err_Ln_B2Ba_SEG3 Err_Ln_B2Ba_SEG4 Err_Ln_B2Ba_SEG5 Err_Ln_B2Ba_SEG6]; Test_E_SCR_2p3_Lerr=Error_Ln_B2Ba_ALL(1:3) % R error Rn_REAL_SEG_A=Rg_2bp1_50Hz; Rn_REAL_SEG_B=Rg_2bp1_50Hz; period=0.9; simulation=9:period:128.7; period_SEG1=18:0.9:30; Rn_B2Ba_SEG1 = interp1(B2Bac.Rn_est.Time,B2Bac.Rn_est.Data(:,(1:1:13000)),period_SEG1); Err_Rn_B2Ba_SEG1=mean((Rn_REAL_SEG_A-mean(Rn_B2Ba_SEG1))./Rn_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Rn_B2Ba_SEG2 = interp1(B2Bac.Rn_est.Time,B2Bac.Rn_est.Data(:,(1:1:13000)),period_SEG2); Err_Rn_B2Ba_SEG2=mean((Rn_REAL_SEG_A-mean(Rn_B2Ba_SEG2))./Rn_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Rn_B2Ba_SEG3 = interp1(B2Bac.Rn_est.Time,B2Bac.Rn_est.Data(:,(1:1:13000)),period_SEG3); Err_Rn_B2Ba_SEG3=mean((Rn_REAL_SEG_A-mean(Rn_B2Ba_SEG3))./Rn_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Rn_B2Ba_SEG4 = interp1(B2Bac.Rn_est.Time,B2Bac.Rn_est.Data(:,(1:1:13000)),period_SEG4); Err_Rn_B2Ba_SEG4=mean((Rn_REAL_SEG_B-mean(Rn_B2Ba_SEG4))./Rn_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Rn_B2Ba_SEG5 = interp1(B2Bac.Rn_est.Time,B2Bac.Rn_est.Data(:,(1:1:13000)),period_SEG5); Err_Rn_B2Ba_SEG5=mean((Rn_REAL_SEG_B-mean(Rn_B2Ba_SEG5))./Rn_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Rn_B2Ba_SEG6 = interp1(B2Bac.Rn_est.Time,B2Bac.Rn_est.Data(:,(1:1:13000)),period_SEG6); Err_Rn_B2Ba_SEG6=mean((Rn_REAL_SEG_B-mean(Rn_B2Ba_SEG6))./Rn_REAL_SEG_B.*100); Error_Rn_B2Ba_ALL=[Err_Rn_B2Ba_SEG1 Err_Rn_B2Ba_SEG2 Err_Rn_B2Ba_SEG3 Err_Rn_B2Ba_SEG4 Err_Rn_B2Ba_SEG5 Err_Rn_B2Ba_SEG6]; Test_E_SCR_2p3_Rerr=Error_Rn_B2Ba_ALL(1:3) %% Test F B3Aaa=load('Batch3a_Rec_noisy_a_v2.mat'); B3Aab=load('Batch3a_Rec_noisy_b_v2.mat'); B3Aac=load('Batch3a_Rec_noisy_c_v2.mat'); B3Aad=load('Batch3a_Rec_noisy_d_v2.mat'); simulation=9:period:128.7; t_checks = 0:0.1:130; Ln_REAL_steps=[0.0165 0.0165 0.0231 0.0231]; % insert real values Timing_imp_step_change=[0 70 70.0001 130]; Impedance_calc=interp1(Timing_imp_step_change,Ln_REAL_steps,t_checks); RN_REAL_steps=[0.5173 0.5173 0.7243 0.7243]; % insert real values Resistance_calc=interp1(Timing_imp_step_change,RN_REAL_steps,t_checks); Ln_est_B3A=interp1(B3Aab.Ln_est.Time,B3Aab.Ln_est.Data(:,(1:1:13000)),simulation);% insert real values Ln_REAL_plot=interp1(t_checks,Impedance_calc,simulation); Rn_est_B3A=interp1(B3Aac.Rn_est.Time,B3Aac.Rn_est.Data(:,(1:1:13000)),simulation);% insert real values Rn_REAL_plot=interp1(t_checks,Resistance_calc,simulation); % L error Ln_REAL_SEG_A=0.0165; Ln_REAL_SEG_B=0.0231; period=0.9; simulation=9:period:128.7; period_SEG1=18:0.9:30; Ln_B3Aa_SEG1 = interp1(B3Aab.Ln_est.Time,B3Aab.Ln_est.Data(:,(1:1:13000)),period_SEG1); Err_Ln_B3Aa_SEG1=mean((Ln_REAL_SEG_A-mean(Ln_B3Aa_SEG1))./Ln_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Ln_B3Aa_SEG2 = interp1(B3Aab.Ln_est.Time,B3Aab.Ln_est.Data(:,(1:1:13000)),period_SEG2); Err_Ln_B3Aa_SEG2=mean((Ln_REAL_SEG_A-mean(Ln_B3Aa_SEG2))./Ln_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Ln_B3Aa_SEG3 = interp1(B3Aab.Ln_est.Time,B3Aab.Ln_est.Data(:,(1:1:13000)),period_SEG3); Err_Ln_B3Aa_SEG3=mean((Ln_REAL_SEG_A-mean(Ln_B3Aa_SEG3))./Ln_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Ln_B3Aa_SEG4 = interp1(B3Aab.Ln_est.Time,B3Aab.Ln_est.Data(:,(1:1:13000)),period_SEG4); Err_Ln_B3Aa_SEG4=mean((Ln_REAL_SEG_B-mean(Ln_B3Aa_SEG4))./Ln_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Ln_B3Aa_SEG5 = interp1(B3Aab.Ln_est.Time,B3Aab.Ln_est.Data(:,(1:1:13000)),period_SEG5); Err_Ln_B3Aa_SEG5=mean((Ln_REAL_SEG_B-mean(Ln_B3Aa_SEG5))./Ln_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Ln_B3Aa_SEG6 = interp1(B3Aab.Ln_est.Time,B3Aab.Ln_est.Data(:,(1:1:13000)),period_SEG6); Err_Ln_B3Aa_SEG6=mean((Ln_REAL_SEG_B-mean(Ln_B3Aa_SEG6))./Ln_REAL_SEG_B.*100); Error_Ln_B3Aa_ALL=[Err_Ln_B3Aa_SEG1 Err_Ln_B3Aa_SEG2 Err_Ln_B3Aa_SEG3 Err_Ln_B3Aa_SEG4 Err_Ln_B3Aa_SEG5 Err_Ln_B3Aa_SEG6]; Test_F_SCR_7_Lerr=Error_Ln_B3Aa_ALL(1:3) Test_F_SCR_5_Lerr=Error_Ln_B3Aa_ALL(4:end) % R error Rn_REAL_SEG_A=0.5173; Rn_REAL_SEG_B=0.7243; period=0.9; simulation=9:period:128.7; period_SEG1=18:0.9:30; Rn_B3Aa_SEG1 = interp1(B3Aac.Rn_est.Time,B3Aac.Rn_est.Data(:,(1:1:13000)),period_SEG1); Err_Rn_B3Aa_SEG1=mean((Rn_REAL_SEG_A-mean(Rn_B3Aa_SEG1))./Rn_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Rn_B3Aa_SEG2 = interp1(B3Aac.Rn_est.Time,B3Aac.Rn_est.Data(:,(1:1:13000)),period_SEG2); Err_Rn_B3Aa_SEG2=mean((Rn_REAL_SEG_A-mean(Rn_B3Aa_SEG2))./Rn_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Rn_B3Aa_SEG3 = interp1(B3Aac.Rn_est.Time,B3Aac.Rn_est.Data(:,(1:1:13000)),period_SEG3); Err_Rn_B3Aa_SEG3=mean((Rn_REAL_SEG_A-mean(Rn_B3Aa_SEG3))./Rn_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Rn_B3Aa_SEG4 = interp1(B3Aac.Rn_est.Time,B3Aac.Rn_est.Data(:,(1:1:13000)),period_SEG4); Err_Rn_B3Aa_SEG4=mean((Rn_REAL_SEG_B-mean(Rn_B3Aa_SEG4))./Rn_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Rn_B3Aa_SEG5 = interp1(B3Aac.Rn_est.Time,B3Aac.Rn_est.Data(:,(1:1:13000)),period_SEG5); Err_Rn_B3Aa_SEG5=mean((Rn_REAL_SEG_B-mean(Rn_B3Aa_SEG5))./Rn_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Rn_B3Aa_SEG6 = interp1(B3Aac.Rn_est.Time,B3Aac.Rn_est.Data(:,(1:1:13000)),period_SEG6); Err_Rn_B3Aa_SEG6=mean((Rn_REAL_SEG_B-mean(Rn_B3Aa_SEG6))./Rn_REAL_SEG_B.*100); Error_Rn_B3Aa_ALL=[Err_Rn_B3Aa_SEG1 Err_Rn_B3Aa_SEG2 Err_Rn_B3Aa_SEG3 Err_Rn_B3Aa_SEG4 Err_Rn_B3Aa_SEG5 Err_Rn_B3Aa_SEG6]; Test_F_SCR_7_Rerr=Error_Rn_B3Aa_ALL(1:3) Test_F_SCR_5_Rerr=Error_Rn_B3Aa_ALL(4:end) % This batch is for RLS+RLCRL+noise+SCR3and1.5 ==== index 8 B3Baa=load('Batch3b_Rec_noisy_a_v2.mat'); B3Bab=load('Batch3b_Rec_noisy_b_v2.mat'); B3Bac=load('Batch3b_Rec_noisy_c_v2.mat'); B3Bad=load('Batch3b_Rec_noisy_d_v2.mat'); simulation=9:period:128.7; t_checks = 0:0.1:130; Ln_REAL_steps=[0.0384 0.0384 0.0768 0.0768]; % insert real values Timing_imp_step_change=[0 70 70.0001 130]; Impedance_calc=interp1(Timing_imp_step_change,Ln_REAL_steps,t_checks); RN_REAL_steps=[1.2071 1.2071 2.4143 2.4143]; % insert real values Resistance_calc=interp1(Timing_imp_step_change,RN_REAL_steps,t_checks); Ln_est_B3B=interp1(B3Bab.Ln_est.Time,B3Bab.Ln_est.Data(:,(1:1:13000)),simulation);% insert real values Ln_REAL_plot=interp1(t_checks,Impedance_calc,simulation); Rn_est_B3B=interp1(B3Bac.Rn_est.Time,B3Bac.Rn_est.Data(:,(1:1:13000)),simulation);% insert real values Rn_REAL_plot=interp1(t_checks,Resistance_calc,simulation); % L error Ln_REAL_SEG_A=0.0384; Ln_REAL_SEG_B=0.0768; period=0.9; simulation=9:period:128.7; period_SEG1=18:0.9:30; Ln_B3Ba_SEG1 = interp1(B3Bab.Ln_est.Time,B3Bab.Ln_est.Data(:,(1:1:13000)),period_SEG1); Err_Ln_B3Ba_SEG1=mean((Ln_REAL_SEG_A-mean(Ln_B3Ba_SEG1))./Ln_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Ln_B3Ba_SEG2 = interp1(B3Bab.Ln_est.Time,B3Bab.Ln_est.Data(:,(1:1:13000)),period_SEG2); Err_Ln_B3Ba_SEG2=mean((Ln_REAL_SEG_A-mean(Ln_B3Ba_SEG2))./Ln_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Ln_B3Ba_SEG3 = interp1(B3Bab.Ln_est.Time,B3Bab.Ln_est.Data(:,(1:1:13000)),period_SEG3); Err_Ln_B3Ba_SEG3=mean((Ln_REAL_SEG_A-mean(Ln_B3Ba_SEG3))./Ln_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Ln_B3Ba_SEG4 = interp1(B3Bab.Ln_est.Time,B3Bab.Ln_est.Data(:,(1:1:13000)),period_SEG4); Err_Ln_B3Ba_SEG4=mean((Ln_REAL_SEG_B-mean(Ln_B3Ba_SEG4))./Ln_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Ln_B3Ba_SEG5 = interp1(B3Bab.Ln_est.Time,B3Bab.Ln_est.Data(:,(1:1:13000)),period_SEG5); Err_Ln_B3Ba_SEG5=mean((Ln_REAL_SEG_B-mean(Ln_B3Ba_SEG5))./Ln_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Ln_B3Ba_SEG6 = interp1(B3Bab.Ln_est.Time,B3Bab.Ln_est.Data(:,(1:1:13000)),period_SEG6); Err_Ln_B3Ba_SEG6=mean((Ln_REAL_SEG_B-mean(Ln_B3Ba_SEG6))./Ln_REAL_SEG_B.*100); Error_Ln_B3Ba_ALL=[Err_Ln_B3Ba_SEG1 Err_Ln_B3Ba_SEG2 Err_Ln_B3Ba_SEG3 Err_Ln_B3Ba_SEG4 Err_Ln_B3Ba_SEG5 Err_Ln_B3Ba_SEG6]; Test_F_SCR_3_Lerr=Error_Ln_B3Ba_ALL(1:3) Test_F_SCR_1p5_Lerr=Error_Ln_B3Ba_ALL(4:end) % R error Rn_REAL_SEG_A=1.2071; Rn_REAL_SEG_B=2.4143; period=0.9; simulation=9:period:128.7; period_SEG1=18:0.9:30; Rn_B3Ba_SEG1 = interp1(B3Bac.Rn_est.Time,B3Bac.Rn_est.Data(:,(1:1:13000)),period_SEG1); Err_Rn_B3Ba_SEG1=mean((Rn_REAL_SEG_A-mean(Rn_B3Ba_SEG1))./Rn_REAL_SEG_A.*100); period_SEG2=(18:0.9:30)+20; Rn_B3Ba_SEG2 = interp1(B3Bac.Rn_est.Time,B3Bac.Rn_est.Data(:,(1:1:13000)),period_SEG2); Err_Rn_B3Ba_SEG2=mean((Rn_REAL_SEG_A-mean(Rn_B3Ba_SEG2))./Rn_REAL_SEG_A.*100); period_SEG3=(18:0.9:30)+40; Rn_B3Ba_SEG3 = interp1(B3Bac.Rn_est.Time,B3Bac.Rn_est.Data(:,(1:1:13000)),period_SEG3); Err_Rn_B3Ba_SEG3=mean((Rn_REAL_SEG_A-mean(Rn_B3Ba_SEG3))./Rn_REAL_SEG_A.*100); period_SEG4=(18:0.9:30)+60; Rn_B3Ba_SEG4 = interp1(B3Bac.Rn_est.Time,B3Bac.Rn_est.Data(:,(1:1:13000)),period_SEG4); Err_Rn_B3Ba_SEG4=mean((Rn_REAL_SEG_B-mean(Rn_B3Ba_SEG4))./Rn_REAL_SEG_B.*100); period_SEG5=(18:0.9:30)+80; Rn_B3Ba_SEG5 = interp1(B3Bac.Rn_est.Time,B3Bac.Rn_est.Data(:,(1:1:13000)),period_SEG5); Err_Rn_B3Ba_SEG5=mean((Rn_REAL_SEG_B-mean(Rn_B3Ba_SEG5))./Rn_REAL_SEG_B.*100); period_SEG6=(18:0.9:30)+98.5; Rn_B3Ba_SEG6 = interp1(B3Bac.Rn_est.Time,B3Bac.Rn_est.Data(:,(1:1:13000)),period_SEG6); Err_Rn_B3Ba_SEG6=mean((Rn_REAL_SEG_B-mean(Rn_B3Ba_SEG6))./Rn_REAL_SEG_B.*100); Error_Rn_B3Ba_ALL=[Err_Rn_B3Ba_SEG1 Err_Rn_B3Ba_SEG2 Err_Rn_B3Ba_SEG3 Err_Rn_B3Ba_SEG4 Err_Rn_B3Ba_SEG5 Err_Rn_B3Ba_SEG6]; Test_F_SCR_3_Rerr=Error_Rn_B3Ba_ALL(1:3) Test_F_SCR_1p5_Rerr=Error_Rn_B3Ba_ALL(4:end)