快捷導(dǎo)航

利用Python?Matlab繪制曲線圖的簡(jiǎn)單實(shí)例

更新時(shí)間：2021年12月16日 11:17:55 作者：hhh江月

們經(jīng)常會(huì)遇到這種情況,有一個(gè)數(shù)學(xué)函數(shù),我們希望了解他的圖像,這個(gè)時(shí)候使用python 的matplotlib就可以幫助我們,下面這篇文章主要介紹了利用Python?Matlab繪制曲線圖的相關(guān)資料,需要的朋友可以參考下

一、簡(jiǎn)介

我們?cè)谶@里采用Python中的matplotlib來(lái)實(shí)現(xiàn)曲線圖形的繪制。matplotlib是著名的python繪圖庫(kù)，它提供了一整套繪圖API，十分適合交互式繪圖。

二、繪制圖形

1、第一個(gè)曲線圖

代碼：

具體的繪制的代碼如下所示：

import matplotlib.pyplot as plt
import numpy as np

r = np.array([2072.54, 2076.84, 2085.51, 2103.01, 2129.93, 2162.16, 2200.22, 2242.15,
              2285.71, 2328.29, 2350.18, 2364.01, 2364.01, 2343.29, 2300.17, 2252.25,
              2208.72, 2166.85, 2132.19, 2103.01, 2085.51, 2075.77, 2072.54])
b_ = np.array([30.159, 27.143, 24.127, 21.111, 18.096, 15.080, 12.064, 9.048,
              6.032, 3.016, 1.508, 0, -1.508, -3.016, -6.032, -9.048, -12.064,
              -15.080, -18.096, -21.111, -24.127, -27.143, -30.159])
b = b_ * pow(10, -4)


plt.plot(b, r)
plt.xlabel("B/T")
plt.ylabel("R/Ω")
plt.title("GMB R-B (decreasing B)")
plt.show()

效果：

2、第二個(gè)圖形

代碼：

代碼與上一個(gè)的代碼其實(shí)是比較相似的：

import matplotlib.pyplot as plt
import numpy as np

r = np.array([2072.53, 2076.81, 2085.47, 2103.00, 2129.90, 2162.11, 2200.20, 2242.06,
              2285.66, 2328.24, 2350.13, 2364.00, 2363.96, 2343.19, 2300.20, 2252.29,
              2208.76, 2166.89, 2132.20, 2103.05, 2085.50, 2075.81, 2072.56])
b_ = np.array([30.159, 27.143, 24.127, 21.111, 18.096, 15.080, 12.064, 9.048,
              6.032, 3.016, 1.508, 0, -1.508, -3.016, -6.032, -9.048, -12.064,
              -15.080, -18.096, -21.111, -24.127, -27.143, -30.159])
b = b_ * pow(10, -4)


plt.plot(b, r)
plt.xlabel("B/T")
plt.ylabel("R/Ω")
plt.title("GMB R-B (increasing B)")
plt.show()

效果：

3、第三個(gè)圖形

代碼：

代碼基本是形同的啦：

import matplotlib.pyplot as plt
import numpy as np

v = np.array([274, 270, 261, 243, 219, 189, 155, 118, 81, 48, 34, 21])
b_ = np.array([30.159, 27.143, 24.127, 21.111, 18.096, 15.080, 12.064, 9.048,
              6.032, 3.016, 1.508, 0])
b = b_ * pow(10, -4)


plt.plot(b, v)
plt.xlabel("B/T")
plt.ylabel("V/mV")
plt.title("GMB V-B")
plt.show()

效果：

4、第四個(gè)圖形

代碼：

代碼其實(shí)都是基本一樣的，只不過(guò)主要是更換了數(shù)據(jù)啦：

import matplotlib.pyplot as plt
import numpy as np

w = np.array([43.5, 44, 47, 50, 53, 56, 59, 62, 65, 68, 71, 74, 77, 80, 83, 86,
              89, 92, 95, 98, 101, 104])
v = np.array([0, 5.7, 35.0, 53.8, 45.9, 7.7, -45.7, -51.9, -32.6, -1.8, 34.5, 53.1,
              39.2, -10.1, -47.9, -51.4, -29.5, 5.6, 34.4, 52.4, 40.9, -5.2])


plt.plot(w, v)
plt.xlabel("θ/rad")
plt.ylabel("V/mV")
plt.title("GMB V-θ")
plt.show()

效果：

5.畫(huà)出指定區(qū)間的一個(gè)多項(xiàng)式函數(shù)：

import numpy as np

import matplotlib.pyplot as plt

X = np.linspace(-4, 4, 1024)

Y = .25 * (X + 4.) * (X + 1.) * (X - 2.)

plt.title('$f(x)=\\frac{1}{4}(x+4)(x+1)(x-2)$')

plt.plot(X, Y, c = 'g')

plt.show()