Python3.0科学计算学习之绘图(二)

2023-12-27 07:38

本文主要是介绍Python3.0科学计算学习之绘图(二),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

(1) np.mashgrid()函数:-----生成网络点坐标矩阵,可以是二维网络矩阵,也可以是三维网络矩阵。其中,每个交叉点就是网络点,描述这些网络点的矩阵就是坐标矩阵(横坐标矩阵X中的每个元素与纵坐标矩阵Y中对应位置元素,共同构成一个点的完整坐标)。

背景示例:网络点与坐标矩阵的解释如下:

import numpy as npimport matplotlib.pyplot as pltx=np.array([[0,1,2],[0,1,2]])        #最简单的方法是,可以把横纵坐标矩阵X,Y写出来,生成坐标点y=np.array([[0,0,0],[1,1,1]])plt.plot(x,y,color='red',marker='.',markersize=10,linestyle='-.')      #点的形状为原点,点设置大一些,线性为电划线。plt.grid(True)plt.show()   

注意:按照矩阵给坐标点信息,matplotlib会把横坐标矩阵中,每一列的点当做同一条线。正如上例中把plot的linestyle=’ ‘改为linestyle=’-.’就会发现A-D,B-E,C-F是连接的。
在这里插入图片描述
对于很多网络点的情况,可用如下meshgrid()函数方法:

因注意到坐标矩阵其中有大量的重复—X的每一行都一样,Y的每一列都一样。故基于此规律性,numpy提供的np.meshgrid()函数可以快速生成坐标矩阵X,Y

x=np.linspace(-0.5,2.,10)    y=np.linspace(-1.5,4.,10)X,Y=np.meshgrid(x,y)     #输入的x,y是网络点的横纵坐标列向量(非矩阵),输出的X,Y就是坐标矩阵。plt.plot(X,Y,color='limegreen',marker='.',linestyle='')plt.grid(True)plt.show()

在这里插入图片描述
(2) Python可视化库matplotlib.pyplot里contour()与contourf()函数

区别:contour()和counterf() 函数功能相同,都是画三维等高线图的,不同点在于contourf会根据不同的水平值用不同的颜色来填充绘图(即会对等高线间的区域进行填充)。

contour()和counterf() 都是绘制三维图,其中前两个参数x和y为两个等长一维数组,第三个参数z为二维数组(表示平面点xi和yi映射的函数值)。

matplotlib等高线的绘制步骤:

等高线是三维图像在二维空间的投影

l 首先准备三维函数及待投影平面的网格坐标。

import numpy as npimport matplotlib.pyplot as pltdef height(x,y):return (1-x/2+x**5+y**3)*np.exp(-x**2-y**2)x=np.linspace(-3,3,300)y=np.linspace(-3,3,300)X,Y=np.meshgrid(x,y)

l 绘制等高线

C=plt.contour(X,Y,height(X,Y),10,colors='black')   # 10代表等高线的密集程度。如果是0,则图像被一分为二plt.contourf(X,Y,height(X,Y),10,alpha=0.75,cmap=plt.cm.hot)
#为等高线填充颜色 10表示按照高度分成10层 alpha表示透明度 cmap表示渐变标准

l 添加高度与数字

plt.clabel(C,inline=True,fontsize=10)     #inline=True表示高度写在等高线上plt.xticks(())     plt.yticks(())#去掉坐标轴刻度plt.show()#显示图片

l 最后结果图为:
在这里插入图片描述

实例:寻找某一函数的最小值

