下面主要看106k、212k、424k、848k速率下基带信号的波形及其频谱分析
import numpy as np
import matplotlib.pyplot as plt
import scipy.signal as signal
import random
fs=13560000
bit_x=[1]*64+[-1]*32+[1]*32
bit_y=[1]*128
bit_z=[-1]*32+[1]*96
CEA_SOF=[1]*512+bit_z
cea_frame_length=16
cea_frame_bit=[]
cea_frame_pattern=[]
for i in range(cea_frame_length):
cea_frame_bit.append(random.choice([0,1]))
cea_frame_pattern=cea_frame_pattern+CEA_SOF
pre_bit0=1
for i in range(cea_frame_length):
if cea_frame_bit[i]==0 and pre_bit0==1:
pre_bit0=1
cea_frame_pattern=cea_frame_pattern+bit_z
elif cea_frame_bit[i]==0 and pre_bit0==0:
pre_bit0=1
cea_frame_pattern=cea_frame_pattern+bit_y
else:
pre_bit0=0
cea_frame_pattern=cea_frame_pattern+bit_x
if cea_frame_bit[-1]==0: #EOF
cea_frame_pattern=cea_frame_pattern+bit_z+bit_y
else:
cea_frame_pattern=cea_frame_pattern+bit_y+bit_y
#106k analysis
# FFT analysis
fft_result=np.fft.fft(np.array(cea_frame_pattern))
freqs=np.fft.fftfreq(len(fft_result),1/fs)
length_fft=len(fft_result)//2
# find 95% energy frequency distribution
energy_threshold = 0.95 * np.sum(np.abs(fft_result[1:length_fft])**2)
cumulative_energy = 0
freq_range = []
for i in range(len(fft_result)//2):
if i==0:
continue
if cumulative_energy < energy_threshold:
freq_range.append(freqs[i])
cumulative_energy += np.abs(fft_result[i])**2
plt.figure(figsize=(12,6))
plt.subplot(2,1,1)
plt.plot(cea_frame_pattern,marker='.')
plt.title("CEA 106k signal")
plt.subplot(2,1,2)
plt.plot(freqs[:length_fft],np.abs(fft_result)[:length_fft])
plt.axvline(freq_range[0], color='r', linestyle='--', label='95% Energy Range')
plt.axvline(freq_range[-1], color='r', linestyle='--')
plt.title("Fequency Spectrum")
plt.xlabel("Fequency(Hz))")
plt.tight_layout()
#212k analysis
ceb_frame_pattern_dwon_sample=[]
for i in range(len(cea_frame_pattern)):
if i%2==0:
ceb_frame_pattern_dwon_sample.append(cea_frame_pattern[i])
# FFT analysis
fft_result=np.fft.fft(np.array(ceb_frame_pattern_dwon_sample))
freqs=np.fft.fftfreq(len(fft_result),1/fs)
length_fft=len(fft_result)//2
# find 95% energy frequency distribution
energy_threshold = 0.95 * np.sum(np.abs(fft_result[1:length_fft])**2)
cumulative_energy = 0
freq_range = []
for i in range(len(fft_result)//2):
if i==0:
continue
if cumulative_energy < energy_threshold:
freq_range.append(freqs[i])
cumulative_energy += np.abs(fft_result[i])**2
plt.figure(figsize=(12,6))
plt.subplot(2,1,1)
plt.plot(ceb_frame_pattern_dwon_sample,marker='.')
plt.title("CEA 212k signal")
plt.subplot(2,1,2)
plt.plot(freqs[:length_fft],np.abs(fft_result)[:length_fft])
plt.axvline(freq_range[0], color='r', linestyle='--', label='95% Energy Range')
plt.axvline(freq_range[-1], color='r', linestyle='--')
plt.title("Fequency Spectrum")
plt.xlabel("Fequency(Hz))")
plt.tight_layout()
#424k analysis
ceb_frame_pattern_dwon_sample=[]
for i in range(len(cea_frame_pattern)):
if i%4==0:
ceb_frame_pattern_dwon_sample.append(cea_frame_pattern[i])
# FFT analysis
fft_result=np.fft.fft(np.array(ceb_frame_pattern_dwon_sample))
freqs=np.fft.fftfreq(len(fft_result),1/fs)
length_fft=len(fft_result)//2
# find 95% energy frequency distribution
energy_threshold = 0.95 * np.sum(np.abs(fft_result[1:length_fft])**2)
cumulative_energy = 0
freq_range = []
for i in range(len(fft_result)//2):
if i==0:
continue
if cumulative_energy < energy_threshold:
freq_range.append(freqs[i])
cumulative_energy += np.abs(fft_result[i])**2
plt.figure(figsize=(12,6))
plt.subplot(2,1,1)
plt.plot(ceb_frame_pattern_dwon_sample,marker='.')
plt.title("CEA 424k signal")
plt.subplot(2,1,2)
plt.plot(freqs[:length_fft],np.abs(fft_result)[:length_fft])
plt.axvline(freq_range[0], color='r', linestyle='--', label='95% Energy Range')
plt.axvline(freq_range[-1], color='r', linestyle='--')
plt.title("Fequency Spectrum")
plt.xlabel("Fequency(Hz))")
plt.tight_layout()
#848k analysis
ceb_frame_pattern_dwon_sample=[]
for i in range(len(cea_frame_pattern)):
if i%4==0:
ceb_frame_pattern_dwon_sample.append(cea_frame_pattern[i])
# FFT analysis
fft_result=np.fft.fft(np.array(ceb_frame_pattern_dwon_sample))
freqs=np.fft.fftfreq(len(fft_result),1/fs)
length_fft=len(fft_result)//2
# find 95% energy frequency distribution
energy_threshold = 0.95 * np.sum(np.abs(fft_result[1:length_fft])**2)
cumulative_energy = 0
freq_range = []
for i in range(len(fft_result)//2):
if i==0:
continue
if cumulative_energy < energy_threshold:
freq_range.append(freqs[i])
cumulative_energy += np.abs(fft_result[i])**2
plt.figure(figsize=(12,6))
plt.subplot(2,1,1)
plt.plot(ceb_frame_pattern_dwon_sample,marker='.')
plt.title("CEA 848k signal")
plt.subplot(2,1,2)
plt.plot(freqs[:length_fft],np.abs(fft_result)[:length_fft])
plt.axvline(freq_range[0], color='r', linestyle='--', label='95% Energy Range')
plt.axvline(freq_range[-1], color='r', linestyle='--')
plt.title("Fequency Spectrum")
plt.xlabel("Fequency(Hz))")
plt.show()