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

1、太陽(yáng)黑子的活動(dòng)周期研究 班級(jí)：電控研14-2 姓名：慶春太陽(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。 因?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)品和電器造成損

2、害。 太陽(yáng)黑子雖然顏色較深，但是在觀測(cè)情況下，與太陽(yáng)耀斑同樣清晰顯眼。天文學(xué)家把太陽(yáng)黑子最多的年份稱為“太陽(yáng)活動(dòng)峰年”，太陽(yáng)黑子最少的年份稱為“太陽(yáng)活動(dòng)谷年如何觀測(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)黑子

3、，東邊較大的黑子叫后隨黑子。就太陽(yáng)整體而言，黑子活動(dòng)有大約11年的周期變化。 2光斑是在太陽(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è)。投影觀測(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度） ，也可以用白而薄的紙（紙上

4、有直徑10厘米的圓）（3）天頂棱鏡：若使用的太陽(yáng)投影板平面與鏡筒平行，需使用天頂棱鏡改變光路。（4）目鏡：不能使用膠合目鏡，用惠更斯目鏡較為合適。不能使用高倍目鏡。（5）濾光片：ND中性濾光片，有ND2、ND4、ND8多種。（6）天文普及年歷：觀測(cè)當(dāng)天的B0、L0、P值都需從當(dāng)年的天文年歷中查得。（7）常用的工具：鉛筆、手表、紅黃彩筆、尺、圓規(guī)、投影儀、小磁鐵或夾子。（8）帶上前一天的描點(diǎn)圖。 1、將望遠(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度、云量、描

5、圖開始時(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í)畫出東西線即可。 6、 確定日面坐標(biāo)方向 若投影板與鏡筒

6、平行，由于描圖紙上的像是經(jīng)過(guò)折射鏡和棱鏡兩次改變光路得到的，故太陽(yáng)黑子視運(yùn)動(dòng)方向恰與日面坐標(biāo)中的東西走向大體一致，再按“上北下南，左西右東”，日面坐標(biāo)方向就確定了。 若投影板與鏡筒垂直，由于只經(jīng)過(guò)一次改變光路，方向完全相反：按上南下北，左東右西。即可確定日面坐標(biāo)方向。 7、把投影相和描圖紙的圓對(duì)齊并開始跟蹤，在日面上找一顆形狀較為規(guī)則的大黑子，適當(dāng)調(diào)整望遠(yuǎn)鏡的赤經(jīng)、赤緯旋鈕，當(dāng)投影日面與描圖紙上的圓重合時(shí)迅速描下該黑子的位置，其他黑子就可以以它為主參照黑子定位。另外邊看邊記黑子群數(shù)、個(gè)數(shù)、形態(tài)、位置等。 此時(shí)要對(duì)照前一天的黑子描點(diǎn)圖，并確定出黑子群的類型（采用蘇黎式分類法） 8、描圖時(shí)要迅速，

7、手不要壓屏，頭不要碰目鏡。 9、其他黑子的描圖需按以下步驟進(jìn)行： 黑子需自西向東描。 先描黑子半影，后描本影。 先描前導(dǎo)黑子，后描后隨黑子。 10、描出光斑。 11、記錄描圖結(jié)束的時(shí)間，以描圖開始結(jié)束時(shí)間的中間時(shí)刻，記為描圖時(shí)間。數(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ī)律。數(shù)據(jù)處理 離散時(shí)間信號(hào)的傅里葉變換 Discrete Time Fourier Transform, DTFT DFS是周期序列的離散傅里葉級(jí)數(shù) DTFT是非周期序列的傅里葉變換，稱離散時(shí)間

8、傅里葉變換，其頻譜 是連續(xù)的函數(shù) DFT是有限長(zhǎng)序列的離散傅里葉變換，是對(duì)其DTFT的等間隔抽樣，是離散的頻譜 DFT是DFS的主值序列，是非周期的。而DFS是DTFT的頻域內(nèi)的抽樣。 FFT是DFT的一種高效快速算法，也稱作快速傅里葉變換。離散傅立葉變換（DFT） 公式 這一對(duì)式子，左、右兩邊都是離散的，有限長(zhǎng)，因此可方便地用來(lái)實(shí)現(xiàn)頻譜分析。 但使用時(shí)，一定要想到，它們均來(lái)自DFS, 即 x（n）和 X（k）都是周期的！ FFT是一種DFT的高效算法，稱為快速傅立葉變換（fast Fourier transform）。FFT算法可分為按時(shí)間抽取算法和按頻率抽取算法，先簡(jiǎn)要介紹FFT的基本原理

9、。從DFT運(yùn)算開始，說(shuō)明FFT的基本原理。實(shí)現(xiàn)步驟及內(nèi)容 (1) 讀取數(shù)據(jù)并繪制Wolfer圖，觀察太陽(yáng)黑子的活動(dòng)現(xiàn)象。程序如下： load sunsb.dat %讀取數(shù)據(jù) year=sunb(:,1); %讀取年份信息 wolfer=sunb(:,2); %讀取黑子活動(dòng)數(shù)據(jù) figure%新建圖像 plot(year,wolfer) %畫出時(shí)域圖 xlabel(Years); ylabel( Sunb Data ); title(Sunb Data) %標(biāo)注橫縱坐標(biāo)和標(biāo)題 pause 為清楚起見(jiàn)，畫出部分時(shí)間段(如實(shí)驗(yàn)中采取最近50年)的Wolfer 圖進(jìn)行觀察。程序如下： figure p

10、lot(year(261:311),wolfer(261:311),b.-); %用進(jìn)50年數(shù)據(jù)做圖 xlabel(Years);ylabel( Sunb Data ); title(At the last 50 years) %標(biāo)注橫縱坐標(biāo)和標(biāo)題(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)tit

11、le(Fourier Coefficients in the Complex Plane);xlabel(Real Axis); %標(biāo)注橫坐標(biāo)為實(shí)部ylabel(Imaginary Axis); %標(biāo)注縱坐標(biāo)為虛部pause(3) 功率與頻率的關(guān)系曲線則被定義為周期圖。用周期圖法對(duì)Wolfer數(shù)進(jìn)行功率譜密度分析。程序：figure n=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è)采樣

12、頻率為Fs，采樣點(diǎn)數(shù)為N，做FFT之%后，某一點(diǎn)n（n從1開始）表示的頻率為：%Fn=(n-1)*Fs/N；該點(diǎn)的模值除以N/2就是對(duì)應(yīng)。真實(shí)頻率plot(freq,power) %畫周期圖xlabel(cycles/year);title(Periodogram)pause 為清楚起見(jiàn)，對(duì)功率和頻率的前50個(gè)分量作它的周期圖。程序如下：figureplot(freq(1:50),power(1:50) %為觀察方便xlabel(cycles/year)pause(4) 確定出太陽(yáng)黑子的活動(dòng)周期。為清楚起見(jiàn)，畫出功率與周期（頻率的倒數(shù)）的關(guān)系曲線圖。程序如下：figureperiod=1./freq; %將橫坐標(biāo)單位轉(zhuǎn)換為年/周期plot(period,power);%畫周期圖axis(0 50 0 2e+7);ylabel(Power);xlabel(Period (Years/Cycle);%單位為年/周期pause(5) 在功率與周期關(guān)系曲線圖中標(biāo)出功率的最高點(diǎn)，該位置對(duì)應(yīng)的周期則為太陽(yáng)黑子活動(dòng)周期。hold on;index=find(power=max(pow