- 博客(30)
- 资源 (1)
- 收藏
- 关注
原创 计算机视觉(Python语言描述)(2)————Matplotlib库
目录1. 绘制图像2. 直方图3.交互式标注1. 绘制图像from PIL import Imagefrom matplotlib import pylab as pl# 读取图像到数组中im = pl.array(Image.open('empire.png'))# 绘制图像pl.imshow(im)pl.title('Plotting: "empire.jpg"')pl.axis('off')pl.show()2. 直方图from PIL import Imagefrom
2020-11-28 21:47:54 87
原创 计算机视觉(Python语言描述)(1)————PIL库
目录1.读取和保存2.改变文件格式3.创建缩略图4.复制和粘贴图像区域5.调整尺寸6.旋转1.读取和保存from PIL import Imageimport osprint(filelist)pil_im = Image.open('empire.jpg')pil_im.show()print(pil_im)pil_im = Image.open('empire.jpg').convert('L')pil_im.save("empirehuidu.jpg")pil_im.show()
2020-11-28 21:47:33 162
原创 Python类与对象
class input_net_args(object): def __init__(self, arg_line, arg_node, arg_init_val): self.Line = arg_line self.Node = arg_node self._Init_val = arg_init_val self.Init_val = [] self.Node_infos = [] self._l
2020-11-27 23:21:19 79
原创 Python的深复制和浅复制
def abc(x): x[0] = 1/x[0] x[1] = 1/x[1] return x[0], x[1]Line_arg = [[1, 2, 0.37, 0.404, 0, 0.12264], [2, 3, 0.23, 0.309, 0, 0.16046], [2, 4, 0.585, 0.445, 0, 1.61600], [4, 5, 0.6, 0.170, 0, 0.53379],
2020-11-27 22:59:11 65
原创 Leetcode题目(1)————顺时针打印矩阵
import numpy as npa=np.array([[1,2,3,4,100], [5,6,7,8,23], [10,1,12,13,26], [10,12,13,15,45], [3,56,34,12,13]])upper=0bottom=a.shape[0]-1right=a.shape[1]-1left=0while a.sum()>0: for i in range(l
2020-11-24 20:46:39 78
原创 潮流计算程序————稀疏技术(3)
用类来封装前面的实现的稀疏矩阵存储方式import numpy as npimport randomfrom matplotlib import pyplot as pltplt.rcParams['font.sans-serif']=['SimHei']plt.rcParams['axes.unicode_minus'] = False# matplotlib画图中中文显示会有问题,需要这两行设置默认字体class sparse: def __init__(self, A):
2020-11-24 20:34:48 419 2
原创 潮流计算程序————稀疏技术(1)
前言稀疏矩阵的存储特点是排零存储,即只存储其中的非零元和有关的检索信息。节省内存方便地检索和存取方便地对存储的信息加以修改1.散居格式import numpy as npimport randomfrom matplotlib import pyplot as pltplt.rcParams['font.sans-serif']=['SimHei']plt.rcParams['axes.unicode_minus'] = False# matplotlib画图中中文显示会有问题,需
2020-11-24 15:51:22 469
原创 潮流计算程序————支路功率计算与输出程序
print('迭代求解最终误差为{}'.format(error))print('各节点的实际电压标幺值E(节点号从小到大排):')for i in range(n): print('{:14}'.format(i+1), end=' ')print('')for i in range(n): if f[0, i] >= 0: print('{:.5f}+{:.5f}j'.format(e[0, i], f[0, i]), end=' ') else:
2020-11-24 09:52:21 1915
原创 潮流计算程序————求解雅克比矩阵
e,f是节点电压;n是节点数;G,B是导纳矩阵的实部和虚部;isb是平衡节点的标号;B2是节点分类标号;P,Q,V是各节点初始值;def jabic(e, f, n, G, B, isb, B2, P, Q, V): N0 = 2 * n N = N0 + 1 C = [0 for i in range(n)] D = [0 for i in range(n)] J = np.mat(np.zeros((2*n, N))) for i in ran
2020-11-24 08:20:01 4079
原创 潮流计算程序————牛顿拉夫逊法潮流计算
主函数import numpy as npfrom page000066 import fgaussfrom page000066 import jabicfrom page000066 import fanshufrom matplotlib import pyplot as pltplt.rcParams['font.sans-serif'] = ['SimHei']plt.rcParams['axes.unicode_minus'] = Falsen = 5
2020-11-23 17:04:40 3660 2
原创 C++学习笔记
目录(1)元素逆置(1)元素逆置#include<iostream>using namespace std;int main(){ int a[] = { 6,7,5,15,36,1,0,8}; int start = 0; int end = sizeof(a) / sizeof(a[0])-1; int end1 = end; int temp; while (end>start) { temp = a[end]; a[end] = a[start];
2020-11-19 23:33:24 47
原创 潮流计算程序————节点导纳矩阵及其算法(Python语言描述)
前言目录前言形成节点导纳矩阵的原始数据程序框图代码例子形成节点导纳矩阵的原始数据程序框图代码例子
2020-11-17 13:30:45 2856 2
原创 Python神经网络(1)————训练简单的分类器
代码:from matplotlib import pyplot as pltx1 = [(3, 1), (1, 3)]k = 0.25 # 初始值item = 0L = 0.5 # 学习率error = 1 # 误差值k_record = []while error > 0.0001:
2020-11-16 22:08:45 174
原创 Python数值积分与数值微分(1)————龙贝格积分法
#龙贝格积分法函数import numpy as npfrom matplotlib import pyplot as pltfrom mpl_toolkits.axes_grid1.inset_locator import inset_axesdef romberg(fun, a, b): sum = 0 t1 = ((b - a)*(fun(a)+fun(b)))/2 for i in range(1, 2): sum = sum + fun(a+(
2020-11-15 09:04:20 721 1
原创 Python函数最优逼近(1)————最小二乘拟合函数
import numpy as npfrom matplotlib import pyplot as pltdef fgauss(A, b): n = A.shape[0] zengguan = np.hstack((A, b.T)) ra = np.linalg.matrix_rank(A) rz = np.linalg.matrix_rank(zengguan) temp1 = rz - ra if temp1 > 0: pr
2020-11-13 23:23:45 902
原创 Python插值法(6)————几种插值方法的比较
这里对三次样条插值函数,拉格朗日插值多项式,牛顿插值多项式这三种插值方法进行比较代码import numpy as npfrom sympy import *from matplotlib import pyplot as pltdef lagrange(x, y): p, la = symbols('p la') n = len(x) s = 0 for k in range(n): la = y[k] for j in ra
2020-11-11 23:34:25 2995 1
原创 Python插值法(5)————三次样条插值函数
代码import numpy as npfrom sympy import *from matplotlib import pyplot as pltdef spline(x, y, t): n = len(x) h = [x[i]-x[i-1] for i in range(1, n)] A = np.eye(n-2) A = A*2 u1 = [h[i]/(h[i]+h[i+1]) for i in range(n-2)] la1 = [1-u
2020-11-11 23:09:11 2142 1
原创 Python插值法(3)————hermite插值多项式(拉格朗日型)
代码import numpy as npfrom sympy import *from matplotlib import pyplot as pltdef hermite(x, y, dy): p, la = symbols('p la') f = 0 n = len(x) for i in range(n): la = 1 lp = 0 for j in range(n): if j !
2020-11-10 23:25:34 1975
原创 Python插值法(2)————牛顿插值多项式
代码import numpy as npfrom sympy import *from matplotlib import pyplot as pltdef newton(x, y): n = len(x) c = np.zeros((n, n)) print(c) for i in range(n): c[i, 0] = y[i] for i in range(1, n): for j in range(i, n):
2020-11-10 21:58:49 1529
原创 sympy符号计算库
引言SymPy是用于符号数学的 Python 库。详细教程见目录引言加载sympy库定义符号,替换代码结果代码结果展开代码结果加载sympy库from sympy import *定义符号,替换代码x = symbols('x ')str_expr = 'x**2 + 2*x + 1'expr = sympify(str_expr)y = expr*expr+1print(expand(y))print(y.subs(x, 1))结果代码x = symbols('x ')
2020-11-10 19:24:31 138
原创 潮流计算程序————电力网络的节点编号优化(Tinney-1编号)(Python语言描述)
Tinney-1编号import numpy as npdef tinney_1(id): it = np.mat(np.zeros((1, len(id)))) flag = np.mat(np.zeros((1, len(id)))) k = 1 newnum = 0 while 1: have = 0 for i in range(len(id)): if flag[0, i] == 0:
2020-11-09 14:25:20 1920 4
原创 数据结构与算法系列(3)————递归解法(找零问题)
递归解法def recDC(coinValueList, change, knownResults): minCoins = change if change in coinValueList: knownResults[change] = 1 return 1, knownResults elif knownResults[change] > 0: return knownResults[change], knownResul
2020-11-04 23:02:45 304
原创 数据结构与算法系列(4)————动态规划算法(找零问题)
动态规划算法import numpy as npdef recDC(coinValueList, change, knownResults): minCoins = change if change in coinValueList: knownResults[change] = 1 return 1, knownResults elif knownResults[change] > 0: return knownResu
2020-11-04 23:01:10 233
原创 数据结构与算法系列(2)————贪心算法(找零问题)
代码:# Greedy Algorithmimport numpy as npdef change(big, real): money = [100, 50, 20, 10, 5, 1, 0.5, 0.1] mchange = 0 t = big-real amount = np.zeros((1, len(money)),dtype=int) i = 0 while t >= mchange: if (t-mchange)
2020-11-04 15:15:44 1258
原创 数据结构与算法系列(1)————线性结构栈
栈class Stack: def __init__(self): self.items = [] def isEmpty(self): return self.items == [] def push(self, item): self.items.append(item) def peek(self): return self.items[len(self.items)-1] def siz
2020-11-03 16:26:56 131
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人