python計算分段函數(shù)值的方法
發(fā)布日期：2022-01-29 18:26 | 文章來源：gibhub
本博文源于python科學(xué)計算，旨在解析如何使用python進(jìn)行計算分段函數(shù)值。下面就以復(fù)雜的二元函數(shù)進(jìn)行演示。
題目再現(xiàn)
求解思路
首先先生成x1,x2，然后我們將其映射到網(wǎng)格里，最后寫出表達(dá)式，進(jìn)行求解。題目效果，再來張圖片，就更加好了！
代碼再現(xiàn)
# -*- coding:UTF-8 -*-
# 開發(fā)時間：2021/9/1 13:27
import numpy as np
import matplotlib.pyplot as plt
x = np.arange(-2,2,0.04)
y = np.arange(-2,2,0.04)
# 根據(jù)x和y生成網(wǎng)格點(diǎn)
[X,Y] = np.meshgrid(x,y)
Pxy = 0.5457*np.exp(-0.75*Y**2 - 3.75*X**2 - 1.5*X) * (X+Y>1) + \
0.7575*np.exp(-Y**2-6*X**2)*((X+Y<=1) & (X+Y>-1)) + \
0.5457*np.exp(-0.75*Y**2 - 3.75*X**2 + 1.5*X)*(X+Y<=-1)
# 計算分段函數(shù)值
print(Pxy)
# 繪制出結(jié)果
fig = plt.figure()
ax = plt.gca(projection="3d")
ax.plot_surface(X,Y,Pxy,cmap="rainbow")
ax.set_xlabel("X")
ax.set_ylabel("Y")
ax.set_zlabel("Z")
plt.show()
計算函數(shù)值結(jié)果
[[4.13780323e-10 7.95789904e-10 1.51222176e-09 ... 8.55819561e-12 3.43799817e-12 1.35484814e-12] [4.65976503e-10 8.96174556e-10 1.70298047e-09 ... 1.00270794e-11 4.02807814e-12 1.58738717e-12] [5.23499048e-10 1.00680297e-09 1.91320516e-09 ... 1.17105410e-11 4.70435829e-12 1.85389601e-12] ... [8.34343054e-13 2.15823201e-12 5.47662686e-12 ... 4.02602659e-09 2.14422875e-09 1.12837657e-09] [7.16690905e-13 1.85389601e-12 4.70435829e-12 ... 3.59225426e-09 1.91320516e-09 1.00680297e-09] [6.13662224e-13 1.58738717e-12 4.02807814e-12 ... 3.19753414e-09 1.70298047e-09 8.96174556e-10]]
繪制圖形
總結(jié)
分段函數(shù)值配合定義域是這段代碼的亮點(diǎn)，以后寫代碼時更應(yīng)該多多揣摩。
以上就是本文的全部內(nèi)容，希望對大家的學(xué)習(xí)有所幫助，也希望大家多多支持本站。
版權(quán)聲明：本站文章來源標(biāo)注為YINGSOO的內(nèi)容版權(quán)均為本站所有，歡迎引用、轉(zhuǎn)載，請保持原文完整并注明來源及原文鏈接。禁止復(fù)制或仿造本網(wǎng)站，禁止在非www.sddonglingsh.com所屬的服務(wù)器上建立鏡像，否則將依法追究法律責(zé)任。本站部分內(nèi)容來源于網(wǎng)友推薦、互聯(lián)網(wǎng)收集整理而來，僅供學(xué)習(xí)參考，不代表本站立場，如有內(nèi)容涉嫌侵權(quán)，請聯(lián)系alex-e#qq.com處理。
相關(guān)文章
關(guān)注官方微信