- 博客(32)
- 收藏
- 关注
RSS订阅原创 如何利用matlab神经网络进行水量预测
数据收集百度搜索EPS数据库登陆，我们选择进入城市数据库我们的目的是要进行某年水量预测，我们的目的是知道了某一年的土地、人口等信息，就可以预测出该年的用水量。因此我们需要搜集一些信息用来训练。可以多选几年如2003-2017，然后点击下载即可。代码的运行这里给出了下载的链接，把压缩包解压成文件夹并添加到路径即可运行首先是matlab读取文件A=xlsread('区域经济-中国城市数据库.xls','区域经济-中国城市数据库','B4:F20');%2019年需水量未知，用0替代了，后面可以
2021-08-25 00:39:38
1797
3
原创 数学建模安排
有些数学建模队伍讲究分工明确，但是去年做了尝试发现效率太低了。都是一个人干活两个人看，因此还是三个人都有交集比较好。分析去年的问题，A题是微分方程与数值分析，涉及到研究生课程，一般工科生不能做；B题是图论与优化，相对比较好做一些，但是很多人没接触过图论；C题是综合评判与预测类，相对好做一些，有理有据即可。我们组的安排是先看《实用大众线性代数》入门matlab，然后跟着吴恩达的《机器学习》来一边学习算法，一边锻炼对矩阵的理解能力，做一遍作业，即可学会Matlab基本操作。...
2021-07-16 18:39:07
856
原创 学习笔记10 微分方程的matlab符号求解方法
开始学微分方程，用汤家凤的基础课配他的辅导讲义。然后再来看代码。我要国二！！！！代码全打一遍是防止思维跳跃，切忌形式主义尽量做到，学一次，用一年。求符号解定义符号变量调用dsolve函数[y1,...,yN]=dsolve(eqns,conds,Name,Value)eqns为符号微分方程或符号微分方程组；conds为初值条件或边界条件；Name和Value为可选的成对参数。求x2+y+(x−2y)y′=0x^2+y+\left( x-2y \right) y'=0x2+y+(x
2021-04-05 19:20:50
2300
原创 学习笔记9 微分方程的安排
我们学习数学建模的最直接目的是为了拿奖。我认为数学建模当中微分方程是极其重要的，或者说微分方程的建立就是为了建立模型。预计花费一个星期时间或者到月底学习微分方程。学完微分方程以后要达到几个标准。能够对模型进行归纳和类比，方便比赛时能够快速想到（比如2020美赛用房室模型去类比海洋垃圾）学会熟练使用matlab工具求解微分方程符号解和数值解学习微分方程基础知识，融会贯通，会根据实际情况修改模型。符号运算要非常熟练，专有名词要会使用，方便队友沟通前面线性规划不管做得如何都放一边吧，当前的主要阶段是
2021-04-05 19:20:42
188
原创 学习笔记8 数学建模的准备方法（3）
最近三天平均花十个小时左右学习，平均每人每天3小时，对lingo的简单编程有了一定的了解，还学会了一些整数规划、分段函数的技巧。我想将来遇到了线性规划类问题我们有能力相互讨论，也有能力对模型进行优化。但是，任务并没有完成，队友还在研究《数学模型》这本书以及比这本书要求更高的lingo语法。当然，只有Lingo是不够的，还需要进一步学习matlab和python来解决各种算法。预计延长2天，22号前完成，根据队友的情况看看要先学习非线性规划还是微分方程，然后8月考虑看优秀论文以及分工去学习时间序列、数理统计等
2021-04-05 19:20:34
318
原创 学习笔记7 matlab矩阵的运算练习
可以转化为线性规划的问题学到这里，应该要对线性规划问题有一个印象，了解一些将其他问题转化为线性规划问题的技巧。1.2投资的收益和风险模型一matlab代码clc,cleara=0;hold onwhile a<0.05 c=[-0.05,-0.27,-0.19,-0.185,-0.185]; A=[zeros(4,1),diag([0.025,0.015,0.055,0.026])]; b=a*ones(4,1); Aeq=[1,1.01,1.02,1.04
2021-04-05 19:20:25
384
原创 学习笔记6 《数学模型》例题完结
今天看完《数学模型》例题，然后转战《数学建模算法与应用》。课后习题放明后天集中解决，争取按时完成任务。《数学建模算法与应用》这本书代码比较多，难度比较大，我争取把代码都打一遍并存在博客里，更清晰。数学建模算法与应用线性规划问题（Linear Programming）这本书以matlab为主，我们指导老师建议我们minxfTx,\underset{x}{\min}\boldsymbol{f}^T\boldsymbol{x},xminfTx,s.t.{A⋅x⩽b,Aeq⋅x=beq,lb⩽x
2021-04-05 19:20:14
318
原创 学习笔记5 练习
两天没学建模了，后面要加把劲，因为最后都是用电脑打的论文，平时多写一些有好处，博客可以存一些代码+图文（其实jupyter也行，但是我暂时不想用python），往年优秀论文还没开始看。今天主要是练习。线性规划课后习题习题1例题看多了应该知道怎么做了，把模型往上套就行，然后根据Lingo代码分析问题。假设证券名称A,B,C,D,E分别对应 iii = 1，2，3，4，5；信用等级aia_iai，到期年限bib_ibi，到期税前收益cic_ici，总收益zzz。目标自然是赚钱要最多。z=m
2021-04-05 19:20:05
150
原创 学习笔记4 数学建模的准备方法（2）
现已获得代表学校参加9月份全国大学生数学建模竞赛的资格，目标二等奖（3%）。由于这几天队员有其他事分心，还未能完全投入。3天三人总投入时间约862分钟。人均每日投入95分钟，状态不好。指导老师暂定一周讲一次课。除了完成原先的安排，还需赶上老师的进度。老师负责扫盲和指导，修行靠个人。本次课讲的是层次分析法，使用的教材同学习笔记2，同时老师给了我们一些课件用于学习python在数学建模中的应用，并且了解到司守奎老师出了新书。总之，接下来还是依照原来的办法，视频课+《数学模型》+《数学建模算法与应用》，学
2021-04-05 19:19:55
180
原创 学习笔记3 （线性规划）
学习笔记3 （线性规划）选址问题本篇主要注重于lingo代码的认识和使用线性规划课件数学模型视频（清华大学p13-18）依题意列出公式min∑j=12∑i=16cij[(xj−ai)2+(yj−bi)2]12\min \sum_{j=1}^2{\sum_{i=1}^6{c_{ij}\left[ \left( x_j-a_i \right) ^2+\left( y_j-b_i \right) ^2 \right]}}^{\frac{1}{2}}minj=1∑2i=1∑6cij[(xj
2021-04-05 19:19:46
191
原创 学习笔记2 数学建模的准备方法
由于校内可能有选拔赛，所以需要采用速成的办法，即先学习线性规划+lingo。队员基础薄弱，直接看二十几页论文看不懂，因此采用小项目分析+练习的方法。教材《数学模型》，《数学建模算法与应用》，以及一本线性代数教材。配套视频，b站搜索“数学模型”（清华大学）。软件lingo9.0，matlab2017及以后的版本使用方法：暂定11日至19日学习线性规划，观看视频p13-18，先学习《数学模型》视频中的奶制品生产与销售，然后结合视频了解lingo。接着学习《数学模型》的其它例题与习题，结合项目学习ling
2021-04-05 19:19:35
498
原创 学习笔记1（线性规划）
学习笔记1（线性规划）奶制品的生产与销售通过学习奶制品的生产与销售案例入门线性规划（Linear Programming,简称LP），教材《数学模型》（姜启源，谢金星，叶俊）名词记忆：决策变量、目标函数、可行域约束条件一般的表述形式Minzx=f(x)\underset{x}{\text{Min}z}=f\left( \boldsymbol{x} \right) xMinz=f(x)s.t.gi(x)⩽0,i=1,2,⋯ ,m\text{s}.\text{t}.g_i\left( \bo
2021-04-05 19:19:19
474
原创 学习笔记25马氏链模型
前段时间担任协会负责人，在搞matlab培训，又报了一堆竞赛。最近几个星期有空学习和记录。上次国赛全国二等奖的目标没达到，拿了个省一等奖，继续努力。下面给出图中没有显示的代码，给出注释。通过写注释理解代码以及算法clc,clearformat rat %数据格式是有理分数fid=fopen('msdata.txt','r'); %以只读方式打开a=[];while (~feof(fid)) %表示若未读到文件末尾则继续循环 a=[a fgetl(fid)];end%上面
2020-11-29 11:51:20
369
原创 2.1 函数绘图作业
1.立方抛物线x=0:0.1:5;y=x.^(1/3);plot(x,y)2.高斯曲线x=0:0.1:5;y=exp(-x.^2);plot(x,y)3.笛卡尔曲线a为参数，可自行取值t=0:0.1:5;x1=3*t./(1+t.^3);y1=3*(t.^3)./(1+t.^3);subplot(1,2,1)plot(x1,y1)4.蔓叶线在这里插入代码片...
2020-11-11 19:41:55
261
原创 极限与方程组
syms xlimit(F,x,a)limit(F,x,inf)limit(F,x,a,'right')limit(F,x,a,'left')例2.2.1clear;x=-1:0.01:-0.01;x3=10:500;y1=sin(1./x)y2=sin(x)./xy3=(1+1./x3).^x3;y4=2.71828;figure(1),plot(x,y1)figure(2),plot(x,y2)figure(3),plot(x3,y3,x3,y4)例2.2.3sym
2020-11-11 17:33:09
190
1
原创 2.1 一元微积分之曲线绘图
11月共计四个星期，针对大一的学生进行培训，主要让他们熟悉matlab基本语法，同时结合他们目前所学的微积分开展相关的教学。以《数学实验》2.1 曲线绘图
2020-11-03 18:38:44
757
原创 学习笔记24 加强阶段总结与未来规划
总结在过去一个月的加强阶段，我们相互讲评了论文，效果不佳。后面采取多种方法学习，并没有一个固定的模式，主要是积累模型，我们积累了一些微分方程模型、评价模型、预测模型。在9月的模拟赛中，我认为我们才算成功地入门数学建模，能够独立地完成一篇论文，基于我所在的二本学校学校炮灰众多，我想省一等奖是一定有的。但是我的目标是全国二等奖，获奖率3%。我们选择了华数杯C，评价类的题目，从头到尾就是文件处理、层次分析、灰色关联，效果不佳。而且后面踩了两个坑，模型评价没写，还有附件太大，下次要注意。未来三天任务到时候比
2020-09-08 00:42:34
117
原创 学习笔记23 神经网络
clc, cleara=load('jingliu.txt'); %把表中第2列到第6列的数据保存到纯文本文件jingliu.txta=a'; %注意神经网络的数据格式，不要把矩阵搞转置了。P=a([1:4],[1:end-1]); [PN,PS1]=mapminmax(P); %自变量数据规格化到[-1,1]T=a(5,[1:end-1]); [TN,PS2]=mapminmax(T); %因变量数据规格化到[-1,1]net1=newrb(PN,TN) %训练RBF网络x=a([1:4]
2020-09-06 22:19:37
210
原创 学习笔记22 灰色关联
clc,cleara=[0.83 0.9 0.99 0.92 0.87 0.95 326 295 340 287 310 303 21 38 25 19 27 10 3.2 2.4 2.2 2.0 0.9 1.7 0.2 0.25 0.12 0.33 0.2 0.09 0.15 0.2 0.14 0.09 0.15 0.17 250 180 300 200 150 175 0.23 0.15 0.27 0.3 0.18 0.26 0.87 0
2020-09-06 10:57:27
165
原创 管网入门1
Combined System图1 雨污合流的下水道这是最常见的管道系统缺点：降水超过一定比例后送去污水厂，污水厂没有办法处理那么多，就流到水体了。Separate System图2 雨污分流的下水道避免了雨污合流缺点：要两套管子，距离长Modified Separate System图3 改进版雨污分流的下水道节约了成本...
2020-09-02 21:34:55
113
原创 学习笔记21 有限差分法求解热传导方程的代码
距离9月份全国大学生数学建模竞赛只剩10天，然而，进度缓慢。有限差分法这部分我可能卡了一个星期，不过好在现在并不是毫无头绪，理解了算法，现在阅读代码更加容易了。未来的准备计划。2018A搞定以后进入冲刺阶段。在加强阶段我们学习了不少模型，也阅读了一部分论文，现在积累了一些评价模型、预测模型、微分方程模型，这样已经有实力解决部分问题了。9月4日至6日要找一个赛题练习一下，查缺补漏。然后再看看自己缺什么。废话不多说，下面贴上matlab代码(同样是在b站上看到的，把它打下来，跟学习笔记20的是一套的)。c
2020-09-01 14:14:46
6870
7
原创 学习笔记20 热传导方程
之前没学过数值分析，对于误差分析、偏微分方程求解什么的一窍不通，但是看队友微分方程建模学得比较好，解模的能力需要跟上。但是B站上面有人专门做了视频讲解。有点长，在此以文字形式保存下来。公式介绍符号意义uuu温度xxx传播的位移ttt时间jjjx离散化后的序号jjj=1,2,…,N+1nnnt离散化后的序号，n=1,2,...n=1,2,...n=1,2,...fff热源带来的温度影响热传导方程∂u∂t(x,t)=α2∂2∂x2u(x,
2020-08-31 18:45:39
4420
1
原创 学习笔记19 关于数学建模的思考与总结
今天是8月27日，距离开赛还有14天。原本计划11日到24日没人看7篇文章，但是困难重重，所以只能放慢进度，慢慢消化。当前队内A讲解了2012B太阳能小屋，B讲解了2019A高压油管（但是效果不佳，这题很难，队友不能理解），C讲解了2010B上海世博会。之前由于三人实力相差不大，就一起学习，由于队友提出要分工，我认为有道理且能增加效率，所以现在A主要负责论文撰写和排版，B主要负责建模，C主要负责模型的求解和代码实现。很显然完全分工是做不到的，遇见模型要求每个人都能理解才能更好地合作。下面引用一段杨振宁
2020-08-27 15:34:13
289
原创 学习笔记18 灰色预测
之前李林老师给我们讲了灰色预测的原理与案例。但是作为队内主要负责代码求解的人，从模型的建立到求解都要熟练，或者最起码比赛的时候能在博客里直接使用从前写过、注释过的代码。谢金星老师的《数学模型》模型求解较少，所以今后会少看。近期主要看司守奎老师的《数学建模算法与应用》和往年优秀论文。预测模型clc,clearx0=[71.1 72.4 72.4 71.4 72.0 71.6]';n=length(x0);lamda=x0(1:n-1)./x0(2:n)%错位相除计算级比range=minmax(
2020-08-27 14:15:38
277
原创 学习笔记17 差分方程
仅仅看论文是不够的，之前的书也要看。之前华数杯的时候选A题犯了难，因为不会解偏微分方程，这次选了2018A来看别人是怎么列差分方程求解偏微分方程的，但是昨天看了3小时也没看懂，所以退而求其次，看看司守奎老师的书，顺便积累一些代码练习。今天主要看《数学建模算法与应用》的差分方程，p409。代码练习线性方程假设我们有一些离散的（xxx,yyy）已知 xxx和yyy的几个数值，知道他们成线性关系y=ax+by=ax+by=ax+b(极特殊的情况，条件十分苛刻，因为现实中很难全部的点都在线上)，求aaa，b
2020-08-25 18:37:13
1658
原创 学习笔记16 数学建模的正确学习方法（进阶的层次分析法）
感觉前面走了不少弯路，低估了自己的智商。现在已经有实力去看论文了。小组内做了合理的分工，共学算法，我侧重于编程这一块。我既要会python，也要会matlab，现在看来，相对于上个月，编程其实也不难，有所进步。好在还有20天，抓紧时间，还能把十年论文过一遍。不懂的算法可以再学一遍。作为一名成功的参赛选手，把问题简化且告诉队友是必要的，也就是说既要看懂很难的专业术语，也要能把问题简化教给别人。以学习笔记14为例。论文链接在那我们聊聊如何用层次分析法和模糊数学。题目一句话概括，让我们从感兴趣的方面建立
2020-08-21 00:40:35
505
原创 学习笔记15 python绘图进阶（一个失败的学习方案）
箱线图# -*- coding: utf-8 -*-import numpy as npimport matplotlib.pyplot as pltplt.style.use('ggplot')N=500normal=np.random.normal(loc=0.0,scale=1.0,size=N)lognormal=np.random.lognormal(mean=0.0,sigma=1.0,size=N)index_value=np.random.random_integers(lo
2020-08-17 14:47:56
150
2
原创 学习笔记14 上海世博会的评价
原计划11号开始，结果到16号我才看完了一篇看了上海世博会的评价。重点看算法和代码。这篇文章写得通俗易懂，没有很难的公式推导，尤其是第二个模型。在这里保存打出来的代码，以及我自己加的注释论文需要的自取提取码ygh3程序一clc,clearb=[1 2 4 1/2 1/2 1,2,1/3 1/4 1/2 1 1/5 2 3 5 1][x,lumda]=eig(b)%求特征值和特征向量，A的特征值构成对角矩阵lumda,按特征值绝对值大小排序；%A的特征向量构成矩阵x
2020-08-16 18:08:11
164
原创 学习笔记12 python绘图
今天把基本绘图给学会了条形图# -*- coding: utf-8 -*-import matplotlib.pyplot as pltplt.style.use('ggplot')#使用ggplot样式customers=['ABC','DEF','GHI','JKL','MNO']customers_index=range(len(customers))获取索引sale_amounts=[127,90,201,111,232]#以上几行为准备数据fig=plt.figure()#这是一个
2020-08-06 15:40:44
156
原创 学习笔记11 Python之进阶
电脑已经安装好anaconda，最近先用anaconda，后面再学习jupyter。学习教材《python数据分析基础》时间的代码from math import exp,log,sqrtimport refrom datetime import date,time,datetime,timedeltatoday=date.today()print("Output #41:today:{0!s}".format(today))print("Output #42:today:{0!s}".f
2020-08-05 17:13:18
89
原创 Spyder编译器的数据存储
Spyder编译器的数据存储Spyder的界面右上角的白色框就是当前Ipython操作的文件夹，可以根据框内按钮前进或后退切换文件夹，Ipython的命令行如添加或读取文件将会在这个文件夹内执行。然后在左下角黑框内我们可以操作文件。命令行式修改路径import os os.getcwd() #获取当前路径 os.chdir("D:\\py作业") #跳到目标路径下的D盘文件夹...
2020-03-20 21:18:07
5314
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人