急求助：帮我分析这个程序问题出在哪

para

程序编出来了，调试怎么也过不了，希望哪位大虾帮忙调试一下！
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)=3e5;
        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

spano

u v 要定义，E不用

xjzuo

请把问题讲清楚一下.
还有: v(t,t)=6e4; ------是对角元,还是(t,t)这个矩阵元,还是整个v被赋值为6e4?

para

原帖由 spano 于 2007-5-20 11:26 发表

                               
登录/注册后可看大图

u v 要定义，E不用

该怎么定义一下，我不太会，谢谢！

para

原帖由 xjzuo 于 2007-5-20 14:49 发表

                               
登录/注册后可看大图

请把问题讲清楚一下.
还有: v(t,t)=6e4; ------是对角元,还是(t,t)这个矩阵元,还是整个v被赋值为6e4?

v(t,t)=6e4是(t,t)这个矩阵元，相当于边界条件。

xjzuo

"边界条件"就一个点? 其它矩阵元呢? ----至少也应该有个初值...

para

原帖由 xjzuo 于 2007-5-20 17:11 发表

                               
登录/注册后可看大图

"边界条件"就一个点? 其它矩阵元呢? ----至少也应该有个初值...

u(1,:)=0;u(num,:)=0;u(:,num)=0;u(t,t)=6e4;这些都是边界条件啊！
今天用估计的u值带进去了，结果却是NANI，问题在哪？

xjzuo

你的程序能运行?  都是些明显的错误...
看来我不应该过问这种问题完全没有讲清楚的帖子...

eight

原帖由 xjzuo 于 2007-5-21 17:41 发表

                               
登录/注册后可看大图

你的程序能运行?  都是些明显的错误...
看来我不应该过问这种问题完全没有讲清楚的帖子...

呵呵，我一般遇到这种情况也是比较无奈，只能给出“请阅读matlab基础书”的回复

para

原帖由 xjzuo 于 2007-5-21 17:41 发表

                               
登录/注册后可看大图

你的程序能运行?  都是些明显的错误...
看来我不应该过问这种问题完全没有讲清楚的帖子...

基本思路是用u,ed相互迭代。
调试通过了，运行结果不收敛

[ 本帖最后由 para 于 2007-5-21 21:16 编辑 ]