RST 6상 그래프 R plot 그리기

   RST 6상 그래프를 구글링해서 찾으러고 하니,  RST  3상만 나왔다.  RST 6상은 특수한 설비에서 작동되는 것이니,  나올 수가 없다.  교육용 자료로 만들려 하다 보니, 찾을 수가 없어 직접 그렸다. 

RST 3상은 360에서 3을 나누니,  120의 파형이 그래프로 보여지지만,   RST 6상은 360에서 6을 나누니 60도의 파형 그래프로 그려 졌다.   아래의 소스코드를 보고 그림을 그릴 수 있다. 

# 라이브러리 로딩  

  library(dplyr)

  library(ggplot2)

  library(tidyr)

# 주파수의 라디안 값을 정한다.  파라메터   

  x <- seq(0, 2*pi, by=pi/179.6)

  amp <- 1

  wav <- 1

# sin 함수

  sin_vector <- function(amp, wav,x){

    signal <- amp*sin(wav*x)

  

  return(signal)

  }

  

# cos 함수

  cos_vector <- function(amp, wav,x){

    

    signal <- amp*cos(wav*x)

    return(signal)

  }

  

  r <-sin_vector(amp, wav,x)

  plot(r, type='l', col=2, ylim=c(-1,1)); abline(h=0, lty=3)

  

  

  par(new = T)

  t <-  sin_vector(amp, wav, x)

  t <- c(t[120:360], t[1:119])

  plot(t, type='l', col=1, ylim=c(-1,1)); abline(h=0, lty=3)

  

  par(new = T)

  s = cos_vector(amp, wav, x)

  s <- c(s[150:360], s[1:149])  

  plot(s, type='l', col=4, ylim=c(-1,1)); abline(h=0, lty=3)

 

  

  # r 하고 60도 차이나게 그리기

  par(new = T)

  t1 = cos_vector(amp, wav, x)

  t_prime <- c(t1[210:360], t1[1:209])  

  plot(t_prime, type='l',lty=3, col=1, ylim=c(-1,1)); abline(h=0, lty=3)

  

  

  r_prime <-  -sin_vector(amp, wav, x)

  par(new = T)

  plot(r_prime, type='l', lty=3, col=2, ylim=c(-1,1)); abline(h=0, lty=3)

  

  # 마지막 하나 각도 찾기   

  par(new = T)

  r1 <-  sin_vector(amp, wav, x)

  s_prime <- c(r1[60:360], r1[1:59])

  plot(s_prime, type='l',lty=3, col=4, ylim=c(-1,1)); abline(h=0, lty=3)

  

  

RST 6상