椭圆型方程matlab程序

椭圆方程求解 各种方法 李荣华版 课后作业

function [xkq]=CG（AbeN）
%e为精度要求
%N为最大迭代次数
%x0为迭代初值
n=length（A）;
x=sparse（zeros（n1））;
k=0;
r=b-A*x;
p=r;
L1=（norm（r2））^2;
if L1    q=1;%表示迭代成功

end
for k=1:N
    q=A*p;
    alfa=L1/dot（pq）;%dot表示向量内积函数
    x=x+alfa*p;
    r=r-alfa*q;
    L2=（norm（r2））^2;
    if L2        q=1;
        break;
    end
    beita=L2/L1;
    L1=L2;
    p=r+beita*p;
end
if k>N
q=0;%表示迭代失败
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     706680  2012-03-24 22:17  solve_possion-Mixed\1.1

     文件     706680  2012-03-24 22:18  solve_possion-Mixed\1.2

     文件     706680  2012-03-24 22:16  solve_possion-Mixed\2.tif

     文件     706680  2012-03-24 22:18  solve_possion-Mixed\3.tif

     文件        454  2012-03-24 12:13  solve_possion-Mixed\CG.m

     文件        245  2012-03-27 18:46  solve_possion-Mixed\data.m

     文件        531  2012-03-24 11:45  solve_possion-Mixed\GaussS.m

     文件        848  2012-03-24 12:19  solve_possion-Mixed\gseid.m

     文件        465  2012-03-24 22:15  solve_possion-Mixed\Jacobi1.m

     文件        320  2012-03-27 18:43  solve_possion-Mixed\QU.m

     文件        104  2012-03-16 17:07  solve_possion-Mixed\rhs_f.m

     文件       3022  2012-03-27 18:37  solve_possion-Mixed\solve_poisson.asv

     文件       3053  2012-03-27 18:47  solve_possion-Mixed\solve_poisson.m

     文件       3310  2012-03-24 22:09  solve_possion-Mixed\solve_possion1.asv

    ..A.SH.      5632  2012-03-26 22:13  solve_possion-Mixed\Thumbs.db

     文件        104  2012-03-16 17:06  solve_possion-Mixed\true_solution.m

     文件        202  2011-05-09 16:36  solve_possion-Mixed\zuixiao.m

     文件      13824  2012-03-27 18:23  solve_possion-Mixed\偏微分方程实验报告.doc

     文件      25088  2012-03-16 17:06  solve_possion-Mixed\问题描述.doc

     目录          0  2012-03-27 18:35  solve_possion-Mixed

----------- ---------  ---------- -----  ----

              2883922                    20