CHƯƠNG2:ĐẠISỐTUYẾNTÍNH  §1.CÁCPHƯƠNGTRÌNHĐẠISỐTUYẾNTÍNH 1.Hệphươngtrìnhđầyđủ:TaxéthệphươngtrìnhAx=B.Đểtìmnghiệmcủa hệtadùnglệnhMATLAB:  x=inv(A)*B hay:  x=A\B  2.Hệphươngtrìnhcóít phươngtrìnhhơnsốẩn(underdetermined):Khigiải hệtrêntađãdùngnghịchđảomatrận.Nhưvậy tachỉnhậnđượckếtquảkhi matrậnAvuông(sốphương trìnhbằngsốẩnsốvàđịnhthứcc ủaAphảikhác không).HệcósốphươngtrìnhíthơnsốẩnhayđịnhthứccủamatrậnAcủa hệđầyđủbằng0gọilàhệunderdetermined.Mộthệnhưvậycóthểcóvôs ố nghiệmvớimộthaynhiềubiếnphụthuộcvàocácbiếncònlại.Vớimộthệnhư vậyphươngphápCramerhayphươngphápmatrậnnghịchđảokhôngdùng được.Khisốphương trìnhnhiềuhơnsốẩnphươngphápchiatráicũngcho nghiệm với một vàiẩn số đượcchobằng0.Mộtvídụđơngiản là phương trìnhx+3y=6.Phươngtrìnhnàycór ấtnhiềunghiệm trongđócómộtnghiệm làx=6vày=0:  a=[13];  b=6;  x=a\b  x= 6 0 Sốnghiệmvôhạncóthểtồn tại ngaycảkhisốphươngtrìnhbằngsốẩn.Điều này xảy ra khi det(A) = 0. Với hệ này ta không dùngđược phương pháp Cramervàphươngpháp matrậnnghịchđảo vàphươngphápchiatráicho thông báolàma trậnAsuybiến.Trongtrườnghợpnhưvậytacóthểdùng phươngphápgiảnghịchđảođểtìmđượcmộtnghiệmgọilànghiệmchuẩn minimum. Vídụ:Chohệphươngtrình   x+2y+z=8  0x+y+0z=2  x+y+z=6 29 Khidùngphépchiatráitanhậnđược:  y=a\b  Warning:Matrixissingulartoworkingprecision. y=   Inf   Inf   Inf Nếutadùngphươngphápgiảnghịchđảothìcó:  a=[121;010;111] b=[8;2;6]  x=pinv(a)*b  x=   2.00000000000000   2.00000000000000   2.00000000000000  Mộthệcũngcóthểcóvôsốnghiệmkhicóđủsốphươngtrình.Vídụta cóhệ:  2x‐4y+5z=‐4 ‐4x‐2y+3z=4  2x+6y‐8z =0 Tronghệnàyphươngtrìnhthứ3làtổngcủahaiphương trìnhtrênnênhệthật sựchỉcó2phươngtrình.  Tómlạimộthệmuốncónghiệmduynhấtphảicócácph ươngtrìnhđộc lập.Việcxácđịnhcácphươngtrìnhtronghệcóđộclậphaykhôngkhákhó, nhấtlàđối vớihệcónhiềuphươngtrình.Tađưaramộtphươngphápcho phépxácđịnhhệphươ ngtrìnhcónghiệmvàliệunghiệmđócóduynhấthay không.Phươngphápnàyđòihỏisựhiểubiếtvềhạngcủamatrận.  Taxemxétđịnhthứccủamatrậnsau:   ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎣ ⎡ − − 379 2106 143 Nếutaloạitrừmộthàngvàmộtcộtcủamatrậnchúngta cònlạimatrận2 ×2. Tuỳtheohàngvàcộtbịloạitacó9matrậncon.Địnhthứccủacácmatrậnnày gọilàđịnhthứccon.Vídụnếutabỏhàng1vàcột1tacó:  44 37 210 = −  30 Cácđịnhthứcconcóthểdùngđểxácđịnhhạngcủamatrận.Hạngcủama trậnđượcđịnhnghĩanhưsau: MộtmatrậnAm × ncóhạngr  ≥ 1nếuvàchỉnếu địnhthứccủaAchứamộtđịnh thứcr × rvàmọiđịnhthứcconvuôngcór+1hàng hayhơnbằng0.  Đểxácđịnhhạngcủamatrậntacólệnhrank Vídụ:   a=[3‐41;6102;9‐73];  rank(a)  ans= 2  HệphươngtrìnhAx=Bcómphươngtrìnhvànẩncónghiệmnếuvàchỉnếu rank(A)=rank([AB]).GọihạngcủaAlàr,nếur=nthìnghiệmlàduynhất. Nế ur<nthìhệcóvôsốnghiệmvàrẩncóth ểbiểudiễnnhưlàtổhợptuyến tínhcủan‐rẩncònlạimàgiátrịcóthểchọnbấtkì. Vídụ:Giảihệphươngtrình  3x‐2y+8z=48 ‐6x+5y+z=‐12  9x+4y+2z=24 Taviết:  a=[3‐28;‐651;942];  b=[48;‐12;24];  rank(a)  ans= 3  rank([ab])  ans= 3 Vậyhệcónghiệmduynhất:  x=a\b  x= 2 ‐1 5  Vídụ:Giảihệ  2x‐4y+5z=‐4 ‐6x‐2y+3z=4  2x+6y‐8z=0 31 Taviết:  a=[2‐45;‐6‐23;26‐8];  b=[‐4;4;0];  rank(a)  ans= 2  rank([ab])  ans= 2  Vậyhệcóvôsốnghiệm.Mộttrongcácnghiệmlà: x=pinv(a)*b x= ‐1.21481481481481 0.20740740740741 ‐0.14814814814815  3.Hệphươngtrìnhoverdetermined:Hệphương trình trongđó số phương trìnhđộc lập nhiều hơn số ẩn gọi là hệ overdetermined.Đối với hệ này phươngphápCramervàphươngphápnghịchđảomatrậnkhôngdùngđược. Tuynhiênmột sốhệchonghiệmđúngxácđịnhbằngphépchiatrái.Đốivới các hệ khác không có nghiệm chính xác. Khi r = rank(a)=rank([ab])hệcó nghiệmvànếur=nnghiệmlàduynhất.Khirank(a) ≠rank([ab])hệkhôngcó nghiệm. Vídụ:Giảimạchđiệngồm3nhánhnốisongsong:nhánh1cótổngtrởZ1= 5+2jvànguồne=100sin(314t+30 0 ),nhánh2cótổngtrởZ2=3+4jvànhánh3 cótổngtrở5+6j.Taviếtphương trìnhcủamạchđiệntheodòngnhánh.Sauđó rútramatrậnAvàB.CáclệnhMATLAB: a=[111;5+2*i3+4*i0;0‐(3+4*i)5+6*i] e=100*exp(i*(30*pi/180)) b=[0;e;0]; i=a\b i= 25.25569272231586+19.27124163998603i ‐15.63482777750950‐11.44276084484129i ‐9.62086494480636‐7.82848079514474i    32 §2.NỘISUY 1.Nộisuyhàmmộtbiến :MATLABdùng2kiểunộisuy:nộisuyđathứcvà nộisuytrêncơsởbiếnđổiFourriernhanh,  a.Nộisuyđathức:MATLABdùnghàminterp1(x,y,xi,<phươngpháp>) vớix,làgiátrịcủahàmtạinhữngđiểmđãchovàxilàgiátrịmàtạiđótacần nộisuyragiátrịyi.<phươngpháp>cóthểlà mộttrongcácgiátrịsau:  ‘nearest’‐ phương pháp nàyđặt giá trị nội suy vàogiá trị đã cho gần nhất,Phươngphápnàynhanhnhưngkếtquảkémchínhxácnhất Vídụ: x=[12345]; y=[5.543.1128290.7498.4]; yi=interp1(x,y,1.6,ʹnearestʹ) yi= 43.1000  ‘linear’‐phươngphápnày coiđường congđi qua 2điểmchotrướclà đườngthẳng. Vídụ: yi=interp1(x,y,1.6,ʹlinearʹ) yi= 28.0600   ‘spline”‐dùngphươngphápnộisuyspline Vídụ: yi=interp1(x,y,1.6,ʹsplineʹ) yi= 24.9782   ‘cubic’‐phươngphápnàycoiđườngcongqua2điểmlàđườngcongbậc 3 Vídụ: yi=interp1(x,y,1.6,ʹcubicʹ) yi= 22.3840  b.NộisuyFTT:Hàminterpftth ựchiệnnộisuyhàmmộtbiếnsửdụng phươngphápFFT(FastFourrierTransform).Phươngphápnàytínhtoánbiến đổiFourriermộtvectơchứacácgiátrịcủamộthàmchukì. Nhưvậyphương phápnàytínhbiếnđổiFourrierngượcsửdụngnhiềuđiểm.Dạnghàmlà:  y=interpft(x,n) Vídụ: y=interpft(x,4) y= 1.00002.62363.00005.3764  33 2.Nộisuy hàmhaibiến:Hàm interp2thựchiệnnộisuyhàm2 biến.Dạng hàmtổngquát:  ZI=interp2(X,Y,Z,XI,YI,<phươngpháp>)  Z–matrậnchữnhậtchứagiátrịcủahàm2biến  X,Y–mảngcócùng kíchthước,chứagiátrịx,yđãcho  XI,YI‐mảngchứagiátrịcầnnộisuy Các<phươngpháp>gồm:‘nearest’,’linear’,’cubic’  3.Nộisuyvàmảngnhiềuchiều: interp3 nộisuyhàm3biến interpn nộisuyhàmnhiềubiến  §3.TÍCHPHÂNVÀPHƯƠNGTRÌNHVIPHÂN 1.Tíchphân :Đểtínhtíchphântadùnghàmquad(tínhtíchphântheophương phápSimpson)vàhàmquadl(tínhtíchphânbằngphươngphápLobatto). Vídụ(lưutrongct2_1.m):  f=inline(ʹ1./((x‐0.3).^2+0.01)+1./((x‐0.9).^2+0.04)‐6ʹ);  q=quad(f,0,1)  q=   29.8583  r=quadl(f,0,1)  r=   29.8583  Vídụ(lưutrongct2_2.m)   y=sin(x)  quad(‘sin’,0,pi)  ans=   2.00001659104794   quadl(ʹsinʹ,0,pi)  ans= 1.99999999999989 Tacũngcóthểdùngphươngpháphìnhthanhđểtínhtíchphân: Vídụ(lưutrongct2_3.m):  y=sin(x)  x=[0:pi/100:pi]];  y=sin(x); 34  trapz(x,y)  ans=   1.99983550388744  2.Viphânsố :Đểtínhviphântadùngdiff Vídụ:    a=[1425748]; diff(a) ans= 3‐232‐34  3.Phươngtrìnhviphân :Phươngtrìnhviphâncấpcaoy (n) =f(t,y,y’,..,y (n‐1) ) cóthểđưavềhệphươngtrìnhviphâncấp1bằngcáchđặty 1=y;y2=y’,.., yn=y (n‐1) .Nhưvậy:   )y.,,.y,y,t(fy yy yy n21n 32 21 = ′ = ′ = ′ làhệcónphươngtrìnhviphâncấp1. Vídụ:   y’’’‐3y”‐y’y=0vớiy(0)=0y’(0)=1y”=‐1 đượcbiếnđổithành   1233 32 21 yyy3y yy yy += ′ = ′ = ′ vớiđiềukiệnđầu:y 1(0)=0y2(0)=1y3(0)=‐1 ĐểnhậpphươngtrìnhnàyvàoMATLABtadùngM‐filef.mnhưsau:  functiondy=f(t,y);  dy=[y(2);y(3);3*y(3)+y(3)*y(1)]; vàgiảiphươngtrìnhbằnglệnh:  [t,f]=solver(‘file’,tspan,y0) với “file”làM‐filechứaODE  tspanlàvectơ[t0tfinal]xacđịnhkhoảngtìmnghiệm  y0làvectơgiátrịđiềukiệnđầu.  solver là cách giải, thường dùng phương pháp Runge‐Kutta bậc 2/3(ode23)hay4/5(ode45) 35  [t,y]=ode45(‘f’,[01],[0;1;‐1]) Mỗihàngtrongvectơnghiệmtươngứngvớimộtthờiđiểmtrongvectơcộtt. Như vậy trong ví dụ trên, y(:,1) là nghiệm, y(:,1) làđạo hàm bậc nhất của nghiệmvày(:,2)là đạohàmbậchaicủanghiệm. Vídụ:TìmdòngquađộkhiđóngmạchRCnốitiếpvàonguồnmộtchiềubiết tíchsốRC=0.1,điệnápnguồnlà10Vvàđiệnápbanđầutrêntụlà2V.  Phươngtrìnhcủamạ chlà: e(t)=RC C C u dt du +  Thaysốvàotacó:  0.1u′+u=10  u′=‐10u+100 TacócáclệnhMATLABđểtạohàm:  functionuc=rc(t,u)  uc=‐10*u+100; vàgiảibàitoán:  [t,u]=ode45(ʹrcʹ,[04],2);  plot(t,u,ʹ‐oʹ)  §4.ĐẠISỐMATRẬN  1.PhântíchCholesky:PhươngphápCholeskyphântíchmatrậnAxácđịnh dương thành tích của hai ma trận  A=R’*RvớiRlàmatrậntamgiác trên. MuốnnhậnđượcmatrậnRtadùnghàmchol(A). Vídụ: n=5; X=pascal(n) X= 11111 12345 1361015 14102035 15153570 R=chol(X)  2.PhântíchLU:TaphântíchmatrậnA=L*UtrongđóLlàmatrậntamgiác dướivàUlàmatrậntamgiáctrên.Taviết[L,U]=lu(A). Vídụ(lưutrongct2_4.m): 36 A=[ 123 456 780] [L,U]=lu(A)  3.PhântíchQR:TaphântíchmatrậnA=Q*RvớiQlàmatrậntrựcgiaovàR làmatrậntamgiáctrên. Vídụ(lưutrongct2_5.m): A=[ 123 456 789 101112] [Q,R]=qr(A)  4.Luỹthừa:NếucómatrậnAvuôngvàsốp>0thìA^plàtíchplầncủaA:  Y=A^2  5.Giátrịriêngvàvectơriêng:eig(A)  [d,r]=eig(A)  6.Quaymatrận :b=rot90(a) a=[210;‐25‐1;346] a= 210 ‐25‐1 346  b=rot90(a) b= 0‐16 154 2‐23  7.Đảomatrận:fliplr(a)đảomatrậntừtráisangphải c=fliplr(a) c= 37 012 ‐15‐2 643   flipud(a)đảomatrậntừtrênxuốngdưới d=flipud(a) d= 346 ‐25‐1 210  8.Cáchàmxửlímatrậnkhác: reshape(a,m,n)địnhdạnglạimatrậnavớisốhàngmớimvàsốcộtmới n  a=[123;567;891]; reshape(a,1,9) ans= 158269371   diag(a)lấycácphầntửtrênđườngchéochínhcủamatrậnavàlưuvào mộtvectơ  diag(a,k)chọnđườngchéotuỳtheogiátrịcủak k=0‐chọnđườngchéochính k>0‐chọnđườngchéothứktrênđườngchéochính k<0‐chọnđườngchéothứkdướiđườngchéochính  a=   123   567   891  v=diag(a,1) v= 2   7  a=diag(v)n ếuvlàvectơthìalàmatrậnvuôngvớivlàđườngchéo chính  b=triu(a)tạora matrậnbcùngcỡvớimatrậna,chứacácphầntửcủa ma trận a nằm trênđường chéo chính và phía trênđường chéo chính. Các phầntửkhácbằng0.  a=[123;456;789] 38