![](https://img-blog.csdnimg.cn/20201014180756913.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
科学计算
文章平均质量分 83
yinlung
geek in sci and tech
展开
-
Matlab 中输入希腊字母
举例:fprintf('Kuiken reports %c''(0) = -0.769861.\n', char([952]))输出为:Kuiken reports θ'(0) = -0.769861.1.Matlab支持tex解析,所以如果是在figure里加希腊字母,或者是数学公式,非常简单。matlab中用转义符来输入希腊字母的方法上标用 ^(指数)转载 2012-11-05 10:16:19 · 18915 阅读 · 2 评论 -
使用Chebfun求解Blasius方程(二)
在上一篇文章使用Chebfun求解Blasius方程(一)里,我们使用chebfun求解了Blasius方程。 由于Blasius方程定义在半无界区间,因此我们将区间进行截断以求解,也就是说,在无穷远处的边界条件f’(+∞)=1被f’(infty)=1替代,此处infty是一个比较大的数字,如10,20,100等,但是并非是无穷大。 Chebfun能很好地表示具有可去奇点的函数,利用这一特原创 2012-12-06 18:32:00 · 2980 阅读 · 2 评论 -
阿达(Ada Lovelace)
转自百度百科:http://baike.baidu.com/view/1107373.htm简介: 阿达·奥古斯塔,19世纪诗人拜伦的女儿,数学家。穿孔机程序创始人,建立了循环和子程序概念。为计算程序拟定“算法”,写作的第一份“程序设计流程图”,被珍视为“第一个给计算机写程序的人”。编辑本段生平: 1815年生于伦敦,她是英国著名诗人拜伦转载 2012-12-10 20:25:27 · 7906 阅读 · 1 评论 -
Chebyshev 展开
Chebyshev展开是将有限区间上的光滑函数以Chebyshev多项式为基做展开。和Taylor展开不同的是,它对展开函数的光滑性要求较低,只需连续即可。著名的Chebfun系统基础之一就是Chebyshev展开。下面是Mathematica上的一个简单的Chebyshev展开,展开系数使用Gauss-Chebyshev积分计算,积分的代数精度是2*M+1,这里M是展开的阶数。(******原创 2013-01-13 15:31:19 · 5013 阅读 · 0 评论 -
Clenshaw–Curtis quadrature
参见 http://en.wikipedia.org/wiki/Clenshaw%E2%80%93Curtis_quadrature转载 2013-01-14 11:52:53 · 2053 阅读 · 0 评论 -
OpenCL开发案例学习
OpenCL开发案例学习本文主要为开发者提供一些OpenCL程序的实现方法,帮助开发者更直观地了解如何使用OpenCL编写程序。【案例】基于OpenCL的并行HDR算法实现【案例】基于OpenCL的K-means算法的并行实现【案例】基于OpenCL的Adaboost并行人脸检测算法【案例】基于OpenCL的三维模型凸包生成算转载 2013-01-18 15:00:31 · 1189 阅读 · 1 评论 -
Shanks transformation
转自:http://en.wikipedia.org/wiki/Shanks_transformationMotivationThe Shanks transformation is motivated by the observation that — for larger — the partial sum quite often behaves approximate转载 2013-02-24 16:01:27 · 2582 阅读 · 0 评论 -
BloomFilter(布隆过滤器)
布隆过滤器 (Bloom Filter)是由Burton Howard Bloom于1970年提出,它是一种space efficient的概率型数据结构,用于判断一个元素是否在集合中。在垃圾邮件过滤的黑白名单方法、爬虫(Crawler)的网址判重模块中等等经常被用到。哈希表也能用于判断元素是否在集合中,但是布隆过滤器只需要哈希表的1/8或1/4的空间复杂度就能完成同样的问题。布隆过滤器可以插入元转载 2013-03-09 11:43:41 · 1969 阅读 · 0 评论 -
初窥Applet
下午在查Fourier近似时,在网页上看到一个小程序,给定一些函数,可以将傅里叶近似的结果用图像显示,可以和用户交互,如下图,一目了然(如果大学教育能多多地使用这种方式,那么学生会更能接受)。一时好奇,就想看看源代码,同时想把这个小程序保存在本地。网页的地址是http://www.jhu.edu/signals/fourier2/index.html,通过查看网页的源代码,发现原创 2013-04-16 21:50:26 · 1664 阅读 · 0 评论 -
使用 Matlab 的 bvp4c 求解边值问题
Kuiken 利用相似变换,得到如下非线性微分方程满足如下边界条件其中, 表示对 求导, 为普朗特数. 此方程是耦合的非线性边值问题,在无穷远点具有奇性.当 时,使用Matlab的bvp4c求解如下:将原方程转化为一阶方程组% kuikenode.mfunction df=kuikenode(eta,f)sigma=1; df=[ f(2)原创 2012-11-06 15:18:15 · 13355 阅读 · 0 评论 -
使用Chebfun求解Blasius方程(一)
Chebfun的特点:1. 基于Chebyshev展开,展开项数由机器精度自适应控制;2. 将符号计算和数值计算结合,以处理数值的速度处理函数;3. 在Matlab上实现,将Matlab处理向量和矩阵的命令重载,以处理函数和算子;4. 基于Newton迭代法求解非线性微分方程;5. 使用自动微分技术计算Frechet导数;6. Chebop的实现利用了谱方法和惰性求值的思想原创 2012-12-06 10:20:11 · 4797 阅读 · 0 评论 -
MATLAB 函数句柄的用法
转自:http://blog.sina.com.cn/s/blog_6aaf46c401010fvd.html函数句柄(Function handle)是MATLAB的一种数据类型。引入函数句柄是为了使feval及借助于它的泛函指令工作更可靠;特别在反复调用情况下更显效率;使“函数调用”像“变量调用”一样方便灵活;提高函数调用速度,提高软件重用性,扩大子函数和私用函数的可调用范围;迅速获得转载 2012-11-23 10:47:58 · 3253 阅读 · 0 评论 -
MATLAB中的一些小技巧
1. Ctrl+C 中断正在执行的操作如果程序不小心进入死循环,或者计算时间太长,可以在命令窗口中使用Ctrl+c来中断。MATLAB这时可能正疲于应付,响应会有些滞后。2. figure命令新建一个绘图窗口figure 可以打开一个空的绘图窗口,接下的绘图命令可以将图画在它里面,而不会覆盖以前的绘图窗口。当有多个figure窗口时,在命令窗口中执行如Plot等命令将覆盖当前fig转载 2012-11-07 08:54:02 · 5334 阅读 · 0 评论 -
Lobatto quadrature
Lobatto积分类似于Guass积分,但有如下差异:1. 积分节点包括积分区间的端点.2. 积分具有2n–3次代数精度,其中n是节点数.公式如下:Lobatto quadrature of function f(x) on interval [–1, +1]:Abscissas: is the st zero of . Here Pn(x) ar翻译 2012-11-06 22:44:32 · 4367 阅读 · 0 评论 -
牛顿下山法
转自:http://blog.sina.com.cn/s/blog_49f7186101009e1b.html因牛顿迭代法受初值选取的限制,为防止迭代发散,对迭代过程再附加一项要求:|f(x(k+1))|% Newton.mfunction [x1,n]=Newton(f,x0,emg1,emg2)n=0;u=1;[f0,d0]=feval(f,x0);x转载 2012-11-27 10:22:39 · 9040 阅读 · 0 评论 -
关于Mathematica系统通讯机制MathLink的研究
计算机应用Computer Applications1999年 第19卷 第7期 Vol.19 No.7 1999关于Mathematica系统通讯机制MathLink的研究吕 民 蔡经球 摘 要 使用MathLink通讯机制可以扩展Mathematica系统的功能。本文介绍了MathLink通讯的基本方式和方法,设转载 2012-11-14 10:28:24 · 4039 阅读 · 0 评论 -
如何利用Mathematica调用C编写的函数
如何利用Mathematica调用C编写的函数[使用范围]Windows XP SP2; Mathematica 6.0; 传递变量类型简单的函数[方法]1.用C语言写好函数,如double f(double x,double y){return x*y;}2.找到路径C:\Program Files\Wolfram Research\Mathematica\6.0转载 2012-11-14 10:30:53 · 4587 阅读 · 0 评论 -
Mathematica 函数调用发生异常时停止计算
问题描述:用Mathematica写程序,如果调用Solve[ ]出现某种异常,则中断计算,返回错误信息。代码实现:sol = Check[Solve[EQ, vars], Print["The number of variables is greater than the number of equations."]; Abort[],So原创 2012-11-16 15:31:44 · 2946 阅读 · 0 评论 -
Mathematica中清除一系列符号定义的函数
Mathematica中使用Clear[symbol_1,symbol_2,... ] 清除符号symbol_i的值和定义。但是,有时需要清除一系列符号的值和定义。比如, 有一系列变量a, b, c构成一个集和,称作 Vars,Vars={a,b,c,...};中间计算得到,a=1; b=2; c=3;...现在想清除这一系列变量a, b, c,...的值。虽然可以将每一个原创 2012-11-20 22:46:02 · 8542 阅读 · 1 评论 -
Excel表格乘法函数公式
简介在Excel表格中,我们常常会利用Excel公式来统计一些报表或数据等,这时就少不了要用到加、减、乘、除法,在前面我们已经详细的讲解了Excel求和以及求差公式使用方法。那么我们又如何利用公式来对一些数据进行乘法计算呢?怎样快速而又方便的来算出结果呢?下面Word联盟就来教大家一步一步的使用Excel乘法公式!步骤/方法1、A1*B1=C1的转载 2012-11-21 16:08:06 · 6855 阅读 · 0 评论 -
Physicists Discover a Whopping 13 New Solutions to Three-Body Problem
转自:http://news.sciencemag.org/sciencenow/2013/03/physicists-discover-a-whopping.html物理家发现三体问题的13个新的周期解It's the sort of abstract puzzle that keeps a scientist awake at night: Can you predict ho转载 2013-04-19 09:32:30 · 1285 阅读 · 0 评论