3 雪shan飞hu

尚未进行身份认证

我要认证

暂无相关简介

等级
TA的排名 11w+

求numpy.array数组里有多少个元素(len,shape,size)

import numpy as npz = np.random.randint(10,size=(3,5))print("z = \n", z)print("len(z) = ",len(z))print("z.size = ", z.size)print("z.shape = ", z.shape)运行结果:z = [[8 5 6 7 7] [8 2 7 5 1] [9 3 0 3 3]]len(z) = 3z.size = 15z.shape = (3, 5)

2020-06-04 17:06:19

玩转Python随机数

本文转载自https://www.cnblogs.com/twilight77/p/7675512.htmlPython中的random模块用于生成随机数。下面介绍一下random模块中最常用的几个函数。random.randomrandom.random()用于生成一个0到1的随机符点数: 0 <= n < 1.0random.uniformrandom.uniform的函数原型为:random.uniform(a, b),用于生成一个指定范围内的随机符点数,两个参数其中一个是上限,

2020-05-28 17:30:07

from sklearn.model_selection import train_test_split 报错 ImportError: cannot import name comb

在出现这个报错前,如果你的报错信息显示没有 sklearn.model_selection,先打开cmd输入"conda list"看一下scikit-learn的版本,比如你可能是0.17.1,低版本是没有“model-selection”模块的,所以需要更新。更新后比如到了0.18.1,报错信息如标题,这是因为已经安装了Anaconda,将scikit-learn升级到最新版本可能会导致与已安装的scipy和numpy版本不匹配。 所以可以使用anaconda的安装程序升级所有三个库,但在此之前必须先

2020-05-26 14:34:41

Anaconda报错CondaHTTPError

一 .condarc(conda 配置文件)Configuration — Conda documentation.condarc以点开头,一般表示 conda 应用程序的配置文件,在用户的家目录(windows:C:\users\username\,linux:/home/username/)。但对于.condarc配置文件,是一种可选的(optional)运行期配置文件,其默认情况下是不存在的,但当用户第一次运行 conda config命令时,将会在用户的家目录创建该文件。二. Conda

2020-05-24 20:34:00

数据库综合练习一

主要使用postgreSQL,数据库下载Question 1:Write an SQL query that gives the number of films whose rental rate is greater than the average rental rate of all movies. Your query must generate a single column cou...

2019-12-03 05:12:46

Python统计list列表各个元素出现次数

计数用的函数是count(),列表中会出现重复的项,所以转化成set去除重复元素import numpy as nplist1 = [1,2,3,1,1,2,2,1,3]list2 = list(set(list1))for i in list2: print(i,": ",list1.count(i))...

2020-04-19 18:37:35

美式期权定价的最小二乘蒙特卡洛算法(理论部分)

* 本文数据来源庄洋洋的论文《基于最小二乘蒙特卡洛法的美式巴黎期权定价》最小二乘蒙特卡洛算法步骤介绍(1)进行N次蒙特卡洛模拟,得到一条N+1个点的股票价格样本路径,重复M次,则得到M条股票价格样本路径。生成每个点的股票价格的公式是举个例子,模拟美式看跌期权,取M=10,N=3,随机生成的路径如下(2)从期权到期日开始逆向求解,用最小二乘法计算各点如果继续持有期权的价...

2020-04-15 00:17:55

C++实现前向欧拉法Forward Euler解决偏微分方程

1. 数学原理以简单的热方程为例,其中D为常数:根据导数的定义式有:如果我们将x理解成delta_x*n,则u(x,t)可以看成只是t的函数,n作为参数来控制x,即定义:这样我们就可以把这个偏微分方程看做是一个常微分方程对于常微分方程的C++实现,可以先看看https://blog.csdn.net/weixin_39374967/article/details...

2020-04-06 19:40:57

C++实现Forward Euler前向欧拉法解常微分方程

用C++处理常微分方程时,一种常用的方法是Forward Euler法。我们用一个简单的ODE举例,第一行可以写成显然这是个迭代的算法,用C++实现时,我们假设lamda=1,h=timestep=0.01,y(0)=1,假设我们要求 y(10) 的值,则是需要迭代1000次const int numberOfTimeSteps = 1000;const double ti...

2020-04-05 21:30:21

文献引用Havard(哈佛)常用格式说明

写论文引用时要求使用Havard格式,简单记录一下常用的一、书籍引用(a)文内引用如果已经点明了作者名字,括号里只写年份即可,如:According to Guy (2001) the Zulus faced many grave dangers when confronting the British…三人以上 This was proved by Dym et al. (2009)…...

