plot_roach_coefficents

0001 p_0(:,1) = dlmread('./cbass_s_matlab_scripts/pumba/amp0real.txt');
0002 p_0(:,2) = csvread('./cbass_s_matlab_scripts/pumba/amp0imag.txt');
0003 
0004 p_1(:,1) = csvread('./cbass_s_matlab_scripts/pumba/amp1real.txt');
0005 p_1(:,2) = csvread('./cbass_s_matlab_scripts/pumba/amp1imag.txt');
0006 
0007 p_2(:,1) = csvread('./cbass_s_matlab_scripts/pumba/amp2real.txt');
0008 p_2(:,2) = csvread('./cbass_s_matlab_scripts/pumba/amp2imag.txt');
0009 
0010 p_3(:,1) = csvread('./cbass_s_matlab_scripts/pumba/amp3real.txt');
0011 p_3(:,2) = csvread('./cbass_s_matlab_scripts/pumba/amp3imag.txt');
0012 
0013 p_4(:,1) = csvread('./cbass_s_matlab_scripts/pumba/amp4real.txt');
0014 p_4(:,2) = csvread('./cbass_s_matlab_scripts/pumba/amp4imag.txt');
0015 
0016 p_5(:,1) = csvread('./cbass_s_matlab_scripts/pumba/amp5real.txt');
0017 p_5(:,2) = csvread('./cbass_s_matlab_scripts/pumba/amp5imag.txt');
0018 
0019 p_6(:,1) = csvread('./cbass_s_matlab_scripts/pumba/amp6real.txt');
0020 p_6(:,2) = csvread('./cbass_s_matlab_scripts/pumba/amp6imag.txt');
0021 
0022 p_7(:,1) = csvread('./cbass_s_matlab_scripts/pumba/amp7real.txt');
0023 p_7(:,2) = csvread('./cbass_s_matlab_scripts/pumba/amp7imag.txt');
0024 
0025 %%
0026 % Make into complex arrays
0027 p(1,:) = complex(p_0(:,1),p_0(:,2))
0028 p(2,:) = complex(p_1(:,1),p_1(:,2))
0029 p(3,:) = complex(p_2(:,1),p_2(:,2))
0030 p(4,:) = complex(p_3(:,1),p_3(:,2))
0031 p(5,:) = complex(p_4(:,1),p_4(:,2))
0032 p(6,:) = complex(p_5(:,1),p_5(:,2))
0033 p(7,:) = complex(p_6(:,1),p_6(:,2))
0034 p(8,:) = complex(p_7(:,1),p_7(:,2))
0035 
0036 %% Now make the arrays into channel data 1:4
0037 pumba = zeros(1,64)
0038 pumba(1,1:2:64)=(p(1,:))
0039 pumba(1,2:2:64)=(p(2,:))
0040 pumba(2,1:2:64)=(p(3,:))
0041 pumba(2,2:2:64)=(p(4,:))
0042 pumba(3,1:2:64)=(p(5,:))
0043 pumba(3,2:2:64)=(p(6,:))
0044 pumba(4,1:2:64)=(p(7,:))
0045 pumba(4,2:2:64)=(p(8,:))
0046 %%
0047 % Plot the amplitude and phase of each set of coefficients
0048 figure(1)
0049 for i=1:4
0050     subplot(2,2,i)
0051     plot(abs(pumba(i,:)))
0052     ylabel('amplitude')
0053     xlabel('channel')
0054     title(['Channel: ',num2str(i-1)])
0055     gtitle('Pumba Coeffs: Amplitude')
0056 end
0057 
0058 figure(2)
0059 for i=1:4
0060     subplot(2,2,i)
0061     plot(rad2deg(angle(pumba(i,:))))
0062     ylabel('phase, rad')
0063     xlabel('channel')
0064     title(['Channel: ',num2str(i-1)])
0065     gtitle('Pumba Coeffs: Phase')
0066 end
0067 %%
0068 %Now the same for Timon
0069 t_0(:,1) = csvread('./cbass_s_matlab_scripts/timon/amp0real.txt');
0070 t_0(:,2) = csvread('./cbass_s_matlab_scripts/timon/amp0imag.txt');
0071 
0072 t_1(:,1) = csvread('./cbass_s_matlab_scripts/timon/amp1real.txt');
0073 t_1(:,2) = csvread('./cbass_s_matlab_scripts/timon/amp1imag.txt');
0074 
0075 t_2(:,1) = csvread('./cbass_s_matlab_scripts/timon/amp2real.txt');
0076 t_2(:,2) = csvread('./cbass_s_matlab_scripts/timon/amp2imag.txt');
0077 
0078 t_3(:,1) = csvread('./cbass_s_matlab_scripts/timon/amp3real.txt');
0079 t_3(:,2) = csvread('./cbass_s_matlab_scripts/timon/amp3imag.txt');
0080 
0081 t_4(:,1) = csvread('./cbass_s_matlab_scripts/timon/amp4real.txt');
0082 t_4(:,2) = csvread('./cbass_s_matlab_scripts/timon/amp4imag.txt');
0083 
0084 t_5(:,1) = csvread('./cbass_s_matlab_scripts/timon/amp5real.txt');
0085 t_5(:,2) = csvread('./cbass_s_matlab_scripts/timon/amp5imag.txt');
0086 
0087 t_6(:,1) = csvread('./cbass_s_matlab_scripts/timon/amp6real.txt');
0088 t_6(:,2) = csvread('./cbass_s_matlab_scripts/timon/amp6imag.txt');
0089 
0090 t_7(:,1) = csvread('./cbass_s_matlab_scripts/timon/amp7real.txt');
0091 t_7(:,2) = csvread('./cbass_s_matlab_scripts/timon/amp7imag.txt');
0092 
0093 %%
0094 % Make into complex arrays
0095 t(1,:) = complex(t_0(:,1),t_0(:,2))
0096 t(2,:) = complex(t_1(:,1),t_1(:,2))
0097 t(3,:) = complex(t_2(:,1),t_2(:,2))
0098 t(4,:) = complex(t_3(:,1),t_3(:,2))
0099 t(5,:) = complex(t_4(:,1),t_4(:,2))
0100 t(6,:) = complex(t_5(:,1),t_5(:,2))
0101 t(7,:) = complex(t_6(:,1),t_6(:,2))
0102 t(8,:) = complex(t_7(:,1),t_7(:,2))
0103 %% Now make the arrays into channel data 1:4
0104 timon = zeros(1,64)
0105 timon(1,1:2:64)=(t(1,:))
0106 timon(1,2:2:64)=(t(2,:))
0107 timon(2,1:2:64)=(t(3,:))
0108 timon(2,2:2:64)=(t(4,:))
0109 timon(3,1:2:64)=(t(5,:))
0110 timon(3,2:2:64)=(t(6,:))
0111 timon(4,1:2:64)=(t(7,:))
0112 timon(4,2:2:64)=(t(8,:))
0113 %%
0114 % Plot the amplitude and phase of each set of coefficients
0115 figure(3)
0116 for i=1:4
0117     subplot(2,2,i)
0118     plot(abs(timon(i,:)))
0119     ylabel('amplitude')
0120     xlabel('channel')
0121     title(['Channel: ',num2str(i-1)])
0122     gtitle('Timon Coeffs: Amplitude')
0123 end
0124 
0125 figure(4)
0126 for i=1:4
0127     subplot(2,2,i)
0128     plot(rad2deg(angle(timon(i,:))))
0129     ylabel('phase, rad')
0130     xlabel('channel')
0131     title(['Channel: ',num2str(i-1)])
0132     gtitle('Timon Coeffs: Phase')
0133 end