python计算分段函数值的方法

本博文源于python科学计算,旨在解析如何使用python进行计算分段函数值。下面就以复杂的二元函数进行演示。

题目再现

求解思路

首先先生成x1,x2,然后我们将其映射到网格里,最后写出表达式,进行求解。题目效果,再来张图片,就更加好了!

代码再现

# -*- coding:UTF-8 -*-
# 开发时间: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)
# 根据x和y生成网格点
[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)

# 计算分段函数值
print(Pxy)

# 绘制出结果
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()

计算函数值结果

[[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]]

绘制图形

总结

分段函数值配合定义域是这段代码的亮点,以后写代码时更应该多多揣摩。

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持我们。

(0)

相关推荐

  • python使用Matplotlib绘制分段函数

    本文实例为大家分享了python使用Matplotlib绘制分段函数的具体代码,供大家参考,具体内容如下 环境 Python3 Mac OS 代码 # coding:utf-8 """ Author: roguesir Date: 2017/8/30 GitHub: https://roguesir.github.com Blog: http://blog.csdn.net/roguesir """ import numpy as np impo

  • python计算分段函数值的方法

    本博文源于python科学计算,旨在解析如何使用python进行计算分段函数值.下面就以复杂的二元函数进行演示. 题目再现 求解思路 首先先生成x1,x2,然后我们将其映射到网格里,最后写出表达式,进行求解.题目效果,再来张图片,就更加好了! 代码再现 # -*- coding:UTF-8 -*- # 开发时间:2021/9/1 13:27 import numpy as np import matplotlib.pyplot as plt x = np.arange(-2,2,0.04) y

  • python计算方程式根的方法

    本文实例讲述了python计算方程式根的方法.分享给大家供大家参考.具体实现方法如下: ''' roots = polyRoots(a). Uses Laguerre's method to compute all the roots of a[0] + a[1]*x + a[2]*x^2 +...+ a[n]*x^n = 0. The roots are returned in the array 'roots', ''' from evalPoly import * from numpy i

  • python计算圆周率pi的方法

    本文实例讲述了python计算圆周率pi的方法.分享给大家供大家参考.具体如下: from sys import stdout scale = 10000 maxarr = 2800 arrinit = 2000 carry = 0 arr = [arrinit] * (maxarr + 1) for i in xrange(maxarr, 1, -14): total = 0 for j in xrange(i, 0, -1): total = (total * j) + (scale * a

  • Python计算字符宽度的方法

    本文实例讲述了Python计算字符宽度的方法.分享给大家供大家参考,具体如下: 最近在用python写一个CLI小程序,其中涉及到计算字符宽度,目标是以友好的方式将一个长字符串截取为等宽的片段. 对于unicode字符,python的len函数可以准确的计算其中所包含的字符个数,但是个数并不代表宽度,如: >>>len(u'你好a') 3 因此无法简单的使用这种方式来计算宽度. GBK decode 首先我想到GBK编码,00–7F范围内的字符是一字节编码,其余是双字节编码,正好与字符的

  • python计算阶乘和的方法(1!+2!+3!+...+n!)

    方法一:使用while循环来计算 n = int(input()) jie = 1 sum = 0 i = 1 while n >= i: jie = jie * i sum = sum + jie i = i + 1 print(sum) 方法二:使用递归函数调用阶乘方法求和(其中n的值在1~40之间) def jie(n): if n == 1: return 1 else: return n*jie(n-1) n = int(input()) sum = 0 if n < 1 or n

  • Python计算三角函数之asin()方法的使用

    asin()方法返回x的反正弦,以弧度表示. 语法 以下是asin()方法语法: asin(x) 注意:此函数是无法直接访问的,所以我们需要导入math模块,然后需要用math的静态对象来调用这个函数. 参数 x -- 这必须是在-1到1的数字值的范围内,如果x大于1,则它会产生一个错误. 返回值 此方法返回的x反正弦,以弧度表示. 例子 下面的例子显示了asin()方法的使用. #!/usr/bin/python import math print "asin(0.64) : ",

  • 用Python计算三角函数之acos()方法的使用

    acos()方法返回x的反余弦值,以弧度表示. 语法 以下是acos()方法的语法: acos(x) 注意:此函数是无法直接访问的,所以我们需要导入math模块,然后需要用math的静态对象来调用这个函数. 参数 x -- 这必须是在范围内的数字值-1到1,如果x大于1,则它会产生一个错误. 返回值 此方法返回的X反余弦,以弧度表示. 例子 下面的例子显示acos()方法的使用. #!/usr/bin/python import math print "acos(0.64) : ",

  • 用Python计算三角函数之atan()方法的使用

    atan()方法返回x的反正切值,以弧度表示. Syntax 以下是atan()方法的语法: atan(x) 注意:此函数是无法直接访问的,所以我们需要导入math模块,然后需要用math的静态对象来调用这个函数. 参数 x -- 这必须是一个数值. 返回值 此方法返回 x 的反正切值,以弧度表示. 例子 下面的例子显示atan()方法的使用. #!/usr/bin/python import math print "atan(0.64) : ", math.atan(0.64) pr

  • python计算时间差的方法

    本文实例讲述了python计算时间差的方法.分享给大家供大家参考.具体分析如下: 1.问题: 给定你两个日期,如何计算这两个日期之间间隔几天,几个星期,几个月,几年? 2.解决方法: 标准模块datetime和第三方包dateutil(特别是dateutil的rrule.count方法)能非常简单迅速的帮你解决这个问题. from dateutil import rrule import datetime def weeks_between(start_date, end_date): week

  • python计算文本文件行数的方法

    本文实例讲述了python计算文本文件行数的方法.分享给大家供大家参考.具体实现方法如下: filename = "somefile.txt" myfile = open(filename) lines = len(myfile.readlines()) print "There are %d lines in %s" % (lines, filename) 希望本文所述对大家的Python程序设计有所帮助.

随机推荐