import marimo

__generated_with = "0.25.1"
app = marimo.App()


@app.cell
def _():
    import marimo as mo

    return (mo,)


@app.cell(hide_code=True)
def _(mo):
    mo.md(r"""
    # Digital to Analog conversion using sinc
    """)
    return


@app.cell
def _():
    import numpy as np
    import matplotlib.pyplot as plt
    def reconstruct_with_sinc(ts,fd,t):
        n, = ts.shape
        dt = ts[1] - ts[0]
        fr = []
        for k,ti in enumerate(t):
            # for each time point
            sumf = 0.0
            for i in range(n):
                # for each point in a sampled set
                sumf += fd[i]*np.sin(np.pi*(ti/dt-i))/(ti/dt-i)
            
        
            fr.append((1./np.pi)*sumf)
        
        return np.asarray(fr,dtype='f')

    return np, plt, reconstruct_with_sinc


@app.cell
def _(np, plt, reconstruct_with_sinc):
    # if __name__ == '__main__':
    t = np.arange(0.0,0.6,0.001)
    fa = 1.0*np.sin(2*np.pi*10*t)+0.2*np.sin(2*np.pi*6*t)
    fs = 10 # Hz
    ts = np.arange(0.0,0.6,1./fs) # sampling time
    fd = 1.0*np.sin(2*np.pi*10*ts)+0.2*np.sin(2*np.pi*6*ts) # sampled data


    plt.figure()
    fr = reconstruct_with_sinc(ts,fd,t)
    plt.plot(t,fa,'b-',ts,fd,'ro',t,fr,'g--')
    plt.xlabel('t [sec]')
    plt.ylabel('y [V]')
    plt.legend(('Original','Sampled','Reconstructed'))
    plt.show()
    return


if __name__ == "__main__":
    app.run()