MATLAB太陽(yáng)黑子活動(dòng)周期

太陽(yáng)黑子的活動(dòng)周期研究班級(jí)：電控研14-2姓名：慶春2021/6/271太陽(yáng)黑子（太陽(yáng)的暗黑斑點(diǎn)）太陽(yáng)黑子（sunspot）是在太陽(yáng)的光球?qū)由习l(fā)生的一種太陽(yáng)活動(dòng)，是太陽(yáng)活動(dòng)中最基本、最明顯的。一般認(rèn)為，太陽(yáng)黑子實(shí)際上是太陽(yáng)表面一種熾熱氣體的巨大漩渦，溫度大約為3000-4500℃。2021/6/272因?yàn)槠錅囟缺忍?yáng)的光球?qū)颖砻鏈囟纫?000到2000攝氏度（光球?qū)颖砻鏈囟燃s為6000攝氏度），所以看上去像一些深暗色的斑點(diǎn)。太陽(yáng)黑子很少單獨(dú)活動(dòng)，通常是成群出現(xiàn)。黑子的活動(dòng)周期為11年，活躍時(shí)會(huì)對(duì)地球的磁場(chǎng)產(chǎn)生影響，主要是使地球南北極和赤道的大氣環(huán)流作經(jīng)向流動(dòng)，從而造成惡劣天氣，使氣候轉(zhuǎn)冷。嚴(yán)重時(shí)會(huì)對(duì)各類電子產(chǎn)品和電器造成損害。太陽(yáng)黑子雖然顏色較"深"，但是在觀測(cè)情況下，與太陽(yáng)耀斑同樣清晰顯眼。天文學(xué)家把太陽(yáng)黑子最多的年份稱為“太陽(yáng)活動(dòng)峰年”，太陽(yáng)黑子最少的年份稱為“太陽(yáng)活動(dòng)谷年2021/6/273如何觀測(cè)一般的光學(xué)望遠(yuǎn)鏡主要是觀測(cè)太陽(yáng)光球?qū)拥幕顒?dòng)現(xiàn)象：黑子、耀斑和米粒組織。我們主要觀測(cè)黑子和耀斑。1．太陽(yáng)黑子是光球上最明顯的活動(dòng)，在日面的形成、發(fā)展、消失、形態(tài)、多少、大小和分布等，都有一定的規(guī)律。黑子在日面上的東西分布是不對(duì)稱的——東半球比西半球多，并且集中分布在日面中緯度區(qū)域。較完整的黑子是由本影和半影組成，中間深黑色的部分叫本影，四周暗淡的部分叫半影。另外大多數(shù)黑子都成群出現(xiàn)，多是沿東西方向分布著。一群中，西邊較大的黑子叫前導(dǎo)黑子，東邊較大的黑子叫后隨黑子。就太陽(yáng)整體而言，黑子活動(dòng)有大約11年的周期變化。2021/6/2742．光斑是在太陽(yáng)光球?qū)拥臇|西邊緣，往往出現(xiàn)一些比光球的背景更明亮的區(qū)域，這就是光斑。光斑常伴隨日面邊緣的黑子一起出現(xiàn)，也存在著11年的周期。觀測(cè)的方法：太陽(yáng)觀測(cè)有幾種不同的方法。我們暫時(shí)只選用投影觀測(cè)，以后進(jìn)行照相觀測(cè)。下面的內(nèi)容中將重點(diǎn)介紹一下投影觀測(cè)，簡(jiǎn)要說(shuō)明照相觀測(cè)。2021/6/275投影觀測(cè)一、測(cè)前準(zhǔn)備（1）觀測(cè)紙：印有日本威信公司生產(chǎn)的太陽(yáng)觀測(cè)網(wǎng)格圖，觀測(cè)者根據(jù)當(dāng)日的B0、P值選擇，共有八種（在下文中介紹）。（2）描圖紙：印有與觀測(cè)紙網(wǎng)格圖同樣大小稀疏網(wǎng)格圖（直徑為10厘米，B0=0度，P=0度），也可以用白而薄的紙（紙上有直徑10厘米的圓）（3）天頂棱鏡：若使用的太陽(yáng)投影板平面與鏡筒平行，需使用天頂棱鏡改變光路。（4）目鏡：不能使用膠合目鏡，用惠更斯目鏡較為合適。不能使用高倍目鏡。（5）濾光片：ND中性濾光片，有ND2、ND4、ND8多種。（6）天文普及年歷：觀測(cè)當(dāng)天的B0、L0、P值都需從當(dāng)年的天文年歷中查得。（7）常用的工具：鉛筆、手表、紅黃彩筆、尺、圓規(guī)、投影儀、小磁鐵或夾子。（8）帶上前一天的描點(diǎn)圖。2021/6/2761、將望遠(yuǎn)鏡極軸大致對(duì)準(zhǔn)北極，鏡筒指向太陽(yáng)，觀察鏡筒的影子，當(dāng)影子最小時(shí)，說(shuō)明太陽(yáng)已經(jīng)大致對(duì)好。！?。≌?qǐng)注意，千萬(wàn)不能用眼睛通過(guò)望遠(yuǎn)鏡直接尋找太陽(yáng)，尋星鏡也最好不要加裝。2、記錄透明度、寧?kù)o度、云量、描圖開(kāi)始時(shí)間等觀測(cè)背景數(shù)據(jù)。關(guān)于寧?kù)o度和透明度的規(guī)定是這樣的：范圍是0-5最好是5，最差是0，一般都要在2.5以上才可以觀測(cè)。寧?kù)o度與風(fēng)有關(guān)，一般4級(jí)以下可以定為3。透明度與污染等有關(guān)，取值看經(jīng)驗(yàn)了。3、裝上、目鏡、投影板，調(diào)好焦距，將描圖紙固定在投影板上（我們用小的吸鐵石或夾子），上下移動(dòng)板，使太陽(yáng)的像與描圖紙上的圓面等大。4、不要跟蹤太陽(yáng)，選擇一顆小而清楚的黑子，觀察其在描圖紙上的視運(yùn)動(dòng)軌跡，確定東西線。5、若使用網(wǎng)格圖來(lái)描點(diǎn)，只需旋轉(zhuǎn)描圖紙，使黑子的視運(yùn)動(dòng)軌跡方向（太陽(yáng)的東西線方向）與描圖紙中的緯線平行。使用白紙只需在數(shù)據(jù)分析時(shí)畫(huà)出東西線即可。2021/6/2776、確定日面坐標(biāo)方向①若投影板與鏡筒平行，由于描圖紙上的像是經(jīng)過(guò)折射鏡和棱鏡兩次改變光路得到的，故太陽(yáng)黑子視運(yùn)動(dòng)方向恰與日面坐標(biāo)中的東西走向大體一致，再按“上北下南，左西右東”，日面坐標(biāo)方向就確定了。②若投影板與鏡筒垂直，由于只經(jīng)過(guò)一次改變光路，方向完全相反：按上南下北，左東右西。即可確定日面坐標(biāo)方向。7、把投影相和描圖紙的圓對(duì)齊并開(kāi)始跟蹤，在日面上找一顆形狀較為規(guī)則的大黑子，適當(dāng)調(diào)整望遠(yuǎn)鏡的赤經(jīng)、赤緯旋鈕，當(dāng)投影日面與描圖紙上的圓重合時(shí)迅速描下該黑子的位置，其他黑子就可以以它為主參照黑子定位。另外邊看邊記黑子群數(shù)、個(gè)數(shù)、形態(tài)、位置等。此時(shí)要對(duì)照前一天的黑子描點(diǎn)圖，并確定出黑子群的類型（采用蘇黎式分類法）2021/6/2788、描圖時(shí)要迅速，手不要壓屏，頭不要碰目鏡。9、其他黑子的描圖需按以下步驟進(jìn)行：①黑子需自西向東描。②先描黑子半影，后描本影。③先描前導(dǎo)黑子，后描后隨黑子。10、描出光斑。11、記錄描圖結(jié)束的時(shí)間，以描圖開(kāi)始結(jié)束時(shí)間的中間時(shí)刻，記為描圖時(shí)間。2021/6/279數(shù)據(jù)見(jiàn)附表太陽(yáng)黑子的活動(dòng)呈現(xiàn)周期性變化的，沃爾夫(R.Wolfer)根據(jù)在過(guò)去的288年(1700年~1987年)間每年太陽(yáng)黑子出現(xiàn)的數(shù)量和大小的觀測(cè)數(shù)據(jù)推算出11年的周期規(guī)律。2021/6/2710數(shù)據(jù)處理離散時(shí)間信號(hào)的傅里葉變換DiscreteTimeFourierTransform,DTFTDFS是周期序列的離散傅里葉級(jí)數(shù)DTFT是非周期序列的傅里葉變換，稱離散時(shí)間傅里葉變換，其頻譜是連續(xù)的函數(shù)DFT是有限長(zhǎng)序列的離散傅里葉變換，是對(duì)其DTFT的等間隔抽樣，是離散的頻譜DFT是DFS的主值序列，是非周期的。而DFS是DTFT的頻域內(nèi)的抽樣。FFT是DFT的一種高效快速算法，也稱作快速傅里葉變換。2021/6/2711離散傅立葉變換（DFT）

