蔡军帅

原创 PAT 甲级 1012 The Best Rank (25 分)（结构体排序）

题意：为了评估我们第一年的CS专业学生的表现，我们只考虑他们的三个课程的成绩：C - C编程语言，M - 数学（微积分或线性代数）和E - 英语。同时，我们鼓励学生强调自己的最优秀队伍 - 也就是说，在三个课程和平均成绩的四个职级中，我们打印每个学生的最佳排名。例如，C，M，E和A - 4名学生的成绩如下：StudentID C M E A310101 98 85 88 903...

原创 anconda + python 3.6安装(以前的anconda,anaconda和python版本对应关系)

anconda + python 3.6安装anaconda + python3.6安装安装老版本带python3.6的anaconda安装老版本带python3.6的anacondaanaconda下载.--------------------- 作者：Longriver111 来源：CSDN 原文：https://blog.csdn.net/qq_35148758/artic...

原创 数学建模python matlab 编程(喷泉模拟)

在无风情况下的喷泉模拟我的python代码import numpy as npimport randomimport matplotlibmatplotlib.rcParams['font.sans-serif']=[u'simHei']matplotlib.rcParams['axes.unicode_minus']=Falseimport matplotlib.pypl...

原创 数学建模python matlab 编程(疾病传播模型)

例12：一只游船上有800（1000）人，一名游客不慎患传染病，12（10）小时后有3人发病，由于船上不能及时隔离，问经过60（30）小时，72小时，患此病的人数。(与人口模型和Logistic模型类似)先用python和matlab模拟我的python代码# -*- coding: utf-8 -*-import numpy as npimport randomimpor...

原创 scikit-learn机器学习(四)使用决策树做分类,并画出决策树,随机森林对比

数据来自 UCI 数据集 匹马印第安人糖尿病数据集载入数据# -*- coding: utf-8 -*-import pandas as pdimport matplotlibmatplotlib.rcParams['font.sans-serif']=[u'simHei']matplotlib.rcParams['axes.unicode_minus']=F...

原创 scikit-learn机器学习(四)使用决策树做分类

我们使用决策树来创建一个能屏蔽网页横幅广告的软件。已知图片的数据判断它属于广告还是文章内容。数据来自http://archive.ics.uci.edu/ml/datasets/Internet+Advertisements其中包含3279张图片的数据，该数据集中的类的比例是不均衡的，459张图片是广告，零位2820张图片是文章内容。首先导入数据，数据预处理#...

原创 从最大似然到EM算法浅解

https://blog.csdn.net/zouxy09/article/details/8537620机器学习十大算法之一：EM算法。能评得上十大之一，让人听起来觉得挺NB的。什么是NB啊，我们一般说某个人很NB，是因为他能解决一些别人解决不了的问题。神为什么是神，因为神能做很多人做不了的事。那么EM算法能解决什么问题呢？或者说EM算法是因为什么而来到这个世界上，还吸引了那么多世人的目光。...

原创 scikit-learn机器学习(三)多项式回归(二阶，三阶，九阶)

我们仍然使用披萨直径的价格的数据import matplotlibmatplotlib.rcParams['font.sans-serif']=[u'simHei']matplotlib.rcParams['axes.unicode_minus']=Falseimport numpy as npimport pandas as pdimport matplotlib.p...

原创 scikit-learn机器学习(二)逻辑回归进行二分类(垃圾邮件分类),二分类性能指标，画ROC曲线，计算acc,recall,presicion,f1...

数据来自UCI机器学习仓库中的垃圾信息数据集数据可从http://archive.ics.uci.edu/ml/datasets/sms+spam+collection下载转成csv载入数据import matplotlibmatplotlib.rcParams['font.sans-serif']=[u'simHei']matplotlib.rcParams['ax...

原创 scikit-learn机器学习(一)简单线性回归

# -*- coding: utf-8 -*-import numpy as npimport matplotlib.pyplot as plt## 设置字符集，防止中文乱码import matplotlibmatplotlib.rcParams['font.sans-serif']=[u'simHei']matplotlib.rcParams['axes....

原创 Python描述性统计numpy

import numpy as npimport pandas as pdimport matplotlib.pyplot as pltfrom sklearn import datasets,preprocessingfrom sklearn.model_selection import learning_curvefrom sklearn.model_sel...

原创 Python统计分析可视化库seaborn(相关性图，变量分布图，箱线图等等)

Visualization of seaborn  seaborn[1]是一个建立在matplot之上，可用于制作丰富和非常具有吸引力统计图形的Python库。Seaborn库旨在将可视化作为探索和理解数据的核心部分，有助于帮人们更近距离了解所研究的数据集。无论是在kaggle官网各项算法比赛中，还是互联网公司的实际业务数据挖掘场景中，都有它的身影。   在本次介绍的这个项目中，我们将利用se...

