extended_ortho_des_frank.m

calculation of the autocorrelation peak sidelobe levels and corsscorrelation. CE sequence, frank aut

clear,close all;
L=9;
sqL=3;
s1=frank_seq(sqL,L,1);
s2=frank_seq(sqL,L,1);
S4=[s1,s2,s1,s2;-s2,s1,-s2,s1;s1,s2,-s1,-s2;-s2,s1,s2,-s1];

m=power_integer_2(L*4);
nfft=2^(m);

sr1fft=fft(S4(1,:),nfft);
sr2fft=fft(S4(2,:),nfft);
sr3fft=fft(S4(3,:),nfft);
sr4fft=fft(S4(4,:),nfft);
time=(-(nfft/2)):1:(nfft/2-1);
figure(1);
%1,2
Cc12=fftshift(abs(ifft(sr1fft.* conj(sr2fft)))/(4*L));
subplot(3,2,1);
plot(time,Cc12);
xlabel('sample shift');
ylabel('Crosscorrelation Amplitude');
title('1,2');
%1,3
Cc13=fftshift(abs(ifft(sr1fft.* conj(sr3fft)))/(4*L));
subplot(3,2,2);
plot(time,Cc13);
xlabel('sample shift');
ylabel('Crosscorrelation Amplitude');
title('1,3');
%1,4
Cc14=fftshift(abs(ifft(sr1fft.* conj(sr4fft)))/(4*L));
subplot(3,2,3);
plot(time,Cc14);
xlabel('sample shift');
ylabel('Crosscorrelation Amplitude');
title('1,4');
%2,3
Cc23=fftshift(abs(ifft(sr2fft.* conj(sr3fft)))/(4*L));
subplot(3,2,4);
plot(time,Cc23);
xlabel('sample shift');
ylabel('Crosscorrelation Amplitude');
title('2,3');
%2,4
Cc24=fftshift(abs(ifft(sr2fft.* conj(sr4fft)))/(4*L));
subplot(3,2,5);
plot(time,Cc24);
xlabel('sample shift');
ylabel('Crosscorrelation Amplitude');
title('1,3');
%3,4
Cc34=fftshift(abs(ifft(sr3fft.* conj(sr4fft)))/(4*L));
subplot(3,2,6);
plot(time,Cc34);
xlabel('sample shift');
ylabel('Crosscorrelation Amplitude');
title('3,4');
gtext('Cross Correlation between Polyphase Sequences of length 36 using Orthogonal Matrix Design from Frank Squences');