数值分析实验报告1.docx

文档介绍：
数值分析实验报告1

数值分析上机试验报告
〔注：本试验报告中全部程序均为MATLAB语言程序〕
数 原理：秦久韶算法 程序：
n=input('n='); x=input('x='); a(1)=5; for k=1:n;
a(k+1)=2.*a(k)+3; end
s(n+1)=a(n+1); for i=n:-1:1
s(i)=x.*s(i+1)+a(i); end
Pnx=s(1)
运行结果：
n=101 x=
Pnx =

n=150 x=13
Pnx =
+213











3、设Y0=28，按递推公式Yn= Yn?1-1/101783计算到Y101，Y500。假设取783?(五位有效数字)，试问计算Y101、Y500将有多大的误差。 目的：运用递推公式计算求数列第n项值 原理：用递推公式迭代 程序：
clear
k=input('k=') Y0=28;
Y(1)=Y0-sqrt(783)/101; for n=2:k
Y(n)=Y(n-1)-sqrt(783)/101; end
y(1)=Y0-; for n=2:k
y(n)=y(n-1)-; end
y101=y(101) y500=y(500)
e1=y(101)-Y(101) e2=y(500)-Y(500)
运行结果：
k=500 k =
500
y101 =

y500 =
-+002 e1 =
-004 e2 =
-004











二章
目的：驾驭插值原理和不同插值方法
3、给出f(x)=lnx的数值表，。 x Lnx - - 原理：matlab线性插值函数 线性插值程序：
- - - x1=input('x1='); x=[::];
y=[- -