Hướng dẫn sử dụng R

linear_algebra_calculus_r

Tài liệu ngắn dưới đây tôi và nhóm bạn sinh viên viết phục vụ cho mục đích học tập.

Hướng dẫn sử dụng R (cho học phần Toán cao cấp – Phiên bản 12/2017):

Huong dan su dung R

MỤC LỤC

Chủ đề 0. Download và cài đặt

Chủ đề 1. Tính toán trên trường số thực

Chủ đề 2. Ma trận và định thức

Chủ đề 3. Hệ phương trình tuyến tính

Chủ đề 4. Dạng toàn phương

Chủ đề 5. Đồ thị của hàm số

Chủ đề 6. Giới hạn

Chủ đề 7. Đạo hàm

Chủ đề 8. Giá trị lớn nhất – nhỏ nhất

Chủ đề 9. Tích phân hàm một biến

Chủ đề 10. Phương trình vi phân

Chủ đề 11. Phương trình sai phân



Lưu ý:

1. Các bạn SV gặp vấn đề khi làm Bài tập nhóm thì chịu khó tìm câu trả lời ở phần comment. Tôi sẽ KHÔNG trả lời những câu hỏi mà đã trả lời trước ở phần comment.

2. Tôi không giải đáp hộ các bạn những lỗi kỹ thuật (những lỗi do bạn gõ lệnh sai). Tôi chỉ giải đáp những vấn đề mang tính nguyên lý. Nhóm làm không được thì không có điểm.

3. Để làm được, các bạn cần đọc kỹ + chạy thử các VD trong Tài liệu hướng dẫn, kiểm tra kỹ câu lệnh mình viết có đúng như hướng dẫn không.

4. Nếu gặp lỗi, xem kỹ xem có phải lỗi do TXĐ (hàm số không xác định) [thường hiện chữ NaNs]. Nếu là lỗi do TXĐ thì chỉ cần copy câu lệnh + giải thích vì sao hàm số không xác định.

5. Khi các bạn đặt câu hỏi, phải copy đầy đủ các câu lệnh thì tôi mới giải đáp được (cần bỏ dấu > ở đầu câu lệnh đi).

————–&&————-

 

21 phản hồi

  1. thay oi.trong R thay huong dan có lẹnh
    > ketqua <- ode(y = dieukien, func = ham,times = 2:20, parms = 0, method = "iteration")
    ở phần sai phân í, nhung may luôn báo là không có chứ năng "ode"
    Error: could not find function "ode"
    em khong hiểu lắm thầy xem hộ em với

  2. thầy ơi .cho em hỏi về cách nhân 2 ma trận cái ạ.
    em khai báo được rồi nhưng khi nhân xong nó lai hiện
    Error: unexpected ‘>’ in “>”
    > >B*C
    thầy giải thích cho em cái ạ.
    em cảm ơn thầy nhiều.

  3. thưa thầy lệnh cuối nghĩa là gì ạ
    x y plot(y,type=’1′)
    Error in plot.xy(xy, type, …) : invalid plot type ‘1’

  4. install.packages(“deSolve”)
    — Please select a CRAN mirror for use in this session —
    Warning: failed to download mirrors file (cannot open URL ‘https://cran.r-project.org/CRAN_mirrors.csv’); using local file ‘E:/DATFIL~1/R-34~1.0/doc/CRAN_mirrors.csv’
    Warning: unable to access index for repository https://cran.rediris.es/src/contrib:
    cannot open URL ‘https://cran.rediris.es/src/contrib/PACKAGES’
    Warning: unable to access index for repository https://cran.rediris.es/bin/windows/contrib/3.4:
    cannot open URL ‘https://cran.rediris.es/bin/windows/contrib/3.4/PACKAGES’
    Warning messages:
    1: In download.file(url, destfile = f, quiet = TRUE) :
    InternetOpenUrl failed: ‘An error occurred in the secure channel support’
    2: package ‘deSolve’ is not available (for R version 3.4.0)

    e cài “packega desolve” như hướng dẫn mà chả đc, nó lại hiện ra cái dãy như trên là sao ạ?

  5. DLSODA- Warning..Internal T (=R1) and H (=R2) are
    such that in the machine, T + H = T on the next step
    (H = step size). Solver will continue anyway.
    In above message, R1 = 10, R2 = 8.49671e-21

    DLSODA- Warning..Internal T (=R1) and H (=R2) are
    such that in the machine, T + H = T on the next step
    (H = step size). Solver will continue anyway.
    In above message, R1 = 10, R2 = 8.49671e-21

  6. library(deSolve)
    ham <- function (n, x,parms) {list(c((9^n)*x + 3^(x^2+2*n) ))}
    dieukien <- c(y=2)
    nghiem <- ode(y = dieukien, func = ham,times = 8:15, parms = 0, method = "iteration")
    plot(nghiem, type = "1", which = "y",lwd = 2, xlab = "Truc n", ylab = "Truc x",main = "PTSP")
    "Error in plot.window(…) : need finite 'ylim' values" em thưa thầy, lỗi này là sao ạ, em làm về phần chương 10 ptsp ạ

    • Bạn làm đúng rồi. Chỉ là giá trị y quá lớn nên ko vẽ đc hình (R kí hiệu là Inf, bạn gõ nghiem rồi enter là hiểu ngay)
      Bạn copy câu lệnh + thông báo lỗi là ok.

  7. PTSP trong sách là hàm x, nhưng trong phần hướng dẫn của thầy là hàm y. bài tập thầy cho cũng là điều kiện của y. Vậy e cho x là y hay như thế nào ạ?

  8. Em chào thầy, bài của em bị lỗi như sau
    > B det(B)
    [1] 0
    > solve(B)
    Error in solve.default(B) :
    Lapack routine dgesv: system is exactly singular: U[2,2] = 0
    Thầy xem giúp em với ạ!

  9. thầy ơi cho em hỏi lại ạ [,1] [,2] [,3] [,4]
    [1,] 1 2 1 1
    [2,] 2 3 2 6
    [3,] -783 -1 0 31
    [4,] 1 1 1 5
    det(B)
    [1] 1.738609e-12
    solve(B)
    Error in solve.default(B) :
    system is computationally singular: reciprocal condition number = 2.00541e-19
    bài trên của bạn là định thức bằng k còn của e định thức k bằng 0 thầy giúp em vs ạ

  10. Trường hợp này det(B) = 1.738609e-12 ~ 0 thì chứng tỏ không tồn tại ma trận nghịch đảo. Bạn copy câu lệnh + thông báo lỗi + giải thích lỗi là vì như thế.
    [Kết quả chính xác phải là 0, nhưng R tính toán ko tốt !!!!, có cách hiệu chính nhưng thôi cứ tạm bỏ qua].

  11. Các bạn đọc các lưu ý ở trên rồi hãy đặt câu hỏi.

Trả lời

Mời bạn điền thông tin vào ô dưới đây hoặc kích vào một biểu tượng để đăng nhập:

WordPress.com Logo

Bạn đang bình luận bằng tài khoản WordPress.com Đăng xuất / Thay đổi )

Twitter picture

Bạn đang bình luận bằng tài khoản Twitter Đăng xuất / Thay đổi )

Facebook photo

Bạn đang bình luận bằng tài khoản Facebook Đăng xuất / Thay đổi )

Google+ photo

Bạn đang bình luận bằng tài khoản Google+ Đăng xuất / Thay đổi )

Connecting to %s

%d bloggers like this: