clear
clc
a=[2139425154];
for?t=1:5
????data{t}=xlsread（“county.xlsx“num2str（a（t）））;
????
????conty{t}=xlsread（“county.xlsx“num2str（a（t）*10））;
end

%t為每張表的標號（一個表表示一個州）
for?t=1:5
????k=1;m=1;%設置跑動計數變量的初值
????while?m~=size（data{t}1）?%當跑到每個表結尾時跳出循環
????????for?j=0:7
????????????%按年份次序用NaN補齊缺省年份數據
????????????if?data{t}（m1）==2010+j
????????????????b{t}（k）=data{t}（m4）;
????????????????k=k+1;m=m+1;
????????????????if?m==size（data{t}1）
????????????????????b{t}（k）=data{t}（m4）;
????????????????????break;
????????????????end
????????????else
????????????????b{t}（k）=NaN;
????????????????k=k+1;
????????????end
????????end
????end
end

for?t=1:5
????for?k=1:8
????????%做個行指標為時間，列指標為縣序號，值為查出次數的矩陣
????????c{t}（:k）=b{t}（k:8:end）;
????end
end

%篩選
%%%一些必要的處理%%%
h=c;
%%%
clear?d?d_
%%%
for?t=1:5
????v=1;
????for?i=1:size（c{t}1）
????????for?j=1:size（h{t}2）
????????????h{t}（isnan（h{t}（:j））==1j）=0;
????????end
????????%控制下面if的參數可以控制界限
????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????????if?（max（c{t}（i:））-min（c{t}（i:）））/mean（h{t}（i:））>2
????????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????????????d{t}（v:）=c{t}（i:）;
????????????d_{t}（v）=conty{t}（i）;
????????????v=v+1;
????????end
????end
end


%作圖
a_=‘s-d-h-o-p-‘;
for?t=1:5
????subplot（23t）;
????xlabel（“year“）;
????ylabel（“num“）;
????title（num2str（a（t）‘））;
????axis?square
????hold?on
????if?max（size（d{t}））~=0
????????for?k=1:size（d_{t}2）
????????????f{t}（k:）=num2str（d_{t}（k））;
????????end
????????
????????plot（2010:2017d{t}a_（（2*t-1）:2*t））;
????????hl{t}=legend（f{t}）;
????????set（hl{t}‘Box‘‘off‘）;
????????
????????%set（hl{t}‘Orientation‘‘horizon‘）
????????%legend（‘toggle‘）
????????
????end
end

%%%%%%%%%%
%%%%%%%%%%

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1881??2019-01-28?23:34??1910007\1910007?代碼\county.m

?????文件????6261093??2019-01-28?23:35??1910007\1910007?代碼\county.xlsx

?????文件?????666150??2019-01-28?23:37??1910007\1910007?代碼\map.xlsx

?????文件????1266702??2019-01-28?22:32??1910007\1910007?代碼\MCM_NFLIS_Data.xlsx

?????文件???????4146??2019-01-28?23:31??1910007\1910007?代碼\p1_2.m

?????文件????6562144??2019-01-28?06:29??1910007\1910007?代碼\p1_2.xlsx

?????文件????????529??2019-01-27?18:46??1910007\1910007?代碼\part2.m

?????文件???11026071??2019-01-27?17:55??1910007\1910007?代碼\part2.xlsx

?????文件????????417??2019-01-28?23:22??1910007\1910007?代碼\state.m

?????文件?????530431??2019-01-29?06:37??1910007\1910007.docx

?????文件?????808373??2019-01-29?06:36??1910007\1910007.pdf

?????文件????1994901??2019-01-26?11:16??1910007\19美賽C題程序參考1\19美賽C題程序參考\KY_Features_20181201.xlsx

?????文件????1276468??2019-01-25?19:35??1910007\19美賽C題程序參考1\19美賽C題程序參考\MCM_NFLIS_Data.xlsx

?????文件?????713897??2019-01-26?11:15??1910007\19美賽C題程序參考1\19美賽C題程序參考\OH_Features_20181201.xlsx

?????文件????1360983??2019-01-26?11:07??1910007\19美賽C題程序參考1\19美賽C題程序參考\PA_Features_20181201.xlsx

?????文件???????1473??2019-01-26?16:08??1910007\19美賽C題程序參考1\19美賽C題程序參考\SHOW.m

?????文件????1281413??2019-01-25?19:14??1910007\19美賽C題程序參考1\19美賽C題程序參考\VA_Features_20181201.xlsx

?????文件????1293419??2019-01-26?13:23??1910007\19美賽C題程序參考1\19美賽C題程序參考\讀我.pdf

?????文件??????23506??2019-01-25?12:16??1910007\2019_MCM_Problem_C.docx

?????文件?????409514??2019-01-25?12:16??1910007\2019_MCM_Problem_C.pdf

?????文件?????447980??2019-01-29?02:31??1910007\論文中文版3.0（3）.docx

?????文件?????781635??2019-01-25?12:16??1910007\面對阿片類藥物_?特點_?CDC.pdf

?????目錄??????????0??2019-02-06?12:20??1910007\19美賽C題程序參考1\19美賽C題程序參考

?????目錄??????????0??2019-01-29?01:49??1910007\1910007?代碼

?????目錄??????????0??2019-02-06?12:20??1910007\19美賽C題程序參考1

?????目錄??????????0??2019-02-06?12:20??1910007

-----------?---------??----------?-----??----

?????????????36713126????????????????????26