超定方程求解 c语言,北航数值分析复习资料2018最新整理_假设对函数fx在步长为h的等距点上造表,且f,北航数值分析复习资料2018最新整理-讲义文档类资源...

北航数值分析复习资料最新整理版本，包含网上搜到的所有往年试题，还有复习提纲、习题课笔记，可以直接下载打印使用，2018年12月26号新鲜出炉

5假设对f(x)以等距h在节点上造表,且有(x)≤M,

(1)证明:任意相邻两节点上的线性插值误差不超过Mh2;

(2)设f(x)=sinx间问取多大能使线性插值误差s号×10°

答:5(2),h≤2×10-3

6试由(x)=2的函数表

0512

建立二次插值多项式p2(x用以求2的近似值并估计误差

证两位有效数字

拉o牛

答:n(x)=0252+075+123≈2(03)=12475

293-p2(0.3)3(3+)0.3-0(03-1)-0.0300

Axf"(x)=2(n2)3=0.60

ISx

≤

7、P59,4.

P59,6,8.

8设f(x)=x+x24+3x+1,求1202,…,2和/20,2…,27]

答:(①),(2)0

9设T(x)是次切比雪夫多项式,证明

(1 TIT(x)]=Tmn(r);

(2)Tm(x)+T-(x=2Tm (x),(x);

(3)Tn(2x2-1)=T2n(x)

10、求(x)=x2+3x-5在[-1,1]上的线性最佳一致逼近多项式

解:(x)-(x)=-(x),p(x)=f(x)-172(x)=3x-

1、求f(x)= arcsinx在[-1,1]上的切比雪夫级数(n=7)

解:P1(x)=2+∑a7(x,xE1

其中

2k(x) arcsin x.奇

2k

x

2 l T2k+ (x)arcsinx

k+l

、2c0co[(2k+1)

(sin ede

sing

4

sin(2k+1)pdB

2k+1)r

x(2+1)

Pn(x)=()+-3(x)+T(x)+71(x),x∈[1,

25

49

P,(x)

4「762483288564

x十

x

丌10531517549

+x2x∈[1,1

P115,1,4(2),6,8,13,15,17(1),19,按基本方法即

第4-5章习题课

(数值积分和数值微分,解线性方程组的直接法)

、数值积分与数值微分

求积公式:∫f(xdx=吗

若一个求积公式对于所有次数不超过m的多项式都准确

成立,而对于某一个m+1次的多项式等式不准确成立,则

称该求积公式具有m次代数精度

由拉格朗日插值L2(x)=∑l2(x),得到求积公式

f(x)dxc∑vf,其中v4=["1(x)dx

称为插值型求积公式

余项:R=L/(x)-L1(x)z=f

(m+(2)

(x-I di

=0

定理求积公式∫(x)dx=w至少具有n次代数精度

分它是插值型求积公式

将求积区间[a,做n等分,步长h=“,在等距节点

x=a+hh上的插值型求积公式

I /(ydx(b-a)2ea'5k

称为 Newton- Cotes公式,C)称为 Cotes系数,

作变换x=a+h,则有

CIm

h

-1)7一h

dr

b

k-j nk (n-k)!30 i-o

当n=时,得到梯形公式

f(x)dx7="[f(a)+/(b)

当n=2时,得到抛物线公式,也称为辛普森( Simps n)公式

b-ar

a十

f(x)dr≈S=-[f(a)+4(2

)+f(b),(2.3)

当n=4时,得到柯特斯(tes)公式

C=[7f(x0)+32f(x1)+12(x2)+32f(x)+7f(x4

其中x=a十h,h=二日

2

4

柯特斯系数表.n≥8时C)出现负值,N-C公式不稳定

若/"(x)在[a2b]上连续,则梯形公式的余项为

R门=1-=-b2)r"(m,n∈nb

若(x)在[ab]上连续,则辛普森公式的余项为

尼U门=1-s=-(x:bU(+4

+b

(。)+∫(b

b-a b-a

l80

f(n), nE[a, b

h

h

∑[f(x)+f(x+=[f(a)+2∑f(x)+f(b)

i=0

-T∑[1h3f"(m)=-nfm)≈、b-a

h fn)

i=0

12

12

Sn=2f(a)+4f(x)+2∑f(x)+f(6

h

(2)

∑f1(m,)=

3-a

h“f((n)n∈(a,b

2880

(1)初值1

[f(a)+f(b)]

2

(2)令h

=0,1,2,…),计算

2

T,=1x,+8(x+

(3)求加速值

S=T+(T-T)/3

Cu=S,mn+(S2u-S,)/15,

R。-C2n+(C2n-C,)/63

(4)满足精度要求:否则,转(2)

若一组节点a≤x

广x)()>w/(x),

j=0

具有2n+1次代数精度,则称此组节点为高斯点,并称此

求积公式为高斯求积公式

定理插值型求积公式的节点a≤x

高斯点an1(x)=(x-x)x-x)…(xx)与任何次数

不超过n的多项式P(x带权p(x)交即

P(xo(x)P(x)dx=0

R[f]=

fg(c:a1()(

)dx,∈[a,b]

(2n+2)!

x

0

h

[/(x)-f(x0)]-2f"(5

f(x)=[(x)-f(x)小

h

f"(

h

f()=2h

[-3f(x0)+4f(x1)-f(x2

+f"(5)

f(x1)=,[-f(x)+f(x2)-f"(

f(x2)=2[(x)-4(x)+3(x2)+f

h2

f(x)=2[f(x0)

-2+f(x2)12(5)

3.试用n=5的复伴普森公式求=的近似值,

1+x

并估计误差

dx

准确值

n2=693147.

1+x

解:h=(1-0)==02,x=0+边x=0+(+)h

5

0.21

2×(

1+0

+0.21+0.41+061+08

0.1]+0.31+0.51+071+0.91+1

0.693150

24

2

(4

(+x),1ls。ma义/“(x)=1.330

TEL

I-S3=-2.81944×10

4若用复合梯矜公式计算1=[;,问区间多少等分

1+x

才能保证计算结果有五位有效数字?

解:h=-(1-0)=-,x1;=0+i(=01…,n)

?

maX

J"(x)=f(0)=2

1+x

(1+x)

R

1-0

≤,h2maxf"(x)

12

0≤

6(n

OX

因05<

<1

有一位小数,

J1+

1+x

故只需

±×10-5,n>√105/3=182.574,即n≥183

6(n

思考:试用复合梯形式计算椭圆+y2=的周长,

保证计算结果有五位有效数字

1=V*g+y2d8=4sin 0+cos28d8= \1+3sin28d8

7=23561945,T2=24199207872-7P=00212421

7=242210310,T-T2=0.0002744

7=2422l1206,T2-72=000986

I≈T=2421有五位有效数字

5求形如

f(2)-dx≈af(x)+af(x两点高斯型

求积公式

解法1:72(x)=2x2-1,x=千,

以及

+w1=丌

→

W+2w1=

0

x*心号(2)+受f)

解法2:由高斯一切比雪夫求积公式

兀

2i-1

dx≈

∑

coSI

z

x

n

2n

得到

f(x)

dx≈号f(号)+号f(2)

x

6.根据y=2的数据表

0.5

1.5

2

25

yi141421

200000

2.82843

400000

5.65685

试用二点、三点数值微分公式计算f(15)和”(1.5)

解:取h=05,则

f(1.5)≈

051(.5)-/()】=-1.5686

f05)=2X05U2)-0=2

f"(1.5)≈

f(2)-2f(1.5)+f(]=034314

三、解线性方程组的直接方法

基本内容及基本要求

1.了解求解方程组的两类方法了解矩阵基础知识

2.掌握高斯消去法会矩阵的三角分解。

3.掌握高斯列主元素消去法了解高斯若当消去法。

4.掌握直接三角分解法,了解平方根法会追赶法了解有关结论。

了解向量和矩阵的几种范数

6.了解矩阵和方程组的性态会求其条件数。

7.会初等反射阵和平面旋转阵,了解QR分解了解用正交约化法解超定方程组。

四、练习

分别用顺序Gaus消去法和直接三角分解法(杜利脱尔分解)求解线性方程组

123Tx]「14

F卓

252‖x2|=18

315x」L20

23141「12314112314

解)25218|→01-4-101→01-4-10个x2=2,

315200-5-4-00-24-72x=3

2

解Ly=b得y=14-10,-72)

2)25

3-51

24

解Uk=y得x=(123)