Panda_1875-CSDN博客

原创 pyplot的绘图转换为numpy矩阵

1.核心原理height_default, dpi_default= (4.8, 100) fig = plt.figure()plt.plot(x, y)fig.canvas.draw()fig_str = fig.canvas.tostring_rgb()data = np.frombuffer(fig_str, dtype=np.uint8).reshape((int(height_default*dpi_default), -1, 3)) / 255.0主要步骤为：（1）创

2021-10-03 15:53:41 1788 1

原创避免pyplot出现Fail to allocate bitmap报错的方法

1. 问题复现当重复绘图（主要是重复创建figure）之后，即使只是在函数内部的局部变量有定义，或者已经使用plt.close()关闭绘图，pyplot内依然保留部分内容，但具体是哪一部分暂时没找到。下面是来自Github的bug复现片段：import matplotlib.pyplot as pltimport numpy as npimage = np.random.randn(256, 256, 3)image= np.array(image/image.max()*255, dtype

2021-09-17 14:59:44 1562

原创解决matplotlib在windows下缩放不兼容导致的绘图模糊问题

1. 原因window应用缩放模糊也是祖传的了，一般的解决方法是：[右键点击软件]->[属性]->[兼容性设置]->[更改高DPI行为]->[勾选"替代高DPI缩放行为",并选择"应用程序"]matplotlib的问题在于找到需要修改兼容性设置的"应用程序"。2. 解决方法对"python.exe"执行以上操作即可python.exe的路径可以通过cmd找到：打开该目录，[右键点击对应版本的python.exe]->[属性]->[兼容性设置]->[更改

2021-08-27 11:26:49 709

原创基于python的掺杂介质六面体nastran网格生成脚本

1.效果图1.1 网格Mesh1Mesh_1Mesh1Mesh2Mesh_2Mesh21.2 应力（使用COMSOL求解）missesmissesmisses应力切应力3.程序import numpy as npimport scipy.interpolate as siimport matplotlib.pyplot as pltN_origin= 10 # 初始矩阵尺寸N_inter= 50 # 插值之后矩阵尺寸t_0=

2021-08-01 22:34:30 1007 1

原创 lammps .eam势文件 metal单位制转real单位制的python程序

1.程序# 1 eV/mol = 1.6021766208×10^-19 J# 1 J*6.02214076×10^23 = J/mol# 1KJ/mol = 0.2389029576186153 Kcal/molC= 1.6021766208*6.02214076*10*0.2389029576186153file_name_read = "Cu_u3.eam"file_name_write = "/home/zxy/Documents/lammps/eam_test/Cu_u3_rea

2021-07-25 11:08:58 1215 4

原创基于python与scipy拟合椭圆

1.效果2. 原理椭圆有个性质：椭圆上的点到椭圆两焦点的距离为常数。令：1、椭圆上的点到两焦点的距离之和为 LtargetL_{target}Ltarget2、两焦点坐标为(xfocus,1 ,yfocus,1),(xfocus,2 ,yfocus,2)(x_{focus,1}\,,y_{focus,1}),(x_{focus,2}\,,y_{focus,2})(xfocus,1,yfocus,1),(xfocus,2,yfocus,2)3、nnn个样本点坐

2021-06-23 11:20:56 5606 9

原创查找曲线峰值并使用梯度下降法拟合为幂函数

1.用途用于查看曲线趋势（衰减规律等）2.效果此处拟合为的形式3.脚本PS：直接运行的图可能不太好看，因为需要载入调整过的style文件，即第47行的 # plt.style.use("zxy")import numpy as npimport matplotlib.pyplot as pltdata= [ 0.6510, 0.6583, 0.6630, 0.6702, 0.6782, 0.6870, 0.7...

2021-05-09 11:40:59 1340

原创 python 计算牛顿差商，计算并化简牛顿插值多项式

1. 程序1.1 导入库import numpy as npimport sympy as sy1.2 函数def NewtonInt(x_value, y_value): """ 计算牛顿插值差商表 :param x_value: x序列 :param y_value: y序列 :return: 差商表矩阵, 字符串插值多项式 """ N=len(x_value) res= np.zeros((N, N)) re

2021-03-28 19:10:30 1683 1

原创 python 牛顿迭代法求解非线性方程有限区间内所有实数数值解

1. 引言上一篇文章python 二分法求解非线性方程所有实数解，初步实现了对非线性方程解的查找。由于采用的是二分法，收敛速度较慢，这里改进为牛顿迭代，并对收敛速度进行对比。2.程序2.1 调用库import numpy as npimport matplotlib.pyplot as plt2.2 求解程序def solveFun(fun, min_x=0, max_x=1, n_step=1000, max_lter=64, max_error= 1e-6, plot=Fa

2021-03-04 16:54:39 1516 1

原创 python 二分法求解非线性方程所有实数解

1. 导入库import numpy as npimport matplotlib.pyplot as plt2.求解方法def solveFun(fun, min_x=0, max_x=1, n_step=1000, max_lter=64, plot=False, print_log=False): """ :param fun: 需要求解的函数 :param min_x: 自变量求解区间下限 :param max_x: 自变量求解区间上限 :p

2021-03-03 22:12:29 2591 4

原创华硕 fx-pro(GL552VW) 修改BIOS，关闭超线程

1.修改原因超线程可以提高部分情况下的CPU利用率，但是当CPU处于高负载下时，超线程反而可能导致性能降低，此外部分数值计算软件对于超线程的支持也不太完善。具体可以查看详细测评。但是，ASUS自带BIOS隐藏了该选项，所以需要修改BIOS才能关闭。2.工具序号描述1XTW-2 编程器及其上位机2与bios芯片对应的编程座（夹），此处为SOP83另一台WIN7、关闭签名验证的win8、win10电脑4电烙铁或热风枪及其配套的焊锡等3一些耐心和动

2020-12-20 19:38:20 4570 2

原创每日自动签到的bash脚本（get获取状态，post提交信息）

1.脚本#!/bin/bash# check status and post dataurl_get="https://___________________" #检查签到状态的链接url_post="https://___________________" #提交数据的链接result_get=$(curl -X GET -H "content-type: application/json" \ --cookie "_________COOKIE设置________" --url

2020-10-17 21:32:17 1106

原创某高校课标页面简化脚本

1.动机学校的网站永远是那么的有个性，电脑还好，手机看课标就比较费眼睛。所以写了个脚本。2.功能1）.自动填充账号密码并登陆2）.自动点击进入课表页面3）.课表页面简化3.效果PC端,修改前 PC端,修改后手机端,修改前手机端,修改后 4.脚本1. 课表页面精简// ==UserScript==// @na...

2020-10-04 21:29:55 156

原创 matplotlib.pyplot 绘制曲线的一些参数设置

1. 概述pyplot是替代matlab的一个很好的选择，不过默认的曲线绘图不够美观，这里对一些细节参数进行了修改，以供各位参考。2. 绘图效果对比修改前（右图）与修改后（左图）对比注：图中曲线数据来源是对文献Hu Chengzhi, Tang Dawei, Lv Jizu, et al. Molecular dynamics simulation of frictional properties of Couette flow with striped superhyd

2020-09-19 15:46:49 2283

原创一个提取lammps的log文件中的数据并绘图的python程序

1. 程序及结果1.1 程序# -*- coding:utf-8 -*-import reimport numpy as npimport matplotlib.pyplot as pltfile_name= "log.lammps" # log 文件名jump_step= 20 # 需要跳过的step数，主要用于忽略势能最小化或弛豫产生的输出数据max_lines= 1

2020-08-04 09:20:34 3341 4

原创 Linux 上Vscde配置python的一些问题和解决方法

1.pip下载缓慢1.1 处理思路对于python2：可以使用清华或者aliyun的镜像对于python3：目前清华镜像显示拒绝连接，aliyun镜像则可以正常使用，推荐使用aliyun，以后清华镜像可能会修复。1.2 实际操作方式临时使用：pip install -i https://mirrors.aliyun.com/pypi/simple some-package设为默认之方式1：pip config set global.index-url https://mi

2020-08-04 08:14:17 215

原创基于Matlab的各地阳光强度计算

1.起因寒假回家那几天，经常看到太阳能热水器中的水沸腾然后和蒸汽一起喷射出来。但最近夏天却不怎么出现这个现象了，开始以为是没水，上楼顶看了一下发现事情并不简单。初步考虑有两个方面的原因： 1）夏天阳光对地面是直射的，但太阳能热水器与地面之间有倾角，这使得夏天热水器不能被阳光直射，冬天反而可以。 2）夏天地日距离较长（坐标昆明），这减弱了到达地球的太阳光强度。2.结果及结论从图d中可以看出，太阳能热水器在冬天接收到的太阳光强度确实是最强的，而夏天则比冬...

2020-06-07 17:04:45 1685

原创基于Matlab的六边形细胞自动机

1.简介目前大多数细胞自动机都是矩形的，这个是六边形的。生存规则在倒数13 - 倒数18行。2.演示下面是两种不同规则的动画。3.程序clear;N_col = 20;% 列数设置N_row = round(N_col/0.433*1.5/1.2);Cell = zeros(N_row,N_col, 14,2);% Cell - [细胞状态细胞中点坐标六个边点坐标六个邻近图形序号]% 细胞状态(1)- [当前状态下一步状态]； 0- 死亡，1-.

2020-06-01 17:04:16 1164 4

原创某PC软件去广告日记

软件A是本人常用的软件，但是有个广告多的缺点，今天来去广告试试1.抓包软件A 在每次启动页面都会加载广告。打开Fiddler，启动软件，果然抓到了线索。可以看出，软件从一个名为 apps_must_config.data 的文件获取到了广告内容索引，文件链接为 res.xxxx.com/apps_must_config.data。2.找出存储该链接的文件planA： notepad++全目录搜索，看看能不能在配置文件中找到，结果和预想的一样，没搜到。planB：一些e.

2020-05-15 10:53:38 215

原创简单Android app之一键签到开发日记

致谢：开发此App期间参考了大量文章与帖子，由于数量众多，来不及一一记录。在此一并感谢！如果在本文中发现与您的程序十分雷同的部分（主要为Http请求部分，Stream转String部分，以及主线程与子线程的通信部分），不胜感激。如有无意冒犯之处，请与我联系。附件：包括： java文件、布局文件、资源文件、AndroidManifest文件。需要Gson和另一个库（忘了）。https://download.csdn.net/download/zxylv/124118771..

2020-05-12 14:56:29 2810

原创 BitLocker 操作命令

由于Win10 1909的控制面板中移除了单个分区加密/解密的操作，在控制面板中无法对单个分区进行解密，导致无法安装双系统。表现为：在Linux下强行格式化磁盘0 的加密 NTFS 分区D为 ext4，会导致磁盘0 的加密NTFS 分区C无法读取。这时候就需要单独对D分区进行解密，再在Linux下进行格式化操作。但是在Power Shell中，针对单分区的相关命令依然有效。...

2020-03-24 10:47:51 6457

原创基于ESP8266的温湿度监测、记录平台

1.先看效果2.硬件元件：ESP8266 12-E模块 *1DHT11 温湿度传感器 *1DS3231时钟模块 *1大约200欧电阻 *1LED指示灯*1硬件的连接：（软件中可以修改，默认是如下连接方式）ESP8266 DS3231 DHT11 LED串联电阻 SCL &l...

2020-02-21 14:28:15 5121 3

原创基于GD32F103RET6与Matlab的串口虚拟示波器日记

1.各功能展示测试用的是大约50KHz的方波，由Arduino Nano产生，程序如下：void setup() { pinMode(6,OUTPUT);}void loop() { digitalWrite(6, HIGH); delayMicroseconds(10); digitalWrite(6, LOW); delayMicroseconds(10);...

2020-02-18 20:16:45 1331 1

原创 GD32f103RET6 （某付宝拆机）折腾日记

很久以前 1 元 1个，拆开是gd32f103。最近没事顺便学习以下STM32类单片机。1.接线比较有用的是以上标红的6个点，各个点定义如下：1：串口 0 的RX，用于烧写程序及通信。2：串口0 的TX，用于烧写程序及通信。3：Boot1，进入bootloader需要将此位置 0。4：3.3V Vcc。5：Boot0，进入bootloader需要将此位置 1。...

2020-02-15 01:46:48 3365

原创 GD32f103 8M晶振改12M ，要修改的地方

手里的单片机是gd32f103ret6，晶振和官方库默认的8M不一致，导致串口乱码。网上找了好久全是STM32的例子，不过还是有参考意义的。以下是gd32f10x 的设置方式。1.Keil中的Target设置：PS：这一项好像会自动设置，安全起见检查一下2. "gd32f10x.h" 中的 HXTAL_VALUE 修改。注意：gd32f10x.h默认是...

2020-02-15 00:54:21 5994

原创基于89C51/52单片机的旋转LED制作、程序、及图片取模

这是很久以前学单片机不久挖的坑，硬件做好了，没找到合适的取模软件。最近比较闲，翻出来收尾一下。先看效果图汤姆猫GitHubCSDN硬件部分电路图，用共阳极是因为51单片机下拉能力比上拉能力强得多。由于LED较多，同时旋转部分应尽量轻，所以原件都采用贴片的。走线用的是飞线。如图：背面是飞线：输电部分用的是电刷。本来想用无线输...

2020-02-12 14:28:59 5764 10

原创 Kudrone/Kudrone 2s 手机app不能登录尝试解决方案[已验证]

群友买了Kudrone的手机遥控四轴飞行器，然后Kudrone倒闭了，app无法登录自己的账号。而飞行器的遥控操作需要先登录账号才能进行，这时就尴尬地导致飞行器随着公司的倒闭而失去使用价值。解决思路是a.反编译app。b.修改AndroidManifest.xml 把activity.MAIN为登录之后的activity（如图），跳过登录界面。c.编译app再签名。...

2020-01-11 23:05:04 3669 24

原创基于Matlab dsp工具箱的实时音频采集及频谱显示

题外话这几天在寻找可以使 Matlab 实时采样音频并处理的函数。刚开始找到的是网上很多例子采用的 analoginput( ) ，但是目前在 64位 Windows上无法使用。其次是 audioRecorder( ) , 可以调用但缺点是每次调用都会开关一次麦克风，延迟十分严重，无法实现要求的实时效果。在 Matlab 帮助文档里游荡一会儿之后，发现DSP工具箱里面有个函数 dsp...

2019-10-26 00:26:33 9206 15

原创基于Matlab的傅里叶变换绘制图形

起源之前制作钢琴频谱识别时候，发现傅里叶变换的神奇之处，激起了内心的好奇。前几天在bilibili又发现了用傅里叶变换绘图的视频。决定自己做一个，顺便仔细学习一下傅里叶变换。过程从对傅里叶变换一窍不通，到反复学习推导过程，以及自己边做边发现问题，加上自己的思考和参考视频中的讲解，终于对傅里叶变换有了一个比较深刻的理解。原理对于一个波，对其进行傅里叶变换之后可以得到若干角速度、初...

2019-08-15 12:25:10 12248

原创 matlab 识别钢琴基频演奏音乐

就是把钢琴曲换个音调再演奏出来。初衷是用来生成matlab的曲谱的，先写这么多吧。识别正确率好像不是很高，推荐用一些单纯的钢琴曲效果好一些，比如张宇桦的。换歌曲的话，注释后面带有“***************”的数的值是影响比较大的，需要修改。一个是阈值，一般3~20之间。一个是延时，一般x0.005~x0.1之间。由于本人不是信息、信号相关专业的，写的不专业，欢迎指正。参考资料：一...

2019-02-18 22:02:38 4580 4

原创 Fortran 求矩阵的逆、行列式的值

#2019，10，8 更新：重写部分程序，增加部分注释学Fortran的第一天，就写了这么点东西，分享一下。内容包括：求矩阵的逆、行列式的值其中：求逆的方法是先求伴随矩阵再除以行列式的值，求行列式的值用的是求余子式的迭代法（从matlab里面的det函数获得的启发），需要注意的是，Fortran中数组的存储是先列后行。注释以后有空再加上吧。初学...

2019-02-11 23:14:46 18583 3

原创将博客搬至CSDN_牛仔1875_新浪博客

因csdn对代码块支持较好，现将博客搬至csdn。个人主页：https://blog.csdn.net/zxylv

2019-02-10 17:23:30 90

原创《有限单元法》--王勖成，习题2.12 MATLAB 程序

节点、单元编号只计算到节点位移，应力没计算%% 《有限单元法》--王勖成，习题2.12 MATLAB 程序% 仅供参考function fem() P=[-10 0 10 0 10 0 -10 0]';%结构节点载荷列向量 [K,P]=xQiYi(P,make_K()); solu(K,P); ...

2019-02-10 17:14:59 1543

原创 89C51单片机读取DHT11数据并用数码管显示_牛仔1875_新浪博客

前几天发了arduino读取dht11的文章【arduino读取DHT11温度、湿度数据并串口显示】，里面DHT11的通信协议已经说明了，这里略过，贴个程序吧。测试过可用的。晶振11.0592Mhz。 #include "reg52.h"unsigned char num[10]={0xde,0x82,0xec,0xe6,0xb2,0x76,0x7e,0xc2,0xfe,0x...

2018-08-06 00:42:34 6078 6

原创 arduino读取DHT11温度、湿度数据并串口显示_牛仔1875_新浪博客

图片及说明在新浪博客，就不搬过来了。ps：两边都是我的博客，csdn更好用一些，所以就搬到这边，不存在抄袭。int out=2; //2号IO口作为通信接口void setup() { Serial.begin(115200); //波特率115200 pinMode(out,INPUT); ...

2018-08-04 12:52:46 2772

原创欢迎您在新浪博客安家_牛仔1875_新浪博客

亲爱的朋友：欢迎您在新浪博客安家，您的博客地址是：http://blog.sina.com.cn/u/3957050892您可以用文字、图片、视频记录和展示最真实的自我，与网友交流，与线上好友聊天，还能通过手机发表博文和上传图片，随时随地记录心情和身边趣闻。我们为您提供了丰富的炫酷模板来装点您在网上的家园，强大的音乐播放功能更能陪伴您的网络生活。准备好了吗？现在就开...

2014-05-12 22:16:30 89

Lammps 安装.pdf

一键签到 源码.zip

8266_services_V7.1.zip

空空如也

一键签到源码.zip