AdalmPluto SDR tx rtl sdr rx

from commpy import modulation, filters

from scipy import signal

import numpy as np

import matplotlib.pyplot as plt

from pylab import *

from bindings.python import iio

from rtlsdr import RtlSdr

sdr = RtlSdr()

sdr.sample_rate = 3.0e6

sdr.center_freq = 1.2e9

sdr.gain = 4

# User configurable

TXLO = 1200000000

TXBW = 5000000

TXFS = 3000000

RXLO = TXLO

RXBW = TXBW

RXFS = TXFS

ctx = iio.scan_contexts();

print(ctx)

list = []

for i,j in ctx.items():

    list.append(i)

ctx = iio.Context(_context=list[0])

ctrl = ctx.find_device("ad9361-phy")

txdac = ctx.find_device("cf-ad9361-dds-core-lpc")

rxadc = ctx.find_device("cf-ad9361-lpc")

# Configure transceiver settings

ctrl.channels[0].attrs["frequency"].value = str(int(RXLO))

ctrl.channels[1].attrs["frequency"].value = str(int(TXLO))

ctrl.channels[4].attrs["rf_bandwidth"].value = str(int(RXBW))

ctrl.channels[5].attrs["rf_bandwidth"].value = str(int(TXBW))

ctrl.channels[4].attrs["sampling_frequency"].value = str(int(RXFS))

ctrl.channels[5].attrs["sampling_frequency"].value = str(int(TXFS))

ctrl.channels[5].attrs['gain_control_mode'].value = 'slow_attack'

#ctrl.channels[5].attrs['hardwaregain'].value = '20'

# Enable all IQ channels

rxadc.channels[0].enabled = True

rxadc.channels[1].enabled = True

txdac.channels[4].enabled = True

txdac.channels[5].enabled = True

# Force DAC to use DMA not DDSs

txdac.channels[0].attrs['raw'].value = str(0)

txdac.channels[1].attrs['raw'].value = str(0)

txdac.channels[2].attrs['raw'].value = str(0)

txdac.channels[3].attrs['raw'].value = str(0)

# Create cyclic buffer for TX data

N = 2**15

txbuf = iio.Buffer(txdac, N, False)

# QPSK Modulation

psk = modulation.PSKModem(4)

input_bits = [1, 1, 0, 0, 1, 0, 0, 1, 1, 1,]

psk_val = psk.modulate(input_bits=input_bits)

x = np.linspace(-4,4,2000)

#Filtering the modulated input bits with the sinc function

conv_re = []

conv_im = []

for sym in psk_val:

    conv_re.append(np.convolve(np.sinc(x),real(sym)))

    conv_im.append(np.convolve(np.sinc(x),imag(sym)))

conv_re = 1000000*np.reshape(conv_re,(1,size(conv_re)))

conv_im = 1000000*np.reshape(conv_im,(1,size(conv_im)))

print(psk_val)

"""plt.figure(5)

plt.plot(conv_im[0])

plt.figure(6)

plt.plot(conv_re[0])"""

iq = np.empty((conv_re.size + conv_im.size,), dtype=conv_im.dtype)

iq[0::2] = conv_re

iq[1::2] = conv_im

iq = np.int16(iq)

# Send data to buffer

txbuf.write(iq)

txbuf.push()

samples = sdr.read_samples()

plt.figure(13)

plt.plot(real(samples))

plt.figure(19)

plt.plot(imag(samples))

plt.show()