主要内容

このページ最新ではありませ。をクリックし,英语のの最新版版を参照参照し。

データのピーク検出

Findpeaksを使用,セットから局所的大とその位置をします。

ファイルspots_num.mat1749年1749年年年年までに観测れ太阳黒点数数のの平均がが含まれてていいます。。。。このこのこのこの

最大その発生を検出し。とデータをプロットします。

加载('spots_num.mat')[pks,locs] = findpeaks(avspots);情节(年,Avspots,Year(LOCS),PKS,'或者')xlabel('年')ylabel('数字')轴紧的

图包含一个轴。轴包含2个类型线的对象。

距离非常ピークがいくつかます。それらは定期定期的は缲り返さ缲り返されれないないないない种类种类种类种类种类ののののものものものものものものですですです。。。。。。そのそのそのその10

周期のを正确推定するに,もう一度Findpeaksを使用。,今回はピークピーク间隔をををを年以上に制限します。。最大値値値间の平均平均间隔间隔

[pks,locs] = findpeaks(avspots,“ minpeakdistance',6);情节(年,Avspots,Year(LOCS),PKS,'或者')xlabel('年')ylabel('数字') 标题(“黑子”)轴紧的传奇('数据',,,,“峰”,,,,'地点',,,,'西北'

图包含一个轴。带有标题日光点的轴包含2个类型线的对象。这些对象表示数据,峰值。

循环= 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(卑鄙的)])

图包含一个轴。轴包含3个类型行的对象,文本。

11年た周波周波数数なりなり,なりなりなりなりごとごと年とという予想予想予想がが确认されましましたた。変换変换の最も

[pk,maxfreq] = findpeaks(dbspots,'npeaks',1,'sorttr',,,,“下降”);周期= 1/f(maxfreq)
周期= 10.8936
抓住图(F(Maxfreq),PK,'或者') 抓住离开传奇('傅里叶变换',,,,'1/卑鄙的',,,,'1/ofere'

图包含一个轴。轴包含4个类型行的对象,文本。这些对象代表傅立叶变换,1/卑鄙的1/周期。

2つのつのぴったり一致します。。

参考

|

关连するトピック