精华内容
下载资源
问答
  • 有限差分法MATLAB程序

    万次阅读 多人点赞 2018-04-22 19:05:01
    利用超松弛法程序如下: d=zeros(1,10); h=0; for e=1:0.1:1.9 a=zeros(21,41); a(1,:)=100; b=zeros(19,39); c=eye(19,39); count=1; g=0; while(count==1) m=0; for i=2:1:20 for j=2:1:40 b(i-1,j-1)=a(i,j);...

    设有一个长直接地金属矩形槽,长a=40,宽b=20,其侧壁与底面电位均为零,顶盖电位为100V(相对值),求槽内电位分布。

    利用高斯迭代求解代码如下(相邻两次迭代值最大允许误差为0.001):

    a=zeros(21,41);
    a(1,:)=100;
    b=zeros(19,39);
    c=eye(19,39);
    count=1;
    d=0;
    while(count==1)
        m=0;
        for i=2:1:20
            for j=2:1:40
                b(i-1,j-1)=a(i,j);
                a(i,j)=0.25*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1));
                c(i-1,j-1)=abs(a(i,j)-b(i-1,j-1));
                if(c(i-1,j-1)<0.001)
                    for k=1:1:19
                        for n=1:1:39
                            if(c(k,n)<0.00001)
                                m=m+1;
                            else
                                m=0;
                                break;
                            end
                        end
                        if(m==0)
                            break;
                        end
                        if(m==741)
                            count=0;
                        end
                    end
                end
                if(count==0)
                    break;
                end
            end
            if(count==0)
                break;
            end
        end
        d=d+1;
    end
    d
    a

    利用超松弛法程序如下:

    d=zeros(1,10);
    h=0;
    for e=1:0.1:1.9
    a=zeros(21,41);
    a(1,:)=100;
    b=zeros(19,39);
    c=eye(19,39);
    count=1;
    g=0;
    while(count==1)
        m=0;
        for i=2:1:20
            for j=2:1:40
                b(i-1,j-1)=a(i,j);
                a(i,j)=a(i,j)+e*0.25*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4*a(i,j));
                c(i-1,j-1)=abs(a(i,j)-b(i-1,j-1));
                if(c(i-1,j-1)<0.001)
                    for k=1:1:19
                        for n=1:1:39
                            if(c(k,n)<0.001)
                                m=m+1;
                            else
                                m=0;
                                break;
                            end
                        end
                        if(m==0)
                            break;
                        end
                        if(m==741)
                            count=0;
                        end
                    end
                end
                if(count==0)
                    break;
                end
            end
            if(count==0)
                break;
            end
        end
       g=g+1;
    end
    h=h+1;
    e
    d(1,h)=g;
    end
    d 

    觉得有用就赞一个~有问题可以在下面留言,我尽量回复~

    展开全文
  • MATLAB程序分享求解偏微分方程扩散方程有限差分法-MATLAB求解偏微分方程(扩散方程)有限差分法 源程序代码.rar 程序代码见附件,拿资料请顺便顶个贴~~ 如果下载有问题,请加我 qq 1530497909,给你在线传
  • 精品文档 精品文档 PAGE PAGE # 欢迎下载 有限差分法Matlab 程序(椭圆型方程) function FD_PDE(fun,gun,a,b,c,d) % 用有限差分法求解矩形域上的 Poisson 方程 tol=10*6; % 误差界 N=1000; % 最大迭代次数 n=20; %...
  • 有限差分法(Finite Differential Method, FDM)PDE数值解一类重要方法。 对一类典型的PDE给出了基于MATLAB的数值解算法。
  • MATLAB有限差分法程序

    2018-04-22 19:18:46
    MATLAB有限差分法求解拉普拉斯(Laplace)方程,长直接地金属矩形槽内部电位分布
  • function FD_PDE(fun,gun,a,b,c,d) %用有限差分法求解矩形域上的Poisson方程 tol=10^(-6); % 误差界 N=1000; % 最大迭代次数 n=20; % x轴方向的网格数 m=20; % y轴方向的网格数 h=(b-a)/n; %x轴方向的步长 l=(d-c)/m...
  • 有限差分法MATLAB程序

    千次阅读 2020-03-22 20:04:21
    lx=61;ly=101; v1=zeros(ly,lx); for j=2:lx-1 v1(ly,j)=5; end v2=v1;m=1;t=0;k=0; while(m>0.00001) k=k+1 m=0; for i=2:ly-1 for j=2:lx-1; v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v1(i-1,j)+v1(i,j-1))/4...t=(v2(i,j...
    lx=61;ly=101;
    v1=zeros(ly,lx);
    for j=2:lx-1
    v1(ly,j)=5;
    end
    v2=v1;m=1;t=0;k=0;
    while(m>0.00001)
    k=k+1
    m=0;
    for i=2:ly-1
    for j=2:lx-1;
    v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v1(i-1,j)+v1(i,j-1))/4;
    t=(v2(i,j)-v1(i,j));
    if(t>m)m=t;end
    end
    end
    v1=v2;
    end
    subplot(1,2,1),mesh(v2)
    axis([0,61,0,101,0,5])
    subplot(1,2,2),contour(v2,32)

    展开全文
  • MATLAB时域有限差分法程序

    热门讨论 2010-12-21 09:59:15
    MATLAB时域有限差分法程序,有利于理解fdtd原理,可直接运行的。
  • %*********************************************************************** [color=#FFFFFF]本文来自:[/color] ?%? ? 3-D FDTD code with PEC boundaries [color=#FFFFFF]本文来自:[/color] ?...

    %*********************************************************************** [color=#FFFFFF]本文来自:[/color] ?%? ? 3-D FDTD code with PEC boundaries [color=#FFFFFF]本文来自:[/color] ?%*********************************************************************** [color=#FFFFFF]本文来自:[/color] ?% [color=#FFFFFF]本文来自:[/color] ?%? ? Program author: Susan C. Hagness [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? Department of Electrical and Computer Engineering [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? University of Wisconsin-Madison [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? 1415 Engineering Drive [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? Madison, WI 53706-1691 [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? 608-265-5739 [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ? ? ? ? ? hagness@ [color=#FFFFFF]本文来自:[/color] ?% [color=#FFFFFF]本文来自:[/color] ?%? ? Date of this version:? February 2000 [color=#FFFFFF]本文来自:[/color] ?% [color=#FFFFFF]本文来自:[/color] ?%? ? This MATLAB M-file implements the finite-difference time-domain [color=#FFFFFF]本文来自:[/color] ?%? ? solution of Maxwell's curl equations over a three-dimensional [color=#FFFFFF]本文来自:[/color] ?%? ? Cartesian space lattice comprised of uniform cubic grid cells. [color=#FFFFFF]本文来自:[/color] ?% [color=#FFFFFF]本文来自:[/color] ?%? ? To illustrate the algorithm, an air-filled rectangular cavity [color=#FFFFFF]本文来自:[/color] ?%? ? resonator is modeled.? The length, width, and height of the [color=#FFFFFF]本文来自:[/color] ?%? ? cavity are 10.0 cm (x-direction), 4.8 cm (y-direction), and [color=#FFFFFF]本文来自:[/color] ?%? ? 2.0 cm (z-direction), respectively. [color=#FFFFFF]本文来自:[/color] ?% [color=#FFFFFF]本文来自:[/color] ?%? ? The computational domain is truncated using PEC boundary [color=#FFFFFF]本文来自:[/color] ?%? ? conditions: [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ex(i,j,k)=0 on the j=1, j=jb, k=1, and k=kb planes [color=#FFFFFF]本文来自:[/color] ?%? ? ? ? ? ey(i,j,k)=0 on the i=1, i=ib, k=1, and k=kb plane

    展开全文
  • 有限差分法MATLAB程序!!!!

    万次阅读 多人点赞 2018-10-26 20:48:15
    设矩形金属槽,长为10cm,...这个实验我今天上午电磁场课程的原实验,我课上没做出来,所以下午自己找教程安装了MATLAB 2016A,自己做了一遍各方面还是很有收获。我是一个硬件技术菜鸟,希望能和大家一起持续学习。

    设矩形金属槽,长为10cm,宽为6cm,其侧板和底板接地,盖板电压为5V(相对值)。
    求槽内场(电位)的分布?
    迭代解法:(取最小步长为0.1cm)
    在这里插入图片描述运行效果图如下:在这里插入图片描述

    源代码如下:
    lx=61;ly=101;
    v1=zeros(ly,lx);
    for j=2:lx-1
    v1(ly,j)=5;
    end
    v2=v1;m=1;t=0;k=0;
    while(m>0.00001)
    k=k+1
    m=0;
    for i=2:ly-1
    for j=2:lx-1;
    v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v1(i-1,j)+v1(i,j-1))/4;
    t=(v2(i,j)-v1(i,j));
    if(t>m)m=t;end
    end
    end
    v1=v2;
    end
    subplot(1,2,1),mesh(v2)
    axis([0,61,0,101,0,5])
    subplot(1,2,2),contour(v2,32)

    这个实验我今天上午电磁场课程的原实验,我课上没做出来,所以下午自己找教程安装了MATLAB 2016A,自己做了一遍各方面还是很有收获。我是一个硬件技术菜鸟,希望能和大家一起持续学习。

    展开全文
  • 许多物理现象随着时间而发生变化、如热传导过程、气体扩散过程和波的传播过程都与时间有关。...利用MATLAB有限差分法解这类问题,就是从初始值出发,通过差分格式沿时间增加的方向,逐步求出微分方程的近似解。
  • MATLAB模拟的电磁学时域有限差分法(配套程序代码) 原书配套代码
  • 采用有限差分方程近似地代替偏微分方程对二维场域中的电位函数进行分析,给出了一种新的处理方法,并设计了相应的二维电位函数有限差分法求解的MATLAB应用程序
  • 有限差分法matlab两点边值代码介绍 它是CUDA运行时API的模板库。 开发该库的目的是使用户从执行内存管理,数组大小验证和编写内核函数的常规工作中解放出来,以执行简单的算术运算,因此将精力集中在实现核心算法的...
  • 有限差分法matlab两点边值代码二维椭圆PDE有限差分法可视化 该程序适用于数学软件第四次作业。 要求如下: A和B是学生证中的最大和第二大数字。有限差分法用于求解椭圆PDE方程。 等式在图1中。 该问题属于二维两点...
  • 二维热传导方程有限差分法的分解与计算步骤,最后还附上MATLAB实现程序以及详细解释,是学习偏微分方程以及差分算法的不错的参考材料
  • 3-D FDTD code with PEC boundaries.This MATLAB M-file implements the finite-difference time-domain solution of Maxwell's curl equations over a three-dimensional Cartesian space lattice comprised of ...
  • 采用古典隐式格式对网格进行划分使用matlab实现多介质的有限差分法,包含追赶法,附带10多张程序实现图 二维热传导方程有限差分法的分解与计算步骤,最后还附上MATLAB实现程序以及详细解释,是学习偏微分方程以及差...
  • 本书英文名叫The Finite Difference Time Domain Method for Electromagnetics With MATLAB Simulations,Atef Elsherbeni...中文名《MATLAB模拟的电磁学时域有限差分法》,俞志远翻译。本文件里还有源程序可供参考学习
  • MATLAB模拟的电磁学时域有限差分法》介绍了近年来在电磁领域内发展较快的时域有限差分法(FDTD)的MATLAB语言编程要点,并配有丰富实例。它可作为电类专业高年级本科生或研究生学习时域有限差分法的入门书,也适合对...
  • 精品文档 第五次作业前三题写在作业纸上 一用有限差分方法求解一维非定常热传导方程 初始条件和边界条件见说明 .pdf 文件 热扩散系数 =const 2 T T 2 t x 1. 用 Tylaor 展开推导出 FTCS格式的差分方程 2. 讨论该...
  • MATLAB应用于有限差分法的正演计算中,充分发挥了其强大而方便的功能。通过对二维稳定电流场模型的试算表明,MATLAB在解决实际的工程和数学问题中,与其它计算机程序设计语言C、FORTRAN相比,具有使用更为简便、...
  • 物理过程,都可用椭圆型方程来描述。其中最典型的方程是泊松(Poisson)方程。 本文件给出了使用有限差分法求解泊松方程的两个算例及MATLAB程序
  • 填鬼有限差分法 gFDM 的核心代码是用 C++ 开发的。 一般函数都是用Matlab编写的。 gFDM 实现分为 4 个高级函数和 32 个低级私有函数。 所有函数都以“gfdm_”前缀文本开头。 图 1 说明了高级函数的管道,从 NET ...
  • 有限差分法求电磁场分布-电磁场程序.rar matlab编程求解,有程序清单,要看到话看附件 电磁场程序设计 ——采用有限差分法计算屏蔽微带传输线的电位及其电容 题目 屏蔽微带传输线 当微带线工作在低频时,其...
  • matlab终止以下代码 Vulture:用于电磁仿真的开源FDTD解算器 应用电磁学()应用求解程序Vulture是用于电磁仿真的非均匀结构化网格代码。 它是在上开发的,用于电磁兼容性(),计算电磁学()和的研究。 代码功能 ...
  • - hw1_explict:程序入口 - 网格大小dx,dy,dt设定 - 初值设定 (平板网格编号成一列) - 调用显式求温度矩阵函数 - 调用画图函数 - heatExplict:显式实现温度矩阵求解函数 ----实现方案均来自于教材。简单...
  • 基于MATLAB有限差分法解决二位瞬态导热问题

    千次阅读 多人点赞 2017-01-08 10:05:37
    程序解决的是应对平板的二维导热问题,可以定义...%本文档为有限差分法计算二位导热问题的实例 %为了计算的方便,设定铝板的边界初始为设定初始的值,且温度在传导过程中保持不变 %导热方程为: T_deta/t_deta=a(d(d

空空如也

空空如也

1 2 3 4 5
收藏数 84
精华内容 33
关键字:

有限差分法matlab程序

matlab 订阅