курсовая / выполнение / ДОФМ

import numpy as np
import matplotlib.pyplot as plt
import math

s=[1,1,0,1,0,1,1,0,0]
mf=math.pi/2
ti=2*10**(-3)
u0=3
w=10**4

def OFM(Sa,m,T,U):
t = [0]
Sig = [0]
t[0]=np.linspace(0,T,num=100)
if Sa[0]==1:
Sig[0]=U*np.cos(w*t[0]+m)
elif Sa[0]==0:
Sig[0]=U*np.cos(w*t[0]-m)
for i in range(1,len(Sa)-1):
t.append(np.linspace(i*T,(i+1)*T,100))
if Sa[i]==Sa[i-1]:
Sig.append(U*np.cos(w*t[i])+m)
elif Sa[i]!=Sa[i-1]:
Sig.append(U*np.cos(w*t[i])-m)
Sig=np.array(Sig)
t=np.array(t)
return np.ravel(Sig),np.ravel(t)

Signal,time = OFM(s,mf,ti,u0)

plt.plot(time,Signal)
plt.show()