用超松弛迭代法求解接地金属槽内电位分布

一、实验内容:

试用超松弛迭代法求解接地金属槽内电位的分布。 已知:a 4cm,h a/4 10mm 给定边值如图所示。 给定初值: i(,0j) 0 误差范围: 10 5 计算迭代次数, i,j分布。

二.实验设计原理:有限差分法

有限差分法(Finite Differential Method)是基于差分原理的一种数值计算法。其基本思想:将场域离散为许多小网格,应用差分原理,将求解连续函数 的泊松方程的问题换为求解网格节点上 的差分方程组的问题。

编程时已经考虑到题目要求,所以直接将边值编入到程序中,编写成function的M文件,这样只要调用此M文件,输入变量为迭代因子,即可输出电位矩阵和迭代次数。 迭代时所用公式为

U2(i,j)=U1(i,j)+a*(U1(i,j+1)+U1(i+1,j)+U2(i-1,j)+U2(i,j-1)-4*U1(i,j))/4 其中U2代表k+1,而U1代表k。

以上分析了迭代程序的实现,但是迭代循环如何终止并未说明。题目中的误差范围ε=0.00001,即当两次迭代结果相差不超过ε时停止,这里必须是九点都满足不超过ε,而并不是其中某一点达到即可。当迭代次数过多时,程序会运行很长时间,(本题要求电位点数较少,不会出现迭代次数过多的情况。当然点数越多结果越精确。)当迭代因子a≥2时,迭代不收敛,程序会陷入死循环,因此需要限制循环次数,迭代100000次无结果则退出循环,防止程序崩溃。

这样可以画出流程图如下所示:

0

0

0

10V

启动

输入迭代因子

迭代次数 k=0

k=k+1 开始循环迭代

函数判断相邻二 否

次差值是否小于

给定值

输出k,电位U1

三、程序运行界面及结果

终止

求矩形槽内电位分布matlab,MATLAB超松弛迭代法求解接地金属槽内电位分布相关推荐

  1. 求矩形槽内电位分布matlab,MATLAB求解接地金属槽内电位分布

    基于MATLAB有限差分法中的迭代法求解接地金属槽内电位分布 要求 运用MATLAB求解接地金属槽内点位分布, 精度.行数M.列数N自己定义. 有限差分法 有限差分法是基于差分原理的一种数值计算法.其 ...

  2. matlab里面接地咋找,MATLAB求解接地金属槽内电位分布

    基于MATLAB有限差分法中的迭代法求解接地金属槽内电位分布 要求 运用MATLAB求解接地金属槽内点位分布, 精度.行数M.列数N自己定义. 有限差分法 有限差分法是基于差分原理的一种数值计算法.其 ...

  3. MATLAB用超松弛迭代法求接地金属槽内的电位分布

    工程电磁场仿真实验-用超松弛法求接地金属槽内电位分布 一.实验目的 1.熟悉利用超松弛法描绘出接地金属槽的电位分布情况: 2.学会使用MATLAB进行数值计算,并绘出相应的图形: 二.实验原理 1.有 ...

  4. 使用有限差分法求解无限长直的金属槽内的电位分布

    使用有限差分法求解无限长直的金属槽内的电位分布 实验内容: 实验源码: 实验截图 实验内容: 使用有限差分法求解无限长直的金属槽内的电位分布 实验源码: %控制宽度 Length = 50; %长度 ...

  5. 非线性方程组牛顿迭代法matlab,matlab实现牛顿迭代法求解非线性方程组

    <matlab实现牛顿迭代法求解非线性方程组>由会员分享,可在线阅读,更多相关<matlab实现牛顿迭代法求解非线性方程组(5页珍藏版)>请在人人文库网上搜索. 1.matla ...

  6. matlab牛顿法解非线性方程组,matlab实现牛顿迭代法求解非线性方程组.pdf

    资源描述 matlab 实现牛顿迭代法求解非线性方程组实现牛顿迭代法求解非线性方程组 已知非线性方程组如下 3*x1-cos(x2*x3)-1/2=0 x1^2-81*(x2+0.1)^2+sin(x ...

  7. matlab的超松弛迭代法,基于MATLAB的超松弛迭代法的电位研究

    一.引言 随着计算机技术.实验研究技术的迅速发展,电磁场学科在教研.工程上的应用也在逐步加深.而经典电磁学理论作为电磁技术发展的理论基础,其核心麦克斯韦方程组的重要性不言而喻.因此在工程上求解电磁场的 ...

  8. 超松弛迭代法求解二维电磁场有限差分方程(附Matlab代码)

    二维电磁场泊松方程差分格式 由泰勒公式: 以及: 两式做和,截断于项,得到二阶差商: 二维场域内泊松方程为: ​​​​​​​        ​​​​​​​        ​​​​​​​         ...

  9. matlab用牛顿迭代法求解方程,牛顿迭代法求方程解 程序如下

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 function [x_reality,n_reality] = Newt( f_name,x_start,tolerance,n_limit) %% % ...

  10. Matlab利用牛顿迭代法求解非线性方程组

    我们得首先了解牛顿迭代法的原理是什么: 在这里,我以二阶非线性方程组为例: f1(x,y)=0 f2(x,y)=0,求解x,y 原理 假设方程组的一组近似解为(x0,y0),将方程f1(x,y)=0与 ...

最新文章

  1. 彻底解决Spring mvc中时间的转换和序列化等问题
  2. 神奇的HyperLogLog算法
  3. ConcurrentHashMap的源码分析-put方法第四阶段
  4. linux文件 内存映射 锁,linux – mmap:将映射文件立即加载到内存中吗?
  5. kubectl命令大全
  6. mybatis传入参数类型parameterType详解
  7. Jsp基本指令和动作
  8. html5移动端转盘,基于zepto适合移动端的幸运大转盘插件
  9. 【机械仿真】基于matlab GUI智能电梯仿真系统【含Matlab源码 1087期】
  10. h5居中loading_H5样式与布局 --常用居中方法
  11. 约束最优化方法之最优性条件
  12. 【Ware】有了这些软件,马上告别360全家桶!
  13. android 双向进度条,Android自定义双向进度条
  14. HC32l130单片机的delay函数怎么写
  15. 欧洲篮球冠军联赛网站
  16. EXCEL一个单元格内容分成多个单元格
  17. UR机器人TCP通讯示例 详细例程,手把手教会你
  18. 【Educoder】HTML答案 (HTML入门——基础、HTML入门——基本标签、HTML——表单类的标签)
  19. 赛普拉斯PSoC6正式接入阿里云Link TEE加强物联网应用的安全设计...
  20. Tensorflow应用之LSTM

热门文章

  1. 快递市场定价 到底有何影响
  2. 2021秋软工实践第二次结对编程作业
  3. 易大师接口自动化测试平台如何创建不同协议的接口并进行测试
  4. 一篇小文入门 Python
  5. socket can 编程
  6. oracle rac 各日志,oracle rac 日志体系结构!
  7. markdown生成chm电子书或在线文档
  8. 指针使用入门与 unsafe.Pointer
  9. fio_generate_plots
  10. checkbox如何默认选中