130 likes | 254 Vues
N EVĚŘTE POČÍTAČŮM. Radek Kučera Ostrava 7.6.2006. Jak vyřešit úlohu ?. Nabouch ám to do počítače. Počítač může umět všechno ???. Hilbertova matice. inverzní matice :. Hilbertova matice. n=7, A=hilb(n), B=inv(A), mesh(A*B) n=12. Hilbertova matice.
E N D
NEVĚŘTE POČÍTAČŮM Radek Kučera Ostrava 7.6.2006
Jak vyřešit úlohu ? • Nabouchám to do počítače. • Počítač může umět všechno ???
Hilbertova matice • inverzní matice:
Hilbertova matice n=7, A=hilb(n), B=inv(A), mesh(A*B) n=12
Hilbertova matice n=12, A=hilb(n), B=invhilb(n), mesh(A*B) n=13
Polynomy • poly výpočet koeficientů z kořenů • roots výpočet kořenů z koeficientů
Polynomy k=1:5, a=poly(k), kk=roots(a) 1 1 5.0000 2 -1 4.0000 3 85 3.0000 4 -22 2.0000 5 274 1.0000 -120
1 20.0003 2 18.9981 3 18.0065 4 16.9875 5 16.0139 6 14.9920 7 13.9996 8 13.0046 9 11.9963 10 11.0012 11 10.0001 12 8.9997 13 8.0001 14 7.0000 15 6.0000 16 5.0000 17 4.0000 18 3.0000 19 2.0000 20 1.0000 1 20.9966 2 20.0315 3 18.7863 4 18.3944 5 16.5688 + 0.5914i 6 16.5688 - 0.5914i 7 14.4519 + 0.6258i 8 14.4519 - 0.6258i 9 12.4132 + 0.2544i 10 12.4132 - 0.2544i 11 10.8915 12 10.0399 13 8.9905 14 8.0019 15 6.9997 16 6.0000 17 5.0000 18 4.0000 19 3.0000 20 2.0000 21 1.0000 k=1:21 k=1:20
Grafy x=0:0.01:1; y=sin(x); plot(x,y)
Grafy x=0:0.01:1; y=x.*sin(1./x); plot(x,y)
Grafy k=1:5; a=poly(k); x=1:0.01:5; y=polyval(a,x); plot(x,y)
Grafy k=1:21; a=poly(k); x=1:0.01:5; y=polyval(a,x); plot(x,y)
Grafy k=1:26; a=poly(k); x=1:0.01:5; y=polyval(a,x); plot(x,y)