[原创]MATLAB流函数与势函数的数值计算程序

fisher · 发表于 2014-3-3 14:09:11

是通过求解poisson方程算的吗

fisher · 发表于 2014-3-3 14:09:48

感谢分享！！！

fisher · 发表于 2014-3-3 14:23:11

感谢分享！！！

fisher · 发表于 2014-3-3 14:25:26

感谢分享！！！

沙颖凯 · 发表于 2014-3-3 16:03:35

fisher 发表于 2014-3-3 14:07

登录/注册后可看大图

是通过求解poisson方程算的吗

是，张弛迭代法。

fisher · 发表于 2014-3-3 21:08:41

沙颖凯发表于 2014-3-3 16:03

登录/注册后可看大图

是，张弛迭代法。

楼主，我用NCEP数据的u,v,lon,lat网格试了一下，为啥散度DIVH算出来都是inf和nan,最后导致算出来的chi都是nan呢

fisher · 发表于 2014-3-3 21:31:18

找到问题，lon和lat应该做成73*144的矩阵，我做的是144*73

沙颖凯 · 发表于 2014-3-3 21:35:08

fisher 发表于 2014-3-3 21:31

登录/注册后可看大图

找到问题，lon和lat应该做成73*144的矩阵，我做的是144*73

是吗？我一直用的73*144，144*73不能用？

沙颖凯 · 发表于 2014-3-3 21:50:45

fisher 发表于 2014-3-3 21:31

登录/注册后可看大图

找到问题，lon和lat应该做成73*144的矩阵，我做的是144*73

我知道了，问题是因为差分和矩阵的行列是有关系的，行换成列差分的值就变成0了，所以计算失败，目前想不到什么好的办法自动区分行列对应经度还是纬度。

fisher · 发表于 2014-3-4 00:33:07

嗯，是这样。可以加一个判断既可以了：
lon1=lon(:,1);lat1=lat(:,1);
LLon=length(lon1);
【M，N】=size(lon);

if M==Llon,
lons=lons';lats=lats';
end

		自动登录	找回密码
密码			立即注册

[源程序] [原创]MATLAB流函数与势函数的数值计算程序