tinytemp

Unnamed repository; edit this file 'description' to name the repository.
Log | Files | Refs

tinytemp.r (2017B)

      1 #d0T0 =-(1A*1h01*(0T0-1T1)                    )/1K0
      2 #d1T1 =-(1A*1h12*(1T1-1T2) - 1A*1h01*(1T1-0T0))/1K1
      3 #d1T2 =-(1A*1h23*(1T2-1T3) + 1A*1h12*(1T2-1T1))/1K2
      4 
      5 K = 10000
      6 ref = c('01','02','03','04')
      7 A = c(100, 100, 100,100)
      8 U = c(5.0,5.0,5.0,5.0)
      9 K = c(1000,1000,1000,1000)
     10 
     11 df <- data.frame(ref, A, U, K)
     12 
     13 print(df)
     14 
     15 
     16 y = c(rep(20, nrow(df)*5+2))
     17 y = rbind(y)
     18 vars = vector()
     19 
     20 vars = append(vars, "NTN")
     21 
     22 # initial node temps
     23 for(j in 1:5) {
     24     for(i in 1:nrow(df)) { 
     25         vars=append(vars,paste0("S",i,"T",j))
     26 
     27     }
     28 }
     29 
     30 vars = append(vars, "XTX")
     31 
     32 colnames(y) = vars
     33 
     34 #print("f")
     35 print(y)
     36 ls(all.names = TRUE)
     37 #print(state[1:nrow(df)])
     38 
     39 
     40 #XTX <- read.csv("iso.csv", header=FALSE,skip = 5)
     41 #XTX <- approxfun(
     42 #	x = XTX$V1,
     43 #	y = XTX$V6, 
     44 #	method = "linear", 
     45 #	rule = 2)
     46 
     47 #NTN=20
     48 #XTX=20
     49 Xh=5
     50 Nh=5
     51 
     52 
     53 dNTN = (-(sum(df$A,Nh,y[,'NTN'])-sum(df$A,Nh,y[2:nrow(df)]))/K)
     54 
     55 dS1T1 = (-df$A[1]*((1/df$U[1]/6)*(y[,()]
     56 assign(paste0("dS",i,"T1"), (-df$A[i]*((1/df$U[i])/6)*((get("S1T1")-get("S1T2")) +(Nh*(get("S1T1"))-get("NTN"))/(df$K[1]/3))))
     57 
     58 #for(i in 1:nrow(df)) {
     59 #        assign(paste0("dS",i,"T2"), (-df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T2"))-get(paste0("S",i,"T3"))) +
     60 #        (df$A[i]*((1/df$U[i])/6)*(get(paste0("S",i,"T2"))-get(paste0("S",i,"T1")))/(df$K[i]/3))))
     61 #}
     62 
     63 #for(i in 1:nrow(df)) {
     64 #        assign(paste0("dS",i,"T3"), (-df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T3"))-get(paste0("S",i,"T4"))) +
     65 #        (df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T3"))-get(paste0("S",i,"T2")))/(df$K[i]/3))))
     66 #}
     67 
     68 #for(i in 1:nrow(df)) {
     69 #        assign(paste0("dS",i,"T4"), (-df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T4"))-get(paste0("S",i,"T5"))) +
     70 #        (df$A[i]*((1/df$U[i])/6)*(get(paste0("S",i,"T4"))-get(paste0("S",i,"T3")))/(df$K[i]/3))))
     71 #}
     72 
     73 #for(i in 1:nrow(df)) {
     74 #        assign(paste0("dS",i,"T5"), (-df$A[i]*(Xh*(get(paste0("S",i,"T5"))-get(paste0("XTX"))) +
     75 #        (df$A[i]*((1/df$U[i])/6)*(get(paste0("S",i,"T5"))-get(paste0("S",i,"T4")))))))
     76 #}
     77 
     78 
     79 
     80 
     81 #ls(all.names = TRUE)
     82 #print(dS1T2)