weixin_51020254的博客

import mathimport randomimport numpy as npimport pandas as pdimport matplotlib.pyplot as pltfrom datetime import datetime# 子程序：定义优化问题的目标函数def cal_Energy(X, nVar, mk): # # m(k)：惩罚因子，随迭代次数 k 逐渐增大 p1 = (max(0, 6*X[0]+5*X[1]-60))**2 p2 = (max(0

基于奇异值分解压缩数据 pythondata5 2 1 4 0 0 2 4 0 0 00 0 0 0 0 0 0 0 0 3 01 0 5 2 0 0 3 0 3 0 10 5 0 0 4 0 1 0 0 0 00 0 0 0 0 4 0 0 0 4 00 0 1 0 0 0 1 0 0 5 05 0 2 4 2 1 0 3 0 1 00 4 0 0 5 4 0 0 0 0 50 0 0 0 0 0 4 0 4 5 00 0 0 4 0 0 1 5 0 0 00 0 0 0 4 5

莱斯利Leslie种群模型 pythonsympy# 莱斯利Leslie种群模型import numpy as npimport sympy as spX0 = np.array([500, 1000, 500])L = np.array([[0, 4, 3], [0.5, 0, 0], [0, 0.25, 0]])X1 = L @ X0X2 = L @ X1X3 = L @ X2Ls = sp.Matrix([[0, 4, 3], [sp.Rational(1, 2), 0, 0]

pythonsympy库的使用 案例斐波那契数列''import sympy as sp# 斐波那契数列# 定义变量k、k = sp.var('k', positive=True, integer=True)# 新建符号变量时可以指定其定义域a = sp.Matrix([[0,1], [1,1]])val = a.eigenvals() # 求特征值vec = a.eigenvects() # 求特征向量P,D = a.diagonalize() # 把a相似对角化ak = P @(D**

学习笔记，仅供参考。

from scipy.optimize import fsolvefrom scipy.optimize import root# print(help(root))## fsolve 第一个参数函数表达式 ，第二个参数，解的预测估计，第三个参数，函数中的参数a, bdef fun(i): x = i[0] y = i[1] return [x ** 2 + y ** 2 - 1, x - y]s1 = fsolve(fun, [1, 1])print(

python实现K均值聚类，层次聚类，

s = input("请输入字符串")a = []b = []c = []d = []for i in s: if 'a'<=i<='z': a.append(i) if '1'<=i<='9': print(i) b.append(i) if i == ' ': c.append(i)print("字母个数：",len(a))print("数字个数：",len(b))print

python，有四个数字，1,2,3,4，能组成多少个互不相同且重复数字的三位数？各是多少？for i in range(1, 5): for j in range(1, 5): for k in range(1, 5): if i != j and j != k and i != k: print(i, j, k, end='\t')...

python 打印九九乘法表，要求对任意给定的一个正整数N,输出从11,到NN的部分口句表N = int(input("您想计算多少"))for i in range(1, N+1): for j in range(1, i+1): print("{}*{}={}".format(i, j, i*j), end='\t') print('')...

a = int(input("您想计算前多少项的和："))sum = 0for i in range(1, a+1): sum = sum + 1/iprint(sum)

请输入一周中某天的名称的第一个英文字母来判断是星期几，如果第一个字母一样继续判断第二个字母'''请输入一周中某天的名称的第一个英文字母来判断是星期几，如果第一个字母一样继续判断第二个字母'''day = input("请输入星期几")if day[0] == 'M': print("星期一")elif day[0] == 'T': if day[1] == 'u': print("星期二") else: print("星期四")elif

有五个人坐在一起，问第五个人多少岁，他说比第四个人大2岁。问第四个人多少岁，他说比第三个人大两岁文第三个人多少岁，他说比第二个人大两岁，问第二个人多少岁，他说比第一个人大两岁，最后问第一个人多少岁他说是十岁，问第五个人多少岁'''有五个人坐在一起，问第五个人多少岁，他说比第四个人大2岁。问第四个人多少岁，他说比第三个人大两岁文第三个人多少岁，他说比第二个人大两岁，问第二个人多少岁，他说比第一个人大两岁，最后问第一个人多少岁他说是十岁，问第五个人多少岁'''peo = 10n = 1while

