Hướng dẫn sử dụng phần mềm r

     

“Giải thưởng ứng dụng mã mối cung cấp msinh hoạt giỏi nhất: Bossie Awards 2010″ R là phần mềm phân tích tài liệu được kiến tạo bởi vì Ross Ihaka cùng Robert Gentleman tại The University of Auckl&, New Zealand, tiếp tục được trở nên tân tiến bởi nhóm R Development Vi xử lý Core Team.

Bạn đang xem: Hướng dẫn sử dụng phần mềm r

R là một trong những phần mượt hoàn toàn miễn giá thành. Tuy miễn phí, mà lại tính năng của R không chiến bại kém các ứng dụng tmùi hương mại. Tất cả dẫu vậy cách thức, mô hình mà lại các ứng dụng thương thơm mại hoàn toàn có thể làm cho được thì R cũng có thể làm cho được....


*

Hướng dẫn sử dụngPHẦN MỀM R (Cho học tập phần Tân oán cao cấp) LỜI NÓI ĐẦU “Trao Giải phần mềm mã mối cung cấp msống xuất sắc nhất: Bossie Awards 2010″R là phần mềm đối chiếu dữ liệu được thiết kế bởi Ross Ihaka cùng Robert Gentleman trên TheUniversity of Auckl&, New Zeal&, liên tục được phát triển do nhóm R Development CoreTeam. R là một trong những phần mượt trọn vẹn miễn mức giá. Tuy miễn phí tổn, tuy nhiên công dụng của R khônglose kém nhẹm các ứng dụng tmùi hương mại. Tất cả nhưng mà cách thức, quy mô nhưng những phần mềmtmùi hương mại có thể làm được thì R cũng hoàn toàn có thể làm được. R có ích cố kỉnh là khả năng phân tích biểuvật dụng tuyệt vời và hoàn hảo nhất. Không một trong những phần mượt như thế nào rất có thể sánh cùng với R về phần biểu đồ!Tuy nhiên, R có mẫu bất lợi là sử dụng lệnh (y hệt như Stata hay SAS) chứ không hề sử dụng “menu”nlỗi SPSS. Như vậy tức là đối với fan “lười biếng” thì vẫn thấy R bất tiện. Nhưng vớingười mong muốn cố mang đầy đủ cơ chế căn phiên bản của R, thì vẫn yêu thích ngôn ngữ này ngay lập tức. Đối vớinhững người tốt vi tính với tân oán, R là phần mềm lí tưởng.Tài liệu nđính thêm này Shop chúng tôi viết dựa vào mục Help của R nhằm giúp các bạn sinh viên làm cho quencùng với câu hỏi áp dụng 1 phần mềm toán thù học tập để giải các bài xích toán của học tập phần Toán thù cao cấp. LÊ VĂN TUẤN – VŨ VĂN SONG NGUYỄN THỊ HOA - NGUYỄN THỊ HUYỀN11 GV & SV Khoa Tin học tập Thương mại - Đại học Tmùi hương mại MỤC LỤCChủ đề 0. Download cùng thiết lập đặtChủ đề 1. Tính toán trên trường số thựcChủ đề 2. Ma trận với định thứcChủ đề 3. Giải hệ phương trình con đường tínhChủ đề 4. Vẽ vật dụng thịChủ đề 5. Đạo hàmChủ đề 6. Tích phân hàm một biếnChủ đề 7. Phương thơm trình vi phânChủ đề 8. Phương thơm trình không nên phân Chủ đề 0. Download với tải đặtBạn truy vấn vào trang chủ: http://www.r-project.org/ (bối cảnh nlỗi hình sinh sống dưới), click vào tải về Rbạn mang đến trang CRAN Mirrors, clichồng một link làm sao đó, các bạn sẽ đến trang The Comprehensive R ArchiveNetwork, clichồng vào Download R for Windows, cliông xã tiếp install R for the first time, click tiếp DownloadR 2.15.0 for Windows đang tải về được tệp tin.R-2.15.0-win.exe pháo (mon 4/2012), cài đặt nlỗi các phầnmượt không giống.Sau lúc thiết đặt, Shoutcut nhằm chạy phần mềm đang xuất hiện trên Desktop, bạn click để chạy phầnmềm. Cửa sổ lệnh của phần mềm vẫn nhỏng hình bên dưới, trên lốt nhắc “>” chúng ta cũng có thể gõ câu lệnh và thừa nhận Enter(↵) để trải đời ứng dụng triển khai câu lệnh.Ghi chú:Quý khách hàng hoàn toàn có thể vào mục Help bên trên menu nhằm áp dụng những lí giải của R. Chủ đề 1. Tính toán thù trên ngôi trường số thựcCác phép tân oán bên trên trường số thực là: cùng (+), trừ (-), nhân (*), chia (/), lũy vượt (^)Các hàm thông dụng:(pi trình diễn số )1. Tính 7/3.5> 7/3.5<1> 22. Tính a= (4^5-1/6) ( √ + )> a=(4^5-1/6)*(exp(1/3)+pi);<1> 4645.3423. Tính log3 (4)> log(4)/log(3) (Ta sử dụng bí quyết thay đổi cơ số)<1> 1.261864. Tính arcsin(1/2)> asin(1/2)<1> 0.52359885. Cho f(x)= (sin(x) +x2)/(ex+1), tính f(π/6)> (sin(pi/6)+(pi/6)^2)/(exp(pi/6)+1)<1> 0.2879945 Chủ đề 2. Ma trận với định thức1.

