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

流云洗出天空蓝

金钱-5啥情况

铅笔削你

您好，请问您附的那张图里的风速箭头的大小大约是多少啊？因为看到你的流函数 量级是10的7次方 m2/s，所以我在想风速应该是多少？ 如果按照量纲来推算的话，风速岂不是有10的7次方 m/s吗？ 有点迷惑

次旺平措

感谢分享！！！！！！！！！！！！！！！！

将军

非常赞  谢谢楼主

1297712690

楼主牛逼啊

1297712690

楼主您好，我想请教一下，这个是局地的对吧？全球似乎不能用？算出来结果不太一致

taiyang

为什么我这里会显示输入的参数数目不足呢？

贫道敬孔

感谢楼主分享！！

zhangruijie

请问下其中的sor_index是指什么

haochizhenhaoch

非常感谢楼主的代码！我在分析RWS时需要知道辐散风场，所以正好需要楼主的代码，也让我了解了泊松方程的数值解法！
有两处疑问提出，希望楼主或者有共同问题的同学们有空瞅一眼：
1、楼主的主贴中介绍到u v形式是二维矩阵，它的第一个维度代表纬度，第二个维度代表经度，那么迭代循环中Res的拉普拉斯项的dx和dy是否顺序颠倒了？即：
Res(i, j)=(chi(i+1, j)+chi(i-1, j)-2*chi(i, j))./dx2(i, j)+...
                          (chi(i, j+1)+chi(i, j-1)-2*chi(i, j))./dy2(i, j)+...
                           divh(i, j);
是否应改为：
Res(i, j)=(chi(i+1, j)+chi(i-1, j)-2*chi(i, j))./dy2(i, j)+...
                          (chi(i, j+1)+chi(i, j-1)-2*chi(i, j))./dx2(i, j)+...
                           divh(i, j);
2、在MAX=1000处有注释，对应eps为10^-7，但下面计算里用的10^-5，是否换成-7更好？差别会不会很大？