数学建模与数学
关注
分享
扫一扫
蔡军帅
好好学习，天天向上！
展开
专栏收录文章
- 默认排序
- 最新发布
- 最早发布
- 最多阅读
- 最少阅读
-
(转载)Python在数学建模中的简单应用
1.方程求导from __future__ import print_functionfrom __future__ import divisionimport numpy as npimport scipy as spimport scipy.misc def f(x): return 2*x*x + 3*x + 1print(sp.misc.derivative(f...原创 2020-01-14 08:50:00 · 641 阅读 · 0 评论 -
【数学建模】线性规划各种问题的Python调包方法
关键词：Python、调包、线性规划、指派问题、运输问题、pulp、混合整数线性规划（MILP）注：此文章是线性规划的调包实现，具体步骤原理请搜索具体解法。 本文章的各个问题可能会采用多种调用方法，为什么？因为这些包各有特点，有些语法特别像matlab，只要稍稍改变即可达成代码交换；而有些包利用了python本身的特性，在灵活度与代码的可读性上更高。我认为这些包各有优劣，各位各持...原创 2019-07-12 12:24:00 · 1906 阅读 · 4 评论 -
PCA的数学原理
PCA的数学原理 原帖地址：http://blog.codinglabs.org/articles/pca-tutorial.html PCA（Principal Component Analysis）是一种常用的数据分析方法。PCA通过线性变换将原始数据变换为一组各维度线性无关的表示，可用于提取数据的主要特征分量，常用于高维数据的降维。网上关于PCA的文章有很多，但是大多...原创 2019-01-26 08:36:00 · 147 阅读 · 0 评论 -
()Python在数学建模中的简单应用
1.方程求导from __future__ import print_functionfrom __future__ import divisionimport numpy as npimport scipy as spimport scipy.misc def f(x):...原创 2020-01-14 08:50:00 · 146 阅读 · 0 评论 -
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 -...原创 2019-07-13 15:17:00 · 994 阅读 · 0 评论 -
数学建模python matlab 编程(指派问题)
指派授课问题现有A、B、C、D四门课程，需由甲、乙、丙、丁四人讲授，并且规定：每人只讲且必须讲1门课；每门课必须且只需1人讲。四人分别讲每门课的费用示于表中：课费用人ABCD甲21097乙154...原创 2019-07-12 14:33:00 · 1641 阅读 · 1 评论 -
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...原创 2019-07-08 11:26:00 · 1083 阅读 · 0 评论 -
MATLAB学习（七）求解优化问题：线性规划 非线性规划 拟合与插值 多目标规划...
Minf(x)=-5x1 -4x2 -6x3 x1 -x2 +x3 <=20 3x1 +2x2 +4x3 <=42 3x1 +2x2 <=30 0<=x1,0<=x2,0<=x3>> ...原创 2019-07-09 10:03:00 · 1506 阅读 · 0 评论 -
matlab学习——05插值和拟合(一维二维插值,拟合)
05插值和拟合1.一维插值(1) 机床加工零件，试用分段线性和三次样条两种插值方法计算。并求x=0处的曲线斜率和13<=x<=15范围内y的最小值。x0=[0 3 5 7 9 11 12 13 14 15];y0=[0 1.2 1.7 2 2.1 2.0 1.8 1.2 1.0 1.6];x=0:0.1:15;% interp1现有插值函数，要求x...原创 2019-08-22 19:05:00 · 1233 阅读 · 0 评论 -
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...原创 2019-07-09 14:09:00 · 7688 阅读 · 0 评论 -
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原创 2019-07-12 14:07:00 · 15078 阅读 · 2 评论 -
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 =...原创 2019-07-12 13:46:00 · 290 阅读 · 0 评论 -
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=...原创 2019-07-09 14:31:00 · 273 阅读 · 0 评论 -
数学建模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 ...原创 2019-07-12 11:26:00 · 1479 阅读 · 0 评论 -
MATLAB学习（一）数组、变量、表达式、常用简单运算
>> x=[1 2 3;4 5 6;7 8 9]x = 1 2 3 4 5 6 7 8 9>> y=[1,2,3]y = 1 2 3>> y=[1,2,34,5,67,8,9]y =...原创 2019-07-08 10:02:00 · 1138 阅读 · 0 评论 -
MATLAB学习（九）系统聚类
>> X=rand(100,2);>> Y=pdist(X,'euclidean');>> Z = linkage(Y,'average');>> dendrogram(Z);>> >> X = [randn(100,2)+ones(100,2)...原创 2019-07-09 14:48:00 · 530 阅读 · 0 评论 -
matlab学习——01线性规划
01线性规划format compact;% min fx% Ax<=b% Aeq*x=beq% lb<=x<=ub%% max z=2x1+3x2-5x3% x1+x2+x3=7% 2x1-5x2+x3>=10% x1+3x2+x3<=12% x1,x2,x3>=0f=[-2;3;5];a=...原创 2019-08-22 15:31:00 · 1536 阅读 · 0 评论 -
matlab学习——02整数规划(蒙特卡洛法,指派问题，混合整数规划)
02整数规划蒙特卡洛法(随机取样法)编写文件mengte.m，目标函数f和约束向量gfunction[f,g]=mengte(x);f=x(1)^2+x(2)^2+3*x(3)^2+4*x(4)^2+2*x(5)^2-8*x(1)-2*x(2)-3*x(3)-x(4)-... 2*x(5);g=[ sum(x)-400 x(1)...原创 2019-08-22 15:39:00 · 1984 阅读 · 0 评论 -
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...原创 2019-07-09 15:34:00 · 447 阅读 · 0 评论 -
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 =...原创 2019-07-08 10:47:00 · 283 阅读 · 0 评论 -
数学建模python matlab 编程(椭圆声学原理画图证明，解析几何)
证明，在椭圆形的音乐厅内，从一个椭圆的一个焦点发出声音，则另一个焦点听到的声音是最大的。分析：证明，从椭圆的一个焦点任意发射的直线经过反射后，并经过另一个焦点。 画图，过一个焦点随机画一条直线l1，再做它与椭圆的交点，过该点作椭圆的一条切线l2，在作一条与该切线垂直的线l3，再作一条l1关于l3对称的直线l4，然后计算另一个焦点到l4的距离小于一个很小的数即可...原创 2019-07-12 23:59:00 · 496 阅读 · 0 评论 -
MATLAB常用快捷键总结
MATLAB 命令栏显示处理的常用命令清屏：clc紧凑显示格式：formatcompact宽松显示格式：formatloose数据高精度显示：formatlongG数据低精度显示：format short编辑器窗口（Editor）下的常用快捷键：自动对齐程序（整理缩进）-自动整理代码 用鼠标选中代码行，按Ctrl + I （我一般先 Ctrl+A...原创 2019-08-21 17:38:00 · 1914 阅读 · 0 评论 -
MATLAB学习（二）读写xls文件
>> N=xlsread('DRINK.xls','DRINK','A1:D8')N = 207.2000 3.3000 15.5000 2.8000 36.8000 5.9000 12.9000 3.3000 72.2000 7.3000 8.2000 2.4000 36.7000 ...原创 2019-07-08 10:21:00 · 694 阅读 · 0 评论 -
matlab学习——05插值和拟合(黄河小浪底调水调沙问题)
05插值和拟合黄河小浪底调水调沙问题data3.txt1800 1900 2100 2200 2300 2400 2500 2600 2650 2700 2720 265032 60 75 85 90 98 100 102 108 112 ...原创 2019-08-24 15:19:00 · 9929 阅读 · 2 评论 -
matlab学习——04图与网络(最短路,最小生成树,最大流)
04图与网络1.最短路(1) 自己写的dijstra算法format compact;clc,clear alla=zeros(6);a(1,2)=50;a(1,4)=40;a(1,5)=25;a(1,6)=10;a(2,3)=15;a(2,4)=20;a(2,6)=25;a(3,4)=10;a(3,5)=20;a(4,5)=10;a(4,6...原创 2019-08-22 15:51:00 · 2486 阅读 · 1 评论 -
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...原创 2019-07-13 15:25:00 · 2171 阅读 · 2 评论 -
LaTex 论文排版(1): Win10 下 LaTex所需软件安装 (Tex live 2018 + Tex studio)
LaTex 论文排版(1): Win10 下 LaTex所需软件安装 (Tex live 2018 + Tex studio) 原创 ...转载 2020-01-12 12:09:39 · 706 阅读 · 0 评论 -
数学建模python matlab 编程(疾病传播模型)
例12：一只游船上有800（1000）人，一名游客不慎患传染病，12（10）小时后有3人发病，由于船上不能及时隔离，问经过60（30）小时，72小时，患此病的人数。(与人口模型和Logistic模型类似)先用python和matlab模拟我的python代码# -*- coding: utf-8 -*-import numpy as npimport randomimpor...原创 2019-07-16 10:35:00 · 2985 阅读 · 0 评论 -
数学建模python matlab 编程(喷泉模拟)
在无风情况下的喷泉模拟我的python代码import numpy as npimport randomimport matplotlibmatplotlib.rcParams['font.sans-serif']=[u'simHei']matplotlib.rcParams['axes.unicode_minus']=Falseimport matplotlib.pypl...原创 2019-07-16 15:32:00 · 1774 阅读 · 0 评论 -
微积分基本公式(牛顿——莱布尼茨公式)的几何解释
牛顿-莱布尼茨公式的内容是一个连续函数在区间 [ a，b ] 上的定积分等于它的任意一个原函数在区间[ a，b ]上的增量。定义原创 2019-10-15 10:52:00 · 5991 阅读 · 0 评论 -
高数(数一)知识点自我归纳(思维导图)
原创 2019-10-28 16:27:00 · 8223 阅读 · 0 评论 -
张宇高数 第一讲 极限(思维导图)
原创 2019-10-29 11:07:00 · 4709 阅读 · 0 评论 -
张宇高数 第二讲 一元函数微分学(思维导图)
原创 2019-10-29 11:10:00 · 2511 阅读 · 0 评论 -
高数 第三讲 一元函数积分学(思维导图)
原创 2019-10-31 13:13:00 · 3645 阅读 · 0 评论 -
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)#概况,数据的可视化可以让我们看的更轻松...原创 2019-07-10 09:37:00 · 2080 阅读 · 0 评论 -
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, ...原创 2019-07-10 15:58:00 · 579 阅读 · 0 评论 -
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....原创 2019-07-10 15:48:00 · 1088 阅读 · 0 评论