SAYISAL ÇÖZÜMLEME DERSİ ÇALIŞMA SORULARI(Harita Müh. Böl. Arasınav) S-1) Aşağıdaki nonlinear denklemlerin sıfır yerini sabit nokta(basit ) iterasyonu kullanarak 𝜀 = 10−3 tolerans için bulunuz. a)𝑥 = 𝑥 5 − 0.25 b) 𝑥 = 2𝑠𝑖𝑛𝑥 𝑥0 = 0 𝑏𝑎ş𝑙𝑎𝑛𝑔𝚤ç 𝑛𝑜𝑘𝑡𝑎𝑠𝚤 𝑖𝑙𝑒 𝑥0 = 2 𝑏𝑎ş𝑙𝑎𝑛𝑔𝚤ç 𝑛𝑜𝑘𝑡𝑎𝑠𝚤 𝑖𝑙𝑒 c) 𝑥 = √3𝑥 + 1 𝑥0 = 2 𝑏𝑎ş𝑙𝑎𝑛𝑔𝚤ç 𝑛𝑜𝑘𝑡𝑎𝑠𝚤 𝑖𝑙𝑒 2−𝑒 −𝑥 +𝑥 2 3 𝑥0 = 1 𝑏𝑎ş𝑙𝑎𝑛𝑔𝚤ç 𝑛𝑜𝑘𝑡𝑎𝑠𝚤 𝑖𝑙𝑒 d) 𝑥 = S-2) Aşağıdaki nonlinear denklemlerin sıfır yerini ikiye bölme yöntemi yardımıyla 𝜀 = 10−3 tolerans için bulunuz. a) 𝑓(𝑥) = 𝑐𝑜𝑠(𝑥) − 𝑥 𝑓𝑜𝑛𝑘𝑠𝑖𝑦𝑜𝑛𝑢𝑛𝑢 [0; 2] 𝑎𝑟𝑎𝑙𝚤ğ𝚤 b) 𝑓(𝑥) = 𝑥 3 − 𝑥 − 1 𝑓𝑜𝑛𝑘𝑠𝑖𝑦𝑜𝑛𝑢𝑛𝑢 [−2; 3] 𝑎𝑟𝑎𝑙𝚤ğ𝚤 c) 𝑓(𝑥) = 𝑙𝑛(𝑥 + 1) − 𝑥2 4 𝑓𝑜𝑛𝑘𝑠𝑖𝑦𝑜𝑛𝑢𝑛𝑢 [1; 4] 𝑎𝑟𝑎𝑙𝚤ğ𝚤 d) 𝑓(𝑥) = 5𝑒 −𝑥 − cosh(𝑥) 𝑓𝑜𝑛𝑘𝑠𝑖𝑦𝑜𝑛𝑢𝑛𝑢 [0; 4] 𝑎𝑟𝑎𝑙𝚤ğ𝚤 Derste tanımlanan bisect.m programı ve MATLAB ait fzero programları yardımıyla Soru 2 de verilen fonksiyonların sıfıryerlerini belirleyiniz. fzero aşağıdaki örnekteki gibi sıfır yeri hesaplar; >>f=inline(.log(x)-x+4.) >> fsifir(f,4) ans =0.3183 Başlangıç tahmin için x0 = a(keyfi değeri siz belirleyin) değerlerini alınız. S-3) Üçüncü dereceden bir polinomun kökleri analitik olarak hesaplanabilir. Ömer Hayyam ve Girolamo Cardano.nun önemli katkı sağladığı alanlar olarak bilinen bu tür polinomların köklerini MATLAB ortamında hesaplayan program aşağıda verilmektedir. Polinomlar P(x) = a(1)x3 + a(2)x2 + a(3)x + a(4) olarak tanımlanmaktadır. function X=kubik(a) %x^3+px^2+qx+r polinomunun sıfıryeri a=a/a(1); p=a(2);q=a(3);r=a(4); a=1/3*(3*q-p^2) b=1/27*(2*p^3-9*p*q+27*r); delta=b^2/4+a^3/27; A=(-b/2+sqrt(delta))^(1/3); B=(-b/2-sqrt(delta))^(1/3); x1=A+B; x2=-1/2*(A+B)+i* sqrt(3)/2*(A-B); x3=-1/2*(A+B)-i*sqrt(3)/2*(A-B); X=[x1 x2 x3]-p/3; end 𝑥 3 − 0𝑥 2 − 𝑥 − 1 = 0 polinomun katsayıları girilerek program çalıştırılmış aşağıda sonuçlar elde edilmiştir. >> a=[1 0 -1 -1] >> kubik(a) ans = 1.3247 -0.6624 + 0.5623i -0.6624 - 0.5623i Elde ettiğiniz sonuçlar yukarıdaki formül ile elde edilen sonuçlarla aynı olmalıdır: >> roots(a) ans = 1.3247 -0.6624 + 0.5623i -0.6624 - 0.5623i Aşağıda verilen polinomlar için yukarıda verilen programı çalıştırarak kökleri elde ediniz ve ayrıca roots(a) komutu ile elde edilen sonuçları karşılaştırınız. a) b) c) d) 𝑥 3 − 0𝑥 2 − 𝑥 − 1 = 0 2𝑥 3 − 3𝑥 2 − 6𝑥 + 1 = 0 3𝑥 3 − 5𝑥 2 + 9𝑥 − 7 = 0 −𝑥 3 − 6𝑥 2 + 8𝑥 − 1 = 0 S-4) MATLAB ortamnda Ax = b denklemsistemi ‘\’ operatörü yardımıyla x = A\b komutu ile çözmektedir. Buna göre aşağıdaki denklem sistemlerini çözünüz: 𝑥 − 2𝑥 = 3 2 a) { 1 3𝑥1 − 𝑥2 = 5 𝑥 − 2𝑥 = −1 1 2 b) { 2𝑥 − 𝑥 1 2 =2 𝑥1 − 3𝑥2 + 4𝑥3 = 1 −𝑥 c) { 1 + 2𝑥2 − 𝑥3 = 3 2𝑥1 − 𝑥2 + 5𝑥3 = 4 S-5) Aşağıda verilen iki tabanlı sayıların on tabanlı gösterimlerini belirleyiniz a) b) c) d) e) (110101)2 (10101011)2 (110101000)2 (10011001)2 (0.101101)2 S-6) Aşağıda, ondalık sistemde verilen sayıların ikili sistemdeki karşılıklarını bulunuz. a) (23532100)𝑜𝑛 b) (30002)𝑜𝑛 c) (1663211)𝑜𝑛 ç) (8.88888)𝑜𝑛 d) (0.256321)𝑜𝑛 e) (0.7777777)𝑜𝑛 f) (2.344431)𝑜𝑛 S-7) Aşağıda verilen fonksiyonların (1.0, 2.0, 4.0) noktasında değeri bulunduğunda yapılacak olan mutlak ve bağıl hata büyüklüklerini bulunuz. Burada 𝑥, 𝑦 𝑣𝑒 𝑧 değişkenlerinde yapılan hata büyüklüğünün 0.4𝑥10−3 olduğunu kabul ediniz. 𝑥 a) 𝑓(𝑥, 𝑦, 𝑧) = 2𝑥 + 4𝑦 − 𝑧 b) 𝑓(𝑥, 𝑦, 𝑧) = 𝑦 𝑧 c) 𝑓(𝑥, 𝑦, 𝑧) = S-8) 4 1 −1 𝐴=[ 1 3 1 ] −1 1 2 S-9) −3 1 0 𝐴 = [ 1 3 1] 0 1 3 𝑎 S-10) 𝐴 = [𝑑 0 𝑐 𝑎 𝑑 0 𝑐] 𝑎 S-11) Aşağıda verilen denklem sistemlerini LU ayrışım yöntemiyle çözünüz. 4𝑥 − 2𝑥2 = −6 a) { 1 3𝑥1 − 𝑥2 = −2 𝑥2𝑧 √𝑦 b) c) 2𝑥 − 9𝑥2 = 7 { 1 7𝑥1 − 13𝑥2 = 6 15𝑥1 + 22𝑥2 = 37 { 37𝑥1 − 19𝑥2 = 16 3𝑥1 + 2𝑥2 + 𝑥3 = 14 ç) { 12𝑥1 + 13𝑥2 + 6𝑥3 = 40 −3𝑥1 + 8𝑥2 + 9𝑥3 = −28 3𝑥1 + 9𝑥2 + 5𝑥3 = 17 d) { 21𝑥1 − 13𝑥2 − 6𝑥3 = 2 −3𝑥1 + 18𝑥2 + 9𝑥3 = 24 −13𝑥1 − 2𝑥2 + 4𝑥3 = 11 e) {12𝑥1 − 13𝑥2 + 6𝑥3 = −5 −3𝑥1 + 8𝑥2 − 9𝑥3 = 4 S-12) Aşağıda verilen denklem sistemlerini Cholesky yöntemiyle çözünüz. 3𝑥1 + 2𝑥2 + 𝑥3 = 5 𝑎) {2𝑥1 + 13𝑥2 + 8𝑥3 = 15 𝑥1 + 8𝑥2 + 9𝑥3 = 9 5𝑥1 + 9𝑥2 + 7𝑥3 = 21 b) {9𝑥1 − 13𝑥2 − 6𝑥3 = −10 7𝑥1 − 6𝑥2 + 3𝑥3 = 4 15𝑥1 − 2𝑥2 − 3𝑥3 = 2 c) { 12𝑥1 − 7𝑥2 + 9𝑥3 = 25 −3𝑥1 + 9𝑥2 + 6𝑥3 = 33 S-13) Aşağıda verilen denklem sistemlerini 2𝑥1 − 𝑥2 + 𝑥3 = 1 𝑎) {2𝑥1 + 3𝑥2 − 4𝑥3 = 5 𝑥1 + 4𝑥2 − 2𝑥3 = 5 3𝑥1 + 2𝑥2 − 𝑥3 = 4 𝑏) {2𝑥1 + 5𝑥2 − 8𝑥3 = −1 𝑥1 − 6𝑥2 + 2𝑥3 = −3 S-14) Aşağıdaki denklemlerin sıfır yerini Newton-Raphson yöntemini kullanarak 𝜀 = 10−4 tolerans için bulunuz. a)𝑥 3 − 2𝑥 2 − 4 = 0 [1, 4] b)𝑥 3 + 5𝑥 2 − 1 = 0 [−3, 0] c) 𝑥 + sin(𝑥) − 1 = 0 d) 𝑥 3 − 𝑥 − 1 = 0 𝜋 2 [0, ] [1,2] Derste tanımlanan newton.m programı ve MATLAB ait fzero programları yardımıyla Soru 14 de verilen fonksiyonların sıfıryerlerini belirleyiniz. S-15) Aşağıdaki denklemlerin sıfır yerini sekant yöntemini kullanarak 𝜀 = 10−4 tolerans için bulunuz. a)𝑥 = 2−𝑒 𝑥 +𝑥 2 3 b)3𝑥 3 − 𝑒 𝑥 = 0 [0, 2] [0, 2] c) 𝑒 𝑥 + 2cos(𝑥) − 5 = 0 d) 𝑥 2 + 5 sin(𝑥) − 1 = 0 [0,2] [0,2] Derste tanımlanan sekant.m programı ve MATLAB ait fzero programları yardımıyla Soru 15 de verilen fonksiyonların sıfıryerlerini belirleyiniz. S-16) Aşağıdaki denklemlerin sıfır yerini regula false yöntemini kullanarak 𝜀 = 10−3 tolerans için bulunuz. a)𝑥 = 2−𝑒 𝑥 +𝑥 2 3 b)3𝑥 3 − 𝑒 𝑥 = 0 [0, 2] [0, 2] c) 𝑒 𝑥 + 2cos(𝑥) − 5 = 0 d) 𝑥 2 + 5 sin(𝑥) − 1 = 0 [0,2] [0,2] Derste tanımlanan regulafalse.m programı ve MATLAB ait fzero programları yardımıyla Soru 16 da verilen fonksiyonların sıfıryerlerini belirleyiniz. Doç. Dr. Mehmet MERDAN Matematik Mühendisliği Bölümü
© Copyright 2024 Paperzz