26
pa = pn + [zeros(1,length(pn)-length(pm)) pm]
ps = pn - [zeros(1,length(pn)-length(pm)) pm]
Cụ thể
>> pa=p5+[zeros(1,length(p5)-length(p2)) p2]
pa = 6 1 -12 55 -50 45
>> ps=p5-[zeros(1,length(p5)-length(p2)) p2]
ps = 6 1 -12 49 -42 33
Không điểm hay nghiệm của phương trình đa thức
Hàm
roots(p)
tính cho ta các nghiệm của phương trình đa thức
( )
0
n
p s =
. Các nghiệm tìm được
có thể là thực hoặc phức. Ví dụ khi tìm nghiệm của các phương trình
5
4
3
2
5
( )
6
1
12
52
46
39
0
p s
s
s
s
s
s
=
+
−
+
−
+
=
2
2
( )
3
4
6
0
p s
s
s
=
−
+
=
ta thực hiện và thu được
>> p2 = [3 -4 6];
>> p5 = [6 1 -12 52 -46 39];
>> roots(p2)
ans =
0.6667 + 1.2472i
0.6667 - 1.2472i
>> roots(p5)
ans =
-2.6672
0.7770 + 1.1589i
0.7770 - 1.1589i
0.4733 + 1.0138i
0.4733 - 1.0138i
Xây dựng đa thức từ các không điểm cho trước
Hàm
p = poly(nu)
đưa ra cho ta một đa thức nhận các giá trị trong véctơ nu làm nghiệm.
Ví dụ cần đưa ra một đa thức nhận các giá trị sau là nghiệm
nu = [1-2i 1+2i 3 4]
>> nu = [1-2i 1+2i 3 4] ;
>> p4=poly(nu)
p4 = 1 -9 31 -59 60
Đa thức cần tìm sẽ là
4
3
2
4
( )
9
31
59
60
p s
s
s
s
s
=
−
+
−
+
Ta có thể kiểm tra lại bằng lệnh roots(p)
>> roots(p)
ans =
4.0000
