1270阅读
3回复

求助编程问题，急急急急急急

UID ：18917

注册：2008-10-10
登录：2014-03-27
发帖：134
等级：仿真二级

0楼发表于: 2009-03-17 22:58:36

— 本帖被 destroyer 从【互助速问速答有求速应】移动到本区(2009-11-28) —

clc
clear
%---------------需要用到的常量------------------%
k=120*pi; %真空波阻抗
l=1/k;
V=1;    %外导体加的电压值
%------------------波导几何尺寸------------------%
b1=5e-3; %外导体的高度
a1=2e-3;
m=b1-a1; %内导体的高度
%---------网格划分----------%
h=0.25e-3; %网格尺寸
Nx=2*b1/h+1; %x方向网格节点数
Ny=2*b1/h+1; %y方向网格节点数
N=(b1-a1)/h+1;
u= zeros(Nx,Ny); %定义数组
%------------------赋初值------------------%
i= 2:Nx;
u(i,Ny)= 1;
j= 2:Ny;
u(Nx,j)= 1;
i=2:N;
u(i,N)= 0;
j= 2:N;
u(N,j)= 0;
for i=2:Nx-1
u(i,1)=0.25.*(u(i+1,1)+u(i-1,1)+2.*u(i,2));
u(1,i)=0.25.*(u(1,i+1)+u(1,i-1)+2.*u(2,i));
end
p=1.17; %超松弛因子
w=1000;%迭代次数

%------------高斯赛德尔迭代求电位分布------------%
for n= 1:w
for i=2:Nx-1
for j= 2:Ny-1
a= u( i,j);
b= u( i,j+1);
c= u(i+1,j);
d= u( i-1,j);
e=u(i ,j-1);
f= 0.25.*(b+ c+ d +e);
u(i,j )= a+ p.*(f-a);
end
end
end
u   ; % 第一象限电位分布u
%-----利用对称性求出整个区域的电场分布u4------%
for i=1:Nx
u1(i,:)=u(Nx+1-i,:);
end
u2=[u1
u];
for i=1:Nx
u3(:,i)=u2(:,Nx+1-i);
end
u4=[u3 u2];
u4 ;
u4(:,41)=[];
u4(41,:)=[];
%--------求出各个点的电场和磁场分布--------%
[DX,DY]=gradient(u4);
Ex=-DX;
Ey=-DY;
Hx=l*DY;
Hy=-l*DX;
改程序是求解方同 ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册