求问如何计算高温热浪的频次及天数

飞行＊想念

有一个矩阵中存放的是夏季日最高气温的资料，规定35℃为阈值，日最高气温连续6天及以上大于35℃的记为一次高温热浪过程。先要计算高温热浪的频次及天数。求问该如何使用matlab编程。多谢！

wangmol

clear all
close all
a=[10,40,40,50,40,50,40,45,50,20,25,40,22,20,40,50,40,45,40,50,40,20,12,15];
b=a(:)>35;                  %b中绝对值大于15的逻辑矩阵

ind1=find(diff([0;b])==1);     %大于35的数的开始位置
ind2=find(diff([b;0])==-1);    %大于35的数的结束位置
mask=(ind2-ind1+1>=6);       %长度大于等于6个的连续
ind3=ind1(mask);                 %连续6个大于35的数的开始位置
ind4=ind2(mask);                 %连续6个大于35的数的结束位置
p=size(ind3,1);                    %发生频次
t=ind4-ind3+1;                      %每次持续的天数

wangmol

网上查的代码，稍微改改，楼主看看对不对，注释有的不完全准确。

纸鹤

wangmol 发表于 2015-6-4 21:28
网上查的代码，稍微改改，楼主看看对不对，注释有的不完全准确。

厉害厉害，虽然不会MATLAB，但是提供了思路啊

学而不厌

tufcc

wangmol 发表于 2015-6-4 21:27
clear all
close all
a=[10,40,40,50,40,50,40,45,50,20,25,40,22,20,40,50,40,45,40,50,40,20,12,15];

非常好的代码！解决了我的一个大问题

kingjia

不错，学习了，用得到

Isabellasomnus

需要金钱来下 攒几天

嘟嘟萌萌哒

wangmol 发表于 2015-6-4 21:27
clear all
close all
a=[10,40,40,50,40,50,40,45,50,20,25,40,22,20,40,50,40,45,40,50,40,20,12,15];

有没有计算寒潮频次和天数的代码啊？

wangmol

嘟嘟萌萌哒 发表于 2016-7-18 10:32
有没有计算寒潮频次和天数的代码啊？

不好意思，我没有啊，那个代码也是临时编的。