import numpy as npimport matplotlib.pyplot as pltimport scipy.optimize as sorosenbrockfunction=lambda x,y: (1-x)**2+100*(y-x**2)**2  #使用lambda关键字所表示的匿名函数来使代码变得更加简洁x=np.linspace(-.5,2.,100)y=np.linspace(-1.5,4.,100)X,Y=np.meshgrid(x,y)Z=rosenbrockfunction(X,Y)cs=plt.contour(X,Y,Z,np.logspace(0,3.5,7,base=10),cmap='gray')   #以第四个参数给出的值为标准来绘制等高线并用gray颜色绘图   #plt.contourf(X,Y,Z,10,alpha=0.75,cmap='gray')rosen=lambda x:rosenbrockfunction(x[0],x[1])solution,iterates=so.fmin_powell(rosen,x0=np.array([0,-0.7]),retall=True)#迭代方法so.fmin_powell利用鲍威尔方法来寻找最小值(由一个给定的起始值x0开始,并且当给出选项retall=True时报告所有的迭代,经过16次迭代后得出x=0,y=0)x,y=zip(*iterates)plt.plot(x,y,'ko')plt.plot(x,y,'k:',linewidth=1)plt.title('Steps of powells method to compute a minimum')plt.clabel(cs)#clabel使用该contour创建的一个等值线集对象来注释对应函数值的水平值

补充:np.logspace()常用于创建等比数列,它也有常用的3个参数,第一个参数表示起始点的指数,第二个参数终止点的指数,第三个参数表示数列的个数;最后,通过base参数可以修改底数;
在这里插入图片描述
(4) 图像

Imshow()函数:用于在指定的窗口中显示一幅图像。它负责对图像进行处理,并显示其格式,但是不显示图片,其后跟着plt.show()才能显示出来。

imshow图像显示函数:imshow(f,G) 其中,f是函数,G是显示该图像的灰度级数。

                Imshow(f,[low,high]) 小于或者等于low的值都显示为黑色,大于或等于high的值都显示为白色Imshow(f,[ ]) [ ] 自动将low设为函数最小值,将high设为函数最大值。

Colorbar:增加颜色类标的代码是plt.colorbar()

例如:

import numpy as npimport matplotlib.pyplot as pltX=np.array([[1,2],[3,4],[5,6]])plt.imshow(X)plt.colorbar()plt.show()   

在这里插入图片描述

实例:将数组可视化为图像,

import numpy as npimport matplotlib.pyplot as plt

下面代码演示了如何通过meshgrid生成复数参数值的矩阵:

def mandelbrot(h,w,maxit=20):X,Y=np.meshgrid(np.linspace(-2,0.8,w),np.linspace(-1.4,1.4,h))c=X+Y*1jz=cexceeds=np.zeros(z.shape,dtype=bool)for iteration in range(maxit):z=z**2+cexceeded=np.abs(z) > 4exceeds_now=exceeded & (np.logical_not(exceeds))exceeds[exceeds_now] = Truez[exceeded]=2return (exceeds)   #命令imshow将矩阵展示为图像,所选的颜色贴图显示了序列出现的白色无界的区域,而其他区域显示为黑色。plt.imshow(mandelbrot(400,400),cmap='gray')plt.axis('off')    #使用plt.axis('off') 关闭坐标轴plt.show()

#使用imshow的无插值的图像
在这里插入图片描述

注意:imshow命令默认使用插值使图像看起来更好,这在矩阵较小时可以看的很清楚。

plt.imshow(mandelbrot(400,400),interpolation='nearest',cmap='gray')

Interpolation代表插值运算,nearest表示最近邻近插值法,只是插值方式中的一种。

cmap表示绘图时的样式,这里选择的是ocean主题。

#使用最近邻近插值法所得的图像
在这里插入图片描述

这篇关于Python3.0科学计算学习之绘图(二)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/542358

相关文章

Go学习记录之runtime包深入解析

《Go学习记录之runtime包深入解析》Go语言runtime包管理运行时环境,涵盖goroutine调度、内存分配、垃圾回收、类型信息等核心功能,:本文主要介绍Go学习记录之runtime包的... 目录前言:一、runtime包内容学习1、作用:① Goroutine和并发控制:② 垃圾回收:③ 栈和

Android学习总结之Java和kotlin区别超详细分析

