MoM(五)解矩阵方程

最新推荐文章于 2024-03-12 14:43:51 发布

MinJinFan

最新推荐文章于 2024-03-12 14:43:51 发布

阅读量1.1k

点赞数

分类专栏：计算电磁学

本文链接：https://blog.csdn.net/qq_37264323/article/details/104800211

版权

这篇博客主要介绍了在MoM方法中如何解决矩阵方程的问题，内容包括阻抗矩阵填充后的处理步骤，虽然暂时没有公式展示，但已提供相关代码实现。

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

在阻抗矩阵填充之后，接下来就是解矩阵方程了。后期在补上公式，先上代码。

注意：方程求得解是 RWG 函数的系数。

#pragma once
#include<iostream>
#include<stdio.h>
#include<armadillo>
#include<math.h>
#include<complex>
#include<time.h>
#include<cmath>
#include"FillZmat.h"

using namespace arma;
using namespace std;

extern double k;
extern mat basis_list;
extern cube element_list;
extern int elem_num;
extern int rwg_num;
extern mat Center;
extern vec Area;

extern mat CommonEdgeP;
extern mat CommonEdgeN;
extern vec EdgeLength;  //公共边长度
extern vec Vertex_pos;  // 正三角形的 非公共边点 编号
extern vec Vertex_naga;

extern cube Center_;   //九点积分划分后的小三角形重心

extern mat Pho_Pos;       //对于正element三角形 非公共边点到重心的向量坐标 
extern mat Pho_Naga;      //对于负element三角形 非公共边点到重心的向量坐标
extern cube Pho_Pos_;     //对于正element三角形 九点划分后 非公共边点到 小三角形重心的