3.计算定积分Q(k)
4.然后朝着inf/-inf或者歧义点移动一个步长,然后计算定积分
Q(k+1)
5.判断(abs(Q(k+1)-Q(k))
QuasiMonteCarlo或者Oscillatory]
39)为什么我的Matlab程序这么慢?
:#FangQ(Qianqian.Fang@http://www.doczj.com/doc/3601188984868762caaed57c.html), 2002/6/22. BigGreen/en_Matlab#
我们工学院的收发室的门上贴着一张小纸条,写的是
"Our policy is always blaming the computer"
大多数的人在遇到问题的时候,总是责备计算机如何如何,别人如何
如何,其实,最最主要的因素是在于自己。
一个程序运行快慢,有很多因素决定,最主要的是算法,简炼而优美
的的数学公式胜过100遍的优化。能从算法上改进,才能比别人有根
本的优势。计算机也很重要,以前我总把自己用的PC看成万能的加
以崇拜,对UNIX嗤之以鼻,结果当自己真正开始算起来,才知道差
别有多大。搞大型数值计算的,没有好的工作站或者并行系统,就输
在了起跑线上了。然后是程序的优化,看看变量是否占用太多内存,
看看是否有功能重复的模块或者计算,经常的是用牺牲内存来换取速
度,具体取舍,具体需要来决定。用profile看看哪些语句占用时间最
多,然后把核心部分进行优化。
如果是使用Matlab,使用vectorization和矩阵整体操作的代码要比大量
的for循环快很多,eval/inline函数如果出现在核心循环,也会让速度下
降几时倍的。
40).Matlab中如何作非线性回归?
:#FangQ(Qianqian.Fang@http://www.doczj.com/doc/3601188984868762caaed57c.html), 2002/6/22. BigGreen/en_Matlab#
请参考
http://www.doczj.com/doc/3601188984868762caaed57c.html/support/solutions/data/10652.shtml
matlab默认只提供了多项式拟合的函数polyfit,对于其他稍微简单
一点的拟合,如标准的指数、对数、高阶多项式拟合,都有解析公式,参见:http://www.doczj.com/doc/3601188984868762caaed57c.html/LeastSquaresFitting.html
对于更加复杂的非线性函数,建议使用Mathematica或者DataFit
Mathematica中提供了Fit[],以及
<< Statistics`NonlinearFit`