2.3. Периодограмма в matlab.
Чтобы построить периодограмму для временного ряда солнечного радиоизлучения необходимо написать программу для расчета спектра и получить график.
Программа для расчета спектра:
y=load('y.txt');
ri=load('ri107d.txt');
%
ri=detrend(ri);
%plot(y,ri);
%
ts=365*24*3600;
%
fs=1/ts;
%
[Pxx,f] = periodogram(ri,[],[],fs);%
f1=f*ts;
plot(f1,Pxx);
Эта программа позволяет рассчитать спектр и построить для нее график. Команда y=load('y.txt') и ri=load('ri107d.txt') загружает данные с текстовых файлов, где в текстовом файле «y.txt» находятся годы от 1047 по 2010, а в «ri107d.txt» находится поток радиоизлучения Солнца на длине волны 10,7 см. Команда plot(y, ri) соответствует построению обычной функции, когда одномерный массив y соответствует значениям аргумента, а одномерный массив ri - значениям функции. ts=365*24*3600- здесь выполняется расчет периода, период равен 1 году. Далее используется команда f1=f*ts для расчета частоты. [Pxx,f] = periodogram(ri,[],[],fs)- в этой команде непосредственно строится периодограмма с помощью функции periodogram и выводится результат на график с помощью команды plot(f1,Pxx). В результате получили спектр солнечного радиоизлучения и пик, который и есть одиннадцатилетний период солнечной активности. На этом графике по оси х частота в единицах 1/год, на оси у спектр радиоизлучения.
Заключение.
В этой работе была рассмотрена тема «Периоды колебания в спектре излучения солнца». Была написана программа для расчета спектра солнечного радиоизлучения и получен график со спектром. В этом графике есть пик который соответствует одиннадцатилетнему циклу солнечной активности. С помощью этой разработанной программы мы убедились в этом.
Литература.
Цифровая обработка сигналов / А. Б. Сергиенко — СПб.: Питер, 2002. — 608 с.: ил.
http://www.moveinfo.ru/picture/stablelink?page=0000277