matplotlib 绘制天线辐射函数的极坐标图时出现问题

b5lpy0ml 于 2023-03-03 发布在其他

关注(0)|答案(1)|浏览(177)

我在画这个辐射图

下式：

我可以做一部分：

%matplotlib inline
import matplotlib.ticker
import matplotlib.pyplot as plt
plt.style.use('seaborn-white')
import numpy as np
import pylab as pl

def gain_dip(theta, phi):
    return 1.641*(np.cos(np.pi/2*np.cos(theta))/np.sin(theta))**2

theta = np.arange(-np.pi, np.pi,0.01)
# plot
ax = plt.subplot(111, polar=True)
# set zero west
ax.set_theta_zero_location('W')
ax.set_theta_direction('clockwise')
# let set an azimuth for example, pi
plt.plot(theta, gain_dip(theta, np.pi))

图中画出了仰角（θ），与示例中的仰角相似。但是，它有一条不应该出现的180度角的线。

我也想画出方位角phi，但是当我尝试时，它给出了误差：

phi = np.arange(-np.pi, np.pi,0.01)
ax = plt.subplot(111, polar=True)
# set zero west
ax.set_theta_zero_location('W')
ax.set_theta_direction('clockwise')
plt.plot(phi, gain_dip(np.pi/2, phi))

ValueError: x and y must have same first dimension, but have shapes (629,) and (1,)
而且，有可能像第一个图一样绘制两个Angular 吗？

matplotlib

来源：https://stackoverflow.com/questions/58612132/problem-doing-polar-plot-of-a-antenna-radiation-function

1条答案

按热度按时间

vxbzzdmp1#

考虑有2个numpy向量pattern = [patternH，patternV]，例如.MSI天线数据或函数点。

pattern = [patternH,patternV]
NameFig = "Plot.jpg"
Titolo = "First Picture"

M_PI_O_180 = .01745329251994329576
r = np.arange(0, 360, 1)
theta = r*M_PI_O_180
fig, axs = plt.subplots(1,2,figsize=(11, 6),subplot_kw={"projection" : "polar"})
fig.suptitle('Modello '+Titolo, fontsize=12)

ax = axs[0]
ax.plot(theta,pattern[0],linewidth= 2.5)
ax.set_rmax(0)
ax.set_theta_direction(-1)
ax.set_thetagrids([0,30, 60,90,120,150,180,210,240,270,300,330])
ax.set_rticks([-3, -10,-20],labels= ["-3dB", "-10dB","-20dB"])
ax.set_rlabel_position(270)
ax.grid(True,which="minor",linestyle= ":")
ax.grid(True,which="major",linewidth= 1.2)
ax.minorticks_on()
ax.set_title("Diagramma Orizzontale", {"va":"bottom"})

axs[1].plot(theta,pattern[1],linewidth= 2.5)
axs[1].set_rmax(0)
axs[1].set_theta_direction(-1)
axs[1].set_thetagrids([0,30, 60,90,120,150,180,210,240,270,300,330])
axs[1].set_rticks([-3, -10,-20,-30,-40],labels= ["-3dB", "-10dB","-20dB","-30dB","-40dB"])
axs[1].set_rlabel_position(270)
axs[1].grid(True,which="minor",linestyle= ":")
axs[1].grid(True,which="major",linewidth= 1.2)
axs[1].minorticks_on()
axs[1].set_title("Diagramma Verticale", {"va":"bottom"})
plt.savefig(NameFig)

赞(0）回复(0）举报 2023-03-03

我来回答

matplotlib 绘制天线辐射函数的极坐标图时出现问题

1条答案

相关问题

热门标签

最新问答