wulimmya的博客

原创 错误使用 subplot 需要使用有效的轴句柄作为输入，有点意思

matlab画图subplot(211)表示2行1列第一个图，结果今天程序需要画400多个图，到第210个图画完就出错了，错误使用 subplot (line 171)需要使用有效的轴句柄作为输入。出错 gps_interp (line 37) subplot(211)原来是因为211！！！！这时候就必须写成subplot(2，1，1)拉！！哈哈哈有意思...

原创 小波分析,从傅里叶变换到小波变换

文章目录先说变换是什么再说傅里叶变换有什么局限简述FTFT因无法处理时变信号而裹足不前STFT为了克服FT的缺点（不能处理时变信号）诞生了小波变换为了克服STFT的缺点（单一分辨率）诞生了仔细说说小波变换的真正含义小波是什么（主要讲母小波）母小波伸缩平移得到子小波连续小波变换先说变换是什么变换就是定义域和值域都是函数集而非数值的函数。 它仍然是一种映射，只是从函数集映射到函数集。比如傅里叶变...

原创 多元线性回归 stata

回归的思想，任务，使命可以看到0-1回归的因变量自变量都是0-1变量因变量y可以分为解释y的形成机制就是通过建立x，y的数学关系式感觉相关分析就是回归分析的前两步，数据的分类...

原创 4种聚类（kmeans, kmeans++，层次聚类，DBSCAN密度聚类）

k-meansk-means++系统聚类/层次聚类

原创 数学建模 典型相关分析Canonical Correlation analysis

典型相关分析CCA和主成分分析PCA非常相似，都是通过现有变量的线性组合构造新的综合变量以实现减少变量数目但代表几乎同样多的信息的目的。只是CCA分析的是两组变量，是要把每一组变量用一个线性组合综合变量表示，从而把两组变量的相关转化为两个综合变量（典型变量）的相关；而PCA是针对一个数据的多组变量，比如一张图片的所有像素，一个人的各种特征（身高，体重，肺活量···），是去把这一组变量实现降维就行了...

原创 matlab 画图总结

markersize,用圆圈画每个点，离散的点之间不连线，不用圈圈的话plot默认把每个离散点用直线连起来linewidth，线宽MATLAB是不可以画出连续函数的图像的，只能用一群离散点（划分更细）去在视觉上像连续的。x1=(0:12)*pi/6;y1=cos(3*x1);x2=(0:360)*pi/180;y2=cos(3*x2);figure(1)subplot(221);pl...

原创 遗传算法优化 matlab代码（不用工具箱）

先画出函数图像：用figure1上工具中的数据游标手动大致估计到最大值的坐标是（1.5,185.1）代码：clear all close allclcx=-2:.1:2;f=200*exp(-.05*x).*sin(x);%.*表示两向量对应元素相乘plot(x,f)grid on...

原创 遗传算法及其matlab实现（示例用python库geatpy，matlab遗传工具箱geatbx和lingo三种方式求解）

智能算法是一个边缘交叉学科，有蚁群、粒子群算法、遗传算法、免疫算法、模拟退火等，都乐意算作机器学习的一种类别。文章目录理论生物理论基础历史算法编解码编码解码染色体交配运算产生子代染色体基因突变运算产生子代染色体染色体倒位运算产生子代染色体个体适应度评估复制运算示例遗传算法是模拟达尔文进化论和孟德尔遗传学原理的生物计算模型，是一种模拟自然进化过程以搜索最优解的方法。理论生物理论基础先列...

原创 matlab二维插值绘制地貌图

