请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4459|回复: 4

[其他] 【已解决】求助-ncl傅里叶分析

[复制链接]
回帖奖励 10 金钱 回复本帖可获得 5 金钱奖励! 每人限 1 次

新浪微博达人勋

发表于 2022-12-6 18:00:22 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 是冉冉升起的冉 于 2022-12-7 18:55 编辑

ncl官网中傅里叶分析用的ezfftf和ezfftb这两个函数,这个例子中分解1波和2波时:
https://www.ncl.ucar.edu/Applications/Scripts/fanal_3.ncl


CF = ezfftf (z)                             ; ezfftf works on right dim  
printVarSummary(CF)                         ; [2] x [12] x [121] x [120]  

cf = CF  
cf(:,:,:,1:)  = 0.0                         ; set waves 2 and higher to 0.0


z_wave1 = ezfftb (cf, 0.0)                  ; [12] x [121] x [240]  
copy_VarMeta(z, z_wave1)                    ; [time | 12] x [lat | 121] x [lon | 240]

这一行cf(:,:,:,1:)  = 0.0                         ; set waves 2 and higher to 0.0
求1波时设置2波以及更高的波数为0,那么为什么是设置最后一维1:为0,不应该是第一维,如果是最后一维那不就是经度么?







密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2022-12-7 09:40:20 | 显示全部楼层
顶顶我自己,两种设置我都试过了,值也print出来看过,如果是例子中就是经度,第一维就是波,两种画出来的图差异还挺大的。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-12-7 15:39:00 | 显示全部楼层

回帖奖励 +5 金钱

因为ezfftb这个函数经度那一维展开成了2项,展开后的最后一维并不是经度,你可以看下这个函数的说明
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-12-7 18:52:35 | 显示全部楼层
穿着人字拖 发表于 2022-12-7 15:39
因为ezfftb这个函数经度那一维展开成了2项,展开后的最后一维并不是经度,你可以看下这个函数的说明

感谢!!!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-12-7 18:54:56 | 显示全部楼层
ezfftf是对最右边的维度进行转换,如果需要对指定的维度转换,可以使用ezfftf_n。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表