Cấu trỳc:
>> cond(A) % A là ma trận kết quả trả lại dạng nhƣ sau: a* 10k
; 0 < a < 9
k là số digits khụng tin cậy trong kết quả giải hệ phƣơng trỡnh tuyến tớnh và trong việc nghịch đảo ma trận. Nếu k xấp xỉ 1 thỡ đúlà ma trận cú well -condition
Vớ dụ >>A=[1/2 1/3 1/4 ; 1/3 1/4 1/5; 1/4 1/5 1/6]; >> cond(A)
ans=
1.3533e+003
Ta thấy rằng k= 3 tức là cú 3 số khụng đỏng tin cậy Tổng kết
Định nghĩa :Hạng ma trận
Ar là một ma trận r hàng r cột đƣợc xõy dựng từ A , khụng nhất thiết theo thứ tự trong ma trận A và det(Ar)#0 .Nếu bất kỳ ma trận Ar+1 nào đƣợc xõy dựng từ r+1 hàng và r+1 cột của A, det(Ar+1)=0 thỡ chỳng ta núi rằng Matrận A cú hạng bằng r
Một hệ thống m phƣơng trỡnh tuyến tớnh trong n biến (chƣa biết) a11*x1 + a12*x2+ . . . +a1n*xn=b1
a21*x1 + a22*x2+ . . . +a2n*xn=b2 . .
am1*x1 + am2*x2+ . . . +amn*xn=bm Cú thể viết dƣới dạng form ma trận AX=B
Trong đú A là ma trận hệ số và X là vector kết quả
Điều kiện cú nghiệm Matrận [A B] đƣợc gọi là ma trận mở rộng của hệ. Theo Kronecker- Capelli thỡ hệ phƣơng trỡnh tuyến tớnh cú nghiệm khi và chỉ khi hạng của ma trận A bằng hạng của ma trận bổ xung
Nếu r= n thỡ nghiệm trờn là duy nhất
Nếu r<n thỡ hệ trờn khụng xỏc định và r biến cú thể đƣợc biểu diễn dƣới dạng hàm của n-r biến khỏc ,cỏc biến khỏc này cú thể cho giỏ trị bất kỳ( núi cỏch khỏc hệ vụ số nghiệm)
Nghiệm của hệ phƣơng trỡnh tuyến tớnh đƣợc tớnh trong Matlab bằng toỏn tử ( \ ) .Nếu hệ cú nghiệm duy nhất Matlab sẽ cung cấp cho nú , nếu hệ là khụng xỏc định(r<n) thỡ toỏn tử ( \ ) sẽ cung cấp cho chỳng ta một nghiệm riờng trong đú n-r biến sẽ đƣợc đặt =0.
Một nghiệm , nghiệm này làm thoả món tổng bỡnh phƣơng của cỏc nghiệm bộ nhất Dựng lệnh X= pinv(A)*B
50
Nếu hạng của A # hạng matrận mở rộng thỡ toỏn tử ( \ ) cung cấp một kết quả nhƣng kết quả này khụng phải là nghiệm của hệ
Hệ thuần nhất khi vector B=0. Một hệ thuần nhất cú một nghiệm tầm thƣờng khi det(A) # 0.Nếu det(A)=0 hệ cú nhiều hơn một nghiệm trong trƣờng hợp này Matlab sẽ cảnh bỏo ngƣời dựng :
TểM TẮT NỘI DUNG CỐT LếI
- Áp dụng matlab trong tớnh toỏn mạch điện, giải mạch điện nhanh chúng tiện lợi- ỏp dụng giải phƣơng trỡnhvà hệ phƣơng trỡnh tuyến tớnh.
- Giải cỏc phƣơng trỡnh phi tuyến và PT tham số. - Giải hệ PT vi phõn.
BÀI TẬP ỨNG DỤNG, LIấN HỆ THỰC TẾ Cõu hỏi ụn tập
1. Cỏc cỏch nhập một ma trận ?
2. Điều kiện cú nghiệm của hệ phƣơng trỡnh đại số tuyến tớnh , cỏch tớnh
3. Lập chƣơng trỡnh mà đầu vào là hai ma trận A và b, đầu ra là kết quả thụng bào hệ cú nghiệm hay khụng.
4. Hệ phƣơng trỡnh điều kiện yếu là gỉ, những ảnh hƣởng của nú.
Bài tập
Thực hiện cỏc phộp toỏn sau
1.Nhập hai ma trận a=[ 1 2 3; 4 5 6], b=[5 6 7 ;8 9 10] . -Tớnh Ađab= a+b .
-Tớnh Subsab= a-b; -Tớnh Multab= a*b
-Tớnh Mulba=b*a;
-Tớnh Divab= a/b; và b/a 2. Cho mạch điện sau
Cho thụng số: R1= 10(omh) , R2= 20(omh) , R3= 10(omh)
R1 R2
R3
51
e1= 20(v) , e2= 30(v)
Tớnh dũng điện I1 và I2 và I3 bằng cỏch lập theo dạng hệ phương trỡnh đại số tuyến tớnh ba ẩn số
Bài tập giải hệ phƣơng trỡnh tuyến tớnh sau: A*X= B
Trong đú: A=[1/2 1/3 1/4 ; 1/3 1/4 1/5; 1/4 1/5 1/6]; B=[0.95 0.67 0.52] 1. Giải hệ đó cho
2. Thay đổi B(3)=0.53 rồi giải hệ phƣơng trỡnh, so sỏnh với trƣờng hợp trờn 3. Tớnh điều kiện của ma trận này và đƣa ra nhận xột
Chỳ ý khi giải hệ phƣơng trỡnh tuyến tớnh với ma trận hệ số là ma trận Hilbert (ma trận điều kiện yếu) thỡ ta dựng hàm tớnh nghịch đảo của nú là hàm invhilb(n)
Định nghĩa Ma trận Hilbert là: A=[ 1 1/2 1/3...1/n;1/2 1/3 ...1/(n+1);
HƢỚNG DẪN TỰ Ở NHÀ
52
CHƢƠNG 4 :ĐỒ HOẠ TRONG MATLAB
MỤC TIấU CỦA CHƢƠNG
- Hiểu rừ khỏi niệm về đồ họa trong Matlab
- Vận dụng đƣợc cỏc lệnh vào việc giải cỏc bài toỏn và vẽ đồ thị, bảng biểụ.. - Về thỏi độ: Học sinh, Sinh : nắm đƣợc cỏc cõu lờnh và vận dụng giải cỏc bài toỏn
NỘI DUNG BÀI GIẢNG Lí THUYẾT 4.1 Màn hỡnh đồ thị
Đõy là nơi trỡnh bày mọi hỡnh ảnh , độ thị …. đó đƣợc giải trỡnh từ khung của
sổ Command của Matlab để xử lý theo cỏc lệnh , cụng cụ màn hỡnh. Cú hai cỏch để hiện khung màn hỡnh đồ thị trắng:
Từ khung cửa sổ command kớch File/New và chọn Figure từ menu xổ.
Cũng trong khung cửa sổ Command gừ lệnh figure và ấn enter.
4.2.Cỏc lệnh trờn menu đồ họa:
Để giỳp cỏc bạn cú thể nắm vững và sử dụng chƣơng trỡnh Matlab; phần này giới thiệu cỏclệnh trờn thanh menu cựng với cỏc chức năng và cụng dụng của từng menu con nằm trong cỏc menu chớnh.
4.2.1 File:
Hiện menu xổ chứa cỏc lệnh con cú chức năng tạo, quản lý, điều hành cũng nhƣ
thay đổi cỏc thụng tin cỏc thụng số mặc định của chƣơng trỡnh cho phự hợp với từng cụng việc.
New Figure
New Figure dựng để mở trang màn hỡnh đồ hoạ mớị
Để mở trang màn hỡnh mới chồng lờn màn hỡnh đồ thị cũ trong khi đồ thị cũ vẫn cũn hiện diện trờn màn hỡnh, chọn New Figure từ menu xổ. Một cửa sổ mới sẽ xuất hiện ra nằm chồng lờn màn hỡnh cũ
Open:
53
Cỏc bƣớc thực hiện mở đồ thị đó lƣu:
Kích File/ Open từ menu xổ, xuất hiện màn hình thoại Open .
- Look in: Nơi chứa cỏc tập tin đồ thị của Matlab. Nơi chứa cú thể là ổ đĩa, thƣ mục
hoặc chƣơng trỡnh khỏc.
- File nane: Tờn tập tin muốn mở trong khung màn hỡnh đồ thị .
- Files of type: Thể loại tập tin đồ thị là .fig
Kích đúp vào tên tập tin muốn mở hoặc đánh tên tập tin vào khung File name hoặc kích một lần vào tên tập tin, kích Open. Tập tin đồ thị vừa chọn sẽ hiện lên màn hình.
Close:
Đúng khung màn hỡnh đồ thị để về khung cửa sổ nhập lệnh của Matlab ( Biểu
tƣợng cú chức năng tƣơng đƣơng với lệnh Close trong menu File)
Save:
Lƣu lại những thay đổi trong khung màn hỡnh đồ thị hiện hành. Tuy nhiờn, cú
một điều khỏc biệt là lệnh này lƣu lại ngay những thay đổi trong tập tin mới sau khi đó đƣợc đặt tờn và đang hiện diện trờn màn hỡnh để tiếp tục xử lý.
Nếu bạn mở tập tin cũ với lệnh Open để xử lý và nếu đó cú những thay đổi bất kỳ trong nội dung hiện hành và sau khi kớch lệnh save, màn hỡnh hiện khung thoại save as. Từ khung thoại này bạn cú thể lƣu lại nhữnh thay đổi theo tờn tập tin cũ hoặc với một tờn mớị
54
Save As: Hiện khung thoại Save As để bạn lƣu tập tin đồ thị mới vẽ theo một tập tin
mới hoặc lƣu lại những thay đổi trong nội dung của tập tin cũ đƣợc mở với lệnh Open theo tờn cũ hoặc với tờn mớị Cỏc bƣớc thực hiện nhƣ sau:
Sau khi thay đổi , kớch File / Save as
- Save in Nơi chƣa cỏc tập tin muốn. Lƣụ Nơi chứa cỏc tập tin đồ thị của Matlab. Nơi chứa cú thể là ổ đĩa, thƣ mục hoặc chƣơng trỡnh khỏc.
- File name Tờn tập tin tuỳ chọn
để lƣu cho đồ thị vừa tạo
- Save as type Thể loại tập tin muốn lƣụ Mặc định là .fig đối với cỏc tập tin đồ thị
Sau khi chọn nguồn chứa (nếu cần thiết), đặt tờn mới cho đồ thị, kớch vào
Save để lƣụ Export:
Lƣu lại tập tin đồ thị hiện hành thành một dạng tập tin khỏc để sau này cú thể chuyển sang chƣơng trỡnh ứng dụng khỏc.
Tạo một đồ thị mới hoặc mở tập tin đồ thị cũ lờn màn hỡnh.
Kớch menu File và chọn Export. Màn hỡnh hiện khung thoại Export.
Save in Nơi chƣa cỏc tập tin muốn chuyển.
File name Tờn tập tin muốn lƣu lại để chuyển. Bạn cú thể đặt tờn theo tờn cũ
nhƣng phần mở rộng lại là một tờn khỏc .
Save as type Chọn loại tập tin muốn lƣu lại để chuyển.
Sau khi chọn song, kớch vào Save để ghi lại tập tin theo dạng khỏc.
55
Hiện khung thoại Graphics Property để ngƣời sử dụng thay đổi cỏc khung thuộc
tớnh mặc định cho phự hợp với tỏc vụ.
Preferences: Hiện khung thoại với ba tuỳ chọn để ngƣời sử dụng cú thể thay đổi tham số cho phự hợp nhiệm vụ.
General: Hiện khung thoại ngay khi
kớch chọn lệnh Preferences từ menu File của khung cửa sổ lệnh MATLAB. Kớch chọn cỏc loại tham số muốn thay đổi hoặc gỏn thờm sau đú kớch OK.
Command Windows Font: Hiện khung danh mục font cựng thuộc tớnh để
ngƣời sử dụng thay đổi font mặc định thành font quờn thuộc.
Copy Options : Hiện khung thụng tin để ngƣời sử dụng tuỳ chọn
những thành phần đối tƣợng muốn sao chộp.
Page Setup: Hiện khung thoại
56
cho phự hợp với đồ thị khi vẽ hoặc in.
Orientation: Chọn hƣớng giấy in dọc hoặc in ngang.
Limits: Đỏnh dấu chọn để ỏp dụng hoặc khụng ỏp dụng chế độ vừa khớp với
cỏc trục và gia trị theo thanh đo trờn trục.
Color: Chọn chế độ mầu hoặc đen trắng khi in hoặc vẽ trờn thiết bị.
Size and Position: Cú hai tham số tuỳ chọn, sử dụng theo mặc định hoặc
chỉnh theo tuỳ chọn.
Center In đồ thị ở giữa trang
Fill In đồ thị trọn cả trang
MaxAspect In đồ thị theo tỷ lệ tối đa
Print Setup:
Nếu cú nhiều loại thiết bị in( vẽ ), thỡ nờn
khai bỏo lại trƣớc khi chuyển bản đồ thị sang mỏy in hoặc mỏy vẽ.
Cỏc bƣớc thực hiện vẽ một đồ thị: Kớch File / Open.
Kớch File / Print Setup, màn hỡnh hiện khung thoại Print Setup.
Kớch chọn loại mỏy in, khổ giấy…. Theo ý muốn.
Kớch Ok để in.
Print Preview: Xen lại đồ thị hiện hành trƣớc
khi chuyển sang in
Print:
Hiện khung thoại print để ngƣời sử dụng chọn loại mỏy in hoặc mỏy vẽ cho phự hợp; chọn số lƣợng bạn in với Copies. Nếu cần thiết kớch Properties để chọn thờm cỏc thụng tin bổ sung cần thiết.
4.2.2 Edit: Hiện menu xổ chứa cỏc lệnh con liờn quan đến việc sử lý cỏc đối tƣợng
57
Undo: Huỷ lệnh hoặc chức năng đó đƣợc thi hành trƣớc đú. Sau khi kớch lệnh Undo
để thực hiện việc huỷlệnh, Undo sẽ chuyển đổi thành Redo để cú thể phục hồi những chức năng đó huỷ.
Cut ( Ctrl + x): Cắt đối tƣợng đó chọn trong bản đồ thị hiện hành đƣa vào Cliboard
để sau đú dỏn vào đồ thị khỏc. Sau khi cắt đối tƣợng đó chọn sẽ biến mất tại vị trớ hiện hành. Cỏc bƣớc thực hiện nhƣ sau:
Kớch vào biểu tƣợng Enable Plote Editing trờn thanh cụng
cụ.
Kớch chọn đối tƣợng muốn cắt.
Kớch Edit / Cut. Đối tƣợng sẽ biến mất tại vị trớ hiện hành và lƣu vào Cliboard.
Copy ( Ctrl + c): Sao chộp đối tƣợng đó chọn trong đồ thị hiện hành đƣa vào
Cliboard để sau đú dỏn vào đồ thị khỏc. Đối tƣợng sau khi sao chộp vẫn hiện diện tại vị trớ cũ. Cỏc bƣớc thực hiện nhƣ sau:
Kớch vào biểu tƣợng Enable Plote Editing trờn thanh cụng
cụ.
Kớch chọn đối tƣợng muốn cắt.
Kớch Edit / Copỵ Đối tƣợng sẽ vẫn hiện tại vị trớ hiện hành và lƣu vào Cliboard.
Paste ( Ctrl + v): Dỏn đối tƣợng đó đƣợc Cut hoặc Copy vào Cliboard sang đồ thị
khỏc. Cỏc bƣớc thực hiện nhƣ sau:
Mở đồ thị muốn đƣợc dỏn đối tƣợng đó đƣợc lƣu vào Cliboard.
Ấn Ctrl + v hoặc kớch vào Pastẹ Đối tƣợng sẽ đƣợc dỏn vào đồ thị mới hoặc vị trớ mới của đồ thị cũ.
Để di chuyển đối tƣợng đến vị trớ mời ta kớch vào đối tƣợng muốn di chuyển. Khi con trỏ đổi thành mũi tờn bốn đầu thỡ kộo, giữ đối tƣợng đến vị trớ mới rồi thả nỳt chuột.
Clear: Xoỏ đối tƣợng đƣợc chọn trong đồ thị.
Kớch vào biểu tƣợng Enable Plot Editing trờn thanh cụng
cụ.
Kớch chọn đối tƣợng muốn xoỏ.
Kớch Edit / Clear. Đối tƣợng sẽ xoỏ mất tại vị trớ hiện hành.
Select All( Ctrl + a): Chọn tất cả cỏc đối tƣợng trong đồ thị hiện hành để xử lý theo nhu cầụ
58
Copy Figure: Sao chộp cả đồ thị sang Cliboard để rồi sau đú dỏn sang trang màn hỡnh đồ thị hoặc sang chƣơng trỡnh khỏc. Cỏc bƣớc thực hiện nhƣ sau:
Mở đồ thị muốn đƣợc sao chộp sang đối tƣợng khỏc. Kớch Edit / Copy Figurẹ Đồ thị sẽ đƣợc ghi vào Cliboard. Thoỏt khỏi Matlab, mở chƣơng trỡnh cần Copy đồ thị đến.
Kớch chọn vị trớ muốn dỏn đồ thị đến, sau đú kớch Pastẹ Đồ thị sẽ đƣợc dỏn ngay vị trớ đó chọn.
Copy Option:
Hiện khung thoại Preferences để ngƣời sử dụng thay đổi (Xem File /
Preferences).
4.2.3 Tools ( Ctrl + t) :
Chứa cỏc lệnh cụng cụ hỗ trợ việc thực hiện, xử lý cỏc thuộc tớnh minh hoạ trong đồ thị để ngƣời quan sỏt dễ hiểu hơn.
Show Toolbar:
Cho hiện hoặc giấu cỏc biểu tƣợng cụng cụ Standard trờn màn hỡnh
Enable Plot Editing: Đỏnh dấu chọn hoặc bỏ chọn đối tƣợng.
Axes Properties: Gỏn thuộc tớnh
cho cỏc trục toạ độ cựng với tiờu đề đồ thị hiện hành. Thớ dụ với đồ thị lực kộo sợi dõy khi quay trũn quả cầụ Thay vỡ dựng cỏc thuộc tớnh lệnh để gỏn chung với số liệu, ta cú thể dựng khung thoại Axes Properties để gỏn thuộc tớnh. Cỏc bƣớc nhƣ sau:
Sau khi vẽ đồ thị , kớch chọn Tool / Axes Properties. Màn hỡnh xuất hiện khung thoại Axes Properties.
59
Trong khung thoại Axes Properties cú thể: đặt tiờu đề cho đồ thị ( Title), gỏn nhón cho trục x và trục y ( Label), giới hạn độ lớn cỏc trục ( Limits), tỷ lệ xớch cỏc trục….. . cho đồ thị. Kớch Apply và kớch OK.
Line Properties: Gỏn thuộc tớnh nhƣ cỏc điểm đỉnh dữ liệu, tạo nột đồ thị dầy hay
mỏng cũng nhƣ mầu sắc cho đƣờng đồ thị hiện hành. Cỏc bƣớc thực hiện nhƣ sau:
Kớch Enable Plot Editing và kớch chọn đƣờng đồ thị. Những điểm đỉnh dữ liệu trờn đƣờng đồ thị sẽ đƣợc đỏnh dấu bằng những hỡnh vuụng nhỏ.
Kớch Tool / Line Properties. Màn hỡnh
hiện khung thoại line properties.
Chọn độ dầy mỏng của đƣờng đồ thị
Line Width; kớch cỡ dấu đỉnh dữ liệu
trờn đồ thị Marker Size; Kiểu dỏng đồ
thị Line Style; mầu sắc Color …… Kớch Apply và kớch OK.
Text Properties: Thay đổi cỏc thuộc tớnh cho cỏc ký tự trong đồ thị. Cỏc bƣớc thực hiện nhƣ sau:
60
Kớch menu Tool / Axes Properties. Xoỏ cỏc ký tự cần thay đổi trƣớc đú. Nhập lại cỏc ký tự theo cỏch gừ tiếng việt, đừng quan tõm đến cỏc ký tự hiện ra cú đỳng khụng. Kớch Apply và kớch OK.
Kớch biểu tƣợng Enable Plot Editing, chọn chuỗi ký tự muốn thay đổi thuộc tớnh.
Kớch menu Tool / Text
Properties. Màn hỡnh xuất hiện
khung thoại Edit Font
Properties.
Trong khung thoại, chọn loại font chữ Vn Time; kớch chọn cỏc thuộc tớnh cho font chữ
( kớch thƣớc, kiểu dỏng). Kớch Ok.
Unlock Axes Position: Chỉnh lại tỷ lệ trục toạ độ hoặc di chuyển đồ thị đến vị trớ
khỏc trong khung cửa sổ đổ hoạ.
Kớch chọn biểu tƣợng Enable Plot Editing. Chọn trục toạ độ x (y) nếu là đồ thị
2D hoặc trục z nếu đồ thị là 3D. Kớch menu Tool / Unloock Position.
Muốn chỉnh tỷ lệ cho trục nào, kớch chuột vào điểm chỉnh trong bốn gúc của đồ thị và kộo vào trong hay ra ngoài theo hƣớng tõm. Nếu muốn di chuyển đồ thị đến vị trớ
khỏc, kớch chuột vào điểm chỉnh nằm trờn một cạnh và kộo đến vị trớ khỏc. Sau khi thực hiện xong, kớch lại menu Tool. Lần này lệnh Unlock…… đổi thành Lock Position. Kớch Lock Position để khoỏ lạị
61
Show Legend: Cho hiện hoặc giấu bảng phõn loại số liệu của đồ thị trờn màn hỡnh
Mở đồ thị muốn hiện bảng phõn loại dữ liệu trờn màn hỡnh. Sau khi đồ thị hiện trong khung màn hỡnh đồ hoạ, kớch Tool /