Trong bài toán bình phƣơng nhỏ nhất chuẩn tắc , chúng ta thay thế tổng các giá trị tuyệt đối cho tổng bình phƣơng đƣợc sử dụng trong chuẩn tắc Tikhonov, để có đƣợc:
‖ ‖ ‖ ‖ (2.3)
Trong đó ‖ ‖ ∑ | biểu thị chuẩn của x, và λ là tham số chuẩn tắc. Ta gọi (2.3) là một LSP chuẩn . Bài toán này luôn luôn có một giải pháp, nhƣng nó không cần phải là duy nhất.
Một số tính chất cơ bản của chuẩn tắc đƣợc liệt kê ở đây, chỉ ra điểm tƣơng đồng và khác biệt so với chuẩn tắc .
30
• Tính phi tuyến. Từ (2.2), chúng ta thấy rằng, trong chuẩn tắc Tikhonov, vector x là một hàm tuyến tính của vector quan sát y. Ngƣợc lại, trong chuẩn tắc , vector x không phải là tuyến tính với y.
• Hạn chế hoạt động khi λ → 0. Chuẩn tắc chỉ ra hạn chế hoạt động khác biệt với chuẩn tắc khi λ → 0. Giải pháp cho chuẩn tắc (2.3) hội tụ đến điểm trong tất cả các điểm thỏa mãn khi λ → 0.
• Sự hội tụ hữu hạn đến không khi λ → ∞. Nhƣ trong chuẩn tắc Tikhonov, giải pháp tối ƣu có xu hƣớng tiến đến không khi λ → ∞. Tuy nhiên, với chuẩn tắc , sự hội tụ xảy ra cho một giá trị hữu hạn của λ:
‖ ‖ (2.4)
Trong đó ‖ ‖ | ký hiệu chuẩn của vector . Đối với
, giải pháp tối ƣu là 0. Ngƣợc lại, giải pháp tối ƣu cho bài toán
Tikhonov bằng không chỉ trong giới hạn khi .
• Đƣờng chuẩn tắc. Giải pháp của bài toán Tikhonov thay đổi mƣợt khi tham số chuẩn tắc thay đổi trong khoảng [0, ∞). Ngƣợc lại, đƣờng chuẩn tắc của (2.3), tức là, họ các giải pháp khi thay đổi trên (0, ∞), có tính chất đƣờng giải pháp tuyến tính từng khúc [2.2]: Có các giá trị , với
, đƣờng chuẩn tắc là đƣờng cong tuyến tính từng
khúc trên
Trong đó giải bài toán (2.3) với . (Vì vậy, và
khi ).
Quan trọng hơn, chuẩn tắc đặc biệt thực hiện đƣợc với vector thƣa, tức là có tƣơng đối ít các hệ số khác không. (Khi giảm, nó có xu hƣớng
31
thƣa hơn, nhƣng điều này không cần thiết thiết [2.3], [2.4].) Ngƣợc lại, chuẩn tắc của bài toán Tikhonov có tất cả các hệ số khác không.
Gần đây, ý tƣởng của chuẩn tắc nhận đƣợc rất nhiều sự quan tâm trong lĩnh vực xử lý tín hiệu và thống kê. Trong lĩnh vực xử lí tín hiệu, ý tƣởng của chuẩn tắc xuất hiện trong một vài bối cảnh, bao gồm Basis pursuit denoising [2.5] và phƣơng pháp khôi phục tín hiệu từ các phép đo không hoàn thiện (ví dụ, [2.6], [2.7], [2.8], [2.9], [2.10], [2.11], [2.12]). Trong thống kê, ý tƣởng của chuẩn tắc đƣợc sử dụng trong thuật toán nổi tiếng Lasso [2.4] để lựa chọn tính chất và mở rộng của nó, bao gồm elastic net
[2.13].
Một trong những bài toán này không có dạng chuẩn (2.3) nhƣng nó có một biểu diễn tổng quát hơn:
‖ ‖ ∑ | |
(2.5) trong đó là tham số chuẩn tắc (biến tƣơng ứng với là không đƣợc chuẩn tắc).
Bây giờ, chúng ta xem xét khía cạnh tính toán của chuẩn tắc . Không có một công thức phân tích hoặc biểu diễn nào cho bài toán tối ƣu đối với chuẩn tắc (2.3), tƣơng tự (2.2); giải pháp của nó phải đƣợc tính số. Hàm mục tiêu trong (2.3) là lồi nhƣng không khả vi, do đó, giải quyết bài toán là một thách thức tính toán hơn là giải quyết bài toán (2.1).
Phƣơng pháp tổng quát cho bài toán lồi không khả vi, chẳng hạn phƣơng pháp ellipsoid hay phƣơng pháp subgradient [2.14], [2.15], có thể đƣợc sử dụng để giải quyết bài toán (2.3). Các giải pháp này thƣờng rất chậm.
Chuẩn tắc (2.3) có thể đƣợc chuyển thành bài toán toàn phƣơng lồi (Convex quadratic problem), với những hạn chế không cân bằng tuyến tính.
32
Bài toán toàn phƣơng (QP) tƣơng ứng có thể đƣợc giải bằng phƣơng pháp tối ƣu lồi chuẩn, nhƣ phƣơng pháp điểm nội (interior-point method) [2.16], [2.17], [2.18], [2.19]. Phƣơng pháp điểm nội chuẩn đƣợc thực hiện theo giải pháp chung phổ quát, bao gồm MOSEK [2.20], nó có thể sẵn sàng thực thi với bài toán có kích thƣớc vừa và nhỏ. Các phƣơng pháp chuẩn không thể thực thi các bài toán lớn, trong đó có thuật toán nhanh cho toán tử ma trận- vector với A và Phƣơng pháp điểm nội chuyên biệt, khai thác các thuật toán nhƣ vậy, có thể mở rộng cho các bài toán lớn, nhƣ đƣợc chứng minh ở [2.5], [2.21]. Việc thực thi chất lƣợng cao của phƣơng pháp điểm nội chuyên biệt, bao gồm ll-magic [2.22] và PDC0 [2.23], chúng sử dụng thuật toán lặp, chẳng hạn nhƣ conjugate gradient (CG) hoặc thuật toán LSQR [2.24], để tính toán các bƣớc tìm kiếm. (adsbygoogle = window.adsbygoogle || []).push({});
Gần đây, một số nhà nghiên cứu đã đề xuất phƣơng pháp homotopy và các biến thể để giải quyết bài toán [2.25], [2.26], [2.2], [2.27], [2.28]. Sử dụng thuộc tính tuyến tính từng khúc của đƣờng chuẩn tắc, phƣơng pháp theo đƣờng (path-following method) có thể tính toán hiệu quả toàn bộ đƣờng giải pháp trong bài toán . Khi giải pháp (13) rất thƣa, giải pháp này có thể rất nhanh, vì số lƣợng các nút mà phƣơng pháp cần để tìm là vừa phải [2.25]. Mặt khác, phƣơng pháp theo đƣờng có thể chậm, nó thƣờng sử dụng trong bài toán quy mô lớn. Các phƣơng pháp tính toán khác đƣợc phát triển gần đây cho bài toán , bao gồm coordinate-wise descent method [2.29], một phƣơng pháp cố định điểm tiếp [2.30], Bregman phƣơng pháp dựa trên quy tắc lặp đi lặp lại [2.31], [2.32], phƣơng pháp tối ƣu không gian con tuần tự (sequential subspace optimization method) [2.33], phƣơng pháp tối ƣu đƣờng bao (bound optimization method) [2.34], phƣơng pháp co lặp (iterated shrinkage method) [2.35], [2.36], phƣơng pháp gradient [2.37], và thuật toán chiếu gradient
33
(gradient projection algorithm) [2.38]. Một số trong các phƣơng pháp này,bao gồm cả thuật toán chiếu gradient [2.38] có thể xử lý hiệu quả bài toán lớn.
34
Chƣơng 3: Mô phỏng và kết quả 3.1. Mô hình phƣơng pháp lặp vi phân Born (DBIM).
3.2. Chi tiết tính toán phƣơng pháp lặp vi phân Born (DBIM).
Đối tƣợng cần khảo sát chính là vật thể hình trụ tròn O(r) có kích thƣớc rất nhỏ (môi trƣờng B1) nằm trong môi trƣờng B2 (tƣơng ứng nhƣ khối u ở trong môi trƣờng nào đó). Mục tiêu của chúng ta là dựng đƣợc ảnh của vật thể trụ tròn, đó chính là vùng cần quan tâm ROI (Region Of Interest). Vùng diện tích quan tâm này đƣợc chia thành N×N ô vuông (mỗi ô vuông gọi là một pixel) có kích thƣớc là h. Số lƣợng máy phát là Nt và máy thu là Nr. Theo lí thuyết về sóng âm, hàm mục tiêu O(r) (vật thể hình trụ tròn) đƣợc tính bởi công thức:
Sử dụng hàm Green tìm phƣơng trình liên tục tính 𝑝𝑠𝑐
Rời rạc phƣơng trình sóng, sử dụng phƣơng pháp moment (MoM) Tính đƣợc áp suất tại các điểm bên trong và bên ngoài (tán xạ)
đối tƣợng là 𝑝, 𝑝𝑠𝑐
Tính ma trận B, C
Sử dụng phƣơng pháp xấp xỉ Born để tìm mối quan hệ tuyến
tính giữa 𝛥𝑝𝑠𝑐 và 𝛥𝑂
Sử dụng phƣơng pháp Tikhonov để tìm 𝛥𝑂
Sử dụng phƣơng pháp lặp để tìm đối tƣợng 𝑂
Phƣơng trình sóng trong môi trƣờng không đồng nhất chứa đối
35 Hình 3.1. Cấu hình hệ đo 1 1 2 ω - if r R 2 2 c c 1 0 Ο r = 0 if r > R (3.1)
Với và là tốc độ truyền sóng trong đối tƣợng và tốc độ truyền trong nƣớc, f là tần số sóng siêu âm, là tần số góc ( là bán kính của đ i tư ng.
Sử dụng sơ đồ cấu hình hệ đo nhƣ trong hình 3.1, bằng cách sử dụng DBIM để tái tạo lại độ tƣơng phản âm thanh tán xạ để xác định khối u trong môi trƣờng. Giả sử có một không gian vô hạn chứa môi trƣờng đồng nhất chẳng hạn là nƣớc, số sóng là . Trong môi trƣờng đó có vật với số sóng là
phụ thuộc vào không gian trong vật. Phƣơng trình truyền sóng của hệ thống có thể đƣợc cho nhƣ phƣơng trình (3.2).
( ) (3.2)
Viết lại dƣới dạng tích phân ta có:
(3.3) ∬ ⃗⃗⃗ ⃗⃗⃗ | ⃗ ⃗⃗⃗| ⃗⃗⃗ (3.4) Ở đó psc r là sóng tán xạ, pinc r là sóng tới và G(.) là hàm Green. h Transmitter Meshing area Receivers
36
Hàm mục tiêu cần đƣợc khôi phục từ dữ liệu tán xạ đƣợc xác định bởi:
⃗ ⃗ (3.5)
Bằng phƣơng pháp moment (MoM) áp suất của các điểm bên trong đối tƣợng có thể đƣợc tính nhƣ sau:
̅ ( ̅ ̅ ) (3.6)
Áp suất của các điểm bên ngoài đối tƣợng (áp suất tán xạ) là:
̅ (3.7) (adsbygoogle = window.adsbygoogle || []).push({});
Hai biến chƣa biết là ̅ và trong công thức (3.6) và (3.7), trong trƣờng hợp này áp dụng xấp xỉ Born loại 1 với phƣơng trình (3.4) và theo (3.6), (3.7) ta có:
̅ ̅ (3.8)
Với ̅ ̅ ; Ở đó B là ma trận ứng với hệ số G0(r,r‟) từ
các pixel tới máy thu, C là ma trận ứng với hệ số G0(r,r‟) giữa các pixel, I là ma trận đơn vị, và D(.) là toán tử chéo hóa.
Với mỗi bộ phát và bộ thu, chúng ta có một ma trận ̅ và một giá trị vô hƣớng . Thấy rằng vector chƣa biết có giá trị bằng với số pixel của RIO. Hàm mục tiêu (Object function) có thể đƣợc tính bằng cách lặp:
(3.9)
Với và là giá trị của hàm mục tiêu ở bƣớc hiện tại và bƣớc trƣớc đó. có thể đƣợc tìm bằng quy tắc Tikhonov:
‖ ̅ ‖ ‖ ‖ (3.10)
Trong đó ̅ là ( vector chứa giá trị sai khác gi a kết quả đo và kết quả tiên đoán tín hiệu siêu m tán ạ là ma trận
37
3.3. Các phƣơng pháp toán học khôi phục ảnh và chƣơng trình mô phỏng. phỏng. Phƣơng pháp Công thức Moo-PenPenrose Psedoinverse ‖ ‖ Tikhonov regularization ‖ ‖ ‖ ‖ L1 regularization ‖ ‖ ‖ ‖ Các vấn đề tính toán: - Moo-PenPenrose Psedoinverse: - Tikhonov regularization: - L1 regularization: ⁄
3.4. Chƣơng trình code của các thuật toán.
3.4.1. Chƣơng trình code của phƣơng pháp Moo-PenPenrose Psedoinverse Psedoinverse
%======Khôi phục dữ liệu sử dụng phƣơng pháp "Moore-Penrose inverse"=========
Code Ý nghĩa
SC1 = pinv(Mt)*p_sc_t;
SC1: Ma trận tán xạ
Mt: Ma trận M trong DBIM p_sc_t: Áp suất tán xạ
SC1 = reshape(SC1,N,N); Chuyển ma trận SC1 thành ma trận vuông
NxN. figure;
38
3.4.2. Chƣơng trình code của phƣơng pháp Tikhonov regularization
1. Chương trình khôi phục chính.
Code Ý nghĩa
RRE = 2^(-5); Lỗi nền, RRE càng nhỏ, chất lƣợng
càng tốt
gama = .5*(5.9303e-008)^2*1e-4; Tham số chuẩn tắc Tikhonov, càng
nhỏ, càng chậm delta_sound = test_NCG(Mt,delta_sc_t,1e3,RRE,g ama); Hàm khôi phục “test_NCG” SC2 =
SC1+reshape(delta_sound,N,N); Cập nhật sự tƣơng phản âm thanh (adsbygoogle = window.adsbygoogle || []).push({});
2. Chương trình con “test_NCG”.
function[delta_sound]=test_NCG(Mt,delta_sc_t,ni,RRE,gama)
3.4.3. Chƣơng trình code của phƣơng pháp L1 regularization
1. Chương trình khôi phục chính.
Code
lambda = 1e-15; rel_tol =RRE;
AA=[real(Mt) – imag(Mt);imag(Mt) real(Mt);]; yy=[real(delta_sc_t);imag(delta_sc_t)];
[xx,status]=l1_ls(AA,yy,lambda,rel_tol);
delta_sound=xx(1:length(xx)/2)+j*xx(1+length(xx)/2:end); SC1=SC1+reshape(delta_sound,N,N);
39 l1_ls(AA,yy,lambda,rel_tol);
3.5. Kết quả mô phỏng.
Tham số mô phỏng: Tần số 1MHz, N = 21, Nsum = 8. Đƣờng kính đối tƣợng 7.3mm, Độ tƣơng phản âm 7%, Nhiễu Gaussian 10%, Khoảng cách máy phát , máy thu đến tâm = 100mm.
Hình 3.2: Hàm mục tiêu lí tư ng
Hình 3.2 là hàm mục tiêu lí tƣởng trong kịch bản mô phỏng (N = 21) (trong thực tế, nó chính là khối u lạ trong môi trƣờng đồng nhất). Mục tiêu của chúng ta là khôi phục đƣợc đối tƣợng trên, sử dụng các mô hình, các kĩ thuật xử lí tín hiệu để khôi phục đƣợc ảnh với chất lƣợng tốt nhất và thời gian khôi phục nhanh nhất.
3.5.1.Mô phỏng sử dụng phƣơng pháp Tikhonov.
Thay đổi số MP,
MT Lỗi chuẩn hóa từ vòng lặp 1 đến vòng lặp 8 Thời gian
NP = NT = 8 0.9490 0.9427 0.9424 0.9424 0.9424 0.9424 0.9424 0.9424 61.846577 NP = NT = 10 0.6207 0.5179 0.4872 0.4707 0.4600 0.4524 0.4470 0.4431 87.693626 NP = NT = 12 0.8259 0.7999 0.7943 0.7925 0.7916 0.7912 0.7911 0.7911 94.968234 NP = NT = 14 0.8000 0.7464 0.7281 0.7195 0.7145 0.7111 0.7087 0.7070 130.845553 -2 -1 0 1 2 -2 -1 0 1 2 0 5 10 15 20 25 30
Ideal object function
p e rc e n t o f th e s o u n d c o n tr a s t
40 NP = NT = 16 0.4290 0.2691 0.2382 0.2320 0.2302 0.2294 0.2291 0.2290 130.258239 NP = NT = 18 0.4017 0.2419 0.2081 0.1960 0.1900 0.1862 0.1834 0.1812 152.313615 NP = NT = 20 0.3548 0.2266 0.1822 0.1584 0.1449 0.1363 0.1302 0.1256 168.393346 NP = NT = 22 0.4038 0.1778 0.0949 0.0755 0.0668 0.0615 0.0576 0.0546 192.820992 NP = NT = 24 0.3179 0.1999 0.1457 0.1253 0.1125 0.1031 0.0959 0.0900 233.365111 NP = NT = 26 0.2677 0.0909 0.0425 0.0286 0.0220 0.0178 0.0150 0.0129 269.807922 NP = NT = 28 0.2545 0.1055 0.0401 0.0230 0.0161 0.0124 0.0101 0.0084 343.149005 NP = NT = 30 0.3551 0.1428 0.0667 0.0438 0.0318 0.0244 0.0194 0.0158 385.055147 Nhận xét:
+ Khi tăng số máy phát, máy thu thì lỗi chuẩn hóa giảm, nhƣng thời gian tạo ảnh tăng.
+ Số máy phát, máy thu NP = NT = 28 cho lỗi nhỏ nhất là 0.0084.
3.5.2.Mô phỏng sử dụng phƣơng pháp L1 (N2
= 484).
Thay đổi số MP, MT Lỗi chuẩn hóa từ vòng lặp 1 đến vòng lặp 8 Thời gian
NP = NT = 8 (NpNt = 64) 1.2559 1.2624 1.2624 1.2623 1.2623 1.2623 1.2623 1.2623 121.439561 NP = NT = 10 (NpNt = 100) 1.1440 1.1354 1.1353 1.1353 1.1353 1.1353 1.1353 1.1353 122.498515 NP = NT = 12 (NpNt = 144) 1.0836 1.0787 1.0786 1.0786 1.0786 1.0786 1.0786 1.0786 258.668851 NP = NT = 14 (NpNt = 196) 0.9686 0.9690 0.9681 0.9681 0.9681 0.9681 0.9681 0.9681 298.946565 NP = NT = 16 (NpNt = 256) 0.7962 0.7972 0.7978 0.7977 0.7977 0.7977 0.7977 0.7977 413.842807 NP = NT = 18 (NpNt = 324) 0.6470 0.6211 0.6221 0.6223 0.6223 0.6223 0.6223 0.6223 761.758941 NP = NT = 20 (NpNt = 400) 0.7664 0.6379 0.6375 0.6377 0.6376 0.6376 0.6376 0.6376 697.251724 NP = NT = 22 (NpNt = 484) 2.8511 0.7848 0.4247 0.3971 0.3946 0.3938 0.3934 0.3932 1124.067465 NP = NT = 24 (NpNt = 576) 0.5535 0.0570 0.0102 0.0080 0.0078 0.0077 0.0076 0.0075 474.090601 NP = NT = 26 (NpNt = 676) 0.3239 0.0450 0.0036 0.0017 0.0016 0.0016 0.0016 0.0016 424.798644 NP = NT = 28 (NpNt = 784) 0.8879 0.0157 0.0018 0.0014 0.0013 0.0013 0.0013 0.0013 445.144473 NP = NT = 30 (NpNt = 900) 0.1484 0.0180 0.0011 0.0005 0.0005 0.0005 0.0005 0.0005 562.230397
Thay đổi số MP, MT Lỗi chuẩn hóa từ vòng lặp 1 đến vòng lặp 8 Thời gian
NP = NT = 16 (NpNt = 256) 0.7962 0.7972 0.7978 0.7977 0.7977 0.7977 0.7977 0.7977 413.842807
41 NP = 24; NT = 12 (NpNt = 288) 0.8594 0.8364 0.8366 0.8366 0.8366 0.8366 0.8366 0.8366 469.800523 NP = NT = 22 (NpNt = 484) 2.8511 0.7848 0.4247 0.3971 0.3946 0.3938 0.3934 0.3932 1124.067465 NP = 30; NT = 16 (NpNt = 484) 2.7425 0.6821 0.4455 0.4191 0.4155 0.4148 0.4144 0.4141 1530.689987 NP = 16; NT = 30 (NpNt = 484) 2.5267 0.7306 0.5479 0.5372 0.5357 0.5351 0.5349 0.5342 1876.443748 NP = NT = 26 (NpNt = 676) 0.3239 0.0450 0.0036 0.0017 0.0016 0.0016 0.0016 0.0016 424.798644 NP = 21; NT = 32 (NpNt = 676) 0.5951 0.0462 0.0096 0.0079 0.0077 0.0076 0.0075 0.0074 676.611183 NP = 32; NT = 21 (NpNt = 676) 0.2649 0.0320 0.0025 0.0014 0.0013 0.0012 0.0012 0.0012 492.551895 Nhận xét: (adsbygoogle = window.adsbygoogle || []).push({});
+ Khi tăng số máy phát, máy thu thì lỗi chuẩn hóa giảm, nhƣng thời gian tạo ảnh tăng.
+ Đặc biệt với phƣơng pháp L1, số máy phát, máy thu nhỏ nhƣng vẫn khôi phục ảnh thành công. Trong khi đó, phƣơng pháp Tikhonov lại không thể. Vậy thì với phƣơng pháp L1, thì số máy phát, máy thu nhỏ nhất bằng 8 mà vẫn khôi phục ảnh thành công.
+ Khi sử dụng phƣơng pháp L1 thì số máy phát và máy thu bằng nhau cho kết quả khôi phục tốt (Np = Nt).
3.5.3. Đánh giá kết quả sử dụng phƣơng pháp Tikhonov và phƣơng pháp L1.
Thay đổi số
MP, MT Lỗi chuẩn hóa từ vòng lặp 1 đến vòng lặp 8 Thời gian
NP = NT = 8 Tikhonov 0.9490 0.9427 0.9424 0.9424 0.9424 0.9424 0.9424 0.9424 61.846577 L1 1.2559 1.2624 1.2624 1.2623 1.2623 1.2623 1.2623 1.2623 121.439561 % Sai lệch 24.44% 25.32% 25.34% 25.34% 25.34% 25.34% 25.34% 25.34% 49.10% NP = NT = 10 Tikhonov 0.6207 0.5179 0.4872 0.4707 0.4600 0.4524 0.4470 0.4431 87.693626 L1 1.1440 1.1354 1.1353 1.1353 1.1353 1.1353 1.1353 1.1353 122.4985 15 % Sai lệch 45.74% 54.38% 57.08% 58.54% 59.48% 60.15% 60.62% 60.97% 28.41% NP = NT = 12 Tikhonov 0.8259 0.7999 0.7943 0.7925 0.7916 0.7912 0.7911 0.7911 94.968234 L1 1.0836 1.0787 1.0786 1.0786 1.0786 1.0786 1.0786 1.0786 258.668851 % Sai lệch 23.78% 25.84% 26.35% 26.52% 26.60% 26.64% 26.65% 26.65% 63.28% NP = NT = 14 Tikhonov 0.8000 0.7464 0.7281 0.7195 0.7145 0.7111 0.7087 0.7070 130.845553 L1 0.9686 0.9690 0.9681 0.9681 0.9681 0.9681 0.9681 0.9681 298.946565 % Sai lệch 17.4% 22.97% 24.79% 25.67% 26.19% 26.54% 26.79% 26.97% 56.23% NP = NT = 16 Tikhonov 0.4290 0.2691 0.2382 0.2320 0.2302 0.2294 0.2291 0.2290 130.258239 L1 0.7962 0.7972 0.7978 0.7977 0.7977 0.7977 0.7977 0.7977 413.842807
42 % Sai lệch 46.11% 66.24% 70.14% 70.91% 70.91% 71.24% 71.27% 71.29% 68.52% NP = NT = 18 Tikhonov 0.4017 0.2419 0.2081 0.1960 0.1900 0.1862 0.1834 0.1812 152.313615 L1 0.6470 0.6211 0.6221 0.6223 0.6223 0.6223 0.6223 0.6223 761.758941 % Sai lệch 37.91% 61.05% 66.54% 68.50% 69.46% 70.07% 70.52% 70.88% 80% NP = NT = 20