このページ最新ではありませ。をクリックし,英语のの最新版版を参照参照し。
Findpeaks
を使用,セットから局所的大とその位置をします。
ファイルspots_num.mat
1749年1749年年年年までに観测れ太阳黒点数数のの平均がが含まれてていいます。。。。このこのこのこの
最大その発生を検出し。とデータをプロットします。
加载('spots_num.mat')[pks,locs] = findpeaks(avspots);情节(年,Avspots,Year(LOCS),PKS,'或者')xlabel('年')ylabel('数字')轴紧的
距离非常ピークがいくつかます。それらは定期定期的は缲り返さ缲り返されれないないないない种类种类种类种类种类ののののものものものものものものですですです。。。。。。そのそのそのその10
周期のを正确推定するに,もう一度Findpeaks
を使用。,今回はピークピーク间隔をををを年以上に制限します。。最大値値値间の平均平均间隔间隔
[pks,locs] = findpeaks(avspots,“ minpeakdistance',6);情节(年,Avspots,Year(LOCS),PKS,'或者')xlabel('年')ylabel('数字') 标题(“黑子”)轴紧的传奇('数据',,,,“峰”,,,,'地点',,,,'西北')
循环= diff(locs);平均值=平均值(循环)
平均= 10.8696
11年约约约年に缲り返されいることははててていますますます。。。。フーリエをを使っ使ってて确认确认ししますますます。信号信号ののの変动変动変动に集中できるで测定ているをてくださいナイキスト周波以下の周波数数を使用使用。。。
FS = 1;NF = 512;df = fs/nf;F = 0:DF:FS/2-DF;trspots = fftshift(fft(avspots-mean(avspots),nf));dbspots = 20*log10(abs(trspots(nf/2+1:nf)));yaxis = [20 85];图(f,dbspots,1 ./ [卑鄙的卑鄙卑鄙],yaxis)xlabel('频率(年^{ - 1})')ylabel('|FFT |(D b)')轴([0 1/2 yaxis])文本(1/eycycle + .02,25,['<== 1/'num2str(卑鄙的)])
11年た周波周波数数なりなり,なりなりなりなりごとごと年とという予想予想予想がが确认されましましたた。変换変换の最も
[pk,maxfreq] = findpeaks(dbspots,'npeaks',1,'sorttr',,,,“下降”);周期= 1/f(maxfreq)
周期= 10.8936
抓住在图(F(Maxfreq),PK,'或者') 抓住离开传奇('傅里叶变换',,,,'1/卑鄙的',,,,'1/ofere')
2つのつのぴったり一致します。。