




版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
W01E0140120100W01E0140120100太陽(yáng)黑子周期分析1計(jì)算太陽(yáng)黑子周期1)、選取歷年的太陽(yáng)黑子數(shù)據(jù)本次作業(yè)選取的是1700—1999年的太陽(yáng)黑子數(shù)據(jù)。將數(shù)據(jù)導(dǎo)入matlab中,并繪制太陽(yáng)黑子數(shù)隨年份變化的關(guān)系曲線。如圖 1所示。程序如下:clearloadsunspot.datyear=sunspot(:,1);sunspot=sunspot(:,2);plot(year(1:300),sunspot(1:300),'b.-');xlabel('years');ylabel('sunspotdata');title('1700 —1999年太陽(yáng)黑子數(shù)是隨年份變化的關(guān)系曲線 ');gridon1TO-1999年衣酣黑子數(shù)罡隨年份變直的黃系曲裟200604020□▼-1700 1750 1800 1850 1900 1950 2000years圖1、太陽(yáng)黑子數(shù)隨年份的變化曲線2):利用功率譜密度函數(shù)分析周期1、對(duì)已經(jīng)得到的Wolfer數(shù)進(jìn)行FFT變換分析它的變化規(guī)律,并作功率與頻率的關(guān)系圖。y=fft(sunspot(1:300));y(1)=[];n=length(y);power=abs(y(1:n/2)).A2;q=1/2;f=(1:n/2)/(n/2)*q;plot(f,power);xlabel('周期/年');title('周期圖');運(yùn)行結(jié)果如圖2所示。圖2、太陽(yáng)黑子的功率譜為了清楚起見(jiàn),取功率和頻率的前 50個(gè)分量作它的周期圖,程序如下plot(f(1:50),power(1:50));xlabel('頻率');運(yùn)行結(jié)果如圖3所示。頻率圖3、功率和頻率的前50個(gè)分量的周期圖2、確定太陽(yáng)黑子的活動(dòng)周期,畫(huà)出功率與周期的關(guān)系圖。程序如下:T=1./f;plot(T,power);axis([05007e+6]);%X軸圍是0-50,Y軸圍是0-7*10人6xlabel('周期');ylabel('功率');gridon殆功率與周期的關(guān)系圖上標(biāo)出功率的最高點(diǎn), 該位置對(duì)應(yīng)的周期即為太陽(yáng)黑子活動(dòng)的周期。程序如下:holdonindex=find(power==max(power));m=num2str(T(index));plot(T(index),power(index),'r.','MarkerSize',25);text(T(index)+2,power(index),['T=',m]);holdoff運(yùn)行結(jié)果如圖4所示:77(x(xTx)1xTy (2)x106W15 20 263D35 40 45 50Xt 1Xtx106W15 20 263D35 40 45 50Xt 1Xt1 2Xt2 atat:NID(0,/)(1)1.6161.41.31211圖4、太陽(yáng)黑子周期圖運(yùn)用功率譜方法計(jì)算出太陽(yáng)黑子的活動(dòng)周期為 T=11.0741,這與Wolfer得出的11年的周期規(guī)律基本一致,說(shuō)明實(shí)驗(yàn)方法是正確的。2、利用ARMA模型,預(yù)測(cè)未來(lái)某年的太陽(yáng)黑子數(shù)1)、建立AR模型選用二階自回歸模型AR(2),方程為:采用最小二乘法對(duì)參數(shù) 1、2進(jìn)行估計(jì):模型殘差方差:1xt1xt12Xt2)2(3)計(jì)算參數(shù)程序如下:x=zeros(298,2);fori=2:1:299x(i-1,1)=sunspot(i);endfork=1:1:298x(k,2)=sunspot(k);endy=zeros(298,1);fort=3:1:300y(t-2,1)=sunspot(t);endA=x';B=x'*x;C=inv(B);D=C*A*y運(yùn)行得D=1.4867-0.5981即1 計(jì)算參數(shù)程序如下:x=zeros(298,2);fori=2:1:299x(i-1,1)=sunspot(i);endfork=1:1:298x(k,2)=sunspot(k);endy=zeros(298,1);fort=3:1:300y(t-2,1)=sunspot(t);endA=x';B=x'*x;C=inv(B);D=C*A*y運(yùn)行得D=1.4867-0.5981即1 1.48672 -0.5981帶入公式(3)解得 364.1380求解程序如下:symssm=0;s=sunspot(1:300);fori=3:300;m=m+(s(i)-1.4867*s(i-1)+0.5981*s(i-2))A2;endn=m/298n=364.1380解得364.1380解得364.1380故得到AR(2)模型方程是:xt1.4867xt10.5981xt2at其中at~NID(0,364.1380)2)、用上述AR(2)模型進(jìn)行檢驗(yàn)并預(yù)測(cè)Sunspot(1998)=64.3,Sunspot(1999)=93.3,Sunspot(2000)=119.6利用上述AR(2)模型計(jì)算得:Sunspot(2000)=1.4867*93.3-0.5981*64.3=100.2513誤差率=(119.6-100.2513 )/119.6=16.18%;Sunspot(2004)=40.4,Sunspot(2005)=29.8,Sunspot(2006)=15.2利用上述AR(2)模型計(jì)算得:Sunspot(2006)=1.4867*29.8-0
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 9 短詩(shī)三首 繁星(一五九)教學(xué)設(shè)計(jì)-2023-2024學(xué)年統(tǒng)編版語(yǔ)文四年級(jí)下冊(cè)
- 2024年泰山版小學(xué)信息技術(shù)二年級(jí)上冊(cè)《3管理計(jì)算機(jī)的資源》教學(xué)設(shè)計(jì)
- 2023八年級(jí)物理下冊(cè) 第十一章 小粒子與大宇宙第二節(jié) 看不見(jiàn)的運(yùn)動(dòng)教學(xué)設(shè)計(jì) (新版)滬科版
- 4我對(duì)世界很好奇 教學(xué)設(shè)計(jì)二年級(jí)下冊(cè)心理健康同步備課系列(浙教版)
- 2024新教材高中歷史 第二單元 中古時(shí)期的世界 第5課 古代非洲與美洲教學(xué)設(shè)計(jì) 部編版必修中外歷史綱要下
- 財(cái)務(wù)風(fēng)險(xiǎn)警示教育
- 激勵(lì)行業(yè)心態(tài)
- 2024-2025學(xué)年高中政治 第3單元 第7課 第2框 弘揚(yáng)中華民族精神教學(xué)設(shè)計(jì) 新人教版必修3
- 2024年七年級(jí)語(yǔ)文上冊(cè) 第一單元 成長(zhǎng)紀(jì)事 第4課《傷仲永》教學(xué)設(shè)計(jì)1 滬教版五四制
- 九年級(jí)歷史上冊(cè) 第二單元 亞洲和歐洲的封建社會(huì) 第4課 亞洲封建國(guó)家的建立教學(xué)設(shè)計(jì)2 新人教版
- 血糖監(jiān)測(cè)技術(shù)操作考核評(píng)分標(biāo)準(zhǔn)
- 堤防、淤地壩運(yùn)行一般危險(xiǎn)源風(fēng)險(xiǎn)評(píng)價(jià)方法-風(fēng)險(xiǎn)矩陣法(LS法)介紹
- 鍶90敷貼器等設(shè)備技術(shù)參數(shù)及配置要求
- 國(guó)家開(kāi)放大學(xué)《成本管理》形考任務(wù)(1-4)試題答案解析
- 2023年05月2023年江蘇南京市博物總館招考聘用人員筆試參考題庫(kù)附答案解析
- 信息技術(shù)基礎(chǔ)知到章節(jié)答案智慧樹(shù)2023年海南師范大學(xué)
- 退出中華人民共和國(guó)國(guó)籍申請(qǐng)表
- 甘肅省專業(yè)標(biāo)準(zhǔn)化技術(shù)委員會(huì)考核評(píng)估細(xì)則
- 2023工會(huì)春游活動(dòng)方案(7篇)
- 二年級(jí)音樂(lè)上冊(cè) 《大頭娃娃》教學(xué)課件
- 第四章莖尖分生組織培養(yǎng)
評(píng)論
0/150
提交評(píng)論