121
1
2
,
exp( )
a
p
b
p
=
=
Ví dụ xét bộ số liệu sau
x
1.2
2.8
4.3
5.4
6.8
7.9
y
7.5
16.1
38.9
67.0
146.6
266.2
Nhập số liệu vào Matlab
>> x = [1.2 2.8 4.3 5.4 6.8 7.9];
>> y = [7.5 16.1 38.9 67.0 146.6 266.2];
>> plot(x,y,'ko', x,y, 'k-'), grid on, xlabel('x'),ylabel('y')
1
2
3
4
5
6
7
8
0
50
100
150
200
250
300
x
y
Hình 5-4. Đường nối các điểm dữ liệu
Xấp xỉ bằng hàm e mũ
>> p = polyfit(x,log(y),1)
p =
0.5366 1.3321
>> a = p(1); b=exp(p(2));
>> w = b*exp(a*x);
>> err = w-y;
>> S3= sum(err.*err)
Se = 17.6259
>> x3 = [0:0.1:8];
>> y3 = b*exp(a*x3);
>> plot(x,y,'ko', x3,y3, 'k-'), grid on,xlabel('x'),ylabel('y')
Sau đây ta thử xấp xỉ bằng đa thức bậc 3:
>> p = polyfit(x,y,3) % Xap xi bang duong bac ba
p = 1.4497 -11.2896 34.9691 -21.6503