公式2021/6/2712這一對(duì)式子，左、右兩邊都是離散的，有限長(zhǎng)，因此可方便地用來(lái)實(shí)現(xiàn)頻譜分析。但使用時(shí)，一定要想到，它們均來(lái)自DFS,即x（n）和X（k）都是周期的！FFT是一種DFT的高效算法，稱為快速傅立葉變換（fastFouriertransform）。FFT算法可分為按時(shí)間抽取算法和按頻率抽取算法，先簡(jiǎn)要介紹FFT的基本原理。從DFT運(yùn)算開(kāi)始，說(shuō)明FFT的基本原理。2021/6/2713

2021/6/2714實(shí)現(xiàn)步驟及內(nèi)容

(1)讀取數(shù)據(jù)并繪制Wolfer圖，觀察太陽(yáng)黑子的活動(dòng)現(xiàn)象。程序如下：loadsunsb.dat%讀取數(shù)據(jù)year=sunb(:,1);%讀取年份信息wolfer=sunb(:,2);%讀取黑子活動(dòng)數(shù)據(jù)figure%新建圖像plot(year,wolfer)%畫(huà)出時(shí)域圖xlabel('Years');ylabel('SunbData');title('SunbData')%標(biāo)注橫縱坐標(biāo)和標(biāo)題pause2021/6/27152021/6/2716為清楚起見(jiàn)，畫(huà)出部分時(shí)間段(如實(shí)驗(yàn)中采取最近50年)的Wolfer圖進(jìn)行觀察。程序如下：figureplot(year(261:311),wolfer(261:311),'b.-');%用進(jìn)50年數(shù)據(jù)做圖xlabel('Years');ylabel('SunbData');title('Atthelast50years')%標(biāo)注橫縱坐標(biāo)和標(biāo)題2021/6/27172021/6/2718(2)對(duì)已經(jīng)得到的Wolfer數(shù)應(yīng)用FFT技術(shù)分析它的變化規(guī)律。程序如下：Y=fft(wolfer);%對(duì)全部數(shù)據(jù)做FFTY(1)=[];%由于算法問(wèn)題，Y(1)為所有fft數(shù)值之和，因此需舍棄第一個(gè)點(diǎn)pause觀察Wolfer數(shù)在FFT分析后得到的在復(fù)平面中的分布圖。程序如下：figureplot(Y,'ro')%在復(fù)平面做圖，空心點(diǎn)title('FourierCoefficientsintheComplexPlane');xlabel('RealAxis');%標(biāo)注橫坐標(biāo)為’實(shí)部’ylabel('ImaginaryAxis');%標(biāo)注縱坐標(biāo)為’虛部’pause2021/6/27192021/6/2720(3)功率與頻率的關(guān)系曲線則被定義為周期圖。用周期圖法對(duì)Wolfer數(shù)進(jìn)行功率譜密度分析。程序：figuren=length(Y);%取FFT結(jié)果長(zhǎng)度為npower=abs(Y(1:n/2)).^2;%取前n/2個(gè)數(shù)據(jù)求其功率nyquist=1/2;%取最大頻率為0.5freq=(1:n/2)/(n/2)*nyquist;%將FFT轉(zhuǎn)換為DTFT%假設(shè)采樣頻率為Fs，采樣點(diǎn)數(shù)為N，做FFT之%后，某一點(diǎn)n（n從1開(kāi)始）表示的頻率為：%Fn=(n-1)*Fs/N；該點(diǎn)的模值除以N/2就是對(duì)應(yīng)。真實(shí)頻率plot(freq,power)%畫(huà)周期圖xlabel('cycles/year');title('Periodogram')pause2021/6/27212021/6/2722為清楚起見(jiàn)，對(duì)功率和頻率的前50個(gè)分量作它的周期圖。程序如下：figureplot(freq(1:50),power(1:50))%為觀察方便xlabel('cycles/year')pause2021/6/2723(4)確定出太陽(yáng)黑子的活動(dòng)周期。為清楚起見(jiàn)，畫(huà)出功率與周期（頻率的倒數(shù)）的關(guān)系曲線圖。程序如下：figureperiod=1./freq;%將橫坐標(biāo)單位轉(zhuǎn)換為年/周期plot(period,power);%畫(huà)周期圖axis([05002e+7]);ylabel('Power');xlabel('Period(Years/Cycle)');%單位為年/周期pause2021/6/27242021/6/2725(5)在功率與周期關(guān)系曲線圖中標(biāo)出功率的最高點(diǎn)，該位置對(duì)應(yīng)的周期則為太陽(yáng)黑子活動(dòng)周期。holdon;index=find(power==max(power));%找到頻率最大點(diǎn)，該點(diǎn)橫坐標(biāo)即為