计算机数学软件Matlab在解决方程求根与概率问题的简单应用.docxVIP

计算机数学软件Matlab在解决方程求根与概率问题的简单应用1、迭代法与开方求值迭代法是一种按同一公式重复计算逐次逼近真值的算法，是数值计算普遍使用的重要方法。以开方运算为例，它不是四则运算，因此在计算机上求开方值就要转化为四则运算，使用的就是迭代法。假定，求等价于解方程：.这是方程求根问题，可用迭代法求解。现在用简单方法构造迭代法，先给一个初始值近似，令，是一个校正量，称为增量，于是方程式化为：，即由于是小量，若省略高阶项，则得：，即.于是：这里不是的真值，但它是真值的进一步近似，重复以上过程可得到迭代公式：它可逐次求得随着迭代次数的增大，的值就越来越趋近于的真值，可以作为的近似值。有了以上迭代法的理论基础，对于任意的，可通过软件包编程求解其近似值。首先，使用软件包找出距离最近的整数作为初值，再进行迭代运算。先建一个名为sqrt1.m的function文件，并保存：function r=sqrt1(a,n)%让计算机找出距离根号a最近的整数b，n表示n次迭代for k=1:aif a-k^20a-(k+1)^20 b=k;break; endend%迭代法求出根号a的值x=zeros(n,1);x(1)=b;for j=2:n x(j)=(x(j-1)+a/x(j-1))/2;endr=x(j);若要求经过5次迭代的值，在指令框内输入：sqrt(2,5)按回车键执行，得到结果：ans = 1.414213562374690与真值近似值相对比，5次迭代的误差小于。对于任意的非完全平方数a，的值均可用此种方法求解。2、随机数与概率提供了产生随机数的函数，产生的随机数是由计算机内部随机而给定的。而在概率论中，对某一随机事件的概率估计，可以通过重复大量次数的同一实验，用频率的值来估计概率值。然而人工的独立重复随机试验有时候会有人为误差，而且随着试验次数的增大，满足条件的频数有时候不方便统计。可以通过产生随机数的方法，模拟随机试验，并统计符合条件的频数，来估计概率。当随机试验次数足够大时，可以利用频率的值近似代替概率。3、直线上两点所截线段构成三角形的概率例：在一线段AB上随机地取两个点把线段截为三段，求这三段可以构成三角形的概率。处理时，将AB的长度看为单位1，而AC与AD的长度则是在0~1之间产生随机数，注意，有可能是产生的随机数表示的长度中，D点在C点之前。那么三角形的三边可表示为：边1：，边2：边3:让经的实验判断，求频率，输入以下程序：ticNA=rand(N,1);B=rand(N,1);for t=1:1:N C(t,1)=1-(A(t,1)+B(t,1)+abs(A(t,1)-B(t,1)))/2; D(t,1)=abs(A(t,1)-B(t,1)); E(t,1)=(A(t,1)+B(t,1)-abs(A(t,1)-B(t,1)))/2; s=0;endfor j=1:Nif E(j,1)+ D(j,1)C(j,1)E(j,1)+ C(j,1)D(j,1)C(j,1)+ D(j,1)E(j,1) s=s+1;endendp=s/Ntoc执行程序，可得到结果：p = 0.250033587500000Elapsed time is 241.661775 seconds.与几何概率的值p=1/4基本相符合。4、两人相约见面的概率例：两人相约7点到8点在某地会面，先到者等候另一人20分钟，这时就可以离开，试求这两人可能会面的概率。分析：以分别表示两人到达时刻，则会面的充要条件为：模拟此过程时，假设进行试验，产生0~60之间的随机数，并统计其满足条件的频数即可。输入程序:ticN=800000000;A=60*rand(N,1);B=60*rand(N,1);s=0;for k=1:1:Nif abs(A(k,1)-B(k,1))=20 s=s+1;endendp=s/Ntoc执行结果为：p = 0.555571187500000Elapsed time is 8.641059 seconds.求得p =0.555571187500000，计算机执行此过程需要8.641059秒，与准确概率基本相符。