144
Một ví dụ phức tạp hơn
2
2
9
4
2
s
s
F
s
s
−
−
+
=
+ +
>> F = (–s^2 – 9*s + 4)/(s^2 + s + 2);
>> ilaplace(F)
ans = –dirac(t) –8*exp(–1/2*t)*cos(1/2*7^(1/2)*t)+20/7*7^(1/2)*
exp(–1/2*t)*sin(1/2*7^(1/2)*t)
Hay một hàm đơn giản hơn
2
1
(
7)
F
s
=
+
>> f = ilaplace(1/(s+7)^2)
f = t*exp(–7*t)
Để vẽ đồ thị các hàm gốc theo thời gian ta có thể sử dụng lệnh ezplot theo cú pháp
ezplot(f, [t_start, t_end])
Với hàm f trên ta gõ lệnh
>> ezplot(f,[-2,2])
>> ezplot(f,[0, 2])
và nhận được đồ thị như trên hình 7-1 và 7-2.
-2
- 1
0
1
2
- 10
- 5
0
x 10
4
t
t exp(- 7 t)
0
0.5
1
1.5
2
0
0.02
0.04
t
t exp(- 7 t)
Hì
nh 7-1. Hàm
7
, 2
2
t
te
t
−
− ≤
≤
Hì
nh 7-2. Hàm
7
, 0
2
t
te
t
−
≤
≤
Một ví dụ khác là tìm hàm gốc và vẽ đồ thị của nó, khi biết hàm ảnh
2
2
2
3
(
1) (
3)
s
F
s
s
+
=
+
+
Gọi hàm ilaplace và nhận được hàm gốc như sau:
>> f = ilaplace((2*s+3)/((s+1)^2*(s+3)^2))
