matlab差分法求解ppt,Matlab实现电磁场数值计算.ppt

Matlab实现电磁场数值计算.ppt

* 第 2 章 * Matlab在电磁场中的几个应用 Matlab 实现电磁场数值计算 均匀带电线段的电场 均匀带电半圆环圆心的电场 Matlab直观描述复杂的物理现象 均匀带电线段的电场和电势分布 电偶极子的电场和电势分布 有限差分法 有限元法 边界元法 矩量法 积分方程法 积分法 分离变量法 镜像法、电轴法 微分方程法 保角变换法 计算法 实验法 解析法 数值法 实测法 模拟法 边 值 问 题 计算电磁场问题的方法 数值法 解析法 分离变量法：严格求解偏微分方程的经典方法 变换数学法：严格求解积分方程的方法 解析法的优点是： ①可将解答表示为已知函数的显式，从而计算出精确的数值结果； ②可以作为近似解和数值解的检验标准； ③在解析过程中和在解的显式中可以观察到问题的内在联系和各个参数对数值结果所起的作用。 但解析法也存在缺点，主要是：它仅能解决很少量的问题，事实上，只有在为数不多的坐标系中才能分离变量，而用积分方程法时往往求不出结果，致使分析过程既困难又复杂。 数值法与解析法比较，在许多方面具有独特的优点。 ①普适性强，用户拥有的弹性大。一个特定问题的边界条件、电气结构、激励等特性可以不编入基本程序，而由用户输入，更好的情况是通过图形界面输入。 ②用户不必具备高度专业化的电磁场理论、数学及数值技术方面的知识就能用提供的程序解决实际问题。 数值法的缺点是数据输入量大、计算量大、受硬件条件的限制。 电磁场分布型问题的数值积分法 梯形法 dz 当 均匀带电线段的电场 真空中有一长为L的均匀带电直导线，带电量为Q ,试求中垂线上P 点的电场。 解析法： 解析解 Matlab实现电磁场数值计算 数值法： 数值解 解析解 参数取值： 棒长L=10,棒的总电荷Q=1,N=100 距离ρ 数值解E 精确解E1 1 1765048479.328 1765045216.2437 5 254561623.1613 254558441.2272 10 1896 19 20 6779 2533 30 9863952.2078 9863939.2383 40 5581567.284 5581563.0565 50 3582135.6405 3582133.8848 100 898877.2168 898877.105 600 24999.1321 24999.132 1000 8999.8875 8999.8875 >> 参数取值： 棒长L=10,棒的总电荷Q=1,N=1000 距离ρ 数值解E 精确解E1 1 1765045248.8813 1765045216.2437 5 254558473.047 254558441.2272 10 6299 19 20 8575 2533 30 9863939.368 9863939.2383 40 5581563.0988 5581563.0565 50 3582133.9023 3582133.8848 100 898877.1061 898877.105 600 24999.132 24999.132 1000 8999.8875 8999.8875 N越大，子区间越窄，数值解和解析解越接近； 当距离ρ大于30后，数值解与解析解十分接近，能够满足电磁场工程的一般要求。 Matlab程序设计 NMAX=1000; Q=1; L=10; p=Q/L; aa=[1,5,10,20,30,40,50,100,600,1000,0]; eplot1=aa eplot2=aa E0=(1/36/pi)*1e-9; disp(' ') disp('参数取值：') disp([ ‘ ','棒长L=',num2str(L),',','棒的总电荷Q=',num2str(Q),',','N=',num2str(NMAX)]) disp(' ') disp([ ‘ ','距离a',' ','数值解E',' ','精确解E1']); i=1; while(1); a=aa(i); if a<=