python计算分段函数值的方法
作者:执念斩长河 发布时间:2022-10-08 20:23:34
标签: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]]
绘制图形
总结
分段函数值配合定义域是这段代码的亮点,以后写代码时更应该多多揣摩。
来源:https://blog.csdn.net/m0_37149062/article/details/120038188


猜你喜欢
- 1.关系模型:用二维表格结构表示实体集,外键表示实体间联系的数据模型称为关系模型。关系模型是由若干个关系模式组成的集合。2.关系模式:关系模
- python永久添加搜索路径_Python sys.path永久添加在用户目录下,找到隐藏文件.bashrc 文件然后在末尾添加export
- XML虽然比JSON复杂,在Web中应用也不如以前多了,不过仍有很多地方在用,所以,有必要了解如何操作XML。DOM vs SAX操作XML
- 以下笔记是我在 xue.cn 学习群之数据分析小组所整理分享的心得。相关背景是:我选择中文词频统计案例作为考察大家python基础功掌握程度
- 内容摘要:网页设计师制作网页最常用的设计软件应该就算adobe的产品Photoshop了,当然Photoshop不仅可以设计网页,不过作为网
- YAML语法规则:http://www.ibm.com/developerworks/cn/xml/x-cn-yamlintro/下载PyY
- Dreamweaver出现乱码,大致为两种情况:一是没有标明主页制作所用的文字,这种情况下很简单就可以
- python 网络编程详解网络编程的专利权应该属于Unix,各个平台(如windows、Linux等)、各门语言(C、C++、Python、
- 当你有多个进程或线程访问相同的数据时,竞争条件是一个威胁。本文探讨了在发现竞争条件后如何测试它们。Incrmnt你在一个名为“Incrmnt
- 协程协程(co-routine,又称微线程)是一种多方协同的工作方式。当前执行者在某个时刻主动让出(yield)控制流,并记住自身当前的状态
- 1.选择File -> Settings2.选择 File and Code Templates -> Files ->
- 如下所示:import urllib.requestimport urllib.parseurl = 'https://weibo.
- CSV文件简单来说CSV文件由任意数目的记录组成,记录间以某种换行符分隔;每条记录由字段组成,字段间的分隔符是其它字符或字符串,最常见的是逗
- 在深度学习中,如果我们想获得某一个层上的feature map,就像下面的图这样,怎么做呢?我们的代码是使用keras写的VGG16网络,网
- 参与测试的浏览器:IE6 / IE7 / IE8 / FF3 / OP9.6 / SF3 / Chrome2操作系统:Windows蓝色理想
- 在网络上看到的数字人整合动网论坛的方法都非常不全,站长们都是抄人家的,也不说明可不可用,提供下载的文件也不能下载.现在我提供一些信息。一、整
- 这篇文章主要介绍了python错误调试及单元文档测试过程解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需
- 组件 (Component) 是 Vue.js 最强大的功能之一。组件可以扩展 HTML 元素,封装可重用的代码。在较高层面上,组件是自定义
- 本文实例为大家分享了Python密码强弱判断的具体代码,供大家参考,具体内容如下程序说明:通过获取用户输入,判断密码长度是否大于8,同时判断
- PHP计算字符串用strlen()只能得到字符串长度,不是宽高像素,使用到了php函数ImageTTFBBox(),就可以根据字体的大小和所