作者profyang (prof)
看板MATLAB
标题Re: [讨论] 电势模拟问题
时间Mon Dec 11 11:50:39 2017
※ 引述《knutten (biker)》之铭言：
: 大家好, 小弟是matlab新手, 正在努力学matlab中.
: 我在一本书看到以下的程式, 但是书上的例子是设定上面的电压是100V, 两侧和底部的电压是0V.
: 但是我想要模拟成上面电压1500V, 两侧和底部的电压是700V, 可是好像得不到答案,
: 所以只能将上部取800V, 两侧和底部的电压是0V, 用压差来模拟, 请问有没有什麽方法可以
: 模拟成上面电压1500V, 两侧和底部的电压是700V?
这不就FDM? 重新定义一下边界条件应该就可以了吧?
v0=1500;hx=11;hy=9;
v00=700;
: v1=zeros(hy,hx);
: v1(hy,:)=ones(1,hx)*v0;
: v1(2:hy-1,2:hx-1)=ones(hy-2,hx-2);
%原文上面这行我其实不知道有啥意义 只是代个起始值而已用0和用1应该没差很多
v2=v00*ones(hy,hx);maxt=1;t=0; %这行这样才会定义到v2 700V那些边的边界条件
: v2(hy,:)=v1(hy,:);
: while (maxt>0.1)
: for i=2:hy-1
: for j=2:hx-1
: v2(i,j)=(v1(i,j-1)+v1(i,j+1)+v1(i-1,j)+...
: v1(i+1,j))/4;
: t=v2(i,j)-v1(i,j);
: maxt=0;
: if(t>maxt) maxt=t; end
: end
: end
: v1=v2;
: end
: subplot(1,2,1),surf(v2)
: axis ([0,11,0,9,1,800])
: subplot(1,2,2),
: contour(v2);
: hold on;
: x=1:1:11;y=1:1:9;
: [xx,yy]=meshgrid(x,y);
: [Gx,Gy]=gradient(v2,0.6,0.6);
: quiver (xx,yy,Gx,Gy,'r'),
: hold off
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 114.37.174.96
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1512964242.A.7F7.html
※ 编辑: profyang (114.37.174.96), 12/11/2017 11:51:28
※ 编辑: profyang (114.37.174.96), 12/11/2017 11:52:44
1F：推 knutten: 感谢你的回答，FDM没学过，matlab助教也只教了回圈，其他 12/17 06:52
2F：→ knutten: 要我们自己google, 所以很多不懂 12/17 06:52