| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 422 人关注过本帖
标题:[求助]:分析一下这个循环计算程序问题出在哪?
只看楼主 加入收藏
para
Rank: 1
等 级:新手上路
帖 子:2
专家分:0
注 册:2007-5-20
收藏
 问题点数:0 回复次数:0 
[求助]:分析一下这个循环计算程序问题出在哪?

程序编出来了,调试怎么也过不了,希望哪位大虾帮忙调试一下!
u,v,ed,E矩阵需要预定义吗?怎么预定义?


tic;
clear;clc;
epsilon=8.9e-12; num=21; a=0.2/(num-1); t=(num+1)/2;cons=300;
dp1=0.1;pg=1.025;ki=2.2e-4;
v(t,t)=6e4;v1=u;N=0;
while max(max(abs(v-v1)))>cons;
N=N+1;
v1=v;
for i=2:num-1
for j=2:num-1
Eox=-(v(i-1,j)-v(i+1,j))/(2*a);
Eoy=-(v(i,j-1)-v(i,j+1))/(2*a);
E(i,j)=sqrt(Eox.^2+Eoy.^2);
u(i,j)=0.25*(v(i,j-1)+v(i,j+1)+v(i-1,j)+v(i+1,j)+a^2*ed(i,j)/epsilon);
ed(i,j)=-epsilon/(2*a)*(Eox+Eoy)-0.5*((epsilon/a*(Eox+Eoy))^2
+4*(epsilon/a*(Eox*ed(i+1,j)+Eoy*ed(i,j+1))))^0.5;
u(1,:)=0;u(num,:)=0;u(:,num)=0;u(t,t)=6e4;ed(t,t)=3e-5;
end;
end;
v=u;
end;
d=[0:0.2/(num-1):0.2];
[x,y]=meshgrid(d);
u=abs(u);ed=abs(ed);E=abs(E);
figure(1);surfl(x,y,abs(u));
figure(2);surfl(x,y,abs(ed));
figure(3);surfl(x,y,abs(E));
toc


[此贴子已经被作者于2007-5-20 15:55:23编辑过]

2007-05-20 10:49
快速回复:[求助]:分析一下这个循环计算程序问题出在哪?
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.025560 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved