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

流云洗出天空蓝 · 发表于 2020-5-28 14:24:01

金钱-5啥情况

铅笔削你 · 发表于 2020-6-18 11:37:58

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

次旺平措 · 发表于 2020-6-18 11:53:28

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

将军 · 发表于 2020-9-26 11:27:03

非常赞谢谢楼主

1297712690 · 发表于 2021-3-12 17:47:05

楼主牛逼啊

1297712690 · 发表于 2021-4-25 21:29:25

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

taiyang · 发表于 2021-10-14 17:56:00

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

贫道敬孔 · 发表于 2022-9-2 22:39:37

感谢楼主分享！！

zhangruijie · 发表于 2023-3-14 15:59:43

请问下其中的sor_index是指什么

haochizhenhaoch · 发表于 2024-6-4 20:40:11

非常感谢楼主的代码！我在分析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更好？差别会不会很大？

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

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