自定义博客皮肤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

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

原创 matlab中的归一化方法

在数据处理的过程中,为了标准的统一,我们往往要对数据进行归一化。大多数人觉得这比较简单,不就是形如(x-min)/(max-min)这样么。事实上,这样想是没错的。事实上,我们使用matlab的内置函数,来进行数据的归一化,在处理大量的数据时,会节省不少的时间。有内置函数用的,尽量用内置函数,可以减少不必要的时间消耗。举些例子来说明问题吧。clc;clearP=magic(4)T=[1 2 3

2016-11-17 23:35:26 3459

原创 Latex巧画表格(从Excel复制)

Latex提供了不少绘制表格的宏包,但在latex里画表并不直观,特别是在表格比较大的时候,有时候也需要先用Excel等软件先对数据稍作处理,因而经常需要将Excel等表格转化成Latex代码。而Excel插件excel2latex就能实现这一功能。下面以Excel 2010为例说明下如何使用。自行搜索下载宏工具excel2latex,找个妥帖的地方安放。打开Excel,在文件-选项-信任中心-信

2016-11-16 22:29:40 2981

原创 如何用matlab和R语言画K线图

数据结构如图,分别表示开高低收,成交量和持仓量。 matlab绘制K线图%%filename = 'data.csv';delimiter = ',';formatSpec = '%f%f%f%f%f%f%f%[^\n\r]';fileID = fopen(filename,'r');dataArray = textscan(fileID, formatSpec, 'Delimite

2016-11-12 00:17:34 4865

原创 matlab中fmincon函数的使用

fmincon是用于求解非线性多元函数最小值的matlab函数。其语法格式比较多,参数也是各有各的意思,我就举几个例子,大家可从代码中去体会其使用方法。示例一%% 只是大概写个板式,这里的fun2和fun3要自己去定义,这里就不写出来了,前面%有介绍过clcclearA=[3 2 1;1 2 3; 2 0 0;0 3 0;0 0 2]b=[120;80;96;102;40]lb=zero

2016-11-09 11:32:20 59310 6

原创 从一个椭圆积分看matlab求定积分

一个积分,能求出解析解固然好,但是求不出解析解,求一个数值解基本能够满足实际的需求了。现实生活中,我们遇到的积分大部分是解不出或者很难解出解析解,这时候,就需要我们求其数值解。matlab提供了一个求积分的强大算法。下面以椭圆积分为例,我来说一说,matlab如何求定积分。clcclearsyms theta;a=5;b=3;c=sqrt(a^2-b^2);e=c/a;hp=abs(a

2016-11-09 11:01:11 4258

原创 微信朋友圈如何发已存的小视频

我们知道,微信发朋友圈的时候,只能即录即发,不能够发已经存好的视频,或者转发别人的视频。下面以转发别人的小视频为例,说说如何发送已经存好的小视频。有时候,看到别人发的直升机,钞票啊等等小视频,感觉很拉风,恨不得是自己的。那么如何将别人的小视频据为己有,也风光一把呢?单刀直入主题。下载手机微信6.2.0版本并安装。高版本的微信对于发小视频加了验证,换言之,加密了,操作难度大,故换成低版本的。进朋友圈

2016-11-08 23:52:52 2427

原创 通过一个小函数看fprintf和dlmwrite的使用

fprintf和dlmwrite的使用示例function k=test(x,y)if x+y>0 fprintf('Hello'); dlmwrite('E:\win7桌面\data1.txt',x,'newline','pc')endk=x+yend

2016-11-06 23:13:42 636

原创 地图上做标注

如何使用matlab做标注?这里有一些小技巧。 编写mercatorProjection.m子函数如下:function [x,y] = mercatorProjection(lon, lat, width, height) x = mod((lon+180)*width/360, width) ; y = height/2 - log(tan((lat+90)*pi/360))*

2016-11-06 20:52:10 649

原创 matlab利用fzero求零点

matlab利用fzero求零点f = @(x)x.^5-2*x-5; %注意函数的书写格式z = fzero(f,[-1000 1000]) 主要学习一下书写格式,函数也可以另建文件这样写function y=f(x)y=1/(1+x^2);end

2016-11-06 20:43:41 22353

原创 插值中产生的龙格现象

插值中产生的龙格现象g=f(x);%%这只是一个抽象化的说明,函数还是要自己去定义的。for k=2:2:12a=-8;b=8;xi=linspace(a,b,k)fi=subs(g,xi)A=[xi(:) double(fi(:))]p=Lagrange(A);Xi=a:0.01:bsubplot(2,3,k/2)plot(Xi,double(subs(g,Xi)),'r')

2016-11-06 20:36:02 3651

原创 复化公式收敛公式比较

数值积分中,有梯形公式、辛普森公式和高斯公式等等。为了使结果更加逼近真实结果,我们采用的方式一般是复化。所谓复化,就是等分成很多块,每一块再用这些公式计算面积。那么问题来了,既然数值积分方式,这么多,遇到实际问题时,该选哪一个呢?下面,就编写程序对几种复化公式的收敛速度进行比较一下。clear;clc;for m=1:14n=2*m;%设置划分个数na=0;%积分上限b=1;%积分下限

2016-11-06 20:29:50 2415

原创 赛马理论测试

赛马理论测试 在信息论里面,有一块是关于赛马理论的,大概意思就是要讨论什么情况下能赢钱,怎么做才能稳赢。那是关于最优双倍率的一系列讨论。clcclearn=292;p=[1/4 1/4 1/2];m=length(p);o=[4 5 3];bw=[0.1 0 0.9];Wbp=sum(p.*log2(bw.*o));Sw=2^(Wbp);for i=1:n Swi(i

2016-11-06 19:47:14 1519

原创 几个画图小程序的分享(隐含波动率)

几个画图小程序的分享程序一%%for huangclc clearclose allx=1:0.01:5;y=0.3*(x-3).^2+1;plot(x,y);xlim([0 6]);ylim([0 5]);title('隐含波动率随执行价格变化图');xlabel('执行价格');ylabel('隐含波动率');%% 下面是第二个图x=2:0.01:12;y=log

2016-11-06 19:34:54 1113

原创 matlab模拟飞机延时gif动图程序

clearclc%飞机追踪潜艇航线图t=linspace(0,1,2000); flag=0;p=zeros(2000,2);q=ones(2000,2);a=zeros(8,2);b=zeros(8,2);a(:,1)=0;a(1,2)=6;b(1,2)=0;b(1,1)=pi/2;for i=1:2000 theta=sqrt(3)*log(10*t(i)+1); p

2016-11-06 18:32:42 1925

原创 求函数最小值

f[x_] := (x^4 + 1)/(x^3 - x);Plot[f[x], {x, -9, 9}]Show[%4, PlotLabel -> HoldForm[老师让我们找的极小值在这], LabelStyle -> {GrayLevel[0]}]Show[%3, AxesStyle -> Black]Show[%4, AxesLabel -> {HoldForm[x], HoldForm[

2016-11-06 11:28:59 819

原创 小观插值逼近的龙格现象

可以直观帮忙理解为什么会产生龙格现象这种怪诞。Animate[f = 1/(1 + x^2); r = 5; g1 = Plot[f, {x, -r, r}, PlotStyle -> RGBColor[1, 0, 0]]; p0 = Table[{x0, f /. x -> x0}, {x0, -r, r, 2 r/n}]; Interf = Interpolation[p0, Inte

2016-11-06 11:17:43 981

原创 Mathematica做拟合

打开mathematica键入如下内容:advol={{0.49,100},{0.69,150},{0.88,200},{1.06,250}}; line = Fit[advol, line = Fit[advol, line = Fit[advol, line = Fit[advol, line = Fit[advol, line = Fit[advol, {1,x,x^2},x] 按 Shi

2016-11-06 11:10:07 8030

原创 Mathematica精华

Clear[d, a, n, t];n = 50;a = Graphics[Line[{{n - 1, 0}, {0, 0}, {0, n - 1}}]];d = Array[f, n - 1];t = Table[{{0, i}, {n - 1 - i, 0}}, {i, 1, n - 1}];For[i = 1, i <= n - 1, i++, d[[i]] = Graphics[L

2016-11-06 11:07:41 522

原创 Mathematica画图函数命令大全

二维作图 Plot[f,{x,xmin,xmax}] 一维函数f[x]在区间[xmin,xmax]上的函数曲线 Plot[,f2.{f1.},{x,xmin,xmax}] 在一张图上画几条曲线 ListPlot[{y1,y2,..}] 绘出由离散点对(n,yn)组成的图 ListPlot[{{x1,y1},{x2,y2},..}] 绘出由离散点对(xn,yn)组成的图 Plarametri

2016-11-06 11:05:31 17230

原创 zzu数学 实验十一最速降线

zzu数学 实验十一最速降线 意大利科学家伽利略在1630年提出一个分析学的基本问题──“一个质点在重力作用下,从一个给定点到不在它垂直下方的另一点,如果不计摩擦力,问沿着什么曲线滑下所需时间最短”。他说这曲线是圆,可是这是一个错误的答案。 瑞士数学家约翰.伯努利在1696年再提出这个最速降线的问题(problem of brachistochrone),征求解答。次年已有多位数学家得到正确答案

2016-11-06 10:55:20 2990 2

原创 zzu数学 实验九迭代一方程求解

方程求解的方法有很多啦,最常用,使用最广泛的就是迭代法,原理好理解,效果也不差。一般方程求解g[x_] := x^3 - 12 x + 2;f = x - g[x]/D[g[x], x];x0 = -3;series = {x0};For[i = 1, i <= 5, i++, x0 = f /. x -> x0; AppendTo[series, x0]]series = N[ser

2016-11-06 10:33:11 518

原创 zzu数学 实验八物理现象之模拟电场线

zzu数学 实验八物理现象之模拟电场线 大家都学过高中物理,看到两个点电荷之间的电场线,等势线等,有没想过这些图是怎么画出来的?如果是随手画出来的,为什么会这么精准?今天就带大家走进,电场线的绘制。提供两块代码作为对比参考。代码一 ** Physics ** ** wave **x0 = 0.0; v0 = 9

2016-11-06 10:22:41 1162 1

原创 zzu数学 实验七几何变换

zzu数学 实验七几何变换在如photoshop等的各大图形处理软件中,我们能看到各种图形变换,斜切,对称,翻转等,其妙无穷,我们却对其原理一无所知。其实,并不难,今天我以mathematica为例。做两个简单的示例。示例1 ** Linear Transformation **Clear[f]; t = Pi/6;a1 = Cos[t]; b1 = -Sin[t]; a2 = Si

2016-11-06 10:07:11 846

傅里叶谱方法相关文献和代码汇总.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的粉丝

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