清晰;
关闭所有;
%的纬度
% lat = 1:4
lat = 2;
%的否决一天带通(即。10到50天)
ta = 10;
结核病= 50;
dt = 1;%每日数据
%为空间带通(我的否决。e, 4到12度经度)
xa = 4;
xb = 12;
dx = 0.25;% 0.25度分辨率
%选择start结束
台北= 2000;
year2 = 2020;
k = 0;
年=(台北:year2);
loadstring =['负载C:\Users\laure\Desktop\MEA593\ Project_MEA593 / Combined_mw_v5_atlantic_ int2str(年)“.mat”);
eval (loadstring);
tm = sstmx;
ym = lat_actual;
xm = lon_actual;
%的这部分代码获取值所需的纬度。如果有
%两个纬度同样接近,它选择南方。也
%关注60 w 20 e经度域。
xm = xm (1:273);
iy = min(找到(abs (ym-lat) < 0.25));%找到最近的纬度。
ym = ym (iy);
tm = tm (:, 1:273 iy);
[nt, nx] = (tm)大小;
tm (tm = = -9999) =南;
如果国防部(,4)= = 0,飞跃= 1;别的,飞跃= 0;结束;
tmyear (k + 1: k + nt:,:) = tm;
k = k +元;
结束;
纽约(nt nx] =大小(tmyear);k =元;
% %填写缺失的数据
tmyear = fillnan_new (xm (1: nt),挤压(tmyear));
%的时间带通
tmean = nanmean (tmyear);
tmyear = tmyear-ones (nt, 1) * tmean;
tm50 = get_lpass_bartlett (tmyear,结核,1);
tmbp = tmyear-tm50;
tmbp = get_lpass_bartlett (tmbp助教1);
tmbp = tmbp ';
%带状带通
tmean = nanmean (tmbp);
tmbp = tmbp-ones (nx 1) * tmean;
tmbp (isnan (tmbp) = = 1) = 0;
tm48 = get_lpass_bartlett (tmbp xb / dx 1);
tmxtbp = tmbp-tm48;
tmxtbp = get_lpass_bartlett (tmxtbp xa / dx 1);
tmxtbp = tmxtbp ';
panel1 = tmxtbp;
panel1 (panel1 = = 0) =南;
%抓323 - 355度1 n从xm
panel2 = panel1 (:, 93:226);
%行panel2关联
行= 1:7671
%在panel2每一行,nanmean的朗1的值
%为每一行而不是134
panel3(行,1)= nanvar (panel2(行,:));
结束
%花时间在轮廓panel3代表另一个列
%时间值,情节简单
% panel3 (:, 2) = (0.5: 1: nt - 0.5) / 365.25 +台北;
图(1);情节(panel3“k”,“线宽”,1)
标题(“频率”);
包含(的时间(年)')
ylabel(海温信号)
%设置(gca, XTickLabel, 2000:2.5:2020)
这是现在的样子。甘氨胆酸的使用并没有真正的帮助。我知道最后一个峰对应于图2020.5和将结束前达到2021,因此一路到2020年。
