- 博客(23)
- 资源 (17)
- 收藏
- 关注
RSS订阅
原创 复化公式收敛公式比较
数值积分中，有梯形公式、辛普森公式和高斯公式等等。为了使结果更加逼近真实结果，我们采用的方式一般是复化。所谓复化，就是等分成很多块，每一块再用这些公式计算面积。那么问题来了，既然数值积分方式，这么多，遇到实际问题时，该选哪一个呢？下面，就编写程序对几种复化公式的收敛速度进行比较一下。clear;clc;for m=1:14n=2*m;%设置划分个数na=0;%积分上限b=1;%积分下限
2016-11-06 20:29:50
3676
原创 matlab中的归一化方法
在数据处理的过程中，为了标准的统一，我们往往要对数据进行归一化。大多数人觉得这比较简单，不就是形如(x-min)/(max-min)这样么。事实上，这样想是没错的。事实上，我们使用matlab的内置函数，来进行数据的归一化，在处理大量的数据时，会节省不少的时间。有内置函数用的，尽量用内置函数，可以减少不必要的时间消耗。举些例子来说明问题吧。clc;clearP=magic(4)T=[1 2 3
2016-11-17 23:35:26
6905
原创 Latex巧画表格（从Excel复制）
Latex提供了不少绘制表格的宏包，但在latex里画表并不直观，特别是在表格比较大的时候，有时候也需要先用Excel等软件先对数据稍作处理，因而经常需要将Excel等表格转化成Latex代码。而Excel插件excel2latex就能实现这一功能。下面以Excel 2010为例说明下如何使用。自行搜索下载宏工具excel2latex，找个妥帖的地方安放。打开Excel，在文件-选项-信任中心-信
2016-11-16 22:29:40
3783
原创 如何用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
6040
原创 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
83308
11
原创 从一个椭圆积分看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
6494
原创 微信朋友圈如何发已存的小视频
我们知道，微信发朋友圈的时候，只能即录即发，不能够发已经存好的视频，或者转发别人的视频。下面以转发别人的小视频为例，说说如何发送已经存好的小视频。有时候，看到别人发的直升机，钞票啊等等小视频，感觉很拉风，恨不得是自己的。那么如何将别人的小视频据为己有，也风光一把呢？单刀直入主题。下载手机微信6.2.0版本并安装。高版本的微信对于发小视频加了验证，换言之，加密了，操作难度大，故换成低版本的。进朋友圈
2016-11-08 23:52:52
3238
原创 通过一个小函数看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
980
原创 地图上做标注
如何使用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
1046
原创 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
32823
2
原创 插值中产生的龙格现象
插值中产生的龙格现象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
5762
原创 赛马理论测试
赛马理论测试 在信息论里面，有一块是关于赛马理论的，大概意思就是要讨论什么情况下能赢钱，怎么做才能稳赢。那是关于最优双倍率的一系列讨论。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
2772
原创 几个画图小程序的分享（隐含波动率）
几个画图小程序的分享程序一%%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
1947
原创 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
2436
原创 求函数最小值
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
1802
原创 小观插值逼近的龙格现象
可以直观帮忙理解为什么会产生龙格现象这种怪诞。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
14770
原创 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
14395
原创 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
810
原创 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
36289
原创 zzu数学 实验十一最速降线
zzu数学 实验十一最速降线 意大利科学家伽利略在1630年提出一个分析学的基本问题──“一个质点在重力作用下，从一个给定点到不在它垂直下方的另一点，如果不计摩擦力，问沿着什么曲线滑下所需时间最短”。他说这曲线是圆，可是这是一个错误的答案。 瑞士数学家约翰．伯努利在1696年再提出这个最速降线的问题（problem of brachistochrone），征求解答。次年已有多位数学家得到正确答案
2016-11-06 10:55:20
4588
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
811
原创 zzu数学 实验八物理现象之模拟电场线
zzu数学 实验八物理现象之模拟电场线 大家都学过高中物理，看到两个点电荷之间的电场线，等势线等，有没想过这些图是怎么画出来的？如果是随手画出来的，为什么会这么精准？今天就带大家走进，电场线的绘制。提供两块代码作为对比参考。代码一 ** Physics ** ** wave **x0 = 0.0; v0 = 9
2016-11-06 10:22:41
2296
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
1162
各种课程大作业、项目等 latex 报告模板，有精美的 logo 封面、标题、摘要等。
2022-03-30
MATLAB 曲线提取分离数据生成傻瓜工具
2021-08-16
LaTeX模板.rar
2021-07-05
fast marching 方法程序.rar
2021-06-07
XMtool2.0.rar
2021-04-30
MATLAB微分方程高效解法.rar
2021-03-10
函数型数据主成分分析数据与代码全部.rar
2021-02-03
song_data.csv
2020-03-13
单位分解有限元方法(PUFEM)_PPT
2019-12-14
Fast Marching.rar
2019-12-14
双曲守恒律的ENO格式和WENO格式.zip
2019-12-14
一维激波管问题的数值求解代码.rar
2019-12-14
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人