原创 Introduction to statistical learning：with Applications in R （书，数据，R代码，链接）

http://faculty.marshall.usc.edu/gareth-james/http://faculty.marshall.usc.edu/gareth-james/ISL/

原创 Python-sympy科学计算与数据处理(求极限及其它功能)

极限 其它功能

原创 Python-sympy科学计算与数据处理(方程，微分，微分方程，积分)

方程a,b,c,x = symbols("a b c x")my_eq = Eq(a*x**2+b*x+c,0)solve(my_eq,x)Out[12]: [(-b + sqrt(-4*a*c + b**2))/(2*a), -(b + sqrt(-4*a*c + b**2))/(2*a)]微分t=Derivative(sin(x...

原创 Python-sympy科学计算与数据处理(数学表达式)

数学表达式from sympy import *1/2+1/3S(1)/2+1/S(3)Out[4]: 5/6Rational(5,10)Out[5]: 1/2x,y = symbols("x y")factor(15*x**2+2*y-3*x-10*x*y)Out[7]: (3*x -...

原创 用Python解方程

一元一次方程例题1： 这是北师大版小学六年级上册课本95页的一道解方程练习题：大家可以先口算一下，这道题里面的x的值为200接下来我们用python来实现，代码如下，每一句代码后面都写有解释语：# 一元一次方程x = sy.symbols("x") # 申明未知数"x"a = sy.solve((x+(1/5)*x-240),[x]) # 写入需要解的方程体pri...

原创 数学建模python matlab 编程(椭圆声学原理画图证明，解析几何)

证明，在椭圆形的音乐厅内，从一个椭圆的一个焦点发出声音，则另一个焦点听到的声音是最大的。分析：证明，从椭圆的一个焦点任意发射的直线经过反射后，并经过另一个焦点。 画图，过一个焦点随机画一条直线l1，再做它与椭圆的交点，过该点作椭圆的一条切线l2，在作一条与该切线垂直的线l3，再作一条l1关于l3对称的直线l4，然后计算另一个焦点到l4的距离小于一个很小的数即可...

原创 数学建模python matlab 编程(指派问题)

指派授课问题现有A、B、C、D四门课程，需由甲、乙、丙、丁四人讲授，并且规定：每人只讲且必须讲1门课；每门课必须且只需1人讲。四人分别讲每门课的费用示于表中：课费用人ABCD甲21097乙154...

原创 matlab中如何给一个矩阵中的某几个特定位置赋值

用sub2ind>> a=zeros(5); i = [2;3;4]; j = [1;4;2];>> a(sub2ind(size(a), i, j))=1a =0 0 0 0 01 0 0 0 00 0 0 1 00 1 0 0 00 0 0 0 0

原创 python matlab 带包实现全排列

>> A=[2,5,7];perms(A)ans = 7 5 2 7 2 5 5 7 2 5 2 7 2 7 5 2 5 7>> perms(1:4)ans =...

原创 【数学建模】线性规划各种问题的Python调包方法

关键词：Python、调包、线性规划、指派问题、运输问题、pulp、混合整数线性规划（MILP）注：此文章是线性规划的调包实现，具体步骤原理请搜索具体解法。  本文章的各个问题可能会采用多种调用方法，为什么？因为这些包各有特点，有些语法特别像matlab，只要稍稍改变即可达成代码交换；而有些包利用了python本身的特性，在灵活度与代码的可读性上更高。我认为这些包各有优劣，各位各持...

原创 数学建模python matlab 编程(随机游走问题)

1(1). 随机游走问题。在-10到10的一维线段上，质点以1/5的概率用左运动1单位，以2/5的概率停止不动，以2/5的概率向右运动2单位，且碰到-10时向右运动3单位，碰到10时向左运动4单位。请画出它的轨迹。我原来的代码matlab：s=0;now=0;for i=1:400 x=unifrnd(0,1); if x<=0.2 ...

原创 B站动手学深度学习第十八课：seq2seq（编码器和解码器）和注意力机制

from mxnet import ndh_forward = nd.array([1,2])h_backward = nd.array([3,4])h_bi = nd.concat(h_forward,h_backward,dim=0)print(h_bi)[1. 2. 3. 4.]<NDArray 4 @cpu(0)>...

原创 R语言与概率统计(六) 主成分分析 因子分析

超高维度分析，N*P的矩阵，N为样本个数，P为指标，N<<PPCA:抓住对y对重要的影响因素主要有三种：PCA,因子分析，回归方程+惩罚函数（如LASSO）为了降维，用更少的变量解决问题，如果是二维的，那么就是找到一条线，要使这些点再线上的投影最大，投影最大，就是越分散，就考虑方差最大。&gt...

原创 R语言与概率统计(五) 聚类分析

#########################################0808聚类分析X<-data.frame( x1=c(2959.19, 2459.77, 1495.63, 1046.33, 1303.97, 1730.84, 1561.86, 1410.11, 3712.31, 2207.5...

原创 R语言与概率统计(四) 判别分析(分类)

Fisher就是找一个线L使得组内方差小，组间距离大。即找一个直线使得d最大。####################################1、判别分析，线性判别：2、分层抽样#install.packages('MASS')library(MASS)#install.packages('sampling')library...

原创 R语言与概率统计(三) 多元统计分析（下）广义线性回归

广义线性回归> life<-data.frame(+ X1=c(2.5, 173, 119, 10, 502, 4, 14.4, 2, 40, 6.6, + 21.4, 2.8, 2.5, 6, 3.5, 62.2, 10.8, 21.6, 2, 3.4, + 5.1, 2.4, 1.7, 1.1, 12.8, 1.2, ...

原创 R语言与概率统计(三) 多元统计分析（中）

模型修正#但是，回归分析通常很难一步到位，需要不断修正模型###############################6.9通过牙膏销量模型学习模型修正toothpaste<-data.frame( X1=c(-0.05, 0.25,0.60,0, 0.25,0.20, 0.15,0.05,-0.15, 0.15, 0.20, 0.10,0.40,0.45,0....

原创 R语言与概率统计(三) 多元统计分析（上）

> #############6.2一元线性回归分析> x<-c(0.10,0.11,0.12,0.13,0.14,0.15,0.16,0.17,0.18,0.20,0.21,0.23)> y<-c(42.0,43.5,45.0,45.5,45.0,47.5,49.0,53.0,50.0,55.0,55.0,60.0)> pl...

原创 R语言与概率统计(二) 假设检验

> ####################5.2> X<-c(159, 280, 101, 212, 224, 379, 179, 264,+ 222, 362, 168, 250, 149, 260, 485, 170)> t.test(X,alternative='greater',mu=225,conf.level = 0.95...

原创 R语言常用包简介

原创 R语言与概率统计(一) 描述性统计分析

查看运行结果：#查看已安装的包，查看已载入的包，查看包的介绍########例题3.1#向量的输入方法w<-c(75.0, 64.0, 47.4, 66.9, 62.2, 62.2, 58.7, 63.5, 66.6, 64.0, 57.0, 69.0, 56.9, 50.0, 72.0)plot(w)#概况,数据的可视化可以让我们看的更轻松...

原创 MATLAB学习（十）实现文件、图像读写

t=1:5;s1=sin(t);s2=cos(t);s=[t;s1;s2];fid1=fopen('test.dat','wt');fprintf(fid1,'\nThis is a Formated file\n');fprintf(fid1,'\n%4d %+12.5E %+12.5E\n',s);fclose...

原创 MATLAB学习（九）系统聚类

>> X=rand(100,2);>> Y=pdist(X,'euclidean');>> Z = linkage(Y,'average');>> dendrogram(Z);>> 　　>> X = [randn(100,2)+ones(100,2)...

原创 MATLAB学习（六）绘图图形功能

>> x=0:.1:2*pi;plot(x,sin(x),x,cos(x))　　>> plot(x,sin(x),'pentagramr',x,cos(x),'.b')　　x=linspace(0,2*pi,30);y=...

原创 MATLAB学习（八）神经网络拟合工具箱 Neural Net Fitting使用示例

>> x=-3:0.2:5;y=x.^2-1;xn=-2:0.1:7;>> 　　>> %多元函数(z=sin(x2+y2)/(x2+y2))拟合>> [X,Y]=meshgrid(-2:0.2:2);Z=sin(X.^2+Y.^2)./(X.^2+Y.^2+eps);>&g...

原创 MATLAB学习（七）求解优化问题：线性规划 非线性规划 拟合与插值 多目标规划...

Minf(x)=-5x1 -4x2 -6x3 x1 -x2 +x3 <=20 3x1 +2x2 +4x3 <=42 3x1 +2x2 <=30 0<=x1,0<=x2,0<=x3>> ...

原创 MATLAB学习（四）线性方程求解，多项式运算，函数局部最优解

>> A=[2 1 3;1 3 2];B=[5;10]B = 5 10>> X=lsqnonneg(A,B)X = 0 2.8571 0.7143>> A\Bans = 0 2.8571 0.7143...

原创 MATLAB学习（三）元素访问和常用代数运算

>> A=[1,2;3,4],B=[0,2;4,5]A = 1 2 3 4B = 0 2 4 5>> C=A>=BC = 1 1 0 0>> D=A==BD =...