python操作mysql数据库import pymysql#db = MySQLdb.connect("127.0.0.1","root","123456","test" )dbpymysql.connect(host='127.0.0.1',user='root',password='123456',port= 3306)curror = db.cursor()#使用cursor()方法获取操作游标curror.execute('SELECT VERSION()')# 使用execute方

python批量发送邮件，办公自动化from email.mime.text import MIMETextimport smtplibfrom email.header import Headerdef fuck(): # 登录邮件服务器 smtp_obj = smtplib.SMTP_SSL('smtp.qq.com',465) # 发件人邮箱的smtp服务器，端口是25 smtp_obj.login('',"") # 发件人的密码，邮箱账号 # 设置邮件头信

'''给定一个头字节为head的非空链表， 返回链表的中间节点。 如果有两个中间节点， 则返回第二个中间节点。例如： [1,2,3,4,5]输出： 此列表的节点为3'''a = [1,2,3,4,5]if len(a)%2 == 0: print("此列表的节点为%d" % a[int(len(a)/2)])else: print("此列表的节点为%d" % a[int(len(a)/2)])...

mysql 的指令 这是我的笔记 show databases;drop show databases;database python;create database pythonew charset=utf8;use mysql;select database(); show tables; create table xxxx(id int, name varchar(30));create table yyyy(id int primary key not null auto_incr

爬取简书网 requests ，csv**这里只爬取了一个页面 只解析了一个 **import requestsimport pymysqlfrom lxml import etreeheaders = { 'User-Agent': 'Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/86.0.4240.75 Safari/537.36'}url = 'h

爬取豆瓣音乐 250， python 保存为csv 格式import requestsfrom bs4 import BeautifulSoupimport csvfp = open('../爬虫/webspider/豆瓣音乐.csv', 'a')write = csv.writer(fp)write.writerow(['歌曲名称', '内容', '评分', '评价'])urls = ['https://music.douban.com/top250?start={0}'.format(

python 爬取豆瓣图书， 保存为CSV格式import requestsfrom bs4 import BeautifulSoupimport csvurls = ['https://book.douban.com/top250?start={0}'.format(i) for i in range(10)]for url in urls: headers = {'User-Agent': 'Mozilla/5.0 (Windows NT 10.0; Win64; x64) Appl

爬取代理ip 保存为CSV格式import requestsfrom lxml import etreeimport csvfrom multiprocessing.dummy import Poolfp = open("代理.csv", 'a')write = csv.writer(fp)write.writerow(['ip', '端口', "位置", '运营商', '更新时间'])def get_info(url): try: resp = reques

爬取百度百科词条 requestimport requestsfrom bs4 import BeautifulSoupunique_url_list = []used_url = []c_lkist = []def m(url): headers = {'User-Agent': 'Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/86.0.4240.75

selenium爬取豆瓣电影 保存CSV格式from selenium import webdriverimport csvfrom multiprocessing import Poolfp = open("电影.csv",'a')write = csv.writer(fp)write.writerow(['电影名称','评分','评论人数','代表名言','电影信息'])brower = webdriver.Chrome('C:\Program Files\Google\Chrome\A

multiprocessing 爬取京东内衣数据 json数据 ， 保存为 excelimport requestsimport jsonfrom openpyxl import Workbook,load_workbookfrom multiprocessing import Poolwb = Workbook()sheet = wb.activesheet.title = "neiyii"wb.save("neiyii.xlsx")headers = {"User-Agent":

爬取猫眼电影 selenium 保存为csv格式from selenium import webdriverimport csvfp = open('猫眼电影.csv','a')write = csv.writer(fp)write.writerow(['电影名称','主演','发行时间','评分','电影链接'])brower = webdriver.Chrome('C:\Program Files\Google\Chrome\Application\Chromedriver.exe')

知乎 cookie登录 爬取,requestsimport requestsheaders = {'accept': 'text/html,application/xhtml+xml,application/xml;q=0.9,image/avif,image/webp,image/apng,*/*;q=0.8,application/signed-exchange;v=b3;q=0.9','accept-encoding': 'gzip, deflate, br','accept-langua

多线程爬取 queue Threadfrom threading import Threadimport requestsfrom lxml import etreefrom queue import Queue# 存储HTML class SpiderUrl(Thread): def __init__(self, url_queue, html_queue): self.url_queue = url_queue self.html_queue =

''python’如果一个数恰好等于它的因子之和，这个数就称为完数，例如6 = 1+ 2+ 3,输出100内的所有完数’’’print("思路： 得到一个数的因子，然后因子相加，等于该数 ，")i = 1while i<100: sum = 0 for m in range(1, i): if i % m == 0: sum = sum + m if sum == i: print("{0}数是完数".for

python a+aa+aaa+…+aaa…aaaaa = int(input("请输入数字a："))num = int(input("你想让几个数字相加："))n = 0sum = 0A = awhile n<num: n += 1 sum = sum+a a =A +a*10print(sum)拓展 另一种解法： List item结语 ： 希望这篇文章能对你有所帮助，我也希望能很你共同进步，欢迎留言提问。谢谢！...

题目：python输入第十个斐波那契数列，0,1,1,2,3,5,8,13,21，34思路：斐波那契数列定义 ：'''递归'''def s(x): n = [] if x ==1 or x == 2: return 1 else: return s(x-2) +s(x-1)s = s(6)print(s)拓展 另一种解法：字符串切片，切后两个数字，然后相加，赋予u， 将u添加到num中 List item结语 ： 希望这篇文章

题目：对于一个正整数分解质因数python思路：’’‘找到最小因数k，然后将正整数n与k相除得到余数，然后余数与k相除，如果不能相除，就k+1再除先用数学来算一遍，然后再用编程写、90’’’‘’'90 / 2 = 45余数045 / 3 = 15 余数015 / 3 = 5余数0因数 2 3 3 5''number = int(input("请输入一个正整数："))yinshu = []# **因数拼音**while number != 1: u = [] for i

题目：实例013：所有水仙花数pyt 打印出所有的"水仙花数"，所谓"水仙花数"是指一个三位数，其各位数字立方和等于该数本身。 例如：153是一个"水仙花数"，因为153=1的三次方＋5的三次方＋3的三次方思路：*获得个位 (i%10)3获得百位 (i//100)**3获得十位 [(i//10)%10]**3for i in range(100,1000): if ((i//10)%10)**3+(i//100)**3+(i%10)**3 == i : print(i)

题目：打印九九乘法表python思路：思路 遍历1到9，两次;;# 多种方式打印九九乘法表# 多多种样式打印九九乘法表'''## 样式1for i in range(1,10): for m in range(1,i+1): print("{0}*{1}={2}".format(i,m,i*m),end="\t") print("\n")## 样式二 s= 11for i in range(1,10): s = s-i for m

题目：# 输入三个数，xyz，将这三个数由小到大排列python思路：判断大小，比较xyz三者的大小，使用if语句x = int(input("请输入一个数"))y = int(input("请输入一个数"))z = int(input("请输入一个数"))if x > y: s = y y = x x = sif x > z: s = z z = x x = sif y > z: s = z z = y

题目：输入某年某月某日，判断这一天是这一年的第几天python输入某年某月某日，判断这一天是这一年的第几天python思路：*判断闰年 能被4整除但不能被100整除，年份能被400整除#方法1sum = 0if (year%4==0and year%100!=0) or year%400==0: feb = 29else: feb = 28month_day = [0,31,feb,31,30,31,30,31,31,30,31,30,31]for i in month_d

输出100~1000内的所有素数思路 ： 首先 什么是素数？素数一般指质数。质数是指在大于1的自然数中，除了1和它本身以外不再有其他因数的自然数。然后利用数学思想： 只有1 和它本身，能除它，for i in range(100,1000): for m in range(2,i): if i%m == 0:# 如果1 和它本身之间的数能除它，就不是素数， break else: print("是素数",i)# 否者不是素数

‘’‘对于一个五位数，判断它是不是回文数，如12321 是回文数，十位与万位相同，，个位与千位相同’’’for i in range(10000,100000): if str(i)[0] == str(i)[4] and str(i)[1] == str(i)[3] : print(i)# 代码更加简短s = [i for i in range(10000,100000) if str(i)[0] == str(i)[4] and str(i)[1] == str(i)[