MATLAB滤波之后，位势高度500hpa纬向波除0波外，1，2，3波均出现负值

Yang.hahha

常山赵子龙

http://bbs.06climate.com/forum.php?mod=viewthread&tid=12210
这个貌似和你的工作比较相关的。帮你找了一下，你可以参考看一下，楼主可否借鉴一下你的滤波程序呢？

Yang.hahha

不好意思，这么晚才回复
function [lb0t3]=lb(data)
%0—3波分别在1—4位
a(1:13,1:73)=0;
b(1:13,1:73)=0;
%有限求和
        for k=1:4
            for y=1:73
                for x=1:144
                    a(k,y)=a(k,y)+data(x,y).*cosd((k-1)*2.5*x);
                    b(k,y)=a(k,y)+data(x,y).*sind((k-1)*2.5*x);
                end
            end
        end

%平均部分
a(1:13,1:73)=a(1:13,1:73)./72;
b(1:13,1:73)=b(1:13,1:73)./72;
a(1,1:73)=a(1,1:73)./2;
b(1,1:73)=b(1,1:73)./2;
%单波
        for y=1:73
            for x=1:144
                lb(1,x,y)=a(1,y);
            end
        end
  
        for k=2:4
            for y=1:73
                for x=1:144
                    lb(k,x,y)=a(k,y).*cosd((k-1)*2.5*x)+b(k,y).*sind((k-1)*2.5*x);
                end
            end
        end        
%行星尺度波0-3波
        for y=1:73
            for x=1:144
                %lb0t3(x,y)=lb(1,x,y)+lb(2,x,y)+lb(3,x,y)+lb(4,x,y);
                lb0t3(x,y)=lb(2,x,y)+lb(3,x,y)+lb(4,x,y);
            end
        end