clear all list_constants; %% define frequency points %N=10; fmin = 1; fmax = 100; %frequency = linspace(fmin,fmax,N); %% calling function freq = frequency; % freq = logspace(log10(fmin),log10(fmax),N); would be another interesting option r_m = MirrorRad; m = M4; %Suspension Thermal Noise susTN = total_sus_therm(freq); susTN_strain = susTN / (arm_length/2); susTN_x = x_therm(freq); susTN_z = z_therm(freq); susTN_v = v_therm(freq); susTN_x_strain = susTN_x / (arm_length/2); susTN_z_strain = susTN_z / (arm_length/2); susTN_v_strain = susTN_v / (arm_length/2); %Seismic (longitudinal transmissibility) long_seis_ap = long_seis_approx(freq); %Seismic (vertical transmissibility) cross_coupling=1/1000; ver_seis_ap = vertical_seis_approx(freq); %Seismic BSC = BSC_req(freq); long_trans_BSC_approx_strain = BSC .* long_seis_ap ./ (arm_length/2); ver_trans_BSC_approx_strain = 0*cross_coupling .* BSC .* ver_seis_ap./ (arm_length/2); total_seis_BSC_approx_strain = sqrt(long_trans_BSC_approx_strain.^2 + ver_trans_BSC_approx_strain.^2); %Coating Thermal %total_coating_thermal = total_coating_therm(freq)./ (arm_length/2); total_coating_thermal = total_coating_therm_final(freq);%./ (arm_length/2); %Clipping Loss total_clipping_loss = total_clipping(r_m) %Total total_strain=total_seis_BSC_approx_strain + susTN_strain + total_coating_thermal; %% Save/name mass_string=num2str(m); L4_string=num2str(F_l); tot_len=num2str(total_length); name_string=[mass_string,'_',L4_string,'_',tot_len,'_770MPa_06m']; %% (1) aLIGO % aLIGO.coat_aLIGO=total_coating_thermal; % aLIGO.seis_aLIGO=total_seis_BSC_approx_strain; % aLIGO.seis_BSC_aLIGO=BSC; % aLIGO.sus_aLIGO=susTN_strain; % aLIGO.sus_x_aLIGO=susTN_x_strain; % aLIGO.sus_z_aLIGO=susTN_z_strain; % aLIGO.sus_v_aLIGO=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/test/aLIGO.mat'],'aLIGO'); %% (1-1') aLIGO (1GPa) % aLIGO_1GPa.coat_aLIGO_1GPa=total_coating_thermal; % aLIGO_1GPa.seis_aLIGO_1GPa=total_seis_BSC_approx_strain; % aLIGO_1GPa.sus_aLIGO_1GPa=susTN_strain; % aLIGO_1GPa.sus_x_aLIGO_1GPa=susTN_x_strain; % aLIGO_1GPa.sus_z_aLIGO_1GPa=susTN_z_strain; % aLIGO_1GPa.sus_v_aLIGO_1GPa=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_1GPa.mat'],'aLIGO_1GPa'); % %% (1-2') aLIGO (1.5GPa) % aLIGO_15GPa.coat_aLIGO_15GPa=total_coating_thermal; % aLIGO_15GPa.seis_aLIGO_15GPa=total_seis_BSC_approx_strain; % aLIGO_15GPa.sus_aLIGO_15GPa=susTN_strain; % aLIGO_15GPa.sus_x_aLIGO_15GPa=susTN_x_strain; % aLIGO_15GPa.sus_z_aLIGO_15GPa=susTN_z_strain; % aLIGO_15GPa.sus_v_aLIGO_15GPa=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_15GPa.mat'],'aLIGO_15GPa'); %% (1-3') aLIGO (2GPa) % aLIGO_2GPa.coat_aLIGO_2GPa=total_coating_thermal; % aLIGO_2GPa.seis_aLIGO_2GPa=total_seis_BSC_approx_strain; % aLIGO_2GPa.sus_aLIGO_2GPa=susTN_strain; % aLIGO_2GPa.sus_x_aLIGO_2GPa=susTN_x_strain; % aLIGO_2GPa.sus_z_aLIGO_2GPa=susTN_z_strain; % aLIGO_2GPa.sus_v_aLIGO_2GPa=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_2GPa.mat'],'aLIGO_2GPa'); %% (1-4') aLIGO (2.5GPa) % aLIGO_25GPa.coat_aLIGO_25GPa=total_coating_thermal; % aLIGO_25GPa.seis_aLIGO_25GPa=total_seis_BSC_approx_strain; % aLIGO_25GPa.sus_aLIGO_25GPa=susTN_strain; % aLIGO_25GPa.sus_x_aLIGO_25GPa=susTN_x_strain; % aLIGO_25GPa.sus_z_aLIGO_25GPa=susTN_z_strain; % aLIGO_25GPa.sus_v_aLIGO_25GPa=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_25GPa.mat'],'aLIGO_25GPa'); %% (1-5') aLIGO (3GPa) % aLIGO_3GPa.coat_aLIGO_3GPa=total_coating_thermal; % aLIGO_3GPa.seis_aLIGO_3GPa=total_seis_BSC_approx_strain; % aLIGO_3GPa.sus_aLIGO_3GPa=susTN_strain; % aLIGO_3GPa.sus_x_aLIGO_3GPa=susTN_x_strain; % aLIGO_3GPa.sus_z_aLIGO_3GPa=susTN_z_strain; % aLIGO_3GPa.sus_v_aLIGO_3GPa=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_3GPa.mat'],'aLIGO_3GPa'); %% (1-6') aLIGO (3.5GPa) % aLIGO_35GPa.coat_aLIGO_35GPa=total_coating_thermal; % aLIGO_35GPa.seis_aLIGO_35GPa=total_seis_BSC_approx_strain; % aLIGO_35GPa.sus_aLIGO_35GPa=susTN_strain; % aLIGO_35GPa.sus_x_aLIGO_35GPa=susTN_x_strain; % aLIGO_35GPa.sus_z_aLIGO_35GPa=susTN_z_strain; % aLIGO_35GPa.sus_v_aLIGO_35GPa=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_35GPa.mat'],'aLIGO_35GPa'); %% (1-7') aLIGO (4GPa) % aLIGO_4GPa.coat_aLIGO_4GPa=total_coating_thermal; % aLIGO_4GPa.seis_aLIGO_4GPa=total_seis_BSC_approx_strain; % aLIGO_4GPa.sus_aLIGO_4GPa=susTN_strain; % aLIGO_4GPa.sus_x_aLIGO_4GPa=susTN_x_strain; % aLIGO_4GPa.sus_z_aLIGO_4GPa=susTN_z_strain; % aLIGO_4GPa.sus_v_aLIGO_4GPa=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_4GPa.mat'],'aLIGO_4GPa'); %% (1-8') aLIGO (4.5GPa) % aLIGO_45GPa.coat_aLIGO_45GPa=total_coating_thermal; % aLIGO_45GPa.seis_aLIGO_45GPa=total_seis_BSC_approx_strain; % aLIGO_45GPa.sus_aLIGO_45GPa=susTN_strain; % aLIGO_45GPa.sus_x_aLIGO_45GPa=susTN_x_strain; % aLIGO_45GPa.sus_z_aLIGO_45GPa=susTN_z_strain; % aLIGO_45GPa.sus_v_aLIGO_45GPa=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_45GPa.mat'],'aLIGO_45GPa'); % %% (1-9') aLIGO (5GPa) % aLIGO_5GPa.coat_aLIGO_5GPa=total_coating_thermal; % aLIGO_5GPa.seis_aLIGO_5GPa=total_seis_BSC_approx_strain; % aLIGO_5GPa.sus_aLIGO_5GPa=susTN_strain; % aLIGO_5GPa.sus_x_aLIGO_5GPa=susTN_x_strain; % aLIGO_5GPa.sus_z_aLIGO_5GPa=susTN_z_strain; % aLIGO_5GPa.sus_v_aLIGO_5GPa=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_5GPa.mat'],'aLIGO_5GPa'); %% (1-1) aLIGO + 2.14m aLIGO_LT214.coat_aLIGO_LT214=total_coating_thermal; aLIGO_LT214.seis_aLIGO_LT214=total_seis_BSC_approx_strain; aLIGO_LT214.sus_aLIGO_LT214=susTN_strain; save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_LT214.mat'],'aLIGO_LT214'); %% (2) aAH: aLIGO + 1.2GPa % aLIGO12GPa.coat_aLIGO12GPa=total_coating_thermal; % aLIGO12GPa.seis_aLIGO12GPa=total_seis_BSC_approx_strain; % aLIGO12GPa.sus_aLIGO12GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/paper/aLIGO12GPa.mat'],'aLIGO12GPa'); %% (2-2) aAH: aLIGO + 1.2GPa + 1.1m + 2.14m % aLIGO_opt_11_214_12GPa.coat_aLIGO_opt_11_214_12GPa=total_coating_thermal; % aLIGO_opt_11_214_12GPa.seis_aLIGO_opt_11_214_12GPa=total_seis_BSC_approx_strain; % aLIGO_opt_11_214_12GPa.sus_aLIGO_opt_11_214_12GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/mat/aLIGO_opt_11_214_12GPa.mat'],'aLIGO_opt_11_214_12GPa'); %% (2-1) aAH: aLIGO + 1.2GPa + 2.14m % aLIGO_LT214_12GPa.coat_aLIGO_LT214_12GPa=total_coating_thermal; % aLIGO_LT214_12GPa.seis_aLIGO_LT214_12GPa=total_seis_BSC_approx_strain; % aLIGO_LT214_12GPa.sus_aLIGO_LT214_12GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_LT214_12GPa.mat'],'aLIGO_LT214_12GPa'); %% (3) aBL: 40kg + 770MPa + 1.1m % aLIGO770MPa_11m.coat_aLIGO770MPa_11m=total_coating_thermal; % aLIGO770MPa_11m.seis_aLIGO770MPa_11m=total_seis_BSC_approx_strain; % aLIGO770MPa_11m.sus_aLIGO770MPa_11m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO770MPa_11m.mat'],'aLIGO770MPa_11m'); %% (4) aBH: 40kg + 1.2GPa + 1.1m % aLIGO12GPa_11m.coat_aLIGO12GPa_11m=total_coating_thermal; % aLIGO12GPa_11m.seis_aLIGO12GPa_11m=total_seis_BSC_approx_strain; % aLIGO12GPa_11m.sus_aLIGO12GPa_11m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO12GPa_11m.mat'],'aLIGO12GPa_11m'); %% (1) 1AL: 80kg + 770MPa + 0.6m % m80_770MPa_06m.coat_m80_770MPa_06m=total_coating_thermal; % m80_770MPa_06m.seis_m80_770MPa_06m=total_seis_BSC_approx_strain; % m80_770MPa_06m.sus_m80_770MPa_06m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/m80_770MPa_06m.mat'],'m80_770MPa_06m'); % m80_770MPa_0582m.coat_m80_770MPa_0582m=total_coating_thermal; % m80_770MPa_0582m.seis_m80_770MPa_0582m=total_seis_BSC_approx_strain; % m80_770MPa_0582m.sus_m80_770MPa_0582m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/m80_770MPa_0582m.mat'],'m80_770MPa_0582m'); %% (2) 1AH: 80kg + 1.2GPa + 0.6m % m80_12GPa_06m.coat_m80_12GPa_06m=total_coating_thermal; % m80_12GPa_06m.seis_m80_12GPa_06m=total_seis_BSC_approx_strain; % m80_12GPa_06m.sus_m80_12GPa_06m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/m80_12GPa_06m.mat'],'m80_12GPa_06m'); %% (3) 1BL: 80kg + 770MPa + 1.1m % m80_770MPa_11m.coat_m80_770MPa_11m=total_coating_thermal; % m80_770MPa_11m.seis_m80_770MPa_11m=total_seis_BSC_approx_strain; % m80_770MPa_11m.sus_m80_770MPa_11m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/m80_770MPa_11m.mat'],'m80_770MPa_11m'); %% (4) 1BH: 80kg + 1.2GPa + 1.1m % m80_12GPa_11m.coat_m80_12GPa_11m=total_coating_thermal; % m80_12GPa_11m.seis_m80_12GPa_11m=total_seis_BSC_approx_strain; % m80_12GPa_11m.sus_m80_12GPa_11m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/paper/m80_12GPa_11m.mat'],'m80_12GPa_11m'); %% (1) 2AL: 160kg + 770MPa + 0.6m % m160_770MPa_06m.coat_m160_770MPa_06m=total_coating_thermal; % m160_770MPa_06m.seis_m160_770MPa_06m=total_seis_BSC_approx_strain; % m160_770MPa_06m.sus_m160_770MPa_06m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/m160_770MPa_06m.mat'],'m160_770MPa_06m'); % m160_770MPa_0582m.coat_m160_770MPa_0582m=total_coating_thermal; % m160_770MPa_0582m.seis_m160_770MPa_0582m=total_seis_BSC_approx_strain; % m160_770MPa_0582m.sus_m160_770MPa_0582m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/m160_770MPa_0582m.mat'],'m160_770MPa_0582m'); %% (2) 2AH: 160kg + 1.2GPa + 0.6m % m160_12GPa_06m.coat_m160_12GPa_06m=total_coating_thermal; % m160_12GPa_06m.seis_m160_12GPa_06m=total_seis_BSC_approx_strain; % m160_12GPa_06m.sus_m160_12GPa_06m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/m160_12GPa_06m.mat'],'m160_12GPa_06m'); %% (3) 2BL: 160kg + 770MPa + 1.1m % m160_770MPa_11m.coat_m160_770MPa_11m=total_coating_thermal; % m160_770MPa_11m.seis_m160_770MPa_11m=total_seis_BSC_approx_strain; % m160_770MPa_11m.sus_m160_770MPa_11m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/m160_770MPa_11m.mat'],'m160_770MPa_11m'); %% (4) 2BH: 160kg + 1.2GPa + 1.1m % m160_12GPa_11m.coat_m160_12GPa_11m=total_coating_thermal; % m160_12GPa_11m.seis_m160_12GPa_11m=total_seis_BSC_approx_strain; % m160_12GPa_11m.sus_m160_12GPa_11m=susTN_strain; % save(['/Users/klee/Documents/Matlab/1718_Thesis/paper/m160_12GPa_11m.mat'],'m160_12GPa_11m'); %% aLIGO upgrade #1: 40/40/22/22, 0.6m, 1.6m, 1.54GPa, 5mm % aLIGO_154GPa.coat_aLIGO_154GPa=total_coating_thermal; % aLIGO_154GPa.seis_aLIGO_154GPa=total_seis_BSC_approx_strain; % aLIGO_154GPa.sus_aLIGO_154GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/aLIGO_154GPa.mat'],'aLIGO_154GPa'); %% aLIGO upgrade #2: 40/opt, 0.6m, 2.14m, 1.54GPa, 5mm % aLIGO_opt_06_214_154GPa.coat_aLIGO_opt_06_214_154GPa=total_coating_thermal; % aLIGO_opt_06_214_154GPa.seis_aLIGO_opt_06_214_154GPa=total_seis_BSC_approx_strain; % aLIGO_opt_06_214_154GPa.sus_aLIGO_opt_06_214_154GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/aLIGO_opt_06_214_154GPa.mat'],'aLIGO_opt_06_214_154GPa'); %% aLIGO upgrade #3: 40/opt, 1.1m, 2.14m, 1.54GPa, 5mm % aLIGO_opt_11_214_154GPa.coat_aLIGO_opt_11_214_154GPa=total_coating_thermal; % aLIGO_opt_11_214_154GPa.seis_aLIGO_opt_11_214_154GPa=total_seis_BSC_approx_strain; % aLIGO_opt_11_214_154GPa.sus_aLIGO_opt_11_214_154GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/aLIGO_opt_11_214_154GPa.mat'],'aLIGO_opt_11_214_154GPa'); %% aLIGO upgrade #4: 80/opt, 1.1m, 2.14m, 1.54GPa, 5mm % m80_11_214_154GPa.coat_80_11_214_154GPa=total_coating_thermal; % m80_11_214_154GPa.seis_80_11_214_154GPa=total_seis_BSC_approx_strain; % m80_11_214_154GPa.sus_80_11_214_154GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_11_214_154GPa.mat'],'m80_11_214_154GPa'); %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_06_214_770MPa'); %% aLIGO upgrade #5: 160/opt, 1.1m, 2.14m, 1.54GPa, 5mm % m160_11_214_154GPa.coat_160_11_214_154GPa=total_coating_thermal; % m160_11_214_154GPa.seis_160_11_214_154GPa=total_seis_BSC_approx_strain; % m160_11_214_154GPa.sus_160_11_214_154GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m160_11_214_154GPa.mat'],'m160_11_214_154GPa'); % %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_11_214_770MPa'); %% aLIGO upgrade #4-1: 80/opt, 0.6m, 1.6m, 1.54GPa, 5mm % m80_06_16_154GPa.coat_80_06_16_154GPa=total_coating_thermal; % m80_06_16_154GPa.seis_80_06_16_154GPa=total_seis_BSC_approx_strain; % m80_06_16_154GPa.sus_80_06_16_154GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_06_16_154GPa.mat'],'m80_06_16_154GPa'); % %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_06_16_154GPa'); %% aLIGO upgrade #4-2: 80/opt, 1.1m, 1.6m, 1.54GPa, 5mm % m80_11_16_154GPa.coat_80_11_16_154GPa=total_coating_thermal; % m80_11_16_154GPa.seis_80_11_16_154GPa=total_seis_BSC_approx_strain; % m80_11_16_154GPa.sus_80_11_16_154GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_11_16_154GPa.mat'],'m80_11_16_154GPa'); % %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_11_16_154GPa'); %% aLIGO upgrade #5-1: 80/opt, 0.6m, 2.14m, 1.54GPa, 5mm % m80_06_214_154GPa.coat_80_06_214_154GPa=total_coating_thermal; % m80_06_214_154GPa.seis_80_06_214_154GPa=total_seis_BSC_approx_strain; % m80_06_214_154GPa.sus_80_06_214_154GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_06_214_154GPa.mat'],'m80_06_214_154GPa'); % %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_06_214_154GPa'); %% aLIGO upgrade #5-2: 80/opt, 1.1m, 2.14m, 1.54GPa, 5mm % m80_11_214_154GPa.coat_80_11_214_154GPa=total_coating_thermal; % m80_11_214_154GPa.seis_80_11_214_154GPa=total_seis_BSC_approx_strain; % m80_11_214_154GPa.sus_80_11_214_154GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_11_214_154GPa.mat'],'m80_11_214_154GPa'); % %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_11_214_154GPa'); %% aLIGO upgrade 1.2GPa: 40/40/22/22, 0.6m, 1.6m, 1.2GPa, 5mm % aLIGO_12GPa.coat_aLIGO_12GPa=total_coating_thermal; % aLIGO_12GPa.seis_aLIGO_12GPa=total_seis_BSC_approx_strain; % aLIGO_12GPa.sus_aLIGO_12GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/aLIGO_12GPa.mat'],'aLIGO_12GPa'); %% aLIGO upgrade 1.2GPa: 160/opt, 1.1m, 2.14m, 1.2GPa, 5mm % m160_11_214_154GPa.coat_160_11_214_154GPa=total_coating_thermal; % m160_11_214_154GPa.seis_160_11_214_154GPa=total_seis_BSC_approx_strain; % m160_11_214_12GPa.sus_160_11_214_12GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m160_11_214_12GPa.mat'],'m160_11_214_12GPa'); % %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_11_214_770MPa'); %% aLIGO upgrade 1.2GPa: 80/opt, 1.1m, 2.14m, 1.2GPa, 5mm % m80_11_214_154GPa.coat_80_11_214_154GPa=total_coating_thermal; % m80_11_214_154GPa.seis_80_11_214_154GPa=total_seis_BSC_approx_strain; % m80_11_214_12GPa.sus_80_11_214_12GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_11_214_12GPa.mat'],'m80_11_214_12GPa'); % %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_11_214_154GPa'); %% aLIGO upgrade 1.2GPa: 40/opt, 1.1m, 2.14m, 1.2GPa, 5mm % aLIGO_opt_11_214_154GPa.coat_aLIGO_opt_11_214_154GPa=total_coating_thermal; % aLIGO_opt_11_214_154GPa.seis_aLIGO_opt_11_214_154GPa=total_seis_BSC_approx_strain; % aLIGO_opt_11_214_12GPa.sus_aLIGO_opt_11_214_12GPa=susTN_strain; % save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/aLIGO_opt_11_214_12GPa.mat'],'aLIGO_opt_11_214_12GPa'); %% (L4_test1) aLIGO 0.6, 2.14 % aLIGO_06_214.coat_aLIGO_06_214=total_coating_thermal; % aLIGO_06_214.seis_aLIGO_06_214=total_seis_BSC_approx_strain; % aLIGO_06_214.sus_aLIGO_06_214=susTN_strain; % aLIGO_06_214.sus_x_aLIGO_06_214=susTN_x_strain; % aLIGO_06_214.sus_z_aLIGO_06_214=susTN_z_strain; % aLIGO_06_214.sus_v_aLIGO_06_214=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_06_214.mat'],'aLIGO_06_214'); %% (L4_test2) aLIGO 1.1, 2.14 % aLIGO_11_214.coat_aLIGO_11_214=total_coating_thermal; % aLIGO_11_214.seis_aLIGO_11_214=total_seis_BSC_approx_strain; % aLIGO_11_214.sus_aLIGO_11_214=susTN_strain; % aLIGO_11_214.sus_x_aLIGO_11_214=susTN_x_strain; % aLIGO_11_214.sus_z_aLIGO_11_214=susTN_z_strain; % aLIGO_11_214.sus_v_aLIGO_11_214=susTN_v_strain; % % save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_11_214.mat'],'aLIGO_11_214'); %% =========figure======== %{ %tic figure(1) title(['Frequency vs Strain for M4 = ',num2str(m),' kg']) loglog(freq,total_seis_BSC_approx_strain,'c'); hold on loglog(freq,susTN_strain,'k'); hold on loglog(freq,total_coating_thermal,'r'); hold on grid on xlabel('Frequency (Hz)'); ylabel('Strain (1/rtHz)'); legend('seis-total','sus ther','coating') axis([0 100 1e-30 1e-12]) figure(2) title(['Frequency vs Strain for M4 = ',num2str(m),' kg']) loglog(freq,total_seis_BSC_approx_strain,'b'); hold on loglog(freq,susTN_strain,'g'); hold on loglog(freq,total_coating_thermal,'r'); hold on loglog(freq,total_strain,'k'); hold on grid on xlabel('Frequency (Hz)'); ylabel('Strain (1/rtHz)'); legend('seis-total','sus ther','coating','total') axis([0 100 1e-30 1e-12]) %toc %}