数值分析第三次Project.PDFVIP

数值分析第三次Project 5080309019F0803028韩立 December2,2010 1 实验目的 1.通过FiniteDifferenceMethod近似地求解微分方程(e.g.利用Poisson矩 阵逼近二阶导数) 2.利用Gauss-Seidel和ChasingMethod求解线性方程组 2 实验程序 两种方法都可在Wiki上找到详细描述，故不赘述。其中Gauss-Seidel的程序不 具扩展性，只适用于本次实验的矩阵。程序主要由2部分组成： 2.1 Gauss-Seidel clc; clear; low = 10; high = 100; for n = low : 1: high k = 0; threshold = 1E-10; now = zeros(1, n); b = sin(2 *pi *[1: 1: n] / (n + 1)); b = b . *((2 *pi / (n + 1)) ^ 2); while k = 0 k = k + 1; last = now; for i = 1 : n tmp = 0; if (i 1) tmp = tmp - now(i - 1); end if (i n) 1 tmp = tmp - last(i + 1); end now(i) = (b(i) - tmp) / 2; end if (k 1 norm(now - last, inf) threshold) break; end end iter(n) = k; u = sin(2 *pi *[1:1:n] / (n + 1)); % Error between discrete and continuous err_dc(n) = norm(u - now); end %plot([low: high], iter(low: high)) plot([low: high], err_dc(low: high)) 2.2 ChasingMethod clc; low = 10; high = 100; for n = low: 1: high A = sparse(n, n); A = zeros(n, n); for i = 1 : n A(i, i) = 2; end for i = 1 : n - 1 A(i + 1, i) = -1; A(i, i + 1) = -1; end cond(n) = condest(A); % Ordinary Calculation b = ones(1, n) . *2; a = ones(1, n) . *(-1); c = ones(1, n) . *(-1); d = sin(2 *pi *[1: 1: n] / (n + 1)); d = d . *((2 *pi / (n + 1)) ^ 2);