题目:用下列方法解线性方程组。
并比较计算结果精度(方程组准确解为x1=x2=x3=x4=1)
1) 顺序消元法;
2) 列主元消元法;
实验报告。1、程序设计:
我们的实验目的是分别通过高斯消元法和列主元消元法两种方法求同一线性方程组的解,并比较它们的结果的精度 ,故在此以题目中线性方程组的系数增广矩阵为参数,分别编写两段程序使之实现:
function c = mygauss(a,n) %高斯消元法求解线性方程组。
function c = lzy(a,n) %列主元消元法求解线性方程组。
%a(系数增广矩阵),n(系数矩阵的秩)
2、源程序:
%高斯消元法
function x=mygauss(a,n)
for k=1:n-1
for i=k+1:n
a(i,k)=str2num(sprintf('%4f',a(i,k)/a(k,k)))
for j=k+1:n
a(i,j)=str2num(sprintf('%4f',a(i,j)-a(i,k)*a(k,j)))
enda(i,n+1)=str2num(sprintf('%4f',a(i,n+1)-a(i,k)*a(k,n+1)))
endend得到上三角系数矩阵。
x=zeros(n,1回带求解过程。
x(n)=str2num(sprintf('%4f',a(n,n+1)/a(n,n)))
for k=n-1:-1:1
s=a(k,n+1);
for j=k+1:n
s=s-a(k,j)*x(j);
endx(k)=str2num(sprintf('%4f',s/a(k,k)))
end%列主元消元法:
function x=lzy(a,n)
for k=1:n-1调整行过程。
for b=n-1:-1:1
for a=k:b
if abs(a(a+1,k)) abs(a(a,k))
c=a(a+1,:)a(a+1,:)a(a,:)a(a,:)c;
end;end;
end;for i=k+1:n
a(i,k)=str2num(sprintf('%4f',a(i,k)/a(k,k)))
for j=k+1:n
a(i,j)=str2num(sprintf('%4f',a(i,j)-a(i,k)*a(k,j)))
enda(i,n+1)=str2num(sprintf('%4f',a(i,n+1)-a(i,k)*a(k,n+1)))
endend得到上三角系数矩阵。
x=zeros(n,1回带求解。
x(n)=str2num(sprintf('%4f',a(n,n+1)/a(n,n)))
for k=n-1:-1:1
s=a(k,n+1);
for j=k+1:n
s=s-a(k,j)*x(j);
endx(k)=str2num(sprintf('%4f',s/a(k,k)))
end3、实验结果:
%command window :
>a=[1.1348,3.8326,1.
1651,3.4017,9.5342;0.
5301,1.7875,2.5330,1.
5435,6.3941; 3.4129,4.
9317,8.7643,1.3142,18.
4231;1.2371,4.9998,10.
6721,0.0147,16.9237];
>myguass(a,4
高斯消元法。ans =
列主元消去法。
>lzy(a,4)ans =
%把以上两程序中的“.4f”改为“.3f”运行,则结果变化为如下:
高斯消元法。
>myguass(a,4ans =
列主元消去法。
>lzy(a,4)ans =
%把以上两程序中的“.3f”再改为“.6f”运行,则结果变化为如下:
高斯消元法。
>myguass(a,4) ans =
列主元消去法。
>lzy(a,4)ans =
%分析:观察题目中的系数矩阵,发现a(4,4)为0.0147,很小接近于零,会在高斯消元法进行时引入误差。
由我们得到的结果可以看出,一旦在运算过程中保留的小数位数较少,高斯消元法得到的结果误差很大,而用列主元消元法通过不断调整系数矩阵行排列,把较大的数置于前排,可以很好地规避误差,得到的结果更加精确可靠。
而随着计算过程保留小数位数的增加,精度可以保持到a(4,4)的数量级,此时我们也可以通过高斯消元法得到比较可靠的答案。
计算机网络作业
物理层作业。一 判断题 针对下面的描述,对的打 错的打 x 1.传输速率单位 bps 代表bytes per second2.rs 232c的电气特性规定逻辑 0 的电平电压为0至 5伏。二 单选题 在下面的描述中只有一个符合题意,请用a b c d标示之 1.rs 232c的电气特性规定逻辑 0 ...
计算机网络作业
计算机网络 第2次作业。注意 作业完成后在规定时间内提交给老师,提交时在 文件名称 和 注释 处也要写清楚 学习中心,学号,姓名,入学时间 以便登记你的作业成绩。注 标 的书后有答案。一 单项选择题 本大题共20小题,每小题2分,共40分 1.ip 地址实际上由网络地址和主机地址两部分组成,其中主机...
计算机网络 作业
1 列出5种非专用因特网应用及他们所使用的应用层协议。答 the web http file transfer ftp remote login telnet network news nntp email smtp.2 网络体系结构和应用程序体系结构之间有什么区别?答 网络体系结构是固定的,并为应...