74
(
1)
,
1, 2,...,
k
start
x
x
k
h
k
n
=
+
−
⋅
=
Kết quả của việc chia miền là một véctơ gồm
n
phần tử được gán cho biến
x
. Sau
đó ta tính giá trị của hàm số tại các điểm vừa chia, và cũng nhận được một véctơ
cùng cỡ với biến
x
. Với hàm đơn giản như
sin
y
x
=
thì ta có thể tính véctơ
y
theo cách viết thông thường. Khi đã có hai véctơ
x
và
y
cũng cỡ, ta chuyển sang
bước ba đó là sử dụng lệnh plot(x,y) ñể vẽ đồ thị. Thực hiện ba dòng lệnh sau
>> x=[0:0.1:10];
>> y=sin(x);
>> plot(x,y)
Sau khi gõ ENTER giây lát, Matlab mở ra một cửa sổ mới với tên Figure 1, trong
đó có hiện lên đồ thị hàm số
sin
y
x
=
. Kết quả của ba dòng lệnh trên được thể
hiện như trên hình 4-1.
0
2
4
6
8
10
-1
- 0.5
0
0.5
1
Hì
nh 4-1. ðồ thị hàm
sin
y
x
=
với h = 0.1
Bây giờ ta sẽ vẽ lại hàm trên với bước chia là 1, tăng 10 lần so với trường hợp trên.
Thực hiện lệnh
>> x=[0:0:10];
Và thử vẽ lại bằng lệnh, nhưng Matlab đã báo lỗi
>> plot(x,y)
??? Error using ==> plot
Vectors must be the same lengths.
Chuyện gì đã xảy ra, chúng ta đã định nghĩa hàm
sin
y
x
=
tại sao Matlab không
vẽ mà lại báo lỗi. Ở đây chúng ta đã chia lại miền vẽ, số phần tử của véctơ
x
đã
thay đổi, nhưng ta chưa tính lại các giá trị của
y
, véctơ này vẫn giữ lại các giá trị
của lần tính trước và như vậy hai véctơ x và y có số phần tử khác nhau. Đó là lý do.
Để vẽ được ta cần tính lại với véctơ y, hãy sử dụng phím lên xuống (↑, ↓) ñể chọn
lại dòng lệnh hoặc viết lại.