matlab误差分析,数值分析之MATLAB实验一误差分析

  • Post author:
  • Post category:其他


实验 第一次:

误差分析

对0,1,2,,20n =L ,按照下面两种算法计算定积分105n n

x y dx x =+?. 算法1: 利用递推公式115n n y y n

-=

-(1,2,,20)n =L ,取 1001ln 6ln 50.1823225

y dx x ==-≈+?

. 算法2: 利用递推公式11155

n n y y n -=-(20,19,,1)n =L 注意到2011120200001111126655105

x x dx dx x dx x =≤≤=+???,取 20111()0.00873020105126

y ≈+≈. 思考:从计算结果看,哪个算法是不稳定的,哪个算法是稳定的。

法一程序:

t=log(6.0)-log(5.0);

n=0;

y=zeros(1,21);

y(1)=t;

for k=2:21

y(k)=1/k-5*y(k-1);

n=n+1;

end

y(1:6)

y(7:21)

法二程序:

clear all

clc

y=zeros(21,1);

n=1;

y1=(1/105+1/126)/2;

y(21)=y1;

for k=21:-1:2

y(k-1)=1/(5*(k-1))-y(k)/5;

n=n+1;

end

y