如何在R下进行波德图的相位部分

我想在R下创建一个名为“syst”的“siso”的经典bode图，因此我使用“control”包中的“freqresp”函数，然后继续：

# FREQUENCY LOG SEQUENCE
w = exp(seq(log(`^`(10,-2)), log(`^`(10,2)), length.out = `^`(10,3)))

H <- freqresp(syst, w)
mag = 20*log10(abs(H))
phase <- atan2(Im(H), Re(H)) * 180/pi

#频率日志序列
w=exp（seq（log（`^`（10，-2）），log（`^`（10,2）），length.out=`^`（10,3）））
H