《Android学习总结之Java和kotlin区别超详细分析》Java和Kotlin都是用于Android开发的编程语言,它们各自具有独特的特点和优势,:本文主要介绍Android学习总结之Ja... 目录一、空安全机制真题 1:Kotlin 如何解决 Java 的 NullPointerExceptio

重新对Java的类加载器的学习方式

《重新对Java的类加载器的学习方式》:本文主要介绍重新对Java的类加载器的学习方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1、介绍1.1、简介1.2、符号引用和直接引用1、符号引用2、直接引用3、符号转直接的过程2、加载流程3、类加载的分类3.1、显示

Java学习手册之Filter和Listener使用方法

《Java学习手册之Filter和Listener使用方法》:本文主要介绍Java学习手册之Filter和Listener使用方法的相关资料,Filter是一种拦截器,可以在请求到达Servl... 目录一、Filter(过滤器)1. Filter 的工作原理2. Filter 的配置与使用二、Listen

Java进阶学习之如何开启远程调式

《Java进阶学习之如何开启远程调式》Java开发中的远程调试是一项至关重要的技能,特别是在处理生产环境的问题或者协作开发时,:本文主要介绍Java进阶学习之如何开启远程调式的相关资料,需要的朋友... 目录概述Java远程调试的开启与底层原理开启Java远程调试底层原理JVM参数总结&nbsMbKKXJx

Java深度学习库DJL实现Python的NumPy方式

《Java深度学习库DJL实现Python的NumPy方式》本文介绍了DJL库的背景和基本功能,包括NDArray的创建、数学运算、数据获取和设置等,同时,还展示了如何使用NDArray进行数据预处理... 目录1 NDArray 的背景介绍1.1 架构2 JavaDJL使用2.1 安装DJL2.2 基本操

HarmonyOS学习(七)——UI(五)常用布局总结

自适应布局 1.1、线性布局(LinearLayout) 通过线性容器Row和Column实现线性布局。Column容器内的子组件按照垂直方向排列,Row组件中的子组件按照水平方向排列。 属性说明space通过space参数设置主轴上子组件的间距,达到各子组件在排列上的等间距效果alignItems设置子组件在交叉轴上的对齐方式,且在各类尺寸屏幕上表现一致,其中交叉轴为垂直时,取值为Vert

Ilya-AI分享的他在OpenAI学习到的15个提示工程技巧

Ilya(不是本人,claude AI)在社交媒体上分享了他在OpenAI学习到的15个Prompt撰写技巧。 以下是详细的内容: 提示精确化:在编写提示时,力求表达清晰准确。清楚地阐述任务需求和概念定义至关重要。例:不用"分析文本",而用"判断这段话的情感倾向:积极、消极还是中性"。 快速迭代:善于快速连续调整提示。熟练的提示工程师能够灵活地进行多轮优化。例:从"总结文章"到"用

AI绘图怎么变现?想做点副业的小白必看!

在科技飞速发展的今天,AI绘图作为一种新兴技术,不仅改变了艺术创作的方式,也为创作者提供了多种变现途径。本文将详细探讨几种常见的AI绘图变现方式,帮助创作者更好地利用这一技术实现经济收益。 更多实操教程和AI绘画工具,可以扫描下方,免费获取 定制服务:个性化的创意商机 个性化定制 AI绘图技术能够根据用户需求生成个性化的头像、壁纸、插画等作品。例如,姓氏头像在电商平台上非常受欢迎,

【前端学习】AntV G6-08 深入图形与图形分组、自定义节点、节点动画(下)

【课程链接】 AntV G6:深入图形与图形分组、自定义节点、节点动画(下)_哔哩哔哩_bilibili 本章十吾老师讲解了一个复杂的自定义节点中,应该怎样去计算和绘制图形,如何给一个图形制作不间断的动画,以及在鼠标事件之后产生动画。(有点难,需要好好理解) <!DOCTYPE html><html><head><meta charset="UTF-8"><title>06