内容简介作者简介 - 科学与工程计算系

18.0803
(5)
opts.AbsTol=1e-6;opts.RelTol=1e-3;
s2= bvp4c(@mat4ode,@mat4bc,s1,opts);
Lambda2=s2.parameters
The solution was obtained on a mesh of 39 points.
The maximum residual is 4.915e-004.
There were 1443 calls to the ODE function.
There were 52 calls to the BC function.
Lambda2 =
17.0973
(6)
plot(s0.x,s0.y(1,:),'ks--',s1.x,s1.y(1,:),'bo:',s2.x,s2.y(1,:),'r*-')
legend('\fontname{隶书}\fontsize{16}猜测解','第一近似解','第二近似解',0)
axis([0,pi,-1,1.2]),xlabel('x'),ylabel('solution y')
solution y
1
0.8
0.6
0.4
0.2
0
-0.2
-0.4
-0.6
-0.8
4.15 稀疏矩阵
-1
0 0.5 1 1.5
x
2 2.5 3
4.15.1 稀疏矩阵的存储方式
4.15.2 稀疏矩阵的创建
图 4.14-6
4.15.2.1 稀疏矩阵创建指令:sparse
4.15.2.2 稀疏带状矩阵创建指令:spdiags
41
猜测解
第一近似解
第二近似解
4.15.2.3 外部数据转换为稀疏矩阵的指令:spconvert
【例 4.15.2.3-1】用两种不同方式创建三对角稀疏矩阵。
n=5;SM1=sparse(1:n,1:n,-2*ones(1,n),n,n,n);