?? carrier_x4d.m
字號:
% QPSK Demod Test
clear;
load c:\qpsk\lecroy_1.dir\sc1005.dat;
offset=1;
size=1024*4;
wave=sc1005';
fs=1e9;
% Display wave form
plot(wave(1:100));
grid
title('IF Signal');
%%%%%%pause;
data_fft=10*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
(wave(offset:offset+size-1))))));
data_fft=data_fft-max(data_fft);
nnnn=size;
ff=-nnnn/2:nnnn/2-1;
ff=ff/nnnn*2*fs/2;
plot(ff,data_fft);
title(' IF Data Spectrum Before Filtering ');
axis([400e6, 500e6 -50 0]);
grid;
xlabel('FREQUENCY (HZ)');
ylabel('MAGNITUDE (dB) ');
%pause;
% generate local carrier
f_lo=(480-50)*1e6; % carrier
t=[0:length(wave)-1]/fs;
carrier_a=sin(2*pi*f_lo*t); % carrier sine;
if_data=wave.*carrier_a; % down converting
% % plot spectrum;
data=if_data;
data_fft=20*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
data(offset:offset+size-1)))));
ff=-nnnn/2:nnnn/2-1;
ff=ff/nnnn*2*fs/2;
plot(ff,data_fft);
title(' I Data Spectrum After Down Converting');
grid;
%%%pause;
axis([0 100e6 -10 30]);
%%%pause;
%************************************
% Design A Low-Pass Filter *
%************************************
f_lowpass=[ 0, 65/500, 80/500 , 1. ];
m= [ 1., 1.0 , 1e-5, 1e-5 ];
n1=80;
b=firls(n1, f_lowpass, m);
a=1;
moduc=2^(-23);
%b=round(b/moduc/max(b))*moduc;
[amp,freq]=freqz(b,a,4096,fs);
amp=20*log10( (abs(amp)+1e-8) / (max(abs(amp)) +1e-8) );
plot(freq,amp);
grid
%%%pause
axis([0e6 100e6 -50 0]);
title('N1=100 taps FIR low pass filter.');
%%%pause;
if_data=filter(b,a,if_data);
if_50_data=if_data;
%***************************
% Plot Spectrum *
%***************************
data=if_data;
offset=100;
data_fft=20*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
data(offset:offset+size-1)))));
ff=-nnnn/2:nnnn/2-1;
ff=ff/nnnn*2*fs/2;
plot(ff,data_fft);
title(' IF Data Spectrum After Low-Pass Filter');
grid;
%%%pause;
axis([0e6 100e6 -50 30]);
%%%pause;
%***************************
% Times 2 *
%***************************
if2_data=if_data.*if_data;
%***************************
% Plot Spectrum *
%***************************
offset=100;
data=if2_data;
data_fft=20*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
data(offset:offset+size-1)))));
ff=-nnnn/2:nnnn/2-1;
ff=ff/nnnn*2*fs/2;
plot(ff,data_fft);
axis([0e6 250e6 -50 30]);
title(' IF Data Spectrum After X2');
grid;
%pause;
%***************************
% Times 4 *
%***************************
if4_data=if2_data.*if2_data;
%***************************
% Plot Spectrum *
%***************************
offset=100;
data=if4_data;
data_fft=20*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
data(offset:offset+size-1)))));
ff=-nnnn/2:nnnn/2-1;
ff=ff/nnnn*2*fs/2;
plot(ff,data_fft);
axis([200e6 250e6 -100 -30]);
title(' IF Data Spectrum After X4');
grid;
%pause;
%*********************************
% Design A Low-Pass Filter *
%*********************************
b=fir1(400,[213/500 220/500]);
a=1;
moduc=2^(-23);
%b=round(b/moduc/max(b))*moduc;
[amp,freq]=freqz(b,a,4096,fs);
amp=20*log10( (abs(amp)+1e-8) / (max(abs(amp)) +1e-8) );
plot(freq,amp);
title(' Band-Pass Filter After X2');
axis([200e6 250e6 -50 0]);
grid
%%pause
if4_data=filter(b,a,if4_data);
plot(if4_data(1000:3100));
grid;
title('IF Waveform X4');
%%pause;
%***************************
% sign *
%***************************
if4_data=sign(if4_data);
%***************************
% Plot Spectrum *
%***************************
offset=100;
data=if4_data;
data_fft=20*log10(abs(fftshift(fft(kaiser(size,11.)'.*...
data(offset:offset+size-1)))));
ff=-nnnn/2:nnnn/2-1;
ff=ff/nnnn*2*fs/2;
plot(ff,data_fft);
%axis([200e6 250e6 -100 -30]);
title(' IF Data Spectrum After X4');
grid;
%pause;
if4_data=filter(b,a,if4_data);
plot(if4_data);
grid;
title('IF Waveform X4 After Sign 1');
%pause;
for ii=1:10
step=300;
plot(if4_data(100+(ii-1)*step:ii*step+100));
grid;
title('IF Waveform after 1 sign ');
%%pause;
end;
%***********************
% count cycle
%***********************
carrier_ref=if4_data(1000:5000-100);
length_ref=length(carrier_ref);
cycle_count=0;
for i=1:length_ref-1;
if carrier_ref(i)*carrier_ref(i+1) <=0
cycle_count=cycle_count+1;
end
end
carrier_cycle=(length_ref-1)/fs/cycle_count*2;
carrier_frequency=1/carrier_cycle;
%*******************
% Generate LO *
%*******************
length_if=length(if4_data);
t=[0:length_if-1]/fs;
sin_lo=sin(2*pi*(carrier_frequency+1e3)*t);
for ii=1:20
step=200;
plot(if4_data(500+(ii-1)*step:ii*step+500));
hold on;
aaa=5;
plot(sin_lo(500+aaa+(ii-1)*step:ii*step+500+aaa),'r');
hold off;
grid;
title('Carrier and Signal Waveform ');
%pause;
end;
true_carrier_frequency=carrier_frequency/4;
%*******************
% Generate LO *
%*******************
for i_phase=34-10:34+10;
length_wave=length(if_50_data);
t=[0:length_wave-1]/fs;
phase_offset=pi*(i_phase*5)/180;
df=-0e3;
lo_cos=cos(2*pi*(df+true_carrier_frequency)*t+phase_offset);
lo_sin=sin(2*pi*(df+true_carrier_frequency)*t+phase_offset);
%************************
% Check Carrier Phase
%***********************
%for i=1:10
% step=400;
%plot(if_50_data(1+(i-1)*step:i*step)*8);
%hold on;
%aaa=-1;
%plot(lo_sin(1+aaa+(i-1)*step:i*step+aaa),'r');
%plot(lo_cos(1+aaa+(i-1)*step:i*step+aaa),'g');
%hold off;
%grid;
%title('***** Check Carrier Phase ***** ');
%%pause
%end
%****************
% Down Convert LO
%***************
i_data=if_50_data.*lo_cos;
q_data=if_50_data.*lo_sin;
%*********************************
% Design A Low-Pass Filter *
%*********************************
b=fir1(100,25/500);
a=1;
moduc=2^(-23);
%b=round(b/moduc/max(b))*moduc;
[amp,freq]=freqz(b,a,4096,fs);
amp=20*log10( (abs(amp)+1e-8) / (max(abs(amp)) +1e-8) );
plot(freq,amp);
title(' Low Pass Filter');
%axis([200e6 250e6 -50 0]);
grid
%%%pause
i_data=filter(b,a,i_data);
q_data=filter(b,a,q_data);
for i=1:4
step=500;
plot(i_data(100+(i-1)*step:100+i*step));
hold on;
plot(q_data(100+(i-1)*step:100+i*step),'r');
hold off;
grid;
title(' Demod I Data');
axis([1 500 -.1 .1]);
pause;
end
sample_i_data=reshape(i_data(1+500:4000+500),20,200);
sample_q_data=reshape(q_data(1+500:4000+500),20,200);
for ii=1:20
i_out=sample_i_data(ii,:);
q_out=sample_q_data(ii,:);
iq=i_out+1i*q_out;
plot(iq,'*');
grid;
%pause
end
end %phase
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -