在科學(xué)計(jì)算中的應(yīng)用4PPT學(xué)習(xí)教案

1、會(huì)計(jì)學(xué)1在科學(xué)計(jì)算中的應(yīng)用在科學(xué)計(jì)算中的應(yīng)用4A=rand(n)第1頁(yè)/共93頁(yè)第2頁(yè)/共93頁(yè)第3頁(yè)/共93頁(yè)第4頁(yè)/共93頁(yè)第5頁(yè)/共93頁(yè)第6頁(yè)/共93頁(yè)第7頁(yè)/共93頁(yè)第8頁(yè)/共93頁(yè)第9頁(yè)/共93頁(yè)第10頁(yè)/共93頁(yè)第11頁(yè)/共93頁(yè)第12頁(yè)/共93頁(yè)第13頁(yè)/共93頁(yè)第14頁(yè)/共93頁(yè)N=norm(A，選項(xiàng))選項(xiàng)可為1,2,inf等第15頁(yè)/共93頁(yè)第16頁(yè)/共93頁(yè)第17頁(yè)/共93頁(yè)第18頁(yè)/共93頁(yè)格式：P=sym2poly(f)第19頁(yè)/共93頁(yè)第20頁(yè)/共93頁(yè)第21頁(yè)/共93頁(yè)第22頁(yè)/共93頁(yè)第23頁(yè)/共93頁(yè)第24頁(yè)/共93頁(yè)第25頁(yè)/共93頁(yè)*a1*a4+a2*a3

2、)第26頁(yè)/共93頁(yè)第27頁(yè)/共93頁(yè)第28頁(yè)/共93頁(yè)第29頁(yè)/共93頁(yè)格式：x=Ab第30頁(yè)/共93頁(yè)12341234123412340.40960.12340.36780.29430.40430.22460.38720.40150.11290.11500.36450.19200.37810.06430.42400.17840.40020.27860.39270.2557xxxxxxxxxxxxxxxx第31頁(yè)/共93頁(yè)第32頁(yè)/共93頁(yè)第33頁(yè)/共93頁(yè)第34頁(yè)/共93頁(yè)第35頁(yè)/共93頁(yè)第36頁(yè)/共93頁(yè)123344232.53.51001xxxxxx第37頁(yè)/共93頁(yè)ans = 4.

3、4409e-015121234230.95422.53.50.732810-0.076301-0.2977xxaaxx第38頁(yè)/共93頁(yè)第39頁(yè)/共93頁(yè)第40頁(yè)/共93頁(yè)第41頁(yè)/共93頁(yè)第42頁(yè)/共93頁(yè)第43頁(yè)/共93頁(yè)第44頁(yè)/共93頁(yè)第45頁(yè)/共93頁(yè)123241467A第46頁(yè)/共93頁(yè)第47頁(yè)/共93頁(yè)格式： Q,R = qr(A)第48頁(yè)/共93頁(yè)第49頁(yè)/共93頁(yè)第50頁(yè)/共93頁(yè)格式： D=chol(A)第51頁(yè)/共93頁(yè)例：進(jìn)行Cholesky分解。 A=16 4 8; 4 5 -4; 8 -4 22; D=chol(A)D = 4 1 2 0 2 -3 0 0 3164

4、84548422A第52頁(yè)/共93頁(yè)123123124x2xx23xx2x1011x3x8第53頁(yè)/共93頁(yè)第54頁(yè)/共93頁(yè)第55頁(yè)/共93頁(yè)第56頁(yè)/共93頁(yè)例：對(duì)此矩陣做三種變換。122111221A第57頁(yè)/共93頁(yè)122111221A第58頁(yè)/共93頁(yè)f=D-1b.第59頁(yè)/共93頁(yè)n第60頁(yè)/共93頁(yè)121232310910272106xxxxxxx第61頁(yè)/共93頁(yè)素）.第62頁(yè)/共93頁(yè)n第63頁(yè)/共93頁(yè)G-S法迭代求解,看是否收斂。123122911172216xxx第64頁(yè)/共93頁(yè)第65頁(yè)/共93頁(yè)其中：w最佳值在1, 2）之間，不易計(jì)算得到，因此 w通常有經(jīng)驗(yàn)給出。第

5、66頁(yè)/共93頁(yè)endn第67頁(yè)/共93頁(yè) 0.7916第68頁(yè)/共93頁(yè) x(k+1)= (D-U)-1 Lx(k+1/2) + (D-U)-1 b第69頁(yè)/共93頁(yè)n第70頁(yè)/共93頁(yè)123410120611113252110111031815xxxx第71頁(yè)/共93頁(yè)第72頁(yè)/共93頁(yè)第73頁(yè)/共93頁(yè)第74頁(yè)/共93頁(yè)第75頁(yè)/共93頁(yè)第76頁(yè)/共93頁(yè)第77頁(yè)/共93頁(yè)第78頁(yè)/共93頁(yè)第79頁(yè)/共93頁(yè)4 4 2.005 5 0第80頁(yè)/共93頁(yè)第81頁(yè)/共93頁(yè) 稀疏矩陣的計(jì)算： 同滿矩陣比較，稀疏矩陣在算法上有很大的不同。具體表現(xiàn)在存儲(chǔ)空間減少，計(jì)算時(shí)間減少。例：比較求解下面方程組n1000時(shí)兩種方法的差別。124111411141nxxx第82頁(yè)/共93頁(yè)第83頁(yè)/共93頁(yè)