
python
addddddddda
自律
展开
-
python AttributeError: ‘Colorbar‘ object has no attribute ‘set_clim‘
运行程序显示错误:AttributeError: ‘Colorbar’ object has no attribute ‘set_clim’解决:是matplotlib版本过高所致matplotlib3.3.0版本取消了Colorbar的set_clim方法,因此需要降低版本3.2.0即可解决问题首先pip uninstall matplotlib之后pip install matplotlib==3.2.0...原创 2021-06-02 10:45:52 · 2540 阅读 · 0 评论 -
牛顿法求解非线性方程
基础知识[步骤1] 准备:选定初始近似值x0x_0x0,计算f0=f(x0),f0′=f′(x0)f_0=f(x_0),f_0'=f'(x_0)f0=f(x0),f0′=f′(x0)[步骤2] 迭代:按公式:x1=x0−f0/f0′ x_1=x_0-f_0/f_0'x1=x0−f0/f0′迭代一次,得到新的近似值x1x_1x1,计算f1=f(x1),f1′...原创 2019-12-02 20:27:28 · 1539 阅读 · 0 评论 -
下降搜索法
基础知识求解Ax=bAx=bAx=b,若A对称正定,问题可化为求解f(x)=min(12xTAx−bTx)f(x)=min(\frac{1}{2}x^TAx-b^Tx)f(x)=min(21xTAx−bTx)对于f(x)f(x)f(x),其梯度为Ax−bAx-bAx−bpython代码import numpy as npimport matplotlib.pyplot as plt...原创 2019-11-30 20:03:31 · 510 阅读 · 0 评论 -
jacobi迭代法、高斯赛德尔迭代法python实现
jacobi迭代法代码#程序目的:实现jacobi迭代法import numpy as npimport matplotlib.pyplot as pltdef draw_fuction(): x=np.arange(0,8,0.01) y=np.arange(0,8,0.01) x,y=np.meshgrid(x,y) z=4*x+y-24 w=2*x+5*y-30 pl...原创 2019-11-28 13:50:40 · 1280 阅读 · 0 评论 -
CBP(卷积反投影)实现
CBP公式如下:f(x1,x2)=∫0πf^(r,ϕ)∗H(r)∣r=(x1,x2)⋅ϕdϕf(x_1,x_2)=\int_0^{\pi}\hat{f}(r,\phi)*H(r)|_{r=(x_1,x_2)\cdot{\phi}}d\phif(x1,x2)=∫0πf^(r,ϕ)∗H(r)∣r=(x1,x2)⋅ϕdϕ求f(x1,x2)f(x_1,x_2)f(x1,x2)就简单地...原创 2019-11-20 13:26:00 · 3088 阅读 · 0 评论 -
几种滤波器频域时域表示
因为做卷积反投影,所以需要推导滤波函数在时域的表达。R-L 滤波器HRL(w)={∣w∣∣w∣≤wm0∣w∣>wmH_{RL}(w)=\left\{\begin{array}{c}|w| && |w|\leq w_m\\0 && |w|>w_m\end{array}\right.HRL(w)={∣w∣0∣w∣≤wm∣w∣>w...原创 2019-11-14 14:59:32 · 5540 阅读 · 0 评论 -
python实现龙贝格算法
基础部分代码如下程序完全照书上描述写得:#龙贝格算法实现:import numpyimport math#外推法def Trap(f,a,b,Iold,k): #返回用k阶复化梯形公式估计a、b范围内函数的积分 #定义Iold的必要性可以看书上例子p108例子3 if k == 1: Inew=(f(a)+f(b))*(b-a)/2 ...原创 2019-11-14 09:55:09 · 3941 阅读 · 1 评论 -
使用python构建学生-老师-课程选课系统
# 实现学生选课系统import shelveimport random#读取数据库(excel表)内信息# 学生类# 属性:姓名、学号、电话、所选课程列表# 方法:查看:显示该学生所有课程信息;添加课程:将选好的课程添加到课程列表中class Student: def __init__(self,name,number=0,phone=0,classname=None):#设置N...原创 2019-11-11 14:08:56 · 2297 阅读 · 3 评论 -
计算扇形与直线相交长,模拟扇形正投影
交线长计算公式:代码如下:import mathimport structimport numpy as npimport matplotlib.pyplot as plt #plt用于显示图片import sympy as sy#solve方法只能求解线性方程组from scipy.optimize import root,fsolve#扇形相关程序(扇形参数:角度phi,半径...原创 2019-10-29 16:21:03 · 784 阅读 · 0 评论 -
基函数做最佳平方逼近
使用基函数1,x,x2,…,xn1,x,x^2,\dots,x^n1,x,x2,…,xn做最佳平方逼近原函数f(x)=1+x2f(x)=\sqrt{1+x^2}f(x)=1+x2代码如下:import numpy as npimport mathimport matplotlib.pyplot as plt #plt用于显示图片import sympyfrom scipy impo...原创 2019-10-28 14:37:13 · 2562 阅读 · 1 评论 -
Homework of AdvancedProgramming
文章目录第一次作业题目:读取一个二进制图像文件,并尝试进行图像处理(更改其像素数值)再写回文件代码结果第一次作业题目:读取一个二进制图像文件,并尝试进行图像处理(更改其像素数值)再写回文件代码代码如下:其中image.raw是已存在的二进制图像文件。其中有四种灰度值(像素数值)变换函数。import numpy as npimport structimport matplotlib....原创 2019-10-13 20:07:47 · 241 阅读 · 0 评论