% matlab script to test efficiency of % Euler's method, classical Runge-Kutta, and ode45 % on Arenstorf orbit problem close all clear all % these are variables we would like the right-hand function to "see" % without actually passing them as arg…

MATLAB语言 function y=f(x) y=f(x); %函数f(t)的表达式 i=0; %二分次数记数 a=a; %求根区间左端 b=b; %求根区间右端 fa=f(a); %计算f(a)的值 fb=f(b); %计算f(b)的值 c=(a+b)/2; %计算区间中点 fc=f(c); %计算区间中点f(c) while abs(fc)>=ε; %推断f(c)是否为零点 if fa*fc>=0; %推断左側区间是否有根 fa=fc; a=c; else fb=fc; b=c; en…

定义函数 function y=f(x) y=f(x).%函数f(x)的表达式 end function z=h(x) z=h(x).%函数h(x)的表达式 end 主程序 x=X;%迭代初值 i=0;%迭代次数计算 while i<= 100%迭代次数 x0=X-f(X)/h(X);%牛顿迭代格式 if abs(x0-X)>0.01:%收敛推断 X=x0; else break end i=i+1; end fprintf('\n%s%.4f\t%s%d','X='.X.'i='.i) %产…

在离散数据的基础上补插连续函数,使得这条连续曲线通过全部给定的离散数据点.插值是离散函数逼近的重要方法,利用它可通过函数在有限个点处的取值状况,估算出函数在其他点处的近似值.曲面插值是对三维数据进行离散逼近的方法,MATLAB中的曲面插值函数有Triscatteredinterp,interp2,griddata等.我们以griddata为例讲解曲面插值及其交叉验证的过程. 一.  gridata曲面插值 gridata不仅可以对三维曲面进行插值,还能对四维的超平面进行插值.griddata的调…

(转于它处,仅供参考) 1.. Ctrl+C 中断正在执行的操作 如果程序不小心进入死循环,或者计算时间太长,可以在命令窗口中使用Ctrl+c来中断.MATLAB这时可能正疲于应付,响应会有些滞后. 2. figure命令新建一个绘图窗口 figure 可以打开一个空的绘图窗口,接下的绘图命令可以将图画在它里面,而不会覆盖以前的绘图窗口.当有多个figure窗口时,在命令窗口中执行如Plot等命令将覆盖当前figure窗口中的对象.所谓的当前figure窗口,也就是最后一次查看的窗口(可以用命令…

Euler's method…

http://blog.sina.com.cn/s/blog_67f37e760101bu4e.html 实例结果http://wenku.baidu.com/link?url=SiGsFZIxuS1E1VZWtixqXdjG5Y9SY4tu1W8TXgk147HDOLLCgpffjX8ywDMIH1PYkQSi5rp8gkmnMPsLhH-IUaGjMl8hsRhZQssTPmnBjdy 关于MATLAB 插值(Interpolation)http://blog.sina.com.cn/s/b…

今天用Matlab跑程序,由于数据量太大,又出现 Out of memory. Type HELP MEMORY for your options.的问题.看到这篇文章非常实用,转过来方便查阅~ 用 Matlab 进行大规模科学计算或仿真时,内存是一个需要时常注意的问题.在matlab里运行 >>system_dependent memstats 就可以看到内存的使用状况.当你写的 Matlab 程序跳出“Out of Memory” 时,以下几点措施是需要优先考虑的解决方法: 1. 升级内存…

http://blog.sina.com.cn/s/blog_67f37e760101bu4e.html 实例结果http://wenku.baidu.com/link?url=SiGsFZIxuS1E1VZWtixqXdjG5Y9SY4tu1W8TXgk147HDOLLCgpffjX8ywDMIH1PYkQSi5rp8gkmnMPsLhH-IUaGjMl8hsRhZQssTPmnBjdy 关于MATLAB 插值(Interpolation)http://blog.sina.com.cn/s/b…

使用matlab过程中经常会出现内存不足的问题,这里转载一篇来自http://blog.csdn.net/xiaojidan2011/article/details/8089532 的博文,解决这一问题. 做图像处理的,对matlab这工具软件一定不会陌生,他的高集成度为图像处理提供了很大的方便,但我们在使用时难免会出现内存不足的问题,2.3G的内存还不够处理一次简单图像,这确实很恐怖.但这里的警告并不是说真的物理内存用完了,一般out of memenry存在以下几种情况:1. 变量需要的存储…