Xem thêm: Diễm My 9X Vô Tình Để Lộ Nội Y, Diễm My 9X Mặc Áo Xuyên Thấu Lộ Nội Y Dự Sự Kiện

Knhì báo thay đổi ma trậnVD: Knhị báo ma trận cỡ 1x3 (vec tơ dòng):> y A y BBD=B+C ↵ (ở chỗ này ta sẽ tạo ra thêm đổi thay D = B+C )VD: >B*C; >B^10; >5*B3. Phnghiền nhân nhị ma trậnVD: >B%*%C4. Ma trận gửi vịVD: >t(B)5. Tìm hạng của ma trậnVD: >qr(B)$rank6. Tìm ma trận nghịch đảoVD: > solve(B)7. Tính định thức (của ma trận vuông)VD: >det(B) Chủ đề 3. Giải hệ pmùi hương trình tuyến đường tính + + =6 − = −1; VD1: Tìm nghiệm riêng biệt của hệ PTTT: + +2 = 9 Ta tiến hành nhỏng sau: Trên màn hình sẽ mở ra kết quả là: > A y lm.fit(A,y)$coefficientsGhi chú:Trong ví dụ này ta được nghiệm nhất, x=(1, 2, 3) +− =0 3 − =3VD2: Tìm nghiệm riêng của hệ PTTT: Ta thực hiện như sau: Trên màn hình vẫn xuất hiện thêm tác dụng là: > A y lm.fit(A,y)$coefficientsGhi chú: Trong ví dụ này nghiệm của hệ có một tsi số, nghiệm riêng là:x1 = 1 ; x2 = -1 ; x3 = 0VD3: Tìm nghiệm riêng rẽ của hệ PTTT + − =3 Ta triển khai nhỏng sau: Trên màn hình hiển thị vẫn lộ diện kết quả là: > A y lm.fit(A,y)$coefficientsGhi chú: Trong ví dụ này nghiệm của hệ tất cả 2 tsi số, nghiệm riêng rẽ là: x1 = 3; x2 = 0; x3 = 0 Chủ đề 4. Vẽ đồ gia dụng thị1. Hàm 1 phát triển thành (2D)VD: Vẽ đồ thị hàm số: y = x2 + 1 trên đoạn <-10; 10>Thực hiện nay nlỗi sau:> x y plot(y,type="l")Phần mượt vẫn xuất ra đồ vật thị (trên 1 cửa sổ khác – R Graphic) nlỗi hình bên dưới (bạn có thể vào File -> Copylớn the clipboard (CTRL+C) cùng paste vào word)Ghi chú: Biến x là một dãy số, đổi thay y cũng là một hàng số, chúng ta gõ demo x; rồi gõ y.Ttê mê khảo Chủ đề 1 Lúc phải vẽ những hàm phức hợp.2. Hàm 2 đổi thay (3D)VD: Vẽ vật dụng thị hàm z=sin(x)*y với x ∈ <0; 2π>, y ∈ <0; 5>Thực hiện tại nlỗi sau:> x> y f Chủ đề 5. Đạo hàm1. Hàm một biếnVD 1: Tính đạo hàm cung cấp 1 và cấp cho 2 của hàm số f(x) = x3 + 1Ta triển khai nlỗi sau:> D(expression(x^3+1), "x")3 * x^2> D(D(expression(x^3+1), "x"),"x")3 * (2 * x)VD 2: Tính đạo hàm của hàm số f(x) = x3 + 1 trên x = 2.Ta tiến hành nhỏng sau:> D(expression(x^3+1), "x")3 * x^2> x 3 * x^2<1> 12Ghi chú: Tđam mê khảo Chủ đề 1 lúc bắt buộc tính đạo hàm của những hàm phức tạp.2. Hàm nhiều biếnVD: Tính các đạo hàm riêng rẽ cấp cho 1 cùng cấp cho 2 của hàm số z = x2 + y2> D(expression(x^2+y^2), "x")2*x> D(expression(x^2+y^2), "y")2*y> D(D(expression(x^2+y^2), "x"),"x")<1> 2 Chủ đề 6. Tích phân hàm một biến1. Tích phân thông thườngVD: TínhTa thực hiện nlỗi sau:> f integrate(f, lower=1, upper=10)198238.2 with absolute error f integrate(f, lower=1, upper=Inf)1 with absolute error Chủ đề 7. Phương thơm trình vi phânĐể giải PTVP, trước hết ta nên thiết đặt package deSolve sầu bằng phương pháp thực hiện lệnh trên R (laptop phảiđang liên kết internet):> install.packages(“deSolve”)Phần mềm đang hiện 1 bảng CRAN mirror, bạn chọn 1 loại làm sao đó rồi OKE.Từ số đông lần chạy R sau (để giải PTVP) chúng ta không hẳn setup nữa, tuy thế cần gọi package deSolve(không phải liên kết internet) bằng lệnh:> library(deSolve)1. Phương thơm trình vi phân cấp cho 1Ghi chú: Xét pmùi hương trình vi phân cấp 1: dy/dx = f(x,y).Giả sử ta đề nghị giải PTVPhường dy/dx = x2. Pmùi hương trình này còn có nghiệm bao quát là y = x3/3 + C. Với điều kiênlúc đầu y(0) = 1 ta được nghiệm riêng biệt y = x3/3 + 1. Ta đã vẽ vật dụng thị nghiệm riêng biệt này.(Phần mềm R lưu lại hàm số dưới dạng bảng, đề xuất ta không có được công thức tường minch tuy vậy vẫn vẽđược đồ gia dụng thị).VD1: Vẽ đồ gia dụng thị nghiệm riêng biệt của phương thơm trình vi phân dy/dx = x2 với điều kiện lúc đầu y(0) = 5 trênmiền <0, 100>.Ta thực hiện nlỗi sau:> mê say dieukien ketqua plot(ketqua, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")Ghi chú: times = 0:100 bộc lộ miền nghiệm là <0, 100>, do trong PTVP trở nên x hay là phát triển thành thời gianphải tất cả quy ước chuẩn chỉnh là times.Phần mềm đang xuất ra vật dụng thị (trên 1 cửa sổ không giống – R Graphic) nhỏng hình bên trên (bạn có thể vào File -> Copyto the clipboard (CTRL+C) với paste vào word)VD2: Vẽ đồ vật thị nghiệm riêng của phương trình vi phânvới điều kiện thuở đầu y(4) = 2 bên trên miền <4, 10>.Ta triển khai nlỗi sau:> mê man dieukien nghiem plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc x", ylab = "Truc y",main = "PTVP")Ghi chú: Tsay đắm khảo Chủ đề 1 khi phải thao tác làm việc cùng với các hàm tinh vi.2. Phương thơm trình vi phân cấp 2VD1: Vẽ vật dụng thị nghiệm riêng của phương trình vi phân: y” – (1 – y2)y’ + y = 0, với ĐK ban đầuy(0) = 2, y’(0) = 0 bên trên miền <0, 100>.Đặt y<1> = y, y<2> = y’, ta chuyển PTVP về hệ:y<1>’ = y<2> + 0.y<1>y<2>’ = (1 – y<1>2)y<2> + y<1>(ĐK ban đầu, y1 = 2, y2 = 0).Ta triển khai nhỏng sau:> say đắm dieukien nghiem plot(nghiem, type = "l", which = "y1",lwd = 2, ylab = "Truc y",main = "PTVPhường cấp 2")VD2: Vẽ đồ dùng thị nghiệm riêng của phương thơm trình vi phân: y” – (1 – y2)y’ + y = exsinx, với điều kiện banđầu y(0) = 2, y’(0) = 0 bên trên miền <0, 100>.Ta triển khai như sau:> yêu thích dieukien nghiem plot(nghiem, type = "l", which = "y1",lwd = 2, ylab = "Truc y",main = "PTVP.. cấp 2") Chủ đề 8. Phương thơm trình không nên phânĐể giải PTSPhường, trước nhất ta bắt buộc thiết đặt package deSolve bằng phương pháp tiến hành lệnh trên R (máy vi tính phảivẫn kết nối internet):> install.packages(“deSolve”)Phần mượt vẫn hiện 1 bảng CRAN mirror, các bạn chọn một chiếc làm sao kia rồi OKE.Từ đa số lần chạy R sau (để giải PTSP) bạn chưa phải setup nữa, nhưng lại yêu cầu Hotline package deSolve(ko nên liên kết internet) bằng lệnh:> library(deSolve)Nếu phương trình tất cả đựng giai vượt (!) đề xuất cài đặt thêm package pracma với Gọi package này trước lúc giải.1. Phương thơm trình sai phân cung cấp 1Ghi chú: Xét phương trình không nên phân cấp cho 1: y(n+1) = f(n,y(n)).Giả sử ta buộc phải giải PTSPhường. y(n+1) + 2y(n) = 0. Phương thơm trình này còn có nghiệm tổng quát là y(n) = C(-2)n. Vớiđiều kiên lúc đầu y(2) = 3 ta được nghiệm riêng y(n) = ¾(-2)n. Ta vẫn vẽ thiết bị thị nghiệm riêng này.(Phần mềm R giữ hàm số bên dưới dạng bảng, đề nghị ta không có được phương pháp tường minh tuy thế vẫn vẽđược đồ gia dụng thị).VD1: Vẽ thứ thị nghiệm riêng của phương trình sai phân y(n+1) + 2y(n) = 0 với ĐK thuở đầu y(2) =3 bên trên miền <2, 20>.Ta tiến hành như sau:> đê mê dieukien ketqua plot(ketqua, type = "l", which = "y",lwd = 2, xlab = "Truc n", ylab = "Truc y",main = "PTSP")Để hiển thị các quý giá của hàm y, chúng ta thực hiện lệnh:> ketquaGhi chú: times = 2:20 biểu thị miền nghiệm là <2, 20>, vì vào PTSPhường. biến hóa n thường xuyên là biến thời giancần bao gồm quy ước chuẩn là times.Phần mượt vẫn xuất ra vật thị (trên 1 hành lang cửa số không giống – R Graphic) như hình trên (bạn có thể vào File -> Copylớn the clipboard (CTRL+C) và paste vào word)Giải PTSPhường. chỉ khác PTVP ở: method = "iteration"VD2: Vẽ đồ gia dụng thị nghiệm riêng biệt của phương thơm trình sai phân y(n+1) = (n+1)y(n) + (n+1)!n; với điều kiện banđầu y(4) = 2 trên miền <4, 10>.Ta tiến hành nlỗi sau:#Cần yêu cầu hotline package pracma trước> mê mệt dieukien nghiem plot(nghiem, type = "l", which = "y",lwd = 2, xlab = "Truc n", ylab = "Truc y",main = "PTSP")Ghi chú: Tmê mẩn khảo Chủ đề 1 Khi buộc phải thao tác làm việc cùng với các hàm phức tạp.2. Phương thơm trình không nên phân cấp 2VD1: Vẽ thứ thị nghiệm riêng rẽ của pmùi hương trình không đúng phân y(n+2) - 5y(n+1) + 6y(n) = 0, cùng với ĐK banđầu y(0) = 2, y(1) = 5 trên miền <0, 10>.ta gửi PTSPhường. về hệ:y(n+1) = y(n+1) + 0*y(n)y(n+2) = 5*y(n+1) - 6*y(n)Đặt y<1> = y(n), y<2> = y(n+1), ĐK ban sơ, y1 = 2, y2 = 5.Ta thực hiện như sau:> tê mê dieukien > nghiem plot(nghiem, type = "l", which = "y1",lwd = 2, ylab = "Truc y",main = "PTVPhường. cấp cho 2")VD2: Vẽ đồ gia dụng thị nghiệm riêng rẽ của pmùi hương trình vi phân: y(n+2) - 5y(n+1) + 6y(n) = n2 + 2n + 3, với điềukhiếu nại ban đầu y(0) = 2, y(1) = 5 trên miền <0, 10>.Ta tiến hành như sau:> yêu thích dieukien nghiem plot(nghiem, type = "l", which = "y1",lwd = 2, ylab = "Truc y",main = "PTVPhường cấp 2") ----------&&---------