用递归法实现多重网格法——MATLAB

最新推荐文章于 2024-07-03 21:27:14 发布

sdhver214

最新推荐文章于 2024-07-03 21:27:14 发布

阅读量4.3k

点赞数 2

分类专栏： matlab代码文章标签：多重网格法 Multigrid Multilevel FMG Matlab

本文链接：https://blog.csdn.net/zhanwl1992/article/details/83780881

版权

本文通过MATLAB递归代码详细介绍了如何实现多重网格法（FMG），并附带了运行图，展示了FMG循环的功能。引用了《Multilevel methods in lubrication》作为参考文献。

摘要由CSDN通过智能技术生成

FMG运行图

在这里插入图片描述

MATLAB递归代码

clear; clc;
global f1
f1 = fopen('out1.txt', 'w'); %该文件只是为了测试
k = 5; ks = 1;		% 总的网格层、开始层
gamma = 1; ncy = 1;	% 类型(V or W), 总的循环次数 
nu0 = 6; nu1 =2; nu2 = 1; 
fmg(k, ks, nu0, nu1, nu2, gamma, ncy)
fclose(f1);
a = load('out1.txt');
figure(1)
plot(a, 'bo-', 'LineWidth', 1.5);

% FMG
function fmg(k, ks, nu0, nu1, nu2, gamma, ncy)
global f1
if k == ks
    if ks == 1
        for i = 1 : nu0
			% 迭代代码
            fprintf(f1, '%d\n', k);
        end
    else
        for j = 1 : ncy
            cycle(k, nu0, nu1, nu2, gamma);
        end
    end
else
    if k > ks
        fmg(k - 1, ks, nu0, nu1, nu2, gamma, ncy)		
		% 向高一层插值代码		
        for j = 1 : ncy
            cycle(k, nu0, nu1