七月流水账
7.1-7.3
><
7.4
看了一會(huì)兒李老大給的Caffe的教程
?重配python + qt + vtk，大概是原來裝了一個(gè)xxx（忘記名字），然后位數(shù)有些不對(duì)，搞了好一陣
果然還是卸載重裝好><
然后又看了會(huì)兒Matlab神經(jīng)網(wǎng)絡(luò)工具箱怎么用
還沒看太明白就被抓去搞界面了>w<（終于稍微有點(diǎn)點(diǎn)懂怎么把代碼和GUI的東西結(jié)合起來了，實(shí)現(xiàn)還不太明白）
間隙看了線性篩
7.5
讀取dcm的標(biāo)簽信息，需要先配置好simplevtk
import SimpleITK as sitkinputImage = sitk.ReadImage(str(file_name))image_laterality = inputImage.GetMetaData('0020|0060')print image_laterality View Code需要什么信息先用ImageJ打開圖片看下對(duì)應(yīng)標(biāo)簽的編號(hào)
然后感覺是VTK這一套東西都不懂，做起事情來很茫然..
---------------------------------------------------------------------
晚上和以前一起實(shí)習(xí)的聊了聊，發(fā)現(xiàn)好多好多都沒有學(xué)阿
計(jì)算機(jī)網(wǎng)絡(luò)，操作系統(tǒng)，計(jì)算機(jī)組成原理，編譯原理算基礎(chǔ)類的過一遍
語言類的C++ primer之類的過一遍
linux XXXXX編程
機(jī)器學(xué)習(xí)的常用算法，用過還有會(huì)推導(dǎo)，用到什么框架
leetcode?
（web服務(wù)器xxxxxxxx）
7.6
準(zhǔn)備老王中 & 繼續(xù)完善軟件的功能
創(chuàng)建一個(gè)新的界面
7.7
覺得周圍的人都很厲害，自己菜的不行...
DCNN原來是由多個(gè)CNN組成的模型
區(qū)域生長
7.8
LDA?將高維樣本數(shù)據(jù)采用線性組合的方式投影到一維空間，并找到使得類內(nèi)間距最小，類間間距最大的組合系數(shù)
劃了一天水叭
7.9
難過了好久...不知道該怎么辦于是去實(shí)驗(yàn)室接著老王了（其實(shí)是蹭網(wǎng)蹭空調(diào)蹭水喝><
主體部分大概好了...但是總覺得還差很多TAT
7.10
放射組學(xué)
?------- ?昏割線 -------------------------------------
剪掉頭發(fā)重新好好當(dāng)狗了要加油喔><
-------------------------------------------------------
?
使用PyQt 編寫第一個(gè)Python GUI程序
但是教程給出的示例代碼在我的電腦會(huì)報(bào)錯(cuò)
原因是電腦上裝的Qt建一個(gè)ui的時(shí)候，會(huì)搞出一個(gè)發(fā)送者和接受者的信號(hào)，和槽，但是在代碼里面沒有寫到
所以再這樣加上兩個(gè)函數(shù)就好
1 import sys 2 from PyQt4 import QtCore, QtGui, uic 3 4 qtCreatorFile = "C:\\Users\\lab\\Desktop\\test\\haha.ui" # Enter file here. 5 6 Ui_MainWindow, QtBaseClass = uic.loadUiType(qtCreatorFile) 7 8 class MyApp(QtGui.QMainWindow, Ui_MainWindow): 9 def __init__(self): 10 QtGui.QMainWindow.__init__(self) 11 Ui_MainWindow.__init__(self) 12 self.setupUi(self) 13 14 def accept(self): 15 print "111" 16 17 def reject(self): 18 print "22" 19 20 21 if __name__ == "__main__": 22 app = QtGui.QApplication(sys.argv) 23 window = MyApp() 24 window.show() 25 sys.exit(app.exec_()) View Code?
Python調(diào)用另外的py文件中的類
?
7.11
編寫 UI 和后端的邏輯！
python中的編碼問題 ? 有點(diǎn)長先碼在這里叭
具體現(xiàn)在遇到的問題也很好解決，就是輸出中文字符串亂碼的話，在前面加一個(gè)u轉(zhuǎn)為unicode編碼就能正常顯示了
?
真的好久沒有摔過狗吃屎了哈哈哈哈哈，只是感覺手指現(xiàn)在打代碼有點(diǎn)抬不起qwq
?
RNN?今晚講課有講，但是感覺聽的很恍惚...先碼叭
?
Problem A.?Ambiguous Cipher
?按照奇偶來想( 按照樣例來想...奇數(shù)無解，偶數(shù)可以知道順二和倒二位置的字母，然后從左到右推一遍再從右到左推一遍
wa的不明所以的原因是，字符串每次輸出答案是直接printf 完整個(gè)字符串，但是可能有上一個(gè)case的在后面(所以就打印到長度就好>o<好傻喔
1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 #include <algorithm> 5 using namespace std; 6 typedef long long LL; 7 const int maxn = 1e5+5; 8 9 char s[maxn]; 10 11 int main() 12 { 13 int T; 14 freopen("A-large-practice (1).in","r",stdin); 15 freopen("A.txt","w",stdout); 16 scanf("%d",&T); 17 for (int kase = 1;kase <= T;kase++) 18 { 19 //scanf("%s",s+1); 20 cin >> (s+1); 21 int len = strlen(s+1); 22 printf("Case #%d: ",kase); 23 if (len%2) {puts("AMBIGUOUS");} 24 else 25 { 26 char t[105]; 27 t[2] = s[1]; 28 t[len-1] = s[len]; 29 30 for (int i = 4;i <= len;i += 2) 31 { 32 t[i] = (s[i-1] - 'A' - (t[i-2] - 'A' ) + 26) % 26 + 'A'; 33 } 34 for (int i = len-3;i >= 1;i -= 2) 35 { 36 t[i] = ( (s[i+1] - 'A') - (t[i+2] - 'A') + 26) % 26 + 'A'; 37 } 38 39 for (int i = 1;i <= len;i++) printf("%c",t[i]); 40 puts(" "); 41 } 42 } 43 return 0; 44 } View Code?
7.12 ?>_<
論文中...發(fā)現(xiàn)有好幾天沒看操作系統(tǒng)了....要割的感覺
?
7.13 ?>_<
論文中 & 看了一些相關(guān)論文關(guān)于實(shí)驗(yàn)結(jié)果的評(píng)估標(biāo)準(zhǔn)（感覺自己看東西都很不仔細(xì)...
?
7.14 >_<
Python roc
還是不懂qt的一些控件該怎么用
?
7.15 >_<
大早上興高采烈出門結(jié)果悲劇咯><
中午看了51nod一個(gè)教學(xué)視頻！
然后還是不懂Qt的那些button要怎么搞成一個(gè)數(shù)組一樣- -| (
Project Euler
Multiples of 3 and 5?第一次交PE的題目！紀(jì)念 ！暴力 或者 容斥
?
7.16 >_<
原來py的大樹 算 2^2009這樣 只要 2 ** 2009 就好了呀 ！
py字符串拼接成命令執(zhí)行
已經(jīng)第很多次這個(gè)錯(cuò)誤了：IndentationError: unindent does not match any outer indentation level
可以把代碼粘到notepad++看看（顯示空格和tab鍵）
button數(shù)組可以加到一個(gè)group里面就好，然后界面類寫的差不多之后應(yīng)該是一些細(xì)節(jié)的問題叭><
?
Even Fibonacci numbers?小于4000000的偶數(shù)斐波那契數(shù)的和
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; typedef long long LL; const int maxn = 1e5+5;LL f[100];int main() {f[1] = 1;f[2] = 2;for (int i = 3;i < 100;i++){f[i] = f[i-1] + f[i-2];if (f[i] > 4000000) break;}LL sum = 0LL;for (int i = 1;i < 100;i++){if (f[i] % 2 == 0) sum += f[i];}printf("%lld\n",sum);return 0; } View Code?Largest prime factor?600851475143的最大素因子
1 #include<cstdio> 2 #include<cstring> 3 #include<cmath> 4 #include<iostream> 5 #include<algorithm> 6 #include<vector> 7 #include<set> 8 using namespace std; 9 10 typedef long long LL; 11 const int maxn = 100005; 12 int dp[maxn]; 13 int Max[maxn]; 14 vector<int> p; 15 16 int vis[maxn]; 17 int T,a[105]; 18 int cnt[maxn]; 19 set<long long > s; 20 LL n; 21 22 void solve( long long x){ 23 LL y = x; 24 for(long long j = 2;1LL*j*j <= y;j++){ 25 if(y%j == 0){ 26 while(y%j == 0){ 27 s.insert(j); 28 y = y/j; 29 } 30 } 31 } 32 if(y > 1) s.insert(y); 33 34 for(set<long long >::iterator it = s.begin();it != s.end();++it){ 35 printf("*it = %d\n",*it); 36 } 37 38 } 39 40 41 int main(){ 42 while(scanf("%I64d",&n) != EOF){ 43 s.clear(); 44 solve(n); 45 } 46 return 0; 47 } View Code?
周杰倫日！撒花！
?
7.17 >_<
好像論文字?jǐn)?shù)湊太多（要開始刪一刪和提煉一下，白翻那么久qwq
字符串最小表示法? ?之前有看，但今天才板子敲了一敲
題目
1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 #include <algorithm> 5 using namespace std; 6 typedef long long LL; 7 const int maxn = 5e6+5; 8 9 char s[maxn]; 10 11 int getmin(char s[]) 12 { 13 int i = 0,j = 1,k = 0; 14 int len = strlen(s); 15 while(i < len && j < len && k < len) 16 { 17 int t = s[(i+k) % len] - s[(j+k) % len]; 18 if(t == 0) k++; 19 else 20 { 21 if(t > 0) i += k+1; 22 else j += k+1; 23 if (i == j) j++; 24 k = 0; 25 } 26 } 27 return i < j ? i : j; 28 } 29 30 int main() 31 { 32 int T; 33 scanf("%d",&T); 34 while(T--) 35 { 36 scanf("%s",s); 37 printf("%d\n",getmin(s)); 38 } 39 return 0; 40 } View Code這周的討論課講了深度學(xué)習(xí)...可能是我一點(diǎn)沒學(xué)習(xí)了解過，所以感覺大部分是在講一個(gè)框架...而沒有通過一個(gè)小的demo講工作原理- -?
?
7.18 >_<
早上看到一個(gè)月色真美的結(jié)局的回答，戳中...瞬間淚目qwq
數(shù)據(jù)集不平衡
繼續(xù)界面，看完影像報(bào)告之后發(fā)現(xiàn)還要加好多東西- -
隱藏一個(gè)控件?（百度谷歌中英文......終于搜到了...哭）
?Weakness and Poorness?三分加前綴和，講可以三分那一段想起白薯上的三分例題><
還沒寫..?
?
Problem B.?X Squared
看好幾天..終于明白做法了，主要是那個(gè)矩形的定義那句話沒讀明白所以搞的一頭霧水..
矩形是指，只有四角有'X'的矩形
前提條件是只有一行一列是只有一個(gè)X的，然后其他行其他列都得有兩個(gè)X
可以從最終狀態(tài)來看，是有(n-1)個(gè)矩形的，無論怎么換位置，都不會(huì)破壞掉矩形（所以這個(gè)算充分性？）
必要性如果一個(gè)矩形是破損的話，但是無論怎么換X都是不能破壞矩形的所以矛盾了？........不會(huì)證明....
1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 #include <algorithm> 5 #include <map> 6 #include <vector> 7 using namespace std; 8 typedef long long LL; 9 const int maxn = 1e5+5; 10 11 char g[105][105]; 12 int n,row[105],col[105],vis[105][105]; 13 14 int main() 15 { 16 int T; 17 freopen("B-large-practice.in","r",stdin); 18 freopen("B2.txt","w",stdout); 19 scanf("%d",&T); 20 for (int kase = 1;kase <= T;kase++) 21 { 22 scanf("%d",&n); 23 memset(row,0,sizeof(row)); 24 memset(col,0,sizeof(col)); 25 memset(vis,0,sizeof(vis)); 26 for (int i = 1;i <= n;i++) 27 { 28 scanf("%s",g[i]+1); 29 for (int j = 1;j <= n;j++) 30 { 31 if (g[i][j] == 'X') row[i]++,col[j]++; 32 } 33 } 34 printf("Case #%d: ",kase); 35 int c1 = 0,c2 = 0,d1 = 0,d2 = 0; 36 for (int i = 1;i <= n;i++) 37 { 38 if (row[i] == 1) c1++; 39 if (row[i] == 2) c2++; 40 if (col[i] == 1) d1++; 41 if (col[i] == 2) d2++; 42 } 43 //printf("c1 = %d c2 = %d d1 = %d d2 = %d\n",c1,c2,d1,d2); 44 if (c1 != 1 || c2 != n-1 || d1 != 1 || d2 != n-1) {puts("IMPOSSIBLE");continue;} 45 46 vector<pair<int,int> > h[105]; 47 for (int j = 1;j <= n;j++) 48 { 49 int p1 = 0,p2 = 0; 50 for (int i = 1;i <= n;i++) 51 { 52 if (g[i][j] == 'X') 53 { 54 if (p1 == 0) p1 = i; 55 else p2 = i; 56 } 57 } 58 h[j].push_back(make_pair(p1,p2)); 59 } 60 61 int cnt = 0; 62 for (int i = 1;i <= n;i++) 63 { 64 if (row[i] == 1) continue; 65 int p1 = 0,p2 = 0; 66 for (int j = 1;j <= n;j++) 67 { 68 if (g[i][j] == 'X') 69 { 70 if (p1 == 0) p1 = j; 71 else p2 = j; 72 } 73 } 74 //printf("i = %d p1 = %d p2 = %d\n",i,p1,p2); 75 pair<int,int> l = h[p1][0]; 76 pair<int,int> r = h[p2][0]; 77 if (l.first == r .first && l.second == r.second && !vis[p1][p2]) 78 { 79 cnt++; 80 vis[p1][p2] = 1; 81 } 82 } 83 //printf("cnt = %d\n",cnt); 84 if (cnt == (n-1)/2) puts("POSSIBLE"); 85 else puts("IMPOSSIBLE"); 86 87 } 88 return 0; 89 } View Code?
7.19 >_<
?論文刪刪改改...
然后界面現(xiàn)在遇到的問題是，把一些Button在初始化的時(shí)候連接起來，沒有被連接的Button們就只能選中一個(gè)了，嘗試寫另一個(gè)函數(shù)可是沒有卵用...想哭
D. "Or" Game?
還有點(diǎn)沒想明白...
?
Problem C.?Magical Thinking
小數(shù)據(jù)暴力2^10或者貪心
1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 #include <algorithm> 5 using namespace std; 6 typedef long long LL; 7 const int maxn = 1e5+5; 8 9 char s[105],t[105]; 10 int n,q,x; 11 12 int main() 13 { 14 int T; 15 freopen("C-small-practice.in","r",stdin); 16 freopen("Csmall.txt","w",stdout); 17 scanf("%d",&T); 18 for (int kase = 1;kase <= T;kase++) 19 { 20 scanf("%d %d",&n,&q); 21 scanf("%s",s+1); 22 scanf("%s",t+1); 23 scanf("%d",&x); 24 printf("Case #%d: ",kase); 25 int a = 0,d = 0; 26 for (int i = 1;i <= q;i++) 27 { 28 if (s[i] == t[i]) a++; 29 else d++; 30 } 31 int ans = 0; 32 if (x <= a) ans = x + d; 33 else ans = a + d - (x-a); 34 printf("%d\n",ans); 35 } 36 return 0; 37 } View Code?
7.20 >_<
手還痛.....qwq
論文被批評(píng)寫得不好了....qwq
?
7.21 >_<
接下來要做的是盡快把論文搞成英文 以及 導(dǎo)管的定位
?
前幾天問wkc YITU面試都是什么....幾乎都是算法的題目( 而且不像是leetcode 那樣的...而且聽他說了說題目還都不會(huì)的樣子TAT
所以感覺自己還是差很多...很多...（所以要好好加油呀！XDDDDD
?C - Multi-judge Solving?這是一道沒讀對(duì)題意wa了5發(fā)的水題
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; typedef long long LL; const int maxn = 1e5+5;int n,k,a[maxn];int main() {scanf("%d %d",&n,&k);for (int i = 1;i <= n;i++) scanf("%d",&a[i]);sort(a+1,a+n+1);int ans = 0,pos;for (int i = 1;i <= n;i++){if (2 * k >= a[i]) {k = max(k,a[i]);continue;}else{while(k*2 < a[i]) k = k << 1,ans++;k = max(k,a[i]);}//printf("k = %d ans = %d\n",k,ans); }printf("%d\n",ans);return 0; } View Code?
7.22 >_<
?Problem C.?Magical Thinking
大數(shù)據(jù)的做法，因?yàn)橹挥袃蓚€(gè)人，所以一共四種情況
枚舉這四種情況分別是錯(cuò)的個(gè)數(shù) 為i j k l 在滿足第一個(gè)人和第二個(gè)人得分的情況下更新答案
dp的做法還沒有看
1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 #include <algorithm> 5 using namespace std; 6 typedef long long LL; 7 const int maxn = 1e5+5; 8 9 int n,q,b[105]; 10 char s[105][105],t[105]; 11 12 int main() 13 { 14 int T; 15 freopen("C-large-practice.in","r",stdin); 16 freopen("Cbig.txt","w",stdout); 17 scanf("%d",&T); 18 for (int kase = 1;kase <= T;kase++) 19 { 20 scanf("%d %d",&n,&q); 21 for (int i = 1;i <= n;i++) scanf("%s",s[i] + 1); 22 scanf("%s",t + 1); 23 for (int i = 1;i <=n;i++) scanf("%d",&b[i]); 24 printf("Case #%d: ",kase); 25 int ans = 0; 26 if (n == 1) 27 { 28 int a = 0,d = 0; 29 for (int i = 1;i <= q;i++) 30 { 31 if (s[1][i] == t[i]) a++; 32 else d++; 33 } 34 if (b[1] <= a) ans = b[1] + d; 35 else ans = a + d - (b[1] - a); 36 } 37 else 38 { 39 int A = 0,B = 0,C = 0,D = 0; 40 for (int i = 1;i <= q;i++) 41 { 42 if (s[1][i] == s[2][i] && s[1][i] == t[i]) A++; 43 else if (s[1][i] != s[2][i] && s[1][i] == t[i]) B++; 44 else if (s[1][i] != s[2][i] && s[2][i] == t[i]) C++; 45 else D++; 46 } 47 48 //printf("A = %d B = %d C = %d D = %d\n",A,B,C,D); 49 50 for (int i = 0;i <= A;i++) 51 { 52 for (int j = 0;j <= B;j++) 53 { 54 for (int k = 0;k <= C;k++) 55 { 56 for (int l = 0;l <= D;l++) 57 { 58 if (i + j + (C - k) + (D - l) == q - b[1]) 59 if (i + (B - j) + k + (D - l) == q - b[2]) 60 { 61 ans = max(ans,q - (i + j + k + l)); 62 // printf("i = %d j = %d k = %d l = %d\n",i,j,k,l); 63 64 } 65 } 66 } 67 } 68 } 69 } 70 printf("%d\n",ans); 71 } 72 return 0; 73 } View Code?
找不到基礎(chǔ)(適合我這種紙張)的教程 qaq
?
7.23 >_<
初稿完成！
LSTM的一個(gè)栗子?碼著要看
?
7.24 >_<
又是新的一周，感覺日子過得又慢又快qaq
配置theano 和 keras?和 這一篇
配置 theano 感覺有點(diǎn)問題
import theano 不會(huì)報(bào)錯(cuò)
print theano.config.blas.ldflags 也不會(huì)報(bào)錯(cuò)
但是 theano.test() 會(huì)報(bào)錯(cuò) 找不到 nose.config
import keras 沒有問題.... 害怕qaq
?
7.25 >_<
卸載py2.7 和 昨天搞那一坨了，然后用 Anaconda(py3.5 ) + tensorflow
環(huán)境配的有點(diǎn)怕了...現(xiàn)在看起來好像還沒有什么問題的樣子
Tensorflow入門教程
LSTM預(yù)測股票每日最高票價(jià)的例子?碼著要看
?
7.26 >_<
栗子
改論文，不管是語法還是內(nèi)容.... 都好多問題
粒子濾波
前兩天一直在搜現(xiàn)成的代碼例子看，然后都本地跑不通，就一直在調(diào)bug....感覺這樣不好
于是還是從頭看起....看不懂阿QAAAQ
主要看的這個(gè)教程
?
mnist 實(shí)例，卷積神經(jīng)網(wǎng)絡(luò)
7.27 >_<
跑通一個(gè)例子，菜菜！開心 ！
7.28 >_<
py 的插值
嘗試用 GBDT
7.29 >_<
迷茫
?
7.30 >_<
py創(chuàng)建數(shù)組
自己造了一個(gè)周期性的數(shù)據(jù)驗(yàn)證了下GBDT，可能數(shù)據(jù)太規(guī)律，還蠻準(zhǔn)的..（總覺得不先驗(yàn)證下這個(gè)東西就不敢用的感覺 - ^ -
# -*- coding:UTF-8 -*- import pandas as pd import numpy as np import matplotlib.pyplot as plt import tensorflow as tf import scipy.interpolate as itp from scipy import interpolate import numpy as np import pylab as pl from sklearn import ensemble, cross_validationx = [] y = []for i in range(0,16 * 7):x.append(i)tmp = 0;if (i % 7 == 6):tmp = 12if (i % 7 == 5):tmp = 11if (i % 7 == 4):tmp = 5if (i % 7 == 3):tmp = 4if (i % 7 == 2):tmp = 3if (i % 7 == 1):tmp = 2if (i % 7 == 0):tmp = 1y.append(tmp)print (x) print (y)#pl.plot(x,y,"r-") #pl.show()# 組成 訓(xùn)練集 和 測試集 trainx = [[0 for i in range(8)]for j in range(92 * 90)] trainy = [[0 for i in range(1)]for j in range(92 * 90)] testx = [[0 for i in range(8)]for j in range(93 * 90)]for i in range(7,15 * 7):trainx[i][0] = y[i-1]#print (i,i-2)trainx[i][1] = y[i-2]trainx[i][2] = y[i-3]trainx[i][3] = y[i-4]trainx[i][4] = y[i-5]trainx[i][5] = y[i-6]trainx[i][6] = y[i-7]trainx[i][7] = i % 7trainy[i][0] = y[i]testdraw = [[0 for i in range(1)]for j in range(93 * 90)] for i in range(0,16 * 7):testdraw[i][0] = ifor i in range(7,16 * 7):testx[i][0] = y[i-1]testx[i][1] = y[i-2]testx[i][2] = y[i-3]testx[i][3] = y[i-4]testx[i][4] = y[i-5]testx[i][5] = y[i-6]testx[i][6] = y[i-7]testx[i][7] = i % 7params = {'n_estimators': 500, 'max_depth': 4, 'min_samples_split': 1,'learning_rate': 0.01, 'loss': 'ls'} clf = ensemble.GradientBoostingRegressor(**params)clf.fit(trainx,trainy) pre = clf.predict(testx)#pl.plot(x,y,'k--') #pl.plot(testdraw,pre,'r--') #pl.show()print (pre)x_ = [] y_ = []for i in range(0,16 * 7):x_.append(i)y_.append(pre[i])print (i,pre[i])pl.plot(x,y,'k--') pl.plot(x_,y_,'r--') pl.show() View Code?
現(xiàn)在能跑的一個(gè)代碼
1 # -*- coding:UTF-8 -*- 2 import pandas as pd 3 import numpy as np 4 import matplotlib.pyplot as plt 5 import tensorflow as tf 6 import scipy.interpolate as itp 7 from scipy import interpolate 8 import numpy as np 9 import pylab as pl 10 from sklearn import ensemble, cross_validation 11 12 f = open('buquan.txt') 13 line = f.readline() 14 15 X = [] 16 Y = [] 17 18 while line: 19 x = '' 20 pos = 0 21 for i in range(0,line.__len__()): 22 if (line[i] == ' '): 23 pos = i 24 break; 25 x = x + line[i] 26 27 x = int(x) 28 X.append(x) 29 30 y = '' 31 for i in range(pos+1,line.__len__()): 32 y = y + line[i] 33 y = float(y) 34 Y.append(y) 35 36 line = f.readline() 37 38 f.close() 39 40 41 x_use = [] 42 y_use = [] 43 x_debug = [] 44 y_debug = [] 45 46 47 # 用五周的數(shù)據(jù)來訓(xùn)練 組成 訓(xùn)練集 測試集 48 for i in range(0,24 * 30 * 7 * 6): 49 x_use.append(i) 50 #print (i) 51 y_use.append(Y[i]) 52 53 #pl.plot(x_use,y_use,'r--') 54 #pl.show() 55 56 trainx = [[0 for i in range(30)]for j in range(24 * 30 * 7 * 7)] 57 trainy = [[0 for i in range(1)]for j in range(24 * 30 * 7 * 7)] 58 testx = [[0 for i in range(30)]for j in range(24 * 30 * 7 * 7)] 59 60 for i in range(30,24 * 30 * 7 * 5): 61 id = x_use[i] / 720 % 7 + 2 62 trainx[i][0] = id; 63 x_debug.append(i) 64 y_debug.append(y_use[i]) 65 for j in range(1,30): 66 #print (i,j,i-(30-j)) 67 trainx[i][j] = y_use[i - (30 - j)] 68 trainy[i][0] = y_use[i] 69 70 #print (trainx) 71 #print (trainy) 72 73 for i in range(30,24 * 30 * 7 * 6): 74 id = x_use[i] / 720 % 7 + 2 75 testx[i][0] = id; 76 for j in range(1,30): 77 #print (i,j,i-(30-j)) 78 testx[i][j] = y_use[i - (30 - j)] 79 80 #print (testx) 81 82 83 84 params = {'n_estimators': 500, 'max_depth': 4, 'min_samples_split': 1, 85 'learning_rate': 0.01, 'loss': 'ls'} 86 clf = ensemble.GradientBoostingRegressor(**params) 87 88 clf.fit(trainx,trainy) 89 pre = clf.predict(testx) 90 91 x_ = [] 92 y_ = [] 93 94 for i in range(0,24 * 30 * 7 * 6): 95 x_.append(i) 96 y_.append(pre[i]) 97 #print (i,pre[i]) 98 99 pl.plot(x_,y_,'r--') 100 pl.plot(x_use,y_use,'k--') 101 pl.plot(x_debug,y_debug,'c--') 102 pl.show() View Code?
?7.31 >_<
好像不會(huì)了，要預(yù)測的那個(gè)月的前面的數(shù)據(jù)缺好多，感覺應(yīng)該不是去擬合或者補(bǔ)數(shù)據(jù)..
卡爾曼濾波
卡爾曼濾波
好像要被抓去搞另外一個(gè)軟件了xdddddd
RadiAnt?讀dcm感覺不卡，感覺這一點(diǎn)比ImageJ 好用
?
--------------- 人生中第一塊板子 ！要慢慢修煉了呀 >< 撒花的混割線 ---------------------------------
?
py里面模型的保存
?
?
?
?
轉(zhuǎn)載于:https://www.cnblogs.com/wuyuewoniu/p/7114781.html
總結(jié)
- 上一篇： 使用export/import导出和导入
- 下一篇： 作为一个软件实施工程师的发展方向