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

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()