嵩悟空-CSDN博客

原创 python：稀疏矩阵存储，相乘，转置

稀疏矩阵存储，相乘，转置存储：利用三个列表进行存储相乘：这里只写了稀疏矩阵与全矩阵的乘法，并返回全矩阵转置：返回稀疏矩阵这里默认稀疏矩阵较大例如（5000X5000）而全矩阵很小（5000X1）import sysimport numpy as npclass TSMatrix: def __init__(self, rows, cols, value, shape): self.rows = rows self.cols = co

2022-04-11 16:11:27 4592 1

原创 python使用numpy数组object对象赋值时注意事项

python使用numpy数组object对象赋值时注意事项赋值时不能赋同一变量，这样会导致改变该数组一个元素，赋同变量的元素也跟着改变

2022-04-04 10:55:41 910

原创 Python代码：数字图像处理（DIP）7.1.2子带编码example7.2

Subband Coding例子实现的代码如下：import cv2import numpy as npfrom matplotlib import pyplot as pltdef correl_1D(img, window, direction = 'Row'): m = window.shape[0] img2 = np.zeros(img.shape) if direction == 'Row': img1 = np.zeros([i

2021-11-06 10:41:19 2008

原创 Python代码：数字图像处理（DIP）7.1.1图像金字塔example7.1

Approximation and prediction residual pyramids知识请看课本，这里交流代码如下：import cv2import numpy as npfrom matplotlib import pyplot as plt#二维卷积操作def correl2d(img, window): #window表示滤波模板，img表示原始矩阵 #获取模板的行和列 m = window.shape[0] n = window.shape[1

2021-11-03 16:46:37 335

原创 Python：cv2.imread读入图片后返回NoneType的解决方案

cv2.imread读入图片后返回NoneType的解决方案这里用的是Spyder重点如果文件和图片不在一个文件夹下，工作目录也不一致，按照上面代码也会输出NoneType，这时需要把路径写好：这里用的是Spyder先看这段代码import cv2img = cv2.imread('wangye.jpg')cv2.imshow('wangye', img)这里要注意几点工作目录与文件目录在一块才能像上面这样写如下图读图时名字要加后缀（这点很重要，许多返回NoneType就是应

2021-11-03 11:17:16 7726 1

原创 c++:opencv读图后mat矩阵的基本操作

#include <opencv2/opencv.hpp>#include <iostream>using namespace std;using namespace cv;int main(int argc, char** argv){ //读图 Mat scr = imread("C:/Users/dell/Desktop/图片变化/1024_699_3.jpg"); //此为你电脑上图片的路径 //窗口命名 namedWindow("林允

2021-04-19 20:19:20 501

原创 python:共轭梯度法，特征值聚堆情况下迭代次数讨论

共轭梯度法，特征值聚堆情况下迭代次数讨论输入各种特征值聚堆与分散时的矩阵，并应用共轭梯度法，观察迭代次数与聚堆情况的关系。因为对角矩阵的对角线元素为其特征值，则用对角矩阵讨论较为方便代码import numpy as npdef cg(x0, A, b): r0 = np.dot(A, x0) - b p0 = -r0 rk = r0 pk = p0 xk = x0 t = 0 #记录迭代次数 while np.linalg.norm

2021-04-04 10:19:55 384 2

原创 python:共轭梯度法（以希尔伯特矩阵为例）

共轭梯度法（以希尔伯特矩阵为例）实现共轭梯度算法并且使用它解决对称矩阵——希尔伯特矩阵为系数矩阵A时，右边常数项列向量为单位列向量，初始点为零向量。展示维度分别为5，8，12，20时的迭代次数与结果。迭代终止条件为残差的二范数小于1e-6import numpy as npdef cg(x0, A, b): r0 = np.dot(A, x0) - b p0 = -r0 rk = r0 pk = p0 xk = x0 t = 0 #记录迭代次数

2021-04-01 09:38:13 2151

原创 c++:栈的基本操作+实例：迷宫求解

栈的基本操作+实例：迷宫求解#include<iostream>using namespace std;struct Node { int datai; int dataj; Node* next;};//栈struct Stack { Node* pTop; //顶部指针 Node* pBottom; //底部指针};void Init(Stack* pS);//初始化//void CreateStack(Stack* pS);//

2020-10-24 15:34:33 296

原创 python:拉格朗日插值实现及求解

拉格朗日插值实现及求解输入：插值节点与插值节点的函数值，以及所要求近似值的x值;输出：拉格朗日插值多项式（符号函数），x处的函数近似值.def Lagrange(xx,yy,a): import sympy as sym from fractions import Fraction x=sym.Symbol('x') l=[1]*len(xx) i=0 for x0 in xx: for xk in xx: if x0!=xk:

2020-10-18 22:34:01 873

原创 python:字典，元组

创建一个字典，有26个键从"A"到"Z"，每个键关联的值是从1到26。D={}i=0a='A'while 26-i: i+=1 D[chr(ord(a)+i-1)]=i print(D)写个表达式将元组 (4, 5, 6) 改为 (1, 5, 6)。temp=list((4, 5, 6))temp[0]=1T=tuple(temp)T...

2020-10-13 19:09:09 490

原创 C++：顺序表的基本操作（待完善）

根据命令提示进行顺序表的基本操作（待完善）#include<iostream>#include<algorithm>using namespace std;#define MAXSIZE 1000#define OVERFLOW -2#define ERROR 1typedef int ElemType;typedef int Status;typedef struct { ElemType* elem; int length;}SqList;

2020-10-11 16:16:49 96

原创 C++:录入班级数学成绩，计算最大值、平均值、不及格人数

录入班级数学成绩，计算最大值、平均值、不及格人数#include<iostream>using namespace std;int main() { float a=0; float student[100]; int i, j ,N = 0; float max,failnum=0,sum=0; //录入成绩 cout << "录入学生成绩以0为结束"<<endl; do { student[N

2020-10-09 19:49:44 2409

qq_37789035的博客