开发工具:
文件大小: 2kb
下载次数: 0
上传时间: 2017-08-30
详细说明： 主要是用R语言，进行初期日期数据。 data1<-read.csv(file='C:/Users/li/Desktop/参数优化/a.csv',header = T) data2<-read.csv(file='C:/Users/li/Desktop/参数优化/b.csv',header = T) #提取日期数据 date1<-as.Date(data1$日期) dates<-paste(data2$年,data2$月,data2$日,sep='/') date2<-as.Date(dat es,'%Y/%m/%d') xdate1<-date1[1:28]#tree 1 xdate2<-date1[29:54]#tree 2 xdate3<-date1[55:81]#tree 3 #构造09-01的日期形式 dates<-paste(data1$年-1,'09','01',sep='/') datestan<-as.Date(dates,'%Y/%m/%d') datestan1<-datestan[1:28]#tree 1 datestan2<-datestan[29:54]#tree 2 datestan3<-datestan[55:81]#tree 3 #提取气温数据 temp<-as.numeric(data2$平均气温) ##------------------------------------------------- #C t0 tb Tb Rc=function(c,Tb){ rc=0 ci<-datn1 while(rc!=c|ci<=datn1+140){ rc<-rc+sum(temp[ci]<=Tb) ci=ci+1 } print(ci) } # # Rf=function(f,Tb){ rf=0 fi<-tb while(rf<=f|fi<=tb+200){ if(temp[fi] <= Tb){ rf<-rf+0} else {rf<-rf+temp[fi]-Tb} fi=fi+1 } print(c(fi,rf)) } # #------------------------------ ypre<-datestan1 for(i in 1:28){ dat1<-datestan1[i] datn1<-which(date2==dat1) tb<-Rc(c=100,Tb=10) result<-Rf(f=100,Tb=10) y<-result[1] rf<-result[2] ypre[i]<-dat1+y-datn1 } r2<-mean(ypre-xdate1) ...展开收缩
(系统自动生成,下载前可以参看下载内容)
下载文件列表
相关说明
- 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
- 本站是交换下载平台，提供交流渠道，下载内容来自于网络，除下载问题外，其它问题请自行百度。
- 本站已设置防盗链，请勿用迅雷、QQ旋风等多线程下载软件下载资源，下载后用WinRAR最新版进行解压.
- 如果您发现内容无法下载，请稍后再次尝试；或者到消费记录里找到下载记录反馈给我们.
- 下载后发现下载的内容跟说明不相乎，请到消费记录里找到下载记录反馈给我们，经确认后退回积分.
- 如下载前有疑问，可以通过点击"提供者"的名字，查看对方的联系方式，联系对方咨询.