2020-04-04 18:57:08

矩阵范数

1、常用的矩阵范数主要是p=1,2,∞2、谱半径3、谱半径与矩阵范数的关系

2020-03-28 18:19:03

数值线性代数之cholesky平方根法(MATLAB)

1. 该方法只适用于对称正定矩阵所以我们可以解得L代入原方程组即可2. MATLAB代码function [L]=mycholesky (A)n=length(A);for k=1:n A(k,k)=sqrt(A(k,k)); A(k+1:n,k)=A(k+1:n,k)/A(k,k); for j=k+1:n ...

2020-03-27 21:32:06

数值线性代数之列主元Gauss消去法(MATLAB)

1. 这是啥?在使用Gauss变换的过程中,如果矩阵对角线的元素太小,会干扰计算结果,所以我们尽量把每列较大的元素通过行变换挪到对角线上。2. 第二个红框只是普通的Guass变换,第一个红框就是选择最大列主元的算法。3. 下图中,左边的L代表的是常规的Gauss变换,P代表行变换4. 于是我们要解的方程就从Ax=b变成了LUx=Pb,接下来使用上下三角矩阵的前代法和回代法解...

2020-03-27 21:11:52

数值线性代数之Gauss分解(MATLAB)

1. 什么是Gauss变换对于一个向量x,左乘一个矩阵Lk,使得x变成向量y,该y向量第k个元素起皆为0。 比如x=(2,2,6)', 将其转化为向量y=(2,0,0)'。其中的Lk为向量ek=(0,0,...,0,1,0,...,0),2. 你可能没怎么看懂或者懒得看,我们来举个栗子所以Gauss变换至少要保证求分母不为0,不是所有矩阵都可以用的。3. 上个例子是针...

2020-03-27 17:57:23

数值线性代数之前代法解下三角方程组和回代法解上三角(MATLAB)

下三角方程组形如解得即MATLAB代码:function b=myFM(L,b)n=length(b);for j=1:n-1 b(j)=b(j)/L(j,j); b(j+1:n)=b(j+1:n)-b(j).*L(j+1:n,j);endb(n)=b(n)/L(n,n);同理,对上三角矩阵MATLAB代码:...

2020-03-27 14:11:59

C++ 蒙特卡洛求积分

最简单的蒙特卡洛法求指数函数在[0,1]上的积分,面积法const int numberOfMCSimulations = 1000000;double lowerBound = 0.0;double upperBound = 1.0;// We need to draw a rectangle which, for positive functions, has a range on...

2020-03-10 00:00:39

MATLAB实现矩阵雅可比(Jacobi)迭代

包含三个function,下列代码最后一段是主函数,其他都是function。三个function建议从下往上看。function [x,k,resvec,DD,ID,JD,D,Ab] = jacobis(AA,IA,JA,b,x,tol,kmax)%This function is an implementation of Jacobi’s%iterative method for so...

2020-03-01 20:13:08

MATLAB实现矩阵given变换

一共三个function文件,分别如下function [H,Q,c,s,B]=givensqr(H)[n,m]=size(H);Q=eye(n);for k=1:ma=H(k,k);b=H(k+1,k);[c,s]=givens(a,b);H=gmatmat(H,c,s,k);M=eye(n);M(k,k)=c;M(k,k+1)=-s;M(k+1,k)=s;M(k+1,k...

2020-03-01 20:08:30

什么是beta分布,贝叶斯共轭先验

一篇非常好的介绍beta分布的文章,讲得非常清楚https://blog.csdn.net/a358463121/article/details/52562940

2020-02-29 07:06:31

C++随机数,掷骰子,蒙特卡洛求Pi

几个典型的随机数分布uniform_int_distributionuniform_real_distributionnormal_distributionlognormal_distributionbernoulli_distributionbinomial_distributionpoisson_distributionexponential_distribution…...

2020-02-25 00:27:48

查看更多

勋章 我的勋章
  • 持之以恒
    持之以恒
    授予每个自然月内发布4篇或4篇以上原创或翻译IT博文的用户。不积跬步无以至千里,不积小流无以成江海,程序人生的精彩需要坚持不懈地积累!
  • 勤写标兵Lv2
    勤写标兵Lv2
    授予每个自然周发布4篇到6篇原创IT博文的用户。本勋章将于次周周三上午根据用户上周的博文发布情况由系统自动颁发。