SHARE
TWEET

Untitled

a guest Oct 18th, 2019 101 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. library(wesanderson)
  2. library(ggplot2)
  3. library(RColorBrewer)
  4. library(quantmod)
  5.  
  6. # Daten werden eingelesen
  7.  
  8. Hendel1 <- read.table("Hendel1.txt", header = TRUE, sep = "")
  9. #dBergmann <- read.table("6Bergmann.txt", header = TRUE, sep = "")
  10. #dAnderes <- read.table("Anderes.txt", header = TRUE, sep = "")
  11.  
  12.  
  13. # Spalten 1 der Wellenlänge = z & Spalten der Intensität = i werden ausgegeben
  14.  
  15. z01 <- Hendel1[,1]
  16. i01 <- Hendel1[,2]
  17. i02 <- Hendel1[,3]
  18. i03 <- Hendel1[,4]
  19. i04 <- Hendel1[,5]
  20. i05 <- Hendel1[,6]
  21. i06 <- Hendel1[,7]
  22. i07 <- Hendel1[,8]
  23. i08 <- Hendel1[,9]
  24. i09 <- Hendel1[,10]
  25. i10 <- Hendel1[,11]
  26. i11 <- Hendel1[,12]
  27. i12 <- Hendel1[,13]
  28. i13 <- Hendel1[,14]
  29. i14 <- Hendel1[,15]
  30. i15 <- Hendel1[,16]
  31. i16 <- Hendel1[,17]
  32. i17 <- Hendel1[,18]
  33. i18 <- Hendel1[,19]
  34. i19 <- Hendel1[,20]
  35. i20 <- Hendel1[,21]
  36. i21 <- Hendel1[,22]
  37. i22 <- Hendel1[,23]
  38. i23 <- Hendel1[,24]
  39. i24 <- Hendel1[,25]
  40. i25 <- Hendel1[,26]
  41. i26 <- Hendel1[,27]
  42. i27 <- Hendel1[,28]
  43. i28 <- Hendel1[,29]
  44. i29 <- Hendel1[,30]
  45. i30 <- Hendel1[,31]
  46. i31 <- Hendel1[,32]
  47. i32 <- Hendel1[,33]
  48. i33 <- Hendel1[,34]
  49. i34 <- Hendel1[,35]
  50. i35 <- Hendel1[,36]
  51. i36 <- Hendel1[,37]
  52.  
  53. # Farbpalette von rot (früheste) nach blau (späteste)
  54.  
  55. RdBu <- brewer.pal(n=9, name="RdBu")
  56. Greens <- brewer.pal(n=9, name="Greens")
  57. PuRd <- brewer.pal(n=9, name="PuRd")
  58. Blues <- brewer.pal(n=9, name="Blues")
  59. # plot: x-Achse & y-Achse entsteht
  60. # points: Spektren entstehen
  61.  
  62. plot(z01, i01, type = "l", lwd=1, col=RdBu[1], main = "Proband Hendel", xlab="Wellenlänge in nm", ylab = "Intensität")
  63. points(z01, i02, type = "l", lwd=1, col=RdBu[2])
  64. points(z01, i03, type = "l", lwd=1, col=RdBu[3])
  65. points(z01, i04, type = "l", lwd=1, col=RdBu[4])
  66. points(z01, i05, type = "l", lwd=1, col=RdBu[5])
  67. points(z01, i06, type = "l", lwd=1, col=RdBu[6])
  68. points(z01, i07, type = "l", lwd=1, col=RdBu[7])
  69. points(z01, i08, type = "l", lwd=1, col=RdBu[8])
  70. #points(z01, i09, type = "l", lwd=1, col=RdBu[9])
  71. #points(z01, i10, type = "l", lwd=1, col=Greens[1])
  72. #points(z01, i11, type = "l", lwd=1, col=Greens[2])
  73. #points(z01, i12, type = "l", lwd=1, col=Greens[3])
  74. #points(z01, i13, type = "l", lwd=1, col=Greens[4])
  75. #points(z01, i14, type = "l", lwd=1, col=Greens[5])
  76. #points(z01, i15, type = "l", lwd=1, col=Greens[6])
  77. #points(z01, i16, type = "l", lwd=1, col=Greens[7])
  78. #points(z01, i17, type = "l", lwd=1, col=Greens[8])
  79. #points(z01, i18, type = "l", lwd=1, col=Greens[9])
  80. #points(z01, i19, type = "l", lwd=1, col=PuRd[1])
  81. #points(z01, i20, type = "l", lwd=1, col=PuRd[2])
  82. #points(z01, i21, type = "l", lwd=1, col=PuRd[3])
  83. #points(z01, i22, type = "l", lwd=1, col=PuRd[4])
  84. #points(z01, i23, type = "l", lwd=1, col=PuRd[5])
  85. #points(z01, i24, type = "l", lwd=1, col=PuRd[6])
  86. #points(z01, i25, type = "l", lwd=1, col=PuRd[7])
  87. #points(z01, i26, type = "l", lwd=1, col=PuRd[8])
  88. #points(z01, i27, type = "l", lwd=1, col=PuRd[9])
  89. #p#oints(z01, i28, type = "l", lwd=1, col=Blues[1])
  90. #points(z01, i29, type = "l", lwd=1, col=Blues[2])
  91. #points(z01, i30, type = "l", lwd=1, col=Blues[3])
  92. #points(z01, i31, type = "l", lwd=1, col=Blues[4])
  93. #points(z01, i32, type = "l", lwd=1, col=Blues[5])
  94. #points(z01, i33, type = "l", lwd=1, col=Blues[6])
  95. points(z01, i34, type = "l", lwd=1, col=Blues[7])
  96. points(z01, i35, type = "l", lwd=1, col=Blues[8])
  97. points(z01, i36, type = "l", lwd=1, col=Blues[9])
  98. #p01 <- findPeaks(i06, thresh = 500)
  99. #points(p06/ 250, i06[p06], pch=20, col=RdBu[1], cex = 2)
  100.  
  101. #closeAllConnections()
  102. #rm(list=ls())
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!
 
Top