数值分析上机实践题11-2013
研究生数值分析上机试题及解答
东华大学研究生数值分析试题(上机部分)A 卷2008年12月 时间:60分钟班级 学号 机号 姓名 得分 注意:要求写出M 函数(如果需要)、MATLAB 命令和计算结果。
1. 求下列方程组在0<α, β<1中的解⎩⎨⎧-=+=βαββααsin 2.0cos 7.0cos 2.0sin 7.0 命令fun=inline('[x(1)-0.7*sin(x(1))-0.2*cos(x(2)),x(2)-0.7*cos(x(1))+0.2*sin(x(2))]','x'); [x,f,h]=fsolve(fun,[0.5 0.5]) 结果α=0.5265,β=0.50792命令>> fun=inline('c(1)+c(2)*x.^2','c','x'); >> x=[1.1 1.3 1.4 1.6 1.8]; >> y=[26 22 23 24 25];>> c=lsqcurvefit(fun,[0 0],x,y) 结果 c =23.7256 0.12873.求解下列微分方程组2(0)2013(0)1x x yx t y x yy '=-=⎧<<⎨'=+=⎩(结果只要求写出t =1时的解) 命令>> fun=inline('[y(1)-2*y(2);3*y(1)+y(2)]','t','y'); >> [t,y]=ode45(fun,[0 1], [2 1]) 结果x(1)=-5.6020, y(1)=2.15634.用定步长Gauss 积分法(课本123页)计算积分31e ln(1)x x dx -+⎰的近似值(等分数取4,每段取2个Gauss 点)。
命令fun=inline('exp(-x).*log(1+x)','x'); nagsint(fun,1,3,4,2) 结果 0.30865.矩阵改进平方根分解(课本25页)的计算公式为: d 1=a 11, 对i =2, 3, ⋯, n ,iki k ik ii i j ij ij j k jk ik ij ij l s a d i j d s l l s a s ∑∑-=-=-=-==-=1111,1,,2,1 ,/ ,试编写矩阵改进平方根分解的程序,并求矩阵1111551514A -⎛⎫ ⎪=-- ⎪ ⎪-⎝⎭的改进平方根分解。
数值分析第七章上机题
在命令窗口输入: >>f = inline('[(x1+3)*(x2^3-7)+18;sin(x2*exp(x1)-1)]','x1','x2'); >>g = inline ('[x2^3-7,3*x2^2*(x1+3);x2*exp(x1)*cos(x2*exp(x1)-1),exp(x1)*cos(x2*exp(x1)-1)]','x1','x2');
hist = 0 -0.428571428571429 -0.000000000000101 0 1.557407724654902 1.000000000000127 -0.141348392468100 1.087738055836075 -0.002875590925150 1.001269946612821
数理强化班
数值分析第七章计算机实习题
写一程序实现下面问题的牛顿算法——求解方程组:
3 ( x1 3)( x2 7) 18, x1 sin( x e 1) 0. 2
源程序如下:
function [x,it,hist] = newton2(x0,f,g,maxit,tol) % Newton method for eqation systerm % INPUTS: % x0 % f % g % maxit % tol % OUTPUTS: % x % it % hist format long; if nargin<5, if nargin<4, if nargin<3, end end end flag = 1; x0 = [0;0]; x = x0; hist = x; it = 0; x = x0 - feval(g,x0(1),x0(2))\feval(f,x0(1),x0(2)); if norm(x0-x)>=tol, x0 = x; else fprintf('\nNewton Iteration successes!!\n') return end it = it + 1; for k = 1:maxit, tol = 1e-7; maxit = 100; error('too few input!!'); solution iteration history of iteration initial point function gradient maximum iteration tolerance for convergence
数值分析练习题附答案
1
2-3 对矩阵 A 进行 LDLT 分解和 GGT 分解,求解方程组 Ax=b,其中
16 4 8
1
A=( 4 5 −4) , b=(2)
8 −4 22
3
解:(注:课本 P26 P27 根平方法)
设 L=(l i j ),D=diag(di),对 k=1,2,…,n,
其中������������=������������������-∑������������=−11 ���������2��������� ������������
������31=(������31 − ∑0������=1 ������3������������1������ ������������)/ ������1=186=12
������32=(������32
−
∑1������=1
������3������������2������
������������ )/
6.6667
,得 ������3 = 1.78570
−1 209
������4
0
������4
0.47847
(
56
−1
780 (������5) 209)
(200)
(������5) ( 53.718 )
1 −1
4
1 −4
15
������1
25
������2
6.6667再由1源自− 15561
− 56
209
x (k1) 1
1 5
(12
2 x2( k )
x (k) 3
)
2 5
x (k) 2
数值分析练习题附答案
目录一、绪论------------------------------------------------------------------------------------- 2-2二、线性方程组直接解法列主元高斯LU LDL T GG T-------------------- 3-6二、线性方程组迭代法----------------------------------------------------------------- 7-10 三、四、非线性方程组数值解法二分法不动点迭代---------------------- 11-13五、非线性方程组数值解法牛顿迭代下山弦截法----------------- 14-15六、插值线性插值抛物线插值------------------------------------------------ 16-18七、插值Hermite插值分段线性插值-----------------------------------------19-22八、拟合------------------------------------------------------------------------------------ 23-24九、数值积分----------------------------------------------------------------------------- 25-29十、常微分方程数值解法梯形欧拉改进----------------------------------- 30-32 十一、常微分方程数值解法龙格库塔------------------------------------------ 33-35绪论1-1 下列各数都是经过四舍五入得到的近似值 ,试分别指出它们的绝对误差限,相对误差限和有效数字的位数.X 1 =5.420, X 2 =0.5420, X 3 =0.00542, X 4 =6000, X 5 =0.6×105注:将近似值改写为标准形式X 1 =(5*10-1+4*10-2+2*10-3+0*10-4)*101 即n=4,m=1 绝对误差限|△X 1|=|X *1-X 1|≤ 12×10m-n =12×10-3 相对误差限|△r X 1|= |X∗1−X1||X∗1|≤|X∗1−X1||X1|= 12×10-3/5.4201-2 为了使101/2 的相对误差小于0.01%, 试问应取几位有效数字?1-3 求方程x 2 -56x+1=0的两个根, 使它们至少具有4位有效数字( √783≈27.982)注:原方程可改写为(x-28)2=783线性方程组解法(直接法)2-1用列主元Gauss消元法解方程组解:回代得解:X1=0 X2=-1 X3=12-2对矩阵A进行LU分解,并求解方程组Ax=b,其中解:(注:详细分解请看课本P25)A=(211132122)→(211(1/2)5/23/2(1/2)3/23/2)→(2111/25/23/21/2(3/5)3/5)即A=L×U=(11/211/23/51)×(2115/23/23/5)先用前代法解L y=P b 其中P为单位阵(原因是A矩阵未进行行变换)即L y=P b 等价为(11/211/23/51)(y1y2y3)=(111)(465)解得 y 1=4 y 2=4 y 3=35再用回代解Ux =y ,得到结果x即Ux =y 等价为(2115/23/23/5)(x 1x 2x 3)=(y 1y 2y 3)=(443/5) 解得 x 1=1 x 2=1 x 3=1即方程组Ax=b 的解为x =(111)2-3 对矩阵A 进行LDL T 分解和GG T 分解,求解方程组Ax=b,其中A=(164845−48−422) , b =(123)解:(注:课本 P 26 P 27 根平方法)设L=(l i j ),D=diag(d i ),对k=1,2,…,n,其中d k =a kk -∑l kj 2k−1j=1d jl ik =(a ik −∑l ij l kj k−1j=1d j )/ d k 即d 1=a 11-∑l 1j 20j=1d j =16-0=16因为 l 21=(a 21−∑l 2j l 1j 0j=1d j )/ d 1=a 21/ d 1=416=14 所以d 2=a 22-∑l 2j 21j=1d j =5-(14)2d 1=4同理可得d 3=9 即得 D=(1649)同理l 11=(a 11−∑l ij l 1j 0j=1d j )/ d 1=1616=1=l 22=l 33 l 21=(a 21−∑l 2j l 1j 0j=1d j )/ d 1=416=14 l 31=(a 31−∑l 3j l 1j 0j=1d j )/ d 1=816=12 l 32=(a 32−∑l 3j l 2j 1j=1d j )/ d 2=−4−12×14×164=−64=-32即L=(114112−321) L T=(114121−321) 即LDL T分解为A=(114112−321)(1649)(114121−321)解解:A=(164845−48−422)→(41212−32−33)故得GG T分解:A=(4122−33)(4122−33) LDL T分解为A=(114112−321)(1649)(114121−321) 由(114112−321)(y 1y 2y 3)=(123) ,得(y 1y 2y 3)=(0.250.8751.7083)再由(4122−33)(x 1x 2x 3)=(0.250.8751.7083) ,得(x 1x 2x 3)=(−0.54511.29160.5694)2-4 用追赶法求解方程组:解:(4−1−14−1−14−1−14−1−14)→(4−14−1154−415−15615−1556−120956−56209−1780209)由(4−1154−15615−120956−1780209)(y1y2y3y4y5)=(100200),得(y1y2y3y4y5)=(256.66671.785700.4784753.718)再由(1−141−4151−15561−562091)(x1x2x3x4x5)=(256.66671.785700.4784753.718),得(x1x2x3x4x5)=(27.0518.20525.769314.87253.718)线性方程组解法(迭代法)2-1 设线性方程组{4x 1−x 2+2x 3=1−x 1−5x 2+x 3=22x 1+x 2+6x 3=3(1) 写出Jacobi 法和SOR 法的迭代格式(分量形式) (2) 讨论这两种迭代法的收敛性(3) 取初值x (0)=(0,0,0)T ,若用Jacobi 迭代法计算时,预估误差 ||x*-x (10)||∞ (取三位有效数字)解:(1)Jacobi 法和SOR 法的迭代格式分别为Jacobi 法迭代格式SOR(2)因为A 是严格对角占优矩阵,但不是正定矩阵,故Jacobi 法收敛,SOR 法当0<ω≤1时收敛.⎪⎪⎪⎩⎪⎪⎪⎨⎧+--=-+-=+-=+++216131525151412141)(2)(1)1(3)(3)(1)1(2)(3)(2)1(1k k k k k k k k k x x x x x x xx x ⎪⎪⎪⎩⎪⎪⎪⎨⎧-++-=+-+-=+-+-+=++++++)216131()525151()412141()(3)1(2)1(1)(3)1(3)(3)(2)1(1)(2)1(2)(3)(2)(1)(1)1(1k k k k k k k k k k k k k k k x x x x x x x x x x x x x x x ωωω(3)由(1)可见||B ||∞=3/4,且取x (0)=(0,0,0)T ,经计算可得x (1)=(1/4,-2/5,1/2)T ,于是||x (1)-x (0)||∞=1/2,所以有2-2 设方程组为{5x 1+2x 2+x 3=−12−x 1+4x 2+2x 3=202x 1−3x 2+10x 3=3试写出其Jacobi 分量迭代格式以及相应的迭代矩阵,并求解。
2013研究生数值分析试题参考答案
太原科技大学硕士研究生2013/2014学年第1学期《数值分析》课程试卷公式提示:1、Legendre 多项式)(x p n 的递推关系式:⎪⎪⎩⎪⎪⎨⎧=+-++===-+,2,1)(1)(112)()(,1)(1110n x p n n x xp n n x p x x p x p n n n2、Chebyshev 多项式)(x T n 的递推关系式:⎪⎩⎪⎨⎧=-===-+ ,2,1)()(2)()(,1)(1110n x T x xT x T x x T x T n n n一、填空题(每小题5分,共35分)1、为提高数值计算精度,当x 充分小时,应将x cos 1-改写为___22sin2x___ 2、已知x=0.03056是经过四舍五入得到的近似值,则x 有____5__位有效数字.3、设⎪⎪⎭⎫⎝⎛=4321A ,则=∞)(A Cond ___21___.4、已知()sin 1f x x x =--,则牛顿法的迭代公式是_____1sin 1,0,1,2,...1cos k k k k kx x x x k x +--=-=-__________5、求解非线性方程310x x --=的一个收敛的简单迭代公式为_______10,1,2,...k x k +==________。
6、设,,2,1,0,,53)( ==+=k kh x x x f k 则=++],,[21n n n x x x f _______3h________。
7、若用Gauss-Seidel 迭代法解方程组⎩⎨⎧-=+=+3242121x ax ax x ,其中a 为实数,则Gauss-Seidel迭代法收敛的充要条件是应使a 满足______a <<_________。
二、(本题满分15分)(1)用列主元Gauss 消去法求解下列方程组:⎪⎩⎪⎨⎧=++=++=-+3221522321321321x x x x x x x x x (2)写出用Jacobi 迭代求解上述方程组的迭代公式的分量形式。
硕士研究生《数值分析》试卷2013(A)
硕士研究生《数值分析》试卷2013(A)一、判断题 (下列各题,你认为正确的,请在题后的括号内打“√ ”,错误的打“×”,每题2分,共10分) 1. 近似数*3.200x =关于准确值 3.200678x =有4位有效数字。
( ) 2. 设(0,1,2,3)i x i =是互异的点,()(0,1,2,3)i l x i =是Lagrange 插值基函数,则3224()4i ii x l x x==∑. ( )3. 设73()32f x x x =-+,则差商1234567[2,2,2,2,2,2,2]1f =。
( ) 4. 设A 是n 阶非奇异方阵,则解方程组A =x b 的迭代法收敛的充要条件是A 的谱半径()1A ρ<。
( )5. 解常微分方程初值问题的四阶Runge-Kutta 方法的整体截断误差是4()O h ,其中h 是步长。
( )二、填空题 (每空2分,共16分) 1. 设T(2,1,3,4)=-x ,2543A -⎡⎤=⎢⎥-⎣⎦. 则 1||||x = , Cond()A ∞= .2. 设20()d I f x x =⎰,若用梯形求积公式计算I ,结果是4;用Simpson 求积公式计算I ,结果是2. 则(1)f = .3. 设S 是函数f 在区间[0,3]上满足第一类边界条件的的三次样条:()()22,01,()111,13,2x x S x x a x b x ⎧≤≤⎪=⎨-+-+≤≤⎪⎩ 则a = ,b = ,(3)f '= .4. 设函数(0.8) 1.2,(0.9) 1.4,(1) 1.0,(1.1)0.2,(1.2)0.5f f f f f =-=-=-==, 步长0.2h =,则用三点数值微分公式计算(1)f '的近似值为 .5. 设函数()f x 是最高次项系数为1-的3次多项式,2()p x 是()f x 在节点1,0,1-上的Lagrange 插值多项式, 则余项2()()f x p x -= .三(本题满分8分)的近似值*x 的相对误差限是0.01%,求*x 至少应具有几位有效数字?四(本题满分10分) 对下列方程组分别建立收敛的Jacobi 和Gauss-Seidel 迭代格式,并说明理由。
2011数值分析试题及答案
由于f(x)二si nx的4阶导数在[0,二]上的最大值为:M4=1,所以
5
误差为:|I-S2|::——44=0.006641
2880x24
6.求解初值问题」y=sin(x+2y),0兰x兰2的改进Euler方法是否收敛?为什
.y(0) = 1
么?
解:由于|sin(x 2y)-sin(x 2y)|二| 2cos(x 2 )(y-y) 2 | y-y |
5.设f(x) = 4x33x-5,求差商f[0,1], f[1,2,3,4]和f[1,2,3,4,5]。
f(D…f(0)
解:f[0,1]==2-(-5) = 7
1-0
f [1,2,3,4^4,f[1,2,3,4,5]=0
3.解线性方程组丿X1-2忑=2的Jacobi迭代法是否收敛,为什么?
+9x2=3
即,函数f(x, y)二sin(x•2y)连续,且关于变量y满足Lipschitz条件,所以,改 进Euler方法收敛。
所以,a=0, b=5/6,拟合曲线为:y=5/6x2
3.求满足条件f(0)=1,f(1)=2,f(2) =0,f(1)=0的三次插值多项式Ha(x)
的表达式。
解:设H3(x)二(^2)(ax2bx c),则有:
1213
所以,H3(x) (x-2)(x2x 1) (x-3x-2)。
22
11
4.确定求积公式Jf(x)dx痒三f(-1)+Af(0)+A2f(1)中的待定系数,使其代数精 度尽可能高,并问此公式是不是插值型求积公式.
解:令公式对f(x) = 1,x都精确成立,得:A,・A2= 3/2, A2= 1/2,
o
• • •
数值分析上机题参考答案.docx
如有帮助欢迎下载支持数值分析上机题姓名:陈作添学号: 040816习题 120.(上机题)舍入误差与有效数N11 3 1 1设S N,其精确值为 。
22 2 NN 1j 2j1(1)编制按从大到小的顺序111 ,计算 S 的通用程序。
S N1 321N 21 N22(2)编制按从小到大的顺序111,计算 S 的通用程序。
S N1(N 1)2 122 1NN 2(3)按两种顺序分别计算S 102 , S 104 , S 106 ,并指出有效位数。
(编制程序时用单精度)(4)通过本上机题,你明白了什么?按从大到小的顺序计算 S N 的通用程序为: 按从小到大的顺序计算 S N 的通用程序为:#include<iostream.h> #include<iostream.h> float sum(float N) float sum(float N) {{float j,s,sum=0; float j,s,sum=0; for(j=2;j<=N;j++) for(j=N;j>=2;j--) {{s=1/(j*j-1); s=1/(j*j-1); sum+=s;sum+=s;}}return sum;return sum;}}从大到小的顺序的值从小到大的顺序的值精确值有效位数从大到小从小到大0.7400490.740050.74004965 S 1020.7498520.74990.74994 4S 1040.7498520.7499990.74999936S 106通过本上机题, 看出按两种不同的顺序计算的结果是不相同的,按从大到小的顺序计算的值与精确值有较大的误差, 而按从小到大的顺序计算的值与精确值吻合。
从大到小的顺序计算得到的结果的有效位数少。
计算机在进行数值计算时会出现“大数吃小数”的现象,导致计算结果的精度有所降低,我们在计算机中进行同号数的加法时,采用绝对值较小者先加的算法,其结果的相对误差较小。
数值分析上机实习题
数值分析上机实习题第2章插值法1. 已知函数在下列各点的值为试⽤四次⽜顿插值多项式)(x p 4及三次样条韩式)(S x (⾃然边界条件)对数据进⾏插值。
⽤图给出(){}10,11,1,0,08.02.0,,x i =+=i x y i i ,)(x p 4及)(x S Python 代码import numpy as npimport matplotlib.pyplot as pltfrom matplotlib.font_manager import FontPropertiesfont_set = FontProperties(fname=r"c:\windows\fonts\simsun.ttc",size=12) #求⽜顿n 次均差 def qiujuncha(x,f,n): for i in range(1,n): for j in range(4,i-1,-1):f[j]= (f[j] - f[j-1])/(x[j]-x[j-i]) #根据⽜顿多项式求值 def niudun(x,f,x1): sum = f[0]; tmp = 1;for i in range(1,5): tmp *= (x1-x[i-1]) sum = sum + f[i]*tmp return sum#⽜顿插值画图 def drawPic(x,f):x1 = np.linspace(0.2, 1, 100) plt.plot(x1, niudun(x,f,x1))plt.title(u"⽜顿四次插值",fontproperties=font_set) plt.xlabel(u"x 轴",fontproperties=font_set) plt.ylabel(u"y 轴", fontproperties=font_set) plt.show() def qiu_h(x,h): n = len(x) -1 for i in range(n): print(i)h[i] = x[i+1]-x[i]#⾃然边界条件下的三次样条插值求Mdef qiu_m(h,f,o,u,d):n = len(h)o[0] = 0u[n] = 0d[n] = d[0] = 0a = []for i in range(1,n):u[i] = h[i-1]/(h[i-1]+h[i])for i in range(1,n):o[i] = h[i]/(h[i-1]+h[i])for i in range(1,n-1):d[i] = 6*(f[i+1]-f[i])/(h[i-1]+h[i])t = [0 for i in range(5)]t[0] =2t[1] = o[0]a.append(t)for i in range(1,n):t = [0 for i in range(5)]t[i - 1] = u [i + 1]t[i] = 2t[i + 1] = o [i + 1]a.append(t)t = [0 for i in range(5)]t[n - 1] = u[n]t[n] = 2a.append(t)tmp = np.linalg.solve(np.array(a),np.array(d))m = []for i in range(5):m.append(tmp[i])return m#根据三次条插值函数求值def yangtiao(x1,m,x,y,h,j):returnm[j]*(x[j+1]-x1)**3/(6*h[j])+m[j+1]*(x1-x[j])**3/(6*h[j])+(y[j]-m[j]*h[j]**2/6)*(x[j+1]-x1)/h[j] +(y[j+1]-m[j+1]*h[j]**2/6)*(x1-x[j])/h[j] def main():x = [0.2, 0.4, 0.6, 0.8, 1.0]y = [0.98, 0.92, 0.81, 0.64, 0.38]f = y[:]f1 = y[:]h = [0.2,0.2,0.2,0.2]u = [0 for n in range(5)]d = [0 for n in range(5)]o = [0 for n in range(5)] qiujuncha(x,f,4) qiujuncha(x,f1,2)m = qiu_m(h,f1,o,u,d) x1 = np.linspace(0.2, 0.4, 10)p1= plt.plot(x1, yangtiao(x1,m,x,y,h,0),color='red') x1 = np.linspace(0.4, 0.6, 10)plt.plot(x1, yangtiao(x1, m, x, y, h, 1),color='red') x1 = np.linspace(0.6, 0.8, 10)plt.plot(x1, yangtiao(x1, m, x, y, h, 2),color='red') x1 = np.linspace(0.8, 1.0, 10)plt.plot(x1, yangtiao(x1, m, x, y, h, 3),color='red') x1 = np.linspace(0.2, 1.0, 40)p2 = plt.plot(x1,niudun(x,f,x1),color='green') plt.xlabel(u"x 轴", fontproperties=font_set) plt.ylabel(u"y 轴",fontproperties=font_set) plt.title("三次样条插值和⽜顿插值")plt.legend(labels=[u'三次样条插值',u'⽜顿插值'],prop=font_set,loc="best") plt.show() main()实验结果运⾏结果可得插值函数图(如图1-1),4次⽜顿插值函数)(x p 4和三次样条插值函数)(x S 如下:)6.0(*)4.0(*)2.0(625.0)4.0(*)2.0(*3.098.0)(4-------=x x x x x x x P 98.0)8.0(*)6.0(*)4.0(*)2.0(*20833.0+-----x x x x]4.0,2.0[),2.0(467.4)4.0(9.4)2.0(167.1)(S 3∈-+-+-=x x x x x]6.0,4.0[),4.0(113.4)6.0(6467.4)4.0(575.1)6.0(167.1)(S 33∈-+-+----=x x x x x x ]8.0,6.0[),6.0(2.3)8.0(113.4)6.0(575.1)(S 3∈-+-+--=x x x x x]0.1,8.0[),8.0(9.1)0.1(2.3)(S ∈-+-=x x x x图1-1三次样条插值和⽜顿插值图2.在区间[-1,1]上分别取n = 10,20⽤两组等距节点对龙格函数做多项式插值三次样条插值,对每个n值画出插值函数及图形。
北理工数值分析考试卷-2013级数值分析考试习题
2013——2014年《数值分析》考试习题一、(15分)设2101202a a ⎡⎤⎢⎥A =⎢⎥⎢⎥⎣⎦ ,123x x x x ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ ,331b ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦(1)求出当a 满足何条件时可以对A 进行cholesky 分解(即T LL 分解,L 为非奇异下三角矩阵)(2)若a=1时利用T LL 分解对方程组Ax=b 进行求解(保留两位小数)二、(10分)求证:111a a a a a a ⎡⎤⎢⎥A =⎢⎥⎢⎥⎣⎦,只有当1122a -<< 时用Jacobi 迭代法求解方程组Ax=b 才收敛 三、(15分) (1)已知i x1 3 4 6 ()i f x-75814对其进行3次插值,并求出()2f 的值 (2)给定函数组i x0 1 2 ()i f x 1 4 15 ()'i f x6对其进行3次Hermite 插值多项式3H 并求出余项()()f x H x - 的表达式四、(15分)利用最小二乘法,求形如2y a bx =+ 的拟合i x 1 2 3 4 5 i y0.12.98.114.924.1五、(15分)(1)选择合适的数值积分方法,对210x e dx ⎰ 进行积分,使其具有2位有效数字 (2)根据(1)中选取的数值积分方法,编写Matlab 源程序代码 六、(15分)(1)构造迭代格式求解40x x e --= 的根,使其具有4位有效数字 (2)从理论上证明你所构造的迭代格式是收敛的七、(15分)对常微分方程初值问题()()'00,y f x y y x y ⎧=⎪⎨=⎪⎩ ,证明其二步公式111(58)12n n n n n hy y f f f ++-=++- 是三阶精度的迭代公式,并求出该二步公式的局部截断误差主项。
数值分析(2011)试题A卷 参考答案
装订线年 级 学 号 姓 名 专 业一、填空题(本题40分, 每空4分)1.设),,1,0()(n j x l j =为节点n x x x ,,,10 的n 次基函数,则=)(i j x l 1,0,1,,0i j i j n i j=⎧=⎨≠⎩ 。
2.已知函数1)(2++=x x x f ,则三阶差商]4,3,2,1[f = 0 。
3.当n=3时,牛顿-柯特斯系数83,81)3(2)3(1)3(0===C C C ,则=)3(3C 1/8 。
4.用迭代法解线性方程组Ax=b 时,迭代格式 ,2,1,0,)()1(=+=+k f Bx x k k 收敛的充分必要条件是 ()1B ρ< 。
5.设矩阵⎥⎦⎤⎢⎣⎡=1221A ,则A 的条件数2)(A Cond = 3 。
6.正方形的边长约为100cm ,则正方形的边长误差限不超过 0.005 cm才能使其面积误差不超过12cm 。
(结果保留小数)7.要使求积公式)()0(41)(111x f A f dx x f +≈⎰具有2次代数精确度,则 =1x23 , =1A 34。
8. 用杜利特尔(Doolittle )分解法分解LUA =,⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=1359 45- 279 126 0 945- 0 45 1827- 9 189A 其中,则=L 10002100121023113⎛⎫⎪ ⎪ ⎪-⎪ ⎪- ⎪⎝⎭=U 918927091890281540009-⎛⎫⎪-⎪ ⎪-⎪⎝⎭。
二、计算题(10分)已知由数据(0,0),(0.5,y ),(1,3)和(2,2)构造出的三次插值多项式)(3x P 的3x 的系数是6,试确定数据y 。
2011级数值分析 试题 A 卷 2011 ~ 2012学年,第 1 学期一 二 三 四 五 六 七 八 九 十 总分年 级2011级研究生 份 数 拟题人 王吉波 审核人装 订线年级 学 号 姓 名 专 业三、计算题(15分)试导出计算)0(1>a a的Newton迭代格式,使公式中(对n x )既无开方,又无除法运算,并讨论其收敛性。
——数值分析上机题
.......................课程名称:数值分析上机实习报告姓名:学号:专业:联系电话:目录序言 (3)第1章必做题 (4)1.1必做题第一题 (4)1.1.1题目 (4)1.1.2 分析 (4)1.1.3 计算结果 (4)1.1.3 总结 (6)1.2必做题第二题 (6)1.2.1题目 (6)1.2.2分析 (6)1.2.3计算结果 (6)1.2.4结论 (8)1.1必做题第一题....................................................................... 错误!未定义书签。
1.1.1题目 ............................................................................ 错误!未定义书签。
第2章选做题 (8)2.1选做题第一题 (8)2.1.1题目 (8)2.1.2分析 (8)2.1.3计算结果 (8)附录 (10)附录一:必做题第一题程序 (10)附录二:必做题第二题程序 (11)附录三:选做题第一题的程序 (13)序言本次数值分析上机实习采用Matlab数学软件。
Matlab是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。
在数值分析应用中可以直接调用Matlab软件中已有的函数,同时用户也可以将自己编写的实用程序导入到Matlab函数库中方便自己调用。
基于Matlab数学软件的各种实用性功能与优点,本次数值分析实习决定采用其作为分析计算工具。
1.编程效率高MATLAB是一种面向科学与工程计算的高级语言,允许使用数学形式的语言编写程序,且比BASIC、FORTRAN和C等语言更加接近我们书写计算公式的思维方式,用MATLAB编写程序犹如在演算纸上排列出公式与求解问题。
因此,MATLAB语言也可通俗地称为演算纸式科学算法语言。
数值分析上机题答案
数值分析上机题答案【篇一:数值分析上机试题对应参考答案】么是近似值x* 有效数字?若近似值x*的误差限是某一位的半个单位,该位到x*的第一位非零数字共有n位,就说x*有n位有效数字。
它可表示为2、数值计算应该避免采用不稳定的算法,防止有效数字的损失. 因此,在进行数值运算算法设计过程中主要注意什么?(1)简化计算过程,减少运算次数;(2)避免两个相近的数相减;(3)避免除数的绝对值远小于被除数的绝对值;(4)防止大数“吃掉”小数的现象;(5)使用数值稳定的算法,设法控制误差的传播。
3、写出“n 阶阵a 具有n 个不相等的特征值”的等价条件(至少写3 个)(1)|a|不为零(2)n阶矩阵a的列或行向量组线性无关(3)矩阵a为满秩矩阵(4)n阶矩阵a与n阶可逆矩阵b等价4、迭代法的基本思想是什么?就是用某种极限过程去逐步逼近线性方程组精确解得方法。
其基本思想为:先任取一组近似解初值x0,然后按照某种迭代原则,由x0计算新的近似解x1,以此类推,可计算出x2,x3,…xk,。
,如果{x}收敛,则取为原方程组的解。
5、病态线性方程组的主要判断方法有哪些?(1)系数矩阵的某两行(列)几乎近似相关(2)系数矩阵的行列式的值很小(3)用主元消去法解线性方程组时出现小主元(4)近似解x*已使残差向量r=b-ax*的范数很小,但该近似解仍不符合问题要求。
6、lagrange 插值的前提条件是什么?并写出二次lagrange 插值的基函数。
1,j?i?(x)? 前提条件是:l i ,j?0,1,2?,n.?ij0,j?i?二次lagrange 插值的基函数: (x?x)(x?x)12??lx0(xx)(xx) 0?10?2 (x?x)(x?x)02?? lx1(xx)(xx)1?01?2(x?x)(x?x)01?? lx2(x?x)(x?x)20217、什么是数值积分的代数精度?如果某一个求积公式对于次数不超过m的多项式均能准确地成立,但对于m+1次多项式就不准确成立,则称该求积公式具有m次代数精度(或代数精确度)。
研究生数值分析2012-2013试卷
山东科技大学 2012-2013 学年第一学期《数值分析》考试试卷[][][][]其收敛阶出的牛顿迭代格式,并指的写出求方程为正数,记设六、计算题插值多项式。
的三次写出时,已知当五、计算题差。
多项式,并估计平方误上的一次最佳平方逼近在区间求设函数四、计算题一个复化求积公式利用该求积公式构造等分,并记作)将区间(斯型,并说明理由;)判断该公式是否为高(数精度;代数精度,并指出其代,使其具有尽可能高的试确定求积系数给定求积公式:三、计算题差限。
的绝对误差限与相对误位有效数字,试分析均具有设二、计算题。
计算、设的值。
与计算、设一、计算题**211212121710610360)(,,)(ewton )(,5,2,3,1)(5,3,2,020)(,)(,,,1,0,1,2n 11-32,,)1()1()0()1()(365.3,1.12,,,,723226131,1252222222,13)(1x x f a x a a x x f N x f x f x x f x x f n k i ih x nh C B A Cf Bf Af d x f x x x x A A x x A x L f L f x x x f n n i x F ==-=====+-==++-≈+==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=++=⎰-∞并指出其精度。
写出改进的欧拉公式,,记取正整数题考虑常微分方程初值问八、计算题消去法求方程组的解。
用列主元迭代格式的敛散性;试分析迭代格式。
迭代格式与写出给定线性方程组七、计算题.0,,n ,)(),,(auss )3(eidel -auss )2(eidel -auss acobi )1(215702031-22-1'321n i ih a x n a b h a y b x a y x f y G S G S G J x x x i ≤≤+=-=⎩⎨⎧=≤≤=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡η。
数值分析上机实习题及答案.docx
方詡文金兴:爭[数值分析]2017-2018第二学期上机实习题1:编程计算亍丄,其中C= 4. 4942x10307,给出并观察计算结心C"果,若有问题,分析之。
解:mat lab 编程如下:E) funct ion diy i ti formatlong g;n 二input C 输入ii 值= c= 4.4942E307; sum 二0; s 二 0;E3 for i = l:n s = l/ (c*i);>> diyiti 输入n 值:10 104.6356e-308 >> diyiti输入ri 值:1001004.6356e-308 >> diyiti 输入n 值:1000 10004.6356e-308 >> diyiti揄入n 值* 1000001000004・ 6356e-308 >> diyiti输入n 值;1000000001000000004.6356e-308图二:运行结果Mat lab 中,forma t long g 对双精度,显示15位定点或浮点格式,由上图 可知,当输入较小的n 值5分别取10, 100, 1000, 100000, 100000000)的时候, 结果后面的指数中总是含有e-308,这和题目中的C 值很相似,我认为是由于分 母中的C 值相对于n 值过大,出现了 “大数吃小数”的现彖,这是不符合算法原 则的。
2:利用牛顿法求方程-1^ = 2于区间241的根,考虑不同初值下牛顿法的收敛情况。
解:牛顿法公式为:利用mat lab 编程function di2ti21 3i=l ;2 2.85208156699784 xO 二input ('输入初值x0:‘ );A 二[i x0];3 2.55030468822809 t=x0+ (x0-log (xO) -2) /(1-1/xO) ; %迭代函数4 1.91547247100476 三 while (abs (t _x0)>0.01)i=i+l; 5 0.37867158538991 xO 二 t; 6 0.774964959780275 A = [A;i xO];t =x0+(x0-log(xO)-2)/(1-1/xO): 7 4.11574081641933 cnd| 8 5.04162436446126 disp (A);96.81782826645596当输入初值二3的时候并不能收敛。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
数值分析上机实践题
第 11次上机题目(Guass 型数值求积公式)
注意:要把一般的把区间[a,b]变换成[-1,1],其公式是:
第一组: 组长:李龙宇,组员:杜彦霖,胡朋,黄湘云,雷盛华,李伟元 利用两点、三点、四点Guass-Legendre 求积公式计算下列积分
(1) ⎰+102
4dx x x ,并与准确值比较精度; (2) 求⎰-102
dx e x 的近似值。
第二组: 组长:王宇彬,组员:马泽川,权涛涛,师楠颉,路世伦,仲晓磊 利用两点、三点、四点Guass-Legendre 求积公式计算下列积分
(1) ⎰-6
02sin 4π
dx x ,并与准确值比较精度;
(2) 求⎰-10
2dx e x 的近似值
第三组: 组长: 薛原 ,组员:谢胜权,杨帆,王正奇,肖特,张锡云 利用两点、三点、四点Guass-Legendre 求积公式计算下列积分
(1) ⎰-10
2dx e x π,并与准确值比较精度并估计误差; (2) 求⎰-102
dx e x 的近似值
第四组: 组长:柴春晓 ,组员: 韩静兰,李金慧,刘从,马超群,孟凯悦 利用两点、三点、四点Guass-Legendre 求积公式计算下列积分
(1) ⎰311dx x
,并与准确值比较精度并估计误差; (2) 求⎰-102
dx e x 的近似值。
第五组: 组长:龙纯鹏,组员:代喜,白鑫,鲍亚强,周邦安,张佳伟 利用两点、三点、四点Guass-Legendre 求积公式计算下列积分
(1) ⎰+3
021dx x x ,并与准确值比较精度并估计误差;
(2) 求⎰-10
2dx e x 的近似值。
第六组: 组长:何关瑶 ,组员:纪伟亮,侯佳意,李济言,李振华,马文磊 利用两点、三点、四点Guass-Legendre 求积公式计算下列积分
(1) ⎰π
20sin xdx x ,并与准确值比较精度并估计误差;
(2) 求⎰-10
2dx e x 的近似值 第七组: 组长:杨钦 ,组员: 王凌宇,吴凯杰,薛小龙,袁权炜,师俊峰 利用两点、三点、四点Guass-Legendre 求积公式计算下列积分
(1) ⎰10
2dx e x x ;并与准确值比较精度并估计误差; (2) 求⎰-1
02dx e x 的近似值,。
第八组: 组长:汪芳 ,组员:张学利,周幸茹,李雨珏,张飞
利用两点、三点、四点Guass-Legendre 求积公式计算下列积分
(1)
⎰2
1ln xdx ,并与准确值比较精度并估计误差; 提示: lnx Matlab 表示log(x)
(2) 求⎰-102dx e x 的近似值。
第九组: 组长:刘永鸿 ,组员:黄尚政,李超,郭新磊,何奎奎
利用两点、三点、四点Guass-Legendre 求积公式计算下列积分
(1) ⎰432cos ππxdx ;并与准确值比较精度并估计误差;
(2) 求⎰-10
2dx e x 的近似值。
第十组: 组长:杨吉望 ,组员:龙力,任金雄,王亮,王文强,谢丁波 利用两点、三点、四点Guass-Legendre 求积公式计算下列积分
(1) ⎰5
.11
.1dx e x ,并与准确值比较精度并估计误差; (2) 求⎰-10
2dx e x 的近似值
第十一组: 组长: 张国强,组员: 赵奇,袁硕,郭凯旋,于沛生,鲍宏雷 利用两点、三点、四点Guass-Legendre 求积公式计算下列积分(1) ⎰-312dx e x ,
并与准确值比较精度并估计误差;
(2) 求⎰-102dx e x 的近似值
第十二组:组长:苏映雪,组员:邓晓庆,钟桂平,崔楚轩,高鹏程利用两点、三点、四点Guass-Legendre求积公式计算下列积分
(1) ⎰3.003dx
x,并与准确值比较精度并估计误差;
提示:3x Matlab表示power(x,1/3)
(2) 求⎰-102dx
e x的近似值。