【正文】 61371529614982072971130996298165749299991863003099630139051030215498303204572304174959330529189247306788350307696146308339462363096247973103161563116819053124059034632313239312731427896131566341683161892473 灰色關(guān)聯(lián)度分析程序clear,clcx1=xlsread(39。39。,39。A1:D31639。)。 %列出各數(shù)值，可修改x=x139。m=length(x(:,1))。n=length(x(1,:))。x2=xlsread(39。39。,39。E1:E31639。)。 %參考序列x0=x239。for i=1:navg(i)=0。 %均值初始化endfor i=1:m for j=1:navg(j)=avg(j)+x(i,j)。 endend %求均值序列for i=1:navg(i)=avg(i)/m。endfor j=1:m for i=1:n x(j,i)=x(j,i)/avg(i)。 %均值化 endendfor i=1:nx0(i)=x0(i)/avg(i)。 %參考序列均值化endfor j=1:m for i=1:n delta(j,i)=abs(x(j,i)x0(i))。 %求序列差 endendmax=delta(1,1)。for j=1:m for i=1:n if delta(j,i)max max=delta(j,i)。 end endend %求兩極差min=0。for j=1:m xgd(j)=0。 for i=1:n glxs(j,i)=*max/(*max+delta(j,i))。 %計算關(guān)聯(lián)系數(shù)及相關(guān)度 xgd(j)=xgd(j)+glxs(j,i)。 end xgd(j)=xgd(j)/n。enddisp(39。下面是四種變量與價格的灰色關(guān)聯(lián)度度分別為：39。)disp(39。 品種 釉色 尺寸 成交的年份 39。)disp(xgd)。 BP神經(jīng)網(wǎng)絡(luò)預測模型clear,clc。%導入數(shù)據(jù)x1=xlsread(39。39。,39。A91:C31639。)。 x2=xlsread(39。39。,39。D91:D31639。)。x3=xlsread(39。39。,39。E91:E31639。)。x1=x1*10^(7)。x2=x2*10^(7)。x3=x3*10^(7)。p=[x1 x2]。x4=xlsread(39。39。,39。A1:C9039。)。 x5=xlsread(39。39。,39。D1:D9039。)。x6=xlsread(39。39。,39。E1:E9039。)。x4=x4*10^(7)。x5=x5*10^(7)。x6=x6*10^(7)。t=[x4 x5]。p1=p39。[pn,~]=mapminmax(p1,0,1)。[xn,ps]=mapminmax(x339。,0,1)。%對數(shù)據(jù)進行歸一化%創(chuàng)建一個BP神經(jīng)網(wǎng)絡(luò)net=newff(pn,xn,7,{39。tansig39。 39。purelin39。},39。traingd39。)。%設(shè)置訓練參數(shù)=25。%訓練25次刷新一下訓練圖?=。%學習率=1000。%訓練次數(shù)=1e10。%訓練次數(shù)net=init(net)。%調(diào)用TRAINGDM算法訓練BP網(wǎng)絡(luò)[net,~]=train(net,pn,xn)。%輸入已知數(shù)據(jù)進行預測data=t39。[tn,~]=mapminmax(data,0,1)。%對BP網(wǎng)絡(luò)進行仿真a=sim(net,tn)。x=mapminmax(39。reverse39。,a,ps)。%對數(shù)據(jù)進行反歸一化%計算仿真誤差MSE=mse(x639。x)。disp(39。預測值是:39。)disp(x)disp(39。對應(yīng)的mse是:39。)disp(MSE)y=1:90。plot(y,x6,39。r39。,y,x,39。d39。) 清代與仿清瓷數(shù)據(jù)補全品種哪年估價或成交價格(元)轉(zhuǎn)換后價格(萬元)拍賣公司菊瓣紋罐2006150180萬165 北京翰海棒槌瓶200955萬55 中國嘉德國際拍賣有限公司圣主得賢臣頌筆筒2000750750 北京翰海柳葉瓶201060萬（1500）60 北京保利國際拍賣有限公司筆筒20101416萬15 北京中博國際拍賣有限公司十二花神杯20101000萬以上1000 北京華辰拍賣有限公司大罐201120萬以上20 北京保利國際拍賣有限公司五彩錦地魚藻紋將軍罐200567萬7 北京華辰拍賣有限公司貫耳方瓶201050萬（2500）50 北京翰海拍賣有限公司菊瓣紋盤2011174萬港元166 香港蘇富比雙耳爐201160萬60 保利香港拍賣有限公司過枝梅竹盤20003252萬港幣3097 香港佳士得碗2007403萬港幣384 香港佳士得罐2010150萬元以上150 北京翰海拍賣有限公司綬帶耳葫蘆平底瓶2003241萬港幣230 香港蘇富比黃地花卉魚紋賞瓶2003387萬港幣369 香港蘇富比扁瓶2013150200萬175 北京保利國際拍賣有限公司釉里紅龍紋玉壺春瓶2002342萬元342 香港蘇富比雕海水云龍紋燈籠瓶2002337萬港幣321 香港佳士得橄欖瓶2000197 佳士得（香港）有限公司琮式瓶200780100萬90 佳士得（香港）有限公司百鹿尊2003335萬335 浙江駿成拍賣有限公司花卉紋題詩碗20033093萬3093 香港蘇富比花蝶紋如意耳葫蘆尊20031372港幣1307 香港佳士得萬花錦碗2013130 佳士得香港國際藝術(shù)品拍賣有限公司黃地花卉魚紋賞瓶20081820萬19 云南典藏拍賣集團有限公司福壽纏枝蓮紋靈芝耳扁瓶20122030萬25 北京東正拍賣有限公司雙耳瓶2002121萬元121 北京中貿(mào)圣佳蒜頭瓶200228萬28 北京翰海綠地粉彩福壽紋獸耳瓶200399萬99 北京嘉德礬紅龍紋溫壺200217萬17 北京翰海玉壺春瓶2007810萬9 浙江錢塘拍賣有限公司雙耳瓶2011約6萬6 香港富得拍賣行有限公司盒20039 廣州嘉德天球瓶200215 中貿(mào)圣佳國際拍賣有限公司纏枝八寶紋碗200356萬6 上海信仁拍賣有限公司 關(guān)聯(lián)變權(quán)后的量化數(shù)據(jù)表 4 關(guān)聯(lián)變權(quán)后的量化數(shù)據(jù)品種釉色最大尺寸哪年估價或成交轉(zhuǎn)換后的價格427142714914491456905690651865186573657366396639678867888317831784748474854385439125912596999699117331173311733117331173311832118321288212882128821288213023130231314713147142251422515498154981610216102162781627816295162951629516295167591675917797177971932319323193981939819555195552125521255212552264102264123831238312408924089242082420824256242562425624797247972560625606257640257642576425764257642576425764257642581225812260730260732769627696281912819128340283402851628516287742877428984289843099630996309963099631148311483151531515315483154831548322053220532205322053220532205322053414134141354523545235849358493585035850359553595537195371953864538645386453864538675386753879638796391093910939224392244091440914431554315543852438524432444324456274562746013460134895148951493004930049538495384953849538495934959350071500715091350913515275152751527515275152751527517285172851728517285214652146521465214653459534595475354753564265642658354583546022860228602286206562065625226252263029630296440964409646606466065182651826518265182654366543665964