Advertisement
Guest User

Untitled

a guest
Mar 27th, 2017
46
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
R 2.06 KB | None | 0 0
  1. library(dplyr)
  2. library(downloader)
  3. library(rafalib)
  4.  
  5. #install.packages("zoo")
  6. library(zoo)
  7.  
  8. #install.packages("ggplot2")
  9. library(ggplot2)
  10. soldata1 = read.csv("010615.txt", sep="\t")
  11. soldata2 = read.csv("020615.txt", sep="\t")
  12. soldata3 = read.csv("030615.txt", sep="\t")
  13. soldata4 = read.csv("040615.txt", sep="\t")
  14. soldata5 = read.csv("050615.txt", sep="\t")
  15. soldata6 = read.csv("060615.txt", sep="\t")
  16. soldata7 = read.csv("070615.txt", sep="\t")
  17. #usoldata1 = select(soldata1, Tamb) %>% unlist
  18. #usoldata2 = select(soldata2, Tamb) %>% unlist
  19. #usoldata3 = select(soldata3, Tamb) %>% unlist
  20. #usoldata4 = select(soldata4, Tamb) %>% unlist
  21. #usoldata5 = select(soldata5, Tamb) %>% unlist
  22. #usoldata6 = select(soldata6, Tamb) %>% unlist
  23. #usoldata7 = select(soldata7, Tamb) %>% unlist
  24. dag1tamb = filter (soldata1, DateAndTime>10615090000, DateAndTime<10615200000) %>% select(Tamb) %>% unlist
  25. dag1pmax = filter (soldata1, DateAndTime>10615090000, DateAndTime<10615200000) %>% select(Pmax.a.Si) %>% unlist
  26. dag1back = filter (soldata1, DateAndTime>10615090000, DateAndTime<10615200000) %>% select(T.a.Si) %>% unlist
  27. dag1sol = filter (soldata1, DateAndTime>10615090000, DateAndTime<10615200000) %>% select(SolData) %>% unlist
  28.  
  29. for(i in dag1tamb){
  30.   dag1tamb[i]
  31. }
  32.  
  33. #plot(dag1sol, type="o", col="brown")
  34. plot(dag1pmax, type="o", col="blue")
  35. plot(dag1back, type="o", col="yellow")
  36. #lines(dag1tamb, type="o", col="red")
  37. tempPV = dag1tamb+ (64/1.01)*dag1sol
  38.  
  39. Enkle modeller for PV-temperaturen avhenger av solinnstrĂ¥ling og ambienttemperatur:
  40.   Temp(PV) = Temp(Ambient) + konstant*SolinnstrĂ¥ling
  41. = Solar irradiance incident on the module (POA) (W/m^2)
  42.  
  43. #lines(dag1back, type="o", col="yellow")
  44. title(main="Tamb")
  45. View(soldata1)
  46. #View(soldata2)
  47. #View(soldata3)
  48. #View(soldata4)
  49. #View(soldata5)
  50. #View(soldata6)
  51. #View(soldata7)
  52. tamb1 = select(soldata1, Tamb) %>% unlist
  53. View(tamb1)
  54. summary(tamb1)
  55. mean(tamb1)
  56. max(tamb1)
  57. popsd(tamb1)
  58. qqnorm(tamb1)
  59. qqline(0,1)
  60. pmax1 = select(soldata1, Pmax.a.Si) %>% unlist  #bytter ut  space og - med punktum
  61. mean(pmax1)
  62. max(pmax1)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement