1

我需要绘制以下函数，也就是我所说的理论数据

在此处输入图像描述

对于以下情况，a=z-(L/2) b=z+(L/2)

  • N=75
  • 我=1
  • L=0.16
  • R=0.013

z 在 -0.011 和 0.011 之间。我的猜测是存在下溢问题。代码如下

import numpy as np
import matplotlib.pyplot as plt

class Solenoid():

def __init__(self,N,I,L,R):
    self.N=N
    self.I=I
    self.L=L
    self.R=R

def field(self,z):

    mu=4*(np.pi)*(10 ** -7)

    coef=(mu*self.N*self.I)/(2*self.L)

    a=z+(self.L/2)
    b=z-(self.L/2)

    sum_a=a/((self.R ** 2 + a ** 2) ** 0.5)

    sum_b=b/((self.R ** 2 + b ** 2) ** 0.5)

    return (coef*(sum_a - sum_b))

sol=Solenoid(75,1,0.160,0.013)
x=np.linspace(-0.011, 0.011, 10000)
y=[]

for point in x:
    y.append(sol.field(point))

plt.plot(x,y)
plt.show()

输出是

在此处输入图像描述

这是完全错误的。预期的输出与实验数据一起或多或少是一条平滑的线：

在此处输入图像描述

更新新输出是

在此处输入图像描述

4

1 回答 1

0

这是 x 比例的 10 倍...

在此处输入图像描述

于 2014-04-01T21:11:13.607 回答