117
Theo cách biểu diễn đa thức bằng một véctơ bắt đầu từ hệ số ứng với số hạng có
bậc lũy thừa lớn nhất, ta có
>> a=p(1)
a = 0.0392
>> b=p(2)
b = 3.6267
Giá trị cần tìn tại x = 11 được tính bằng lệnh polyval(p,x)
>> polyval(p,11)
ans = 4.0574
Để thấy được rõ đường thẳng cùng với các số liệu đã cho, ta vẽ chúng trên cùng
một đố thị như trên hình 5-1 bằng các dòng lệnh:
>> x1=[6:0.1:24];
>> y1=a*x1+b;
>> plot(x,y,'ko', x1,y1, 'k-'), grid on,
>> xlabel('x'), ylabel('y')
6
8
10
12
14
16
18
20
22
24
3.6
3.8
4
4.2
4.4
4.6
x
y
Hình 5-1a. Đường thẳng cùng với các điểm dữ liệu
Nếu tính theo đường thẳng vừa tìm thì giá trị tại các điểm
i
x cho trong bảng sẽ là
>> w = a*x+b
w = 3.8616 3.9399 4.0182 4.0965 4.1748 4.2532 4.3315 …
4.4098 4.4881 4.5664
Khi chỉnh lại các trục tọa độ bằng lệnh axis, ta thấy các điểm dữ liệu nằm rất gần
và phân bố hai bên đường thẳng
y
ax
b
=
+
vừa tìm được.
>> axis([6 24 0 5])
