您好,登錄后才能下訂單哦!
這篇文章將為大家詳細(xì)講解有關(guān)python怎么實(shí)現(xiàn)二維插值的三維顯示,小編覺得挺實(shí)用的,因此分享給大家做個(gè)參考,希望大家閱讀完這篇文章后可以有所收獲。
具體內(nèi)容如下
# -*- coding: utf-8 -*- """ 演示二維插值。 """ # -*- coding: utf-8 -*- import numpy as np from mpl_toolkits.mplot3d import Axes3D import matplotlib as mpl from scipy import interpolate import matplotlib.cm as cm import matplotlib.pyplot as plt def func(x, y): return (x + y) * np.exp(-5.0 * (x ** 2 + y ** 2)) # X-Y軸分為20*20的網(wǎng)格 x = np.linspace(-1, 1, 20) y = np.linspace(-1, 1, 20) x, y = np.meshgrid(x, y) # 20*20的網(wǎng)格數(shù)據(jù) fvals = func(x, y) # 計(jì)算每個(gè)網(wǎng)格點(diǎn)上的函數(shù)值 15*15的值 fig = plt.figure(figsize=(9, 6)) #設(shè)置圖的大小 # Draw sub-graph2 ax = plt.subplot(1, 2, 1, projection='3d') #設(shè)置圖的位置 surf = ax.plot_surface(x, y, fvals, rstride=2, cstride=2, cmap=cm.coolwarm, linewidth=0.5, antialiased=True) #第四個(gè)第五個(gè)參數(shù)表示隔多少個(gè)取樣點(diǎn)畫一個(gè)小面,第六個(gè)表示畫圖類型,第七個(gè)是畫圖的線寬,第八個(gè)表示抗鋸齒 ax.set_xlabel('x') ax.set_ylabel('y') ax.set_zlabel('f(x, y)') #標(biāo)簽 plt.colorbar(surf, shrink=0.5, aspect=5) # 標(biāo)注 # 二維插值 newfunc = interpolate.interp2d(x, y, fvals, kind='cubic') # newfunc為一個(gè)函數(shù) # 計(jì)算100*100的網(wǎng)格上的插值 xnew = np.linspace(-1, 1, 100) # x ynew = np.linspace(-1, 1, 100) # y fnew = newfunc(xnew, ynew) # 僅僅是y值 100*100的值 np.shape(fnew) is 100*100 xnew, ynew = np.meshgrid(xnew, ynew) ax2 = plt.subplot(1, 2, 2, projection='3d') surf2 = ax2.plot_surface(xnew, ynew, fnew, rstride=2, cstride=2, cmap=cm.coolwarm, linewidth=0.5, antialiased=True) ax2.set_xlabel('xnew') ax2.set_ylabel('ynew') ax2.set_zlabel('fnew(x, y)') plt.colorbar(surf2, shrink=0.5, aspect=5) # 標(biāo)注 plt.show()
關(guān)于“python怎么實(shí)現(xiàn)二維插值的三維顯示”這篇文章就分享到這里了,希望以上內(nèi)容可以對(duì)大家有一定的幫助,使各位可以學(xué)到更多知識(shí),如果覺得文章不錯(cuò),請(qǐng)把它分享出去讓更多的人看到。
免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如果涉及侵權(quán)請(qǐng)聯(lián)系站長(zhǎng)郵箱:is@yisu.com進(jìn)行舉報(bào),并提供相關(guān)證據(jù),一經(jīng)查實(shí),將立刻刪除涉嫌侵權(quán)內(nèi)容。