实 验 报 告
实验课程:MATLAB数值分析
学生姓名:
学 号: 20103010
学 院:
专 业:
指导老师:
20##年12月28日
一、问题描述
考虑如下Dirichlet问题
其中为正方形区域的边界,类似于模型问题,我们得到差分方程
得到系数矩阵
A=其中B=-I,I为n-1阶单位矩阵S’为对角元为1+次对角元为的n-1阶的对称三对角矩阵。对=sin(xy),=,n=20.用超松弛以及共轭梯度求解差分方程的解。
二、Matlab代码及其结果
Gauss-Seidel方法:
n=20;
h=1/n;
a=h*h/4;
m=n*n;
S=zeros(n-1,n-1);
for i=1:n-1
S(i,i)=1+h*h/4;
end
for i=2:n-1
S(i,i-1)=-1/4;
end
for i=1:n-2
S(i,i+1)=-1/4;
end
B=zeros(n-1,n-1);
for i=1:n-1
B(i,i)=-1/4;
…… …… 余下全文