Scribd Logo
Upload

Welcome to Scribd!

  • 1 views

    DSB SC

    Uploaded by

    rbtemp1149
    NA

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    DSB SC

    Uploaded by

    rbtemp1149
    1 views4 pages
    NA

    Original Title

    dsb-sc

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    NA

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    1 views4 pages

    DSB SC

    Uploaded by

    rbtemp1149
    NA

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 4

    dsb-sc

    October 28, 2024

    [1]: import numpy as np


    import matplotlib.pyplot as plt
    from scipy import signal

    [61]: #setting parameters


    Ac=1
    Am=1
    fm=10
    fc=100
    fs=1000
    t=np.arange(0,1,1/fs)

    [62]: #message signal


    message=Am*(np.cos(2*np.pi*fm*t))
    #carrier signal
    carrier=Ac*(np.cos(2*np.pi*fc*t))
    #frequency axis
    freq=np.fft.fftfreq(len(t),1/fs)
    #message fourier transform
    fft_message=np.fft.fft(message)
    #carrier fourier transform
    fft_carrier=np.fft.fft(carrier)

    #Because of symmetry taking only positive spectrum


    fft_positive_message=fft_message[:int(len(fft_message)/2)]/max(fft_message)␣
    ↪#normalizing the amplitude

    fft_positive_carrier=fft_carrier[:int(len(fft_carrier)/2)]/max(fft_carrier)␣
    ↪#normalizing the amplitude

    positive_freq=freq[:int(len(freq)/2)]

    [63]: #plotting message and carrier


    plt.figure(figsize=(12,6))
    plt.subplot(2,2,1)
    plt.plot(t,message,color="red")
    plt.xlabel("Time")
    plt.ylabel("Amplitude")
    plt.title("Message Signal")

    1
    plt.subplot(2,2,2)
    plt.plot(t,carrier,color="blue")
    plt.xlabel("Time")
    plt.ylabel("Amplitude")
    plt.title("Carrier Signal")
    plt.subplot(2,2,3)
    plt.plot(positive_freq,np.abs(fft_positive_message),color="red")
    plt.xlabel("Frequency")
    plt.ylabel("Amplitude")
    plt.title("Fourier Transform of Message Signal")
    plt.subplot(2,2,4)
    plt.plot(positive_freq,np.abs(fft_positive_carrier),color="blue")
    plt.xlabel("Frequency")
    plt.ylabel("Amplitude")
    plt.title("Fourier Transform of Carrier Signal")
    plt.tight_layout()
    plt.show()

    [98]: #dsc-sc modulation


    dsb_sc_mod=message*carrier
    fft_dsb_sc=np.fft.fft(dsb_sc_mod)

    #Because of symmetry taking only positive spectrum


    fft_positive_dsb_sc=fft_dsb_sc[:int(len(fft_dsb_sc)/2)]/max(fft_dsb_sc)

    plt.figure(figsize=(12,5))
    plt.subplot(2,1,1)
    plt.plot(t,dsb_sc_mod,color="purple")

    2
    plt.xlabel("Time")
    plt.ylabel("Amplitude")
    plt.title("DSB-SC Modulated Signal")
    plt.subplot(2,1,2)
    plt.plot(positive_freq,np.abs(fft_positive_dsb_sc),color="brown")
    plt.xlabel("Frequency")
    plt.ylabel("Amplitude")
    plt.title("Fourier Transform of DSB-SC Modulated Signal")
    plt.tight_layout()
    plt.show()

    [91]: def demodulation(mod_signal,cutt_off,fs,n,type):


    b,a=signal.butter(n,cutt_off,btype=type,fs=fs)
    w, h = signal.freqs(b, a)
    demod_signal=signal.filtfilt(b,a,mod_signal)
    return w,h,demod_signal

    [99]: w,h,demod_dsb_sc=demodulation(dsb_sc_mod*carrier,fm,fs,7,"low")
    plt.figure(figsize=(10,10))
    plt.subplot(2,1,1)
    plt.plot(t,demod_dsb_sc,color="red")
    plt.xlabel("Time")
    plt.ylabel("Amplitude")
    plt.title("Demodulated signal")

    [99]: Text(0.5, 1.0, 'Demodulated signal')

    3
    4

    You might also like