import?numpy?as?np
import?matplotlib.pyplot?as?plt
from?matplotlib.pyplot?import?MultipleLocator

def?main（）:
????f0_List?=?np.arange（100?1001?50）
????#print（f0）
????fs_List?=?np.arange（1000?99?-50）
????#print（fs）
????t?=?np.arange（0?0.0101?0.00001）
????#print（t）
????while?1:
????????a?=?input（“鍵入?1?觀察固定頻率不同采樣率時還原曲線\n\
鍵入?2?觀察固定采樣率不同頻率時還原曲線\n鍵入?quit?退出\n“）
????????if?a?==?“1“:
????????????Change_fs（f0_List[2]?fs_List?t）
????????elif?a?==?“2“:
????????????Change_f0（f0_List?fs_List[0]?t）
????????elif?a?==?“quit“:
????????????break
????????else:
????????????print（“無法識別輸入QAQ“）

def?Change_f0（f0_List?fs?t）:
????Ts?=?1?/?fs
????Points_x?=?np.arange（-1?1.001?Ts）#抽樣點橫坐標x
????plt.ion（）#開啟交互模式
????plt.figure（figsize=（14?6））
????for?f0?in?f0_List:
????????#設(shè)置畫布外觀
????????plt.cla（）
????????plt.grid（True）
????????plt.xlim（0?0.01）
????????plt.ylim（-1?1）
????????plt.ylabel（“Amplitude“）
????????plt.title（“fo?=?%ffs“%（f0/fs））
????????##設(shè)置橫坐標間距
????????x_major_locator=MultipleLocator（0.001）
????????ax=plt.gca（）
????????ax.xaxis.set_major_locator（x_major_locator）
????????#設(shè)置圖像線條
????????F?=?np.cos（2*np.pi*f0*t）#原函數(shù)
????????Points_y?=?np.cos（2*np.pi*f0*Points_x）#抽樣點縱坐標y
????????fa?=?[]
????????for?tt?in?t:
????????????fa.append（np.dot（Points_y?np.sinc（fs*（tt-Points_x））））#抽樣后還原函數(shù)
????????#畫圖
????????plt.plot（t?F?“r“?label=“Orignal?Wave“）
????????plt.plot（Points_x?Points_y?“*g“?label=“Sampling?Points“）
??????