如何在R 4.0.2中实现与Matlab Datenum等效的日期序列生成？

I’ve run into similar compatibility issues with older R versions and niche packages, so let’s break down how to get the exact same datenum values as Matlab’s datenum function for the 1988 calendar year—no problematic packages required.

Method 1: Quick Offset Calculation (Simple & Fast) #

Matlab’s datenum counts days starting from January 0, 0000, while R’s base date functions use January 1, 1970 as the origin. We can use a known reference point to calculate an offset and convert R dates directly.

Step-by-Step Code: #

# 1. Generate the full 1988 date sequence in R
target_year <- 1988
r_dates <- seq(
  as.Date(paste0(target_year, "-01-01")),
  as.Date(paste0(target_year, "-12-31")),
  by = "day"
)

# 2. Calculate the offset using your sample reference: Matlab datenum(1988,1,1) = 726103
matlab_ref <- 726103
r_ref_value <- as.numeric(as.Date("1988-01-01"))
offset <- matlab_ref - r_ref_value

# 3. Convert R dates to Matlab-style datenum
datenum_vect <- as.numeric(r_dates) + offset

# Verify the first few values (matches your sample!)
head(datenum_vect)
# Output: 726103 726104 726105 726106 726107 726108

This method is perfect if you only need dates for a specific year and have a known Matlab datenum reference. It’s lightweight and works with all R versions.

Method 2: Manual Matlab Datnum Implementation (Universal & Package-Free) #

If you need to generate datenum values for any date (not just 1988), you can replicate Matlab’s underlying algorithm directly in R. This uses the astronomical formula Matlab relies on to calculate days since January 0, 0000.

Step-by-Step Code: #

# Define a function that mimics Matlab's datenum logic
matlab_datenum <- function(year, month, day) {
  # Matlab's core datenum calculation formula
  367 * year - 
    floor(7 * (year + floor((month + 9)/12))/4) - 
    floor(3 * (floor((year + floor((month - 9)/7))/100) + 1)/4) + 
    floor(275 * month / 9) + 
    day + 
    1721029
}

# Generate 1988's datenum sequence (1988 is a leap year, so 366 days)
target_year <- 1988
start_datenum <- matlab_datenum(target_year, 1, 1)
datenum_vect <- start_datenum:(start_datenum + 365)

# Check the first few values
head(datenum_vect)
# Output: 726103 726104 726105 726106 726107 726108

This function will return identical results to Matlab’s datenum(y, m, d) for any valid date, and it works flawlessly in R 4.0.2 since it uses only base R functions.

内容的提问来源于stack exchange，提问作者Sss