clearclc[x,y]=meshgrid(1:10);% 测量数据h=[0 .02 -.12 0 -2.09 0 -.58 -.08 0 0; .02 0 0 -2.38 0 -4.96 0 0 0 -.1; 0 .1 1 0 -3.04 0 -.53 0 .1 0 ; 0 0 0 3.52 0 0 0 0 0 0; -.43 -1.98 0 0 ...

原创 数学建模 非线性规划

一个很好的例子，列出优化模型和里面的约束条件求解：示例求解结果(x取（1,1）时获得最小值10)>> mainLocal minimum found that satisfies the constraints.Optimization completed because the objective function is non-decreasing in ...

原创 数学建模 二次规划(普通求解方式 && 罚函数法)

matlab一般求解方式示例可以看出12xTHx=2x12−4x1x2+4x22\frac12x^THx=2x_1^2-4x_1x_2+4x_2^221​xTHx=2x12​−4x1​x2​+4x22​即[x1x2][abcd][x1x2]=ax12+(b+c)x1x2+dx22=4x12−8x1x2+8x22\left[ \begin{matrix} x_1 &...

原创 模拟退火算法理论及其MATLAB实现

像遗传算法，模拟退火，免疫算法，粒子群算法，神经网络，蚁群算法，禁忌搜索这样的智能算法，都是独立于问题的算法。对于TSP这种NP难问题，不可能有解析的确定的算法去求解全局最优解，只能依靠这种随机搜索算法去求解局部最优解。这个 “独立于问题” 描述地特别到位，精辟，道出了本质。...

原创 数学建模 灰色预测

原创 数学建模准备之 微分方程建模和求解

微分方程建模一阶线性齐次常微分方程的通解求解：一阶线性非齐次常微分方程的通解求解：常数变易法:问题一定非常现实，每个都看起来很难很难无法解决，而且很可能是别的学科的知识，关键是抽象，抽象的关键是去找出问题中的变量以及变量之间的关系，这样就可以列出式子了，就建模了。这是卓金武那本书，薄膜渗透率的测定，根据A测增加的质量等于B侧减少的质量建立一个等式，并通过让变化时间趋近于0得...

原创 excel根据数据画饼状图等

把Word的数据粘贴到excel在Excel中的部分数据在一边计算各类病总数选择二维饼图改变标题更改图表类型可以选择其他图三维饼图圆环图 旭日图柱状图条状图折线图面积图雷达图...

原创 数学建模准备 插值（拉格朗日多项式插值，牛顿多项式插值，分段线性插值，分段三次样条插值，分段三次Hermite插值）

文章目录一、基础概念插值是什么拟合是什么插值和拟合的相同点插值和拟合的不同点二、常用的基本插值方法高次多项式插值拉格朗日多项式插值牛顿插值差商矩阵低次多项式插值（不易震荡）分段线性插值Hermite插值三次样条插值一、基础概念插值是什么求出过已知个有限数据点的近似函数。拟合是什么求出一个不要求过已知数据点的近似函数，不要求过那些数据点，只要求在这些点上的总偏差最小。插值和拟合的相同点...

原创 Technique for Order Preference by Similarity to Ideal Solution逼近理想解排序法

文章目录基本过程归一化真是个好东西四种指标以及相互转换的公式excel数据做成MATLAB的mat文件简直超easy超好理解看到一半又对数模有了一点信心基本过程得到正向化矩阵----得到标准化矩阵----计算评分里面的计算确实超简单的，只有加加减减平方啥的这个想法确实不错，和最低成绩的插值占最大差距的比例，这在分段线性插值的理论公式中也有应用举一反三：可以立刻想到另一个评分公式...

原创 数学建模 评价类问题， 层次分析法（AHP）

找评价指标知网找相关文献小组头脑风暴层次分析法的判断矩阵（正互反矩阵）一致矩阵

原创 axmath && mathpix snipping tool配合神器，截图得公式（通过latex代码）

axmath下载地址Mathpix下载地址下载安装完成后，axmath是需要注册的，淘宝购买注册文件.lic才34，终身使用，非常划算的使用mathpix 截图想要的公式（ctrl+alt+m），复制下面的latex代码到注册后的axmath中（打开Word，上面有axmath插件）得到的公式是图形，不容易改动，非常快捷便利...

原创 18年研究生国赛 F题 机场新增卫星厅对中转旅客影响的评估方法

终端厅分为航站楼和卫星厅，航站楼的功能包括出发、到达、出入境和候机，而卫星厅作为航站口的延伸，不能出入境。分配的是转场（航班对），飞机总数并不重要。卫星厅是一种特殊的航站楼，设计初衷是缓解登机口不足带来的分配压力，但不可避免地引入了中转距离和中转时间的增加，增加了中转旅客的换乘紧张程度和换乘压力。根据题目给出的中转最短流程时间表格可以看出，无论是国内中转国内还是国内中转国际，需要的中转航站楼之...

原创 Python学习 字符串拼接，pycharm再探索，记录一个异常，Python高维矩阵切片索引

parts = ['Is', 'Chicago', 'Not', 'Chicago?']' '.join(parts)Out[21]: 'Is Chicago Not Chicago?'','.join(parts)Out[22]: 'Is,Chicago,Not,Chicago?'''.join(parts)Out[23]: 'IsChicagoNotChicago?'平时shi...

原创 anaconda初体验

LINUX和mac系统一般都自带一个Python的命令行版本，如果需要图形界面的版本（IDLE(Python GUI)），需要自己另外安装。但是Windows系统默认不提供任何Python版本，这无所谓，自己装一个就是了，反正是自由下载的开源软件，很简单。Python可以跑的平台很多，mac,类unix系统（比如Linux），Windows，sun solaris, os/2, amiga, ...

原创 pickle模块学习和一次惨痛（sb）的debug经历

pickle模块实现了把Python对象结构进行序列化和反序列化（ serializing and de-serializing a Python object structure）的二进制协议。pickling: 把对象序列化为字节流a byte streamunpickling: 把二进制文件binary file的字节流转换为对象JSON： 文本序列化格式a text seriali...

原创 皮尔逊相关系数，斯皮尔曼等级相关系数，(易错！！)假设检验 ，SPSS

本文参考了维基和minitab网站，只是汇总梳理和强调了相关知识点。Spearman’s rank correlation coefficient，是衡量两个变量的依赖性的 非参数 指标（衡量分级定序变量之间的相关程度），利用单调方程评价两个统计变量的相关性。 如果数据中没有重复值， 并且当两个变量完全单调相关时，斯皮尔曼相关系数则为+1或−1。用途是描述数据间的相关性。可以看到斯皮尔...

原创 隐枚举法求解0-1整数规划

隐枚举法就是在 枚举过程中引入过滤条件。0-1整数规划是整数规划的特殊情形，部分或全部变量只能取0或者1, 这个约束表示为0≤xj≤1,xj∈N,j=1,2,⋯n0\leq x_j\leq1, x_j \in N, j=1,2,\cdots n0≤xj​≤1,xj​∈N,j=1,2,⋯n或者表示为xj(1−xj)=0,j=1,2,⋯nx_j(1-x_j)=0, j=1,2,\cdots ...

原创 matlab 数据滤波 (传递函数)

t=linspace(0,10,100);s=sin(2*pi/5*t);noise=0.2*rand(size(t));%xnoise = x + 2.5*gallery('normaldata',size(t),4);%标准正态分布的随机数据矩阵xs=s+noise;y=zeros(size(t));A=[1 -0.9];B=[.05 .06];y=filter(B,A,x...

原创 texstudio中使用Bibtex插入参考文献

最近用一个期刊的latex模板写论文，遇到插入参考文献的问题，以前一直是个MSworder，第一次用latex，搞了好半天才明白到底咋回事.bst文件，是参考文献的格式/风格说明文件，期刊模板中一并提供.bib文件，参考文献的源文件，自己需要引用的参考文献都要放在这里面，.bib文件需要自己创建，创建方法很简单，新建一个txt文本文档，用文献管理工具比如mendeley，右击想引用的文献，...

原创 latex学习 day3 elsarticle.cls学习

elsarticle.cls要用到的一些包natbib.sty 处理引用 for citation processing，处理各种各样的引用，和hyperref.sty，hypernat.sty一起完美工作geometry.sty 页边距设置 for margin settingsfleqn.clo 左对齐公式 for left aligned equationsgraphicx.st...

原创 latex学习 day2

C:\Users\Administrator>texdoc texdocC:\Users\Administrator>texdoc article第一句代码直接打开了文件texdoc.pdf, 路径D:\texlive\2019\texmf-dist\doc\support\texdoc，这里还有编译生成texdoc.pdf文件的源码texdoc.tex第二句代码直接打开了cl...

原创 texstudio深色主题设置

深色主题配置代码下载地址下面是我搬运的一个示例，我用的就是这个用everything工具（下载地址）查找texstudio.ini,双击打开，翻到最下面，在texstudio未打开使用的情况下，用下面的代码替换下图红框开始的所有代码，然后打开texstudio，就搞定啦[formats]version=1.0data\normal\priority=-1data\normal\bold...

原创 win 10局域网共享文件夹设置完整教程

在教研室用台式机，但有的东西在笔记本上，用u盘考来考去很麻烦，于是弄了个共享文件夹，非常方便先在你喜欢的位置新建一个文件夹，专门用于共享文件，我在F盘下新建的对新建的用于共享的文件夹右键，属性，共享权限改为读取/写入点击上图右下角的共享进入网络和共享中心在此电脑界面空白处右键，添加一个网络位置win+R cmd进入命令行，输入ipconfig，查看ipv4地址，...