科学计算课程设1.docVIP

科学计算课程设计 Riemann积分 对于一在区间[a,b]上之给定非负函数f(x)，我们想要确定f(x)所代表的曲线与X坐标轴所夹图形的面积，我们可以将此记为积分，黎曼积分的核心思想就是试图通过无限逼近来确定这个积分值。 function y=f(x) y=exp(-x^2) m=0; n=100000 //………………………. //(n用100、1000、10000、100000分别代入) for s=1:n x=s/n; m=m+f(x); end t=m/n n 积分结果 1000 0.74650801121977 10000 0.74679252617135 100000 0.74682097220349 1000000 0.74682381675210 复合梯形积分 将 [a, b] 分成 n 等分 [xi , xi+1] ，其中 (i = 0, 1, …, n) 复合梯形公式 余项 function y=f(x) y=exp(-x^2) n=10000; h=1/n; m=0; for k=1:n-1 x=k*h;f=exp(-x^2); m=m+f; end m=m*2; t=1;s=exp(-1); T=h/2*(t+s+m) n 积分结果 1000 0.74682407149918 10000 0.74682413219929 100000 0.74682413280629 1000000 0.74682413281238 复合simpson积分 复合 Simpson 公式 余项 function y=f(x) y=exp(-x^2) n=10000; m=0;t=0; h=1/n; for k=0:n-1 x=k*h+h/2; f=exp(-x^2); m=m+f; end m=m*4; for k=1:n-1 x=k*h; f=exp(-x^2); t=t+f; end t=t*2; S=h/6*(1+m+t+exp(-1)) n 积分结果 1000 0.74682413281243 10000 0.74682413281243 100000 0.74682413281242 1000000 0.74682413281245 Romberg积分 龙贝格求积公式也称为逐次分半加速法。它是在梯形公式、辛卜生公式和柯特斯公式之间的关系的基础上，构造出一种加速计算积分的方法。 作为一种外推算法, 它在不增加计算量的前提下提高了误差的精度. 在等距基点的情况下，用计算机计算积分值通常都采用把区间逐次分半的方法进行。这样，前一次分割得到的函数值在分半以后仍可被利用,且易于编程 。function y=f(x) y=exp(-x^2) s=4; q=0; T(1,1)=1/2*(1+exp(-1)); for k=1:s for i=1:2^(k-1) y=f((2*i-1)/(2^k)); q=q+y; end T(k+1,1)=1/2^k*T(k,1)+q/2^k; for m=1:k T(k+1,m+1)=4^m/(4^m-1)*T(k+1,m)-1/(4^m-1)*T(k,m); end end p=T(s+1,s+1) d=T(s+1,s+1)-T(s,s) 投点式Monte Calor法 蒙特卡罗分析，是一种使用随机抽样统计来估算数学函数的计算方法。它需要一个良好的随机数源。这种方法往往包含一些误差，但是随着随机抽取样本数量的增加，结果也会越来越精确。蒙特卡罗方法在纯数学方面一般用来求解一个函数的定积分。它的计算过程如下：先在一个区间或区域内随机抽取一定数量的独立变量样本，然后求相应的独立因变量的平均值，最后用随机样本所在区间（或区域）的长度（或大小）除以所求出的平均值。它与传统的估算定积分的方法有很大差别，传统方法在区间或区域内抽取样本点时是间隔相等、均匀抽取的。蒙特卡罗方法以其在第二次世界大战时被用于原子弹的设计而闻名于世。现在它也已经被应用于多种领域，如超高速公路的运输流量分析、行星演变模型的建立以及股票市场波动的预测。这种方法同样也可应用于集成电路设计、量子力学和通信工程。function y=f(x) y=exp(-x^2) n=1000; X=rand(n);Y=rand(n); s=0; for i=1:n if Y(i)f(X(i)) s=s+1; end end t=s/n n 积分结果 1000