自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+

LSEC小陆的博客

生于Code,死于Bug

  • 博客(34)
  • 资源 (13)
  • 论坛 (1)
  • 收藏
  • 关注

原创 几个R语言画图小程序分享

三维图x=seq(-5,5,by=0.1) #步长很小时画的图就是黑色的了,因为都是画格子的黑线的颜色y=xx1=dnorm(x,0,1) #dnorm()为正态分布密度函数z=outer(x1,x1)persp(x,y,z,theta =30,phi = 25,expand = 0.5,col = "Blue2")等高线图library(MASS)Sigma <- matrix(c

2017-08-16 14:02:52 4525

原创 R语言基础练习与入门实践

提供几个R语言的入门练习,让同学们更快地了解和掌握R语言编程。练习: 熟练使用R软件实践1:最初几步x=1:100#把1,2,...,100个整数向量赋值到x(x=1:100) #同上, 只不过显示出来sample(x,20) #从1,...,100中随机不放回地抽取20个值作为样本set.seed(0);sample(1:10,3)#先设随机种子再抽样.#从1,...,200000中随机

2017-08-16 13:44:37 10240 1

原创 假设检验与参数估计的R语言实现

假设检验是推论统计中用于检验统计假设的一种方法。而“统计假设”是可通过观察一组随机变量的模型进行检验的科学假说。一旦能估计未知参数,就会希望根据结果对未知的真正参数值做出适当的推论。统计上对参数的假设,就是对一个或多个参数的论述。而其中欲检验其正确性的为零假设(null hypothesis),零假设通常由研究者决定,反应研究者对未知参数的看法。相对于零假设的其他有关参数之论述是备择假设(alt

2017-08-16 13:08:57 4376

原创 曲线拟合

曲线拟合,俗称拉曲线,是一种把现有数据透过数学方法来代入一条数式的表示方式。科学和工程问题可以通过诸如采样、实验等方法获得若干离散的数据,根据这些数据,我们往往希望得到一个连续的函数(也就是曲线)或者更加密集的离散方程与已知数据相吻合,这过程就叫做拟合 (fitting)。交通和aqi拟合w=read.table("日拥堵指数和AQI数据.txt")plot(w$V1,w$V2,pch=20,x

2017-08-16 12:43:34 645

原创 R语言做回归分析

回归分析(英语:Regression Analysis)是一种统计学上分析数据的方法,目的在于了解两个或多个变数间是否相关、相关方向与强度,并建立数学模型以便观察特定变数来预测研究者感兴趣的变数。更具体的来说,回归分析可以帮助人们了解在只有一个自变量变化时因变量的变化量。一般来说,通过回归分析我们可以由给出的自变量估计因变量的条件期望。回归的最早形式是最小二乘法,勒让德和高斯都将该方法应用于从天文

2017-08-16 12:34:09 11405

原创 相关性分析

在概率论和统计学中,相关(Correlation,或称相关系数或关联系数),显示两个随机变量之间线性关系的强度和方向。在统计学中,相关的意义是用来衡量两个变量相对于其相互独立的距离。在这个广义的定义下,有许多根据数据特点而定义的用来衡量数据相关的系数。对于不同测量尺度的变数,有不同的相关系数可用:Pearson相关系数(Pearson’s r):衡量两个等距尺度或等比尺度变数之相关性。是最常见的

2017-08-16 11:57:13 1983

原创 中心极限定理

中心极限定理是概率论中的一组定理。中心极限定理说明,大量相互独立的随机变量,其均值的分布以正态分布为极限。这组定理是数理统计学和误差分析的理论基础,指出了大量随机变量之和近似服从正态分布的条件。中心极限定理有着有趣的历史。这个定理的第一版被法国数学家棣莫弗发现,他在1733年发表的卓越论文中使用正态分布去估计大量抛掷硬币出现正面次数的分布。这个超越时代的成果险些被历史遗忘,所幸著名法国数学家拉普拉

2017-08-16 11:51:28 657

原创 正态分布

神说,要有正态分布,就有了正态分布。 神看正态分布是好的,就让随机误差服从了正态分布。关于正态分布的由来,推荐文章正态分布的前世今生。下面是几个R语言的小demo。二项分布的正态近似set.seed(234)x=rbinom(10000,100,0.6)x=matrix(x,nc=4)m=apply(x,1,mean)summary(m)var(m)sd(m)hist(m,

2017-08-16 11:42:31 643 1

原创 A*算法:Dijkstra改进算法

A*算法:Dijkstra改进算法A算法Dijkstra改进算法寻路模型A算法算法简介算法步骤格点模型格点模型分析改进思路方向导航网络算法步骤及其优缺点A算法启发函数的选取伪代码参考外链寻路模型最优路规划问题处理的核心在于寻路地图的建模以及最优路规划模型。A*算法算法简介A*搜索算法,俗称A星算法。这是一种在有多个节点的路径的图形平面上,求出最低成本的路径算法。A*算法是一种基

2017-08-16 10:50:11 10784 1

原创 函数估计的非参数方法

用于函数估计的非参数方法主要有核密度估计、局部多项式回归估计等。非参的函数估计的优点在于稳健,对模型没有什么特定的假设,只是认为函数光滑,避免了模型选择带来的风险;但是,表达式复杂,难以解释,计算量大是非参的一个很大的毛病。所以说使用非参有风险,选择需谨慎。核密度估计(kernel density estimation)是在概率论中用来估计未知的密度函数,属于非参数检验方法之一,由Rosenblat

2017-08-16 09:34:20 2348 1

原创 一个c程序

作为入门,分享一个c程序,主要功能是输入一个数,判断其为几位数,并且输出各个数位上的数。# include <stdio.h># include <math.h>int main() { int a,axxxx,b,bxxx,c,cxx,d,dx,e,f; scanf("%d",&f); e=f%10; dx=f%100; d=(dx-e)/10;

2017-08-15 21:26:08 462

原创 Dev-C++ 5.10的调试功能

Dev-C++ 5.10的调试功能 安装好Dev-C++之后,要在工具-编译选项-代码生成/优化-连接器-产生调试信息处改成“yes”,才能进行调试。删繁就简,调试步骤如下:写好程序,点击代码窗口左边的数字设置若干断点。 按F5或者点击工具栏上的对号进行调试运行,程序运行到第一个断点处,底部窗口自动打到调试选项卡。 添加查看对变量进行监视,不想打字,有个方便的方法是用鼠标选择变量,邮件点击添

2017-08-15 21:16:41 2790

原创 百度贴吧爬虫python示例

引用一个简单的Python百度贴吧爬虫示例程序,功能为下载对应页码内的所有页面并存储为html文件,算是一个入门吧。# -*- coding: utf-8 -*- #--------------------------------------- # 程序:百度贴吧爬虫 # 版本:0.1 # 作者:joker# 日期:2017-05-14 # 语言:Pyth

2017-08-14 22:15:50 493

原创 python查看如何查看版本信息

问题如题,解决方案为编写脚本如下,运行。#!/usr/bin/pythonimport sysprint sys.version print sys.version_info

2017-08-14 21:07:45 11909 2

原创 Python中re的用法示例

re模块的介绍点这里,不再详述。给出一个简单的re实例,匹配字符串中的hello字符串 ,希望各位看官能从中悟出点什么。# -*- coding: utf-8 -*- #一个简单的re实例,匹配字符串中的hello字符串 #导入re模块 import re # 将正则表达式编译成Pattern对象,注意hello前面的r的意思是“原生字符串” pattern = re.compile

2017-08-14 20:57:55 1290

原创 Python编写的一个简单的石头剪刀布游戏

还记得“石头剪刀布”吗?这项此前被讥为“吃饱了撑的”的研究,入选了“麻省理工学院科技评论2014年度最优”,成为中国首次入选“麻省理工科技评论”的社科领域成果。      这项由浙江大学、浙江工商大学、中国科学院组成的跨学科团队共同完成的成果,叫作“石头剪刀布”中的社会循环与条件响应,解释了“剪刀石头布”的制胜原理。“石头剪刀布”研究采用经济学控制性实验方法—也就是实验经济学的方法,利用“石头剪

2017-08-14 20:46:59 20004 2

原创 模糊综合分析程序

模糊综合评价法是一种基于模糊数学的综合评价方法。该综合评价法根据模糊数学的隶属度理论把定性评价转化为定量评价,即用模糊数学对受到多种因素制约的事物或对象做出一个总体的评价。它具有结果清晰,系统性强的特点,能较好地解决模糊的、难以量化的问题,适合各种非确定性问题的解决。clcclearclose alldisp('请输入判断矩阵A(n阶)'); A=input('A='); [n,n]=s

2017-08-13 22:29:19 4414 1

原创 Matlab算累次积分的一个例子

Matlab符号运算能力并不如Mathematica,但是并不说明其在符号计算领域一无是处,下面给出一个Matlab求累次积分的小例子。clcclearsyms x1 x2 x3;syms f1 f2 f3;L1=1-x1-x2;L2 = x1;L3 = x2;g = ((f1*L1+f2*L2+f3*L3)*L1);int(int(g,x2,0,1-x1),x1,0,1)

2017-08-13 22:23:18 1882

原创 MATLAB做曲线拟合

实际工作中,变量间未必都有线性关系,如服药后血药浓度与时间的关系;疾病疗效与疗程长短的关系;毒物剂量与致死率的关系等常呈曲线关系。曲线拟合(curve fitting)是指选择适当的曲线类型来拟合观测数据,并用拟合的曲线方程分析两变量间的关系。公式推导:MATLAB内置函数做拟合:汽车保有量预测线性拟合程序:clearclcclfx=2005:2014y=[764312 985445 10

2017-08-13 22:15:45 1601

原创 层次分析法MATLAB程序

作为一个数学人,对于层次分析法,我的评价同下: 其主观作用太强,例如,那个矩阵的由来就无根据,仅仅是由若干尺度表主观推出来的,并不是真正定量的严谨的计算推理,类似的还有模糊评价里的专家评定,虽无力吐槽,作为一个知识,还是贴一段代码。clearclca=[1 3 5;1/3 1 3;1/5 1/3 1];b1=[1 1/3;3 1];b2=[1 5;1/5 1];b3=[1 3;1/3

2017-08-13 22:09:40 25214 21

原创 最短路算法M程序

最短路径问题是图论研究中的一个经典算法问题,旨在寻找图(由结点和路径组成的)中两结点之间的最短路径。算法具体的形式包括:1、确定起点的最短路径问题:即已知起始结点,求最短路径的问题。2、确定终点的最短路径问题:与确定起点的问题相反,该问题是已知终结结点,求最短路径的问题。在无向图中该问题与确定起点的问题完全等同,在有向图中该问题等同于把所有路径方向反转的确定起点的问题。3、确定起点终点的最短路径

2017-08-13 22:00:33 459

原创 一个DEA程序模板

data envelopment analysis的缩写,即数据包络分析。《数据包络分析》(DEA)是一本关于数据包络分析(DEA)方法、模型和理论的专著,是作者十几年工作的总结。 DEA有效性与相应的多目标规划问题的pareto有效解(或非支配解)是等价的。数据包络分析(即DEA)可以看作是一种统计分析的新方法。它是根据一组关于输入-输出的观察值来估计有效生产前沿面的。1978年由著名的运筹学家

2017-08-13 21:53:41 2801

原创 BP神经网络代码示例

BP神经网络 BP(back propagation)神经网络是1986年由Rumelhart和McClelland为首的科学家提出的概念,是一种按照误差逆向传播算法训练的多层前馈神经网络,是目前应用最广泛的神经网络。发展背景在人工神经网络的发展历史上,感知机(Multilayer Perceptron,MLP)网络曾对人工神经网络的发展发挥了极大的作用,也被认为是一种真正能够使用的人工神经网

2017-08-13 21:45:17 22514 14

原创 灰度模型与灰色预测

概念灰色系统是相对于黑色系统和白色系统而言的。白色系统:系统的内部特征是完全已知的,即系统的信息是完全充分的。黑色系统:一个系统的内部信息对外界来说是一无所知的,只能通过他与外界的联系来加以观测研究。灰色系统:一部分信息是已知的,另一部分信息是未知的,系统内各因素之间具有不确定关系。其特点是‘少数据建模’,着重研究‘外延明确,内涵不明确’的对象。 灰色系统具有相对性与广泛性。指系统对于不同对象

2017-08-13 21:13:03 10388

原创 markdown编辑器使用技巧

各级标题 “#标题一” “##标题二” “###标题三”尾注 “尾注” 在尾巴写上:[1]: http://www.baidu.com/列表与加粗 “- **短横是列表,双星号是加粗。**”分割线 “——————-”(一簇短横,快捷键是 Ctrl + R)引用与插入链接 “> 这里是引用内容” h5语言式插入链接:<a href="https://zh.wikipedia

2017-08-10 21:22:56 739

原创 径向Kohn-Sham方程的谱有限元方法

径向Kohn-Sham方程的谱有限元方法本文以Kohn-Sham方程为例,介绍了现在已发展成型的几种近似交换相关泛函和密度泛函理论的计算步骤,主要以代码的形式呈现。径向密度泛函理论给出了求解下面一个Kohn-Sham方程的一套思路: 其中,VxcV{xc}是近似交换相关泛函,Vn(r)=ZrV_n(r)=\frac{Z}{r}是来源于原子核库伦引力产生的势能,VhV_h是哈特里势能,哈特里势能

2017-08-03 20:48:27 992

原创 一维声子晶体的谱有限元方法

一维声子晶体的谱有限元方法对于布洛赫周期边界条件的声子问题,通过物理建模,得到它的强形式是: 这里的杨氏模量和密度是分段函数: 对上述问题,做简单的代换,可以得到另一个周期化形式: 这里的E(x)E(x)和ρ(x)\rho(x)定义同上。对于这两个问题,做简单的变分,最后可以得到下面两个弱形式:布洛赫周期对应的广义特征值问题 周期对应的二次特征值问题 下面给出程序相应的程序和子函数

2017-08-03 19:36:36 1408

原创 单位分解有限元方法(PUFEM)

有限单元法是求解波传播问题的主要数值方法之一。虽然它有很多优点,并成功地模拟了很多波传播问题,但同样存在许多不足之处。 比如说,有些问题上,有限元方法的要达到一定的精度,计算量是很大的。另外,很多的问题用单纯的有限方法无法解决。举例来说,短波问题的数值模拟就是有限元法尚未解决的两个主要问题之一。例如,为使结果达到可接受的程度,一般说来低阶的有限单元每个波长需要至少布置10个节点。由此导致计算

2017-08-02 22:16:40 4378 2

原创 PLA(Padé近似)技术在二次特征值问题中的应用

PAL for QEP with low-rank damping一类特殊的二次特征值问题:二次特征值问题(QEP) Q(λ)x=(λ2M+λC+K)x=0)Q(\lambda)x = (\lambda^2M+\lambda C+K)x=0)低秩属性: rank(C)=r<<nrank(C)=r<<nPAL算法:给定一个位移σ\sigma和Padé近似的阶m。计算比例因子。通过某种线

2017-08-02 21:13:18 391

原创 matlab如何提取程序所依赖的所有自定义函数文件到一个文件夹中去

matlab如何提取程序所依赖的所有自定义函数文件到一个文件夹中去 如题,有时候拿到别人的matlab代码,有多个主文件和一堆子函数,这时候想将某个主文件和其用到的一些函数从一堆繁杂的程序中提取出来而不需要一个一个手动挑选,就可以利用matlab的matlab.codetools.requiredFilesAndProducts函数,这个函数在旧版的matlab中名称并不叫这个。下面给出一个简

2017-08-02 20:29:43 1776 1

原创 QEP(二次特征值问题)的SOAR和TOAR解法

QEP(二次特征值问题)的SOAR(二阶Arnoldi)和TOAR(二水平正交Arnoldi)解法关于二次特征值问题,在我的上一篇博文中,有过简单介绍。 参见二次特征值问题的一个解法。 今天要介绍的是,在求解二次特征值问题中SOAR和TOAR等几个方法的实现以及他们相互之间的比较。- 主程序如下:% demo routine for second-order Arnoldi method.%

2017-08-01 22:27:55 1404 13

原创 使用MATLAB自带函数求解二次特征值问题

使用MATLAB自带函数求解二次特征值问题二次特征值问题(QEP): 给定M,D,K∈CnxnM,D,K \in \mathbb{C^{nxn}},找到一个λ∈C\lambda \in \mathbb{C},以及一个非零向量x∈Cnx \in \mathbb{C^n},使得: Q(λ)x=(λ2M+λD+K)x=0Q(\lambda)x=(\lambda ^2M+\lambda D+K)

2017-08-01 22:08:11 1627

原创 牛顿迭代法在求解特征值问题中的应用

牛顿迭代法在求解特征值问题中的应用牛顿迭代法常被应用于求解非线性方程(组),通过几步简单的迭代,就能得到真解。 下面给出几个牛顿迭代法在求解矩阵特征值文中的几个算例:牛顿反迭代% inverse_iterationclcclearn = 2; F = @(z) [exp(1i*z.^2) 1; 1 1];Fp = @(z) [2i*z*exp(1i*z.^2) 0; 0 0];tol

2017-08-01 21:49:47 1689

原创 特征值问题的有限元MATLAB程序(一维)

薛定谔特征值问题和拉普拉斯特征值问题的有限元MATLAB程序(一维)−(p(x)u′(x))′+q(x)u(x)=λω(x)u(x)-(p(x)u'(x))'+q(x)u(x)=\lambda \omega(x)u(x) in Ω=(a,b),u(a)=0,u(b)=0 in \space \Omega=(a,b),u(a)=0,u(b)=0 这是对一般的Sturm-Liouville问题的一个表

2017-08-01 20:26:02 2758 4

傅里叶谱方法相关文献和代码汇总.rar

近20篇关于傅里叶谱方法的文章书籍(包括基于MATLAB工程实现),以及相关的代码。包括Nnmerical Recies(中文版和英文版)、Matlab微分方程高效解法:谱方法原理与实现、Spectral methods in Matlab等等。讲解参考我的博客:https://blog.csdn.net/lusongno1/article/details/84259533、https://blog.csdn.net/lusongno1/article/details/84849545等。

2019-12-14

马建国数学分析课堂讲义手稿

这是马建国老师数学分析上课用的讲义,红蓝两色,收集不易,请勿外传。

2021-01-09

博弈论与机制设计.rar

博弈论与机制设计相关资料,中文,内拉哈里,Narahari,囚徒困境,游戏理论,机制设计,VCG,二价密封拍卖,内容详实,清晰,我写得很清楚。

2020-06-08

beamer简明大方英文模板.rar

这是一份 beamer 的模板,简明大方,别具一格,看起来舒心,用起来顺手。欢迎下载。是在 metro 主题上的一种改进。

2020-05-20

song_data.csv

我的博客 https://blog.csdn.net/lusongno1/article/details/81039022 中使用的数据。Python 音乐推荐的用户数据。包含 [user_id,song_id,listen_count,title,artist,song] 等标签。

2020-03-13

单位分解有限元方法(PUFEM)_PPT

单位分解有限元方法(PUFEM)介绍的PPT,参考文档见:https://blog.csdn.net/lusongno1/article/details/76599281。

2019-12-14

Fast Marching.rar

快速行进算法(fast marching)完整的运行部分和函数输入变量的说明。采用了三种Fast Marching方法,包括传统的一阶fast marching方法、二阶msfm方法、和matlab工具箱方法。文档说明可参考博客:https://blog.csdn.net/lusongno1/article/details/88409735。

2019-12-14

双曲守恒律的ENO格式和WENO格式.zip

双曲守恒律的ENO格式和WENO格式的MATLAB代码,求解的是Burgers方程的初值问题,分别使用有限体积法4阶ENO格式,有限体积法3阶和5阶WENO方法求解。时间方向用三阶TVD Runge-Kutta方法。研究格式在解光滑情况下和有间断的情况下的数值精度并作图。说明文档参考博客https://blog.csdn.net/lusongno1/article/details/79070049。

2019-12-14

一维激波管问题的数值求解代码.rar

一维激波管问题的数值求解的MATLAB程序,包括Lax-Friedrichs, Lax-Wendroff,Roe和ENO等格式。说明文档可参考博客https://blog.csdn.net/lusongno1/article/details/90636623。

2019-12-14

g77安装包for tetview

g77的编译器的安装包。用以支持老版本的一些软件,比如tetview。

2018-11-22

OpenMP例子程序

这是OpenMP的例子程序,只要仔细研读,很快就能入门,不需要看网上所谓的教程。

2018-11-09

matlab入门PPT课件

MATLAB的入门PPT,看完之后,即使是没有任何matlab基础也基本入门了。

2018-11-09

书籍微分方程高效解法:谱方法原理与实现

微分方程高效解法:谱方法原理与实现的书籍,不是代码。

2018-11-09

lsec小陆的留言板

发表于 2020-01-02 最后回复 2020-08-02

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人 TA的粉丝

提示
确定要删除当前文章?
取消 删除