Nghiên Cứu Nâng Cao Hiệu Quả Thu Hồi Năng Lượng Của Hệ Thống Phanh Tái Sinh Trên Ô Tô.pdf

Untitled BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN ÁN TIẾN SĨ NGÀNH Tp Hồ Chí Minh, tháng 2020 S KA 0 0 0 0 1 7 LUẬN ÁN TIẾN SĨ NGÀNH KỸ THUẬT CƠ KHÍ Tp Hồ Chí[.]

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN ÁN TIN S '1*7817ễ1* 1*+,ầ1&81ặ1*&$2+,8487+8+,11*/1* &$+7+1*3+$1+7ẩ,6,1+75ầ1é7é NGNH: 7+87&.+ậ SKA0 0 Tp Hồ Chí Minh, tháng 2020 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH DƯƠNG TUẤN TÙNG NGHIÊN CỨU NÂNG CAO HIỆU QUẢ THU HỒI NĂNG LƯỢNG CỦA HỆ THỐNG PHANH TÁI SINH TRÊN Ô TÔ LUẬN ÁN TIẾN SĨ NGÀNH: KỸ THUẬT CƠ KHÍ Tp Hồ Chí Minh, tháng 11/2020 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH DƯƠNG TUẤN TÙNG NGHIÊN CỨU NÂNG CAO HIỆU QUẢ THU HỒI NĂNG LƯỢNG CỦA HỆ THỐNG PHANH TÁI SINH TRÊN Ơ TƠ NGÀNH: KỸ THUẬT CƠ KHÍ - 9520103 Hướng dẫn khoa học: PGS-TS ĐỖ VĂN DŨNG PGS-TS NGUYỄN TRƯỜNG THỊNH Phản biện 1: TS Nguyễn Chí Thanh Phản biện 2: PGS.TS Cao Hùng Phi Phản biện 3: PGS.TS Dương Việt Dũng LÝ LỊCH KHOA HỌC CỦA NGHIÊN CỨU SINH I LÝ LỊCH SƠ LƯỢC Họ tên: DƯƠNG TUẤN TÙNG Giới tính: Nam Ngày, tháng, năm sinh: 15/07/1980 Nơi sinh: Hà Nam Quê quán: Hà Nam Dân tộc: Kinh Học vị cao nhất: Thạc Sỹ Năm, nước nhận học vị: 2010, Việt Nam Chức vụ: Trưởng ngành CNKT ô tô Khoa Đào tạo Chất lượng cao Đơn vị công tác : Trường Đại học Sư Phạm Kỹ Thuật TP.Hồ Chí Minh Chỗ ở riêng địa liên lạc: 40A, tổ 9, khu phố 3, phường An Bình, TP Biên Hịa, Tỉnh Đồng Nai Điện thoại liên hệ: 0914805623 Email: tungdt@hcmute.edu.vn II QUÁ TRÌNH ĐÀO TẠO Đại học - Hệ đào tạo: Chính qui - Nơi đào tạo: Trường Đại học Sư hạm Kỹ thuật TP HCM - Ngành học: Cơ khí động lực - Nước đào tạo: Việt Nam Năm tốt nghiệp: 2004 Sau đại học - Thạc sĩ chuyên ngành: Cơ khí tơ Năm cấp bằng: 2010 - Nơi đào tạo: trường Đại học Sư phạm Kỹ thuật Tp HCM – Việt Nam - Ngoại ngữ: Tiếng Anh Mức độ sử dụng: B2 III Q TRÌNH CƠNG TÁC - 2004-2005: Làm việc công ty ô tô TOYOTA Biên Hòa - 2005 đến nay: Giảng viên Trường Đại học Sư phạm Kỹ thuật Tp.HCM IV LĨNH VỰC NGHIÊN CỨU/CHUYÊN MƠN - Hệ thống truyền lực tơ - Hệ thống điều khiển chuyển động ô tô - Hệ thống phanh tái sinh ô tô - Kỹ thuật sửa chữa thân vỏ sơn xe i V CÁC CÔNG TRÌNH NGHIÊN CỨU KHOA HỌC 5.1 Các cơng trình cơng bố TT Tên cơng trình Thiết kế, chế tạo m áy phân tích màu sơn tơ An overview of research and proposed experiment model of regenerative braking system based on the conventional vehicle powertrain Nghiên cứu thiết kế mô động lực học thu hồi lượng từ hệ thống phanh ô tô Một nghiên cứu thực nghiệm thu hồi lượng tái tạo phanh áp dụng cho xe ô tô có kiểu hệ thống truyền lực truyền thống Nghiên cứu mơ đặc tính động lực ô tô dựa phần mềm AVL Cruise Research on kinetic energy recovery of conventional vehicle based on regenerative braking system Tên tác giả Dương Tuấn Tùng Duong Tuan Tung Dương Tuấn Tùng Dương Tuấn Tùng Dương Tuấn Tùng Duong Tuan Tung ii Nơi công bố (tên tạp chí, tuyển tập) Tạp chí khoa học cơng nghệ trường Đại học Kỹ thuật số 85-2011, ĐH Bách khoa Hà Nội Năm công bố 2011 The International Conference of Green Technology and Sustainable Development 2nd, HCMUTE, 2014 2014 Hội nghị khoa học cơng nghệ tồn quốc Cơ khí lần thứ Tp.Hồ Chí Minh, ngày 06 tháng 11 năm 2015 2015 Hội nghị khoa học cơng nghệ tồn quốc Cơ khí lần thứ Tp.Hồ Chí Minh, ngày 06 tháng 11 năm 2015 2015 Tạp chí khoa học giáo dục kỹ thuật, Đại học SPKT Tp.HCM 2016 The Fifth International Multi-Conference on Engineering and T echnology Innovation 2016 (IMETI2016), 2016 Taichung Taiwan, November, 2016 Research on using PID algorithm to control the inertial energy recovery of vehicle based on regenerative braking system Research on braking force distribution in regenerative braking system apply to conventional vehicle 10 11 5.2 TT Duong Tuan Tung IEEE International Conference on Systems Science and Engineering, July 2123, 2017, HCMUTE 2017 Duong Tuan Tung IEEE International Conference on Green Technology and Sustainable Development (GTSD) December 2018 2018 Research on controlling of experiment model to Journal of Technical Duong Tuan evaluate of kinetic Education Science, Tung energy recovery system 2018 based on driving cycles Research on Designing Journal of Science & the Regenerative Duong Tuan Technology, Braking System Apply Tung Technical University to Conventional No 135 (2019) Vehicle Research on improving IJSRD - International the fuel consumption of Journal for Scientific conventional Research & Duong Tuan powertrain vehicle by Development| Vol 7, Tung using a kinetic energy Issue 03, 2019 | ISSN recovery system (online): 2321-0613 Các đề tài nghiên cứu khoa học thực Tên đề tài nghiên cứu/ Lĩnh vực ứng dụng Nghiên cứu, thiết kế, chế tạo mơ hình hệ thống sơn tự động Thiết kế phần mềm tra công thức màu sơn tơ Năm hồn thành Đề tài cấp (NN, Bộ, ngành, trường) 2018 2019 2019 Trách nhiệm tham gia đề tài 2006 Nghiên cứu khoa học cấp Trường T2006-64 Chủ trì 2010 Nghiên cứu khoa học cấp Trường Chủ trì iii T2010-22 Nghiên cứu, thiết kế, chế tạo máy phân tích màu sơn tơ Thiết kế mạch điều khiển mơ hình hệ thống phanh ABS máy tính Nghiên cứu đề xuất phương án tích trữ lượng phanh tơ Nghiên cứu, tính tốn mơ hình thử nghiệm bánh đà siêu tốc ứng dụng cho hệ thống phanh tái tạo lượng ô tô Nghiên cứu thiết kế chế tạo thử nghiệm hệ thống phanh tái tạo lượng áp dụng cho xe tơ có kiểu hệ thống truyền lực truyền thống Nghiên cứu điều khiển hệ thống thu hồi lượng quán tính ô tô dựa hệ thống phanh tái tạo lượng Nghiên cứu vấn đề quản lý lượng thu hồi phân phối lực phanh hệ thống phanh tái sinh ô tô 10 Nghiên cứu tối ưu hóa lượng thu hồi từ hệ thống phanh tái sinh ô tô 2012 Nghiên cứu khoa học cấp Trường trọng điểm T2011-12TĐ Chủ trì 2013 Nghiên cứu khoa học cấp Trường T2012-13 Chủ trì 2014 Nghiên cứu khoa học cấp Trường T2014-60 Chủ trì 2015 Nghiên cứu khoa học cấp Trường T2015-63 Chủ trì 2016 Nghiên cứu khoa học cấp Trường trọng điểm T2016-63TĐ Chủ trì 2017 Nghiên cứu khoa học cấp Trường trọng điểm T2017-28TĐ Chủ trì 2018 Nghiên cứu khoa học cấp Trường trọng điểm T2018-97TĐ Chủ trì 2019 Nghiên cứu khoa học cấp Trường trọng điểm T2019-97TĐ Chủ trì Tơi xác nhận thơng tin được ghi ở hồn tồn xác TP, Hồ Chí Minh., ngày tháng năm 20 Xác nhận quan công tác Người khai iv LỜI CAM ĐOAN Tôi cam đoan cơng trình nghiên cứu tơi Các số liệu, kết nêu Luận án trung thực chưa được cơng bố cơng trình khác Tôi xin cam đoan tham khảo cho việc thực luận án được trích dẫn rõ ràng Tp Hồ Chí Minh, ngày 09 tháng 11 năm 2020 (Ký tên ghi rõ họ tên) v 6.1 Xác định số vòng quay trục Công suất trục vào: P01  24  kw  Chọn sơ tỷ số truyền: ih  0.9 n01  2600  v / ph  Số vòng quay trục vào Số vòng qua trục n02  n01.ih  0,9.2600  2893(v / ph) Số vòng quay bánh trung tâm: nZ 31  nZ 32  n01 n  8979  v / ph  ; nZ 12  01  2893 (v / ph) i1 i nZ 21  12759  v / ph  ; nZ 22  5938  v / ph  n01  2600  v / ph  ; n02  nZ 11  6.2 Xác định mô men xoắn trục Từ số vòng quay công suất trục vào ta chọn sơ momen xoắn trục sau: + Mômen xoắn trục vào T01  9,55.106 P01  88017  N mm  n01 + Mômen xoắn bánh mặt trời sau T12  79216  N mm  + Mômen xoắn bánh hành tinh trước T21  55042( N mm 6.3 Tính tốn thơng số bánh hành tinh kép Xác định ứng suất cho phép  Đối với cặp bánh ăn khớp Ứng suất tiếp xúc cho phép H H  lim Z Z K K      V R XH HL  H  Trong đó: 128 +  H0 lim : ứng suất tiếp xúc cho phép: Cấp nhanh: H lim  HRC  70  2.290  700  650MPa Cấp chậm: H lim  17 HRC  200  17.50  200  1050MPa + ZR: hệ số xét đến độ nhám bề mặt làm việc + ZV: hệ số xét đến ảnh hưởng vận tốc vòng + KXH: hệ số xét đến ảnh hưởng kích thước bánh + KHL: hệ số tuổi thọ + H: hệ số an toàn tiếp xúc, tra bảng - Cấp nhanh : H = 1,1 - Cấp chậm : H = 1,2 Khi thiết kế sơ chọn: ZVZRKXH = Cơng thức tính ứng suất tiếp xúc có dạng: H  lim K    H    HL  H   N K HL  mH HO NH Trong đó: NHO: số chu kì thay đổi ứng suất sở thử tiếp xúc - Cấp nhanh: N HO  30 HB - Cấp chậm: N HO1 2,4  30.2902,4  24,37.106  30 HB 2,4  30.4822,4  82,50.106 mH: bậc đường cong mỏi thử tiếp xúc HB: độ cứng Brinen: - Cấp nhanh: HB = 290 - Cấp chậm: HB = 482 (HRC = 50) + NH: Số chu kì thay đổi ứng suất tương đương 129 N H  60ct  n2  n0  N H  60ct  n1  n0  Trong đó: - c: số bánh vệ tinh, c = - t : tổng thời gian làm việc (tải tĩnh) t  24 24 365.10  46720  h  35 Với cấp nhanh: truyền (Z12 - Z22) N H 12  60.4.46720  n02  n12   6,8241.106 N H 12  N H 02  K HL  N H 22  60.46720  n22  n02   8535.106 N H 22  N H 02  K HK  Do với bánh thẳng:  F0 lim K HL1 650   590,91MPa  H    H 12    H 22   1,1 H - Với cấp chậm: truyền (Z11 - Z21) N H 21  60.46720  n21  n01   3240, 49.106 N H 21  N H 01  K HL  N H 11  N H 01  K HL1  Do với bánh thẳng  F0 lim K HL1 1050   875MPa  H    H 21    H 11   F 1, Ứng suất uốn cho phép   F0 lim  F     F   YRYS K XF K FC K FL  130 Trong đó: +  F  : ứng suất uốn cho phép +  F lim  : ứng suất giới hạn, tra bảng 6.2 [2]:  F0 lim  1,5HB  1,5.290  522MPa - Cấp nhanh: Cấp chậm : +  F0 lim  550MPa  F : hệ số an tồn tính uốn, tra bảng 6.2 [2] δF = 1,75 + YR : hệ số xét đến ảnh hưởng độ nhám bề mặt lượn chân + YS : hệ số xét đến ảnh hưởng độ nhạy ứng suất vật liệu tập trung ứng suất + K XF : hệ số xét đến kích thước bánh ảnh hưởng đến độ bền uốn Khi tính sơ : YRYS K XF  + K FC : hệ số xét đến ảnh hưởng đặt tải trọng + K FL : hệ số tuổi thọ K FL  mF N FO NF Trong đó: mF : bậc đường cong mỏi thử tải uốn; N F : chu kì thay đổi ứng suất, N F  4.10 ; N F : số chu kì thay đổi ứng suất tương đương Vì: N F1  N H 1; N F  N H K FL1  K FL  nên : N F1  N F ; N F1  N F => Vậy:   F0 lim   F     K FL K FC   F  Với bánh trung tâm Z11, Z12: K FC  - bánh làm việc phía 131  F11   550  314, 29MPa 1,75  F12   522  298, 29MPa 1,75 Với bánh vệ tinh Z21, Z22: K FC  0,75 - bánh làm việc phía  F 22    F12 .0,75  298, 29.0,75  223,71MPa  F 21    F11 .0,75  314, 29.0,75  235,72MPa Đối với cặp bánh ăn khớp Ứng suất tiếp xúc cho phép bánh cấp nhanh Z32 2,4 N HO  30 H HB  30.2302,4  13,97.106 ; N H 32  N H 22 q  7939,5 106 u22 32 Z 32  4,3 Trong đó: N H 22  8535.10 ; u22  32  ;q=4 Z 22 => K H 32  K HO => K HL32  Tra bảng:  H  1,1  H0 lim  HB  70  2.230  70  530MPa  H0 lim 530.1 K HL  418,8MPa  H 32  1,1 H Với bánh cấp chậm Z31 2,4 N HO  30 H HB  302302,4  13,97.106 ; N H 31  Trong đó: q = 4; N H 21  3240, 49.10 ; u21 31  => K H 31  K HO => K HL31  132 N H 21.q  3812 106 u2131 Z31  3, Z 21 Tra bảng 6.2 [2]:  h  1,1 H lim  HB  70  2.230  70  530MPa  H0 lim 530.1  418,8MPa K HL   H 31  1,1 H Ứng suất uốn cho phép  F0 lim K FC K FL  F 3  F Trong đó: phía)  F0 lim =1,8HB=1,8.230=414=MPa; K FC  (bánh làm việc  F  1,75 Vì: N F 31  N H 31 ; N F 32  N H 32 ; Nên N F 31  N FO ; N F 32  N FO => K FL32  K FL31  ; N F 31  N FO ; N F 32  N FO  F 3   F 31    F 32   414  236,57 MPa 1,75 6.4 Xác định thông số truyền bánh hành tinh kép Do vật liệu bánh ăn khớp bánh ăn khớp khác nhau, số truyền chọn trước thoả mãn yêu cầu tỷ số truyền kích thước hình học nên ta tính thiết kế truyền ăn khớp (Z11 - Z21) (Z12 - Z22) Với truyền ăn khớp (Z21 - Z31) (Z22 - Z32) tính kiểm nghiệm bền Tính truyền bánh ăn khớp ngồi Bộ truyền cấp nhanh (Z12 - Z22) Đường kính vịng lăn bánh nhỏ được xác định theo công thức sau: d w22  K d Tn K h  u  1  H 2 u. bd Trong đó: + d w22 : đường kính bánh hành tinh sau + K d : hệ số phụ thuộc vào vật liệu cặp bánh loại bánh 133 1/3 răng,tra bảng K d  77 MPa + Tn : Mômen xoắn bánh hành tinh Tn  + T12  19804  Nmm   H  : ứng suất tiếp xúc cho phép:  H   590,91 MPa + q :số bánh vệ tinh, q = + bd : hệ số, tra bảng 4.2: bd  bw d wn Với truyền đồng A, tra bảng 6.28 [2]: Chọn  bd  + u : tỷ số truyền: u  bw d wn Z12  2, 29 Z 22 Với Z22, Z12 chọn ở phần trước + K H  : hệ số kể đến phân bố không tải trọng chiều rộng vành cho bánh vệ K H   Kc  K H  1 Với : - KC : hệ số phân bố không tải trọng cho bánh vệ tinh Kc  1,1 (khi sử dụng bánh trung tâm bánh vệ tinh) - K H  : hệ số phân bố không tải trọng chiều rộng vành KH   1, với c =  bd  0.63 Ta có: K H   1,1  1,   1,3 Thay số vào công thức (2.5): 134 tinh: d w22  77 19804.1,3. 2, 29  1 590,91 2, 29.0,63  42,5 (mm) + Chiều rộng vành răng: bw12   bd12 d w12  0,63.42.5  26  mm  Chọn bw  30 Theo 6.8 [2] ta có : m bw12 30     2.5  12  15 12  15 Tra bảng 6.8 [2]: Chọn m = 2,5 + Đường kính vịng lăn bánh vệ tinh Z22 d w22  mZ 22  2,5.17  42,5  mm  + Đường kính vòng lăn bánh bao Z32 d w32  mZ32  2,5 73  182.5  mm  + Đường kính vịng lăn bánh trung tâm Z12 d w12  mZ12  2,5.39  97,5  mm  Tính truyền cấp chậm (Z11 - Z21) Tương tự với truyền (Z11 - Z21) 1/3 Ta có: K d  77 MPa u ; Tn  T11 /  55042  Nmm  ; q = 4; Z11  1, Z 21  bd 31  0,63 ;  K F   1, Như ta có: d w21  77 55042.1, 2.1,  1 8752.1, 4.0,63 Chọn bw  30, m  2,5 135  47,5  mm  + Đường kính vòng lăn bánh vệ tinh cấp chậm Z21 d w21  mZ 22  2,5.19  47,5  mm  + Đường kính vịng lăn bánh trung tâm Z11 d w11  mZ11  2,5.27  67.5  mm  + Đường kính vịng lăn bánh trung tâm Z31 d w31  mZ 22  2,5.65  162,5  mm  Code chương trình Code SPline Ta xây dựng hàm cubicspline () để nội suy: function y = cubicspline (xData, yData, x) %Ham xap xi bang da thuc bac spline %Cu phap: [yi,f] = cubicspline(xData, yData, x) n = length(xData); c = zeros(n‐1, 1); d = ones(n, 1); e = zeros(n‐1, 1); k = zeros(n, 1); c(1:n‐2) = xData(1:n‐2) ‐ xData(2:n‐1); d(2:n‐1) = 2*(xData(1:n‐2) ‐ xData(3:n)); e(2:n‐1) = xData(2:n‐1) ‐ xData(3:n); k(2:n‐1) = 6*(yData(1:n‐2) ‐ yData(2:n‐1)) /(xData(1:n‐2) ‐ xData(2:n‐1)) ‐ 6*(yData(2:n‐1) ‐ yData(3:n)) /(xData(2:n‐1) ‐ xData(3:n)); [c, d, e] = band3(c, d e); k = band3sol(c, d, e, k); i = findseg(xData, x); h = xData(i) ‐ xData(i+1); y = ((x ‐ xData(i+1))^3/h ‐ (x ‐ xData(i+1))*h)*k(i)/6.0 ‐ ((x ‐ xData(i))^3/h ‐ (x ‐ xData(i))*h)*k(i+1)/6.0 + yData(i)*(x ‐ xData(i+1))/h ‐ yData(i+1)*(x ‐ xData(i))/h; Ta có chương trình ctcubicspline.m dùng nội suy: clear all, clc x1 = 0:0.1:5; y1 = (x1+1).^2; while x = input(ʹx = ʹ); if isempty(x) fprintf(ʹKet thucʹ); break 136 end y = cubicspline(xData, yData, x) fprintf(ʹ\nʹ) end Trendline1 moothing spline: f(x) = piecewise polynomial computed from p Smoothing parameter: p = 0.9999952 Goodness of fit: SSE: 1.03e+05 R-square: 0.8158 Adjusted R-square: 0.7154 RMSE: 28.71 Code function [fitresult, gof] = createFit(v, nl) %CREATEFIT(V,NL) % Create a fit % % Data for 'trendline1' fit: % X Input : v % Y Output: nl % Output: % fitresult : a fit object representing the fit % gof : structure with goodness-of fit info % % See also FIT, CFIT, SFIT % Auto-generated by MATLAB on 10-Feb-2020 22:01:24 %% Fit: 'trendline1' [xData, yData] = prepareCurveData( v, nl ); % Set up fittype and options ft = fittype( 'smoothingspline' ); opts = fitoptions( 'Method', 'SmoothingSpline' ); opts.SmoothingParam = 0.9999951973847; % Fit model to data [fitresult, gof] = fit( xData, yData, ft, opts ); % Plot fit with data figure( 'Name', 'trendline1' ); h = plot( fitresult, xData, yData ); legend( h, 'nl vs v', 'trendline1', 'Location', 'NorthEast' ); % Label axes xlabel v ylabel nl 137 grid on ################################################################### Trendline2 Smoothing spline: f(x) = piecewise polynomial computed from p Smoothing parameter: p = 0.99988886 Goodness of fit: SSE: 1.511e+05 R-square: 0.7161 Adjusted R-square: 0.5666 RMSE: 23.75 Code function [fitresult, gof] = createFit(v2, nl2) %CREATEFIT(V2,NL2) % Create a fit % Data for 'trendline2' fit: % X Input : v2 % Y Output: nl2 % Output: % fitresult : a fit object representing the fit % gof : structure with goodness-of fit info % % See also FIT, CFIT, SFIT % Auto-generated by MATLAB on 10-Feb-2020 22:07:34 %% Fit: 'trendline2' [xData, yData] = prepareCurveData( v2, nl2 ); % Set up fittype and options ft = fittype( 'smoothingspline' ); opts = fitoptions( 'Method', 'SmoothingSpline' ); opts.SmoothingParam = 0.999888858650065; % Fit model to data [fitresult, gof] = fit( xData, yData, ft, opts ); % Plot fit with data figure( 'Name', 'trendline2' ); h = plot( fitresult, xData, yData ); legend( h, 'nl2 vs v2', 'trendline2', 'Location', 'NorthEast' ); % Label axes xlabel v2 ylabel nl2 grid on ################################################################### 138 Trendline3 Smoothing spline: f(x) = piecewise polynomial computed from p Smoothing parameter: p = 0.9999997 Goodness of fit: SSE: 1.439e+06 R-square: 0.6752 Adjusted R-square: 0.3238 RMSE: 69.4 function [fitresult, gof] = createFit(v3, nl3) %CREATEFIT(V3,NL3) % Create a fit % Data for 'trendline3' fit: % X Input : v3 % Y Output: nl3 % Output: % fitresult : a fit object representing the fit % gof : structure with goodness-of fit info % % See also FIT, CFIT, SFIT % Auto-generated by MATLAB on 10-Feb-2020 22:14:21 %% Fit: 'trendline3' [xData, yData] = prepareCurveData( v3, nl3 ); % Set up fittype and options ft = fittype( 'smoothingspline' ); opts = fitoptions( 'Method', 'SmoothingSpline' ); opts.SmoothingParam = 0.999999695868438; % Fit model to data [fitresult, gof] = fit( xData, yData, ft, opts ); % Plot fit with data figure( 'Name', 'trendline3' ); h = plot( fitresult, xData, yData ); legend( h, 'nl3 vs v3', 'trendline3', 'Location', 'NorthEast' ); % Label axes xlabel v3 ylabel nl3 grid on ########################################################### Smoothing spline: f(x) = piecewise polynomial computed from p 139 Smoothing parameter: p = 0.99996173 Goodness of fit: SSE: 2.815e+06 R-square: 0.511 Adjusted R-square: 0.3288 RMSE: 55.83 Code function [fitresult, gof] = createFit(v4, nl4) %CREATEFIT(V4,NL4) % Create a fit % Data for 'trendline4nl' fit: % X Input : v4 % Y Output: nl4 % Output: % fitresult : a fit object representing the fit % gof : structure with goodness-of fit info % See also FIT, CFIT, SFIT % Auto-generated by MATLAB on 10-Feb-2020 22:21:35 %% Fit: 'trendline4nl' [xData, yData] = prepareCurveData( v4, nl4 ); % Set up fittype and options ft = fittype( 'smoothingspline' ); opts = fitoptions( 'Method', 'SmoothingSpline' ); opts.SmoothingParam = 0.999961729799589; % Fit model to data [fitresult, gof] = fit( xData, yData, ft, opts ); % Plot fit with data figure( 'Name', 'trendline4nl' ); h = plot( fitresult, xData, yData ); legend( h, 'nl4 vs v4', 'trendline4nl', 'Location', 'NorthEast' ); % Label axes xlabel v4 ylabel nl4 140