原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值实验实验二
数值实验一
实验1 电阻网络问题的求解
1、实验目的：Gauss消去法是求解低阶线性方程组的基本方法。理解Gauss消去法的原理和解题步骤，掌握该方法是计算实践的第一步；解对数据误差的敏感性也是计算时间中必须注意的问题。
2、实验步骤和要求
（1）、写出方程组的具体形式，注意其稀疏特点，非零元集中在主对角线的附近；
方程组的具体形式如下：
A矩阵为：
向量V为：
V= [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 5 0 5]’
（2）、采用高斯消元法的程序如果
function x=gauss(A,b);
na=size(A,2);
nb=size(b,1);
C=[A,b];
if na~=nb
error(A 和 b必须有一致的维数 );
end
%消元过程
for i=1:1:na-1
for j=i+1:1:na
C(j,:)=C(j,:)-C(j,i)/C(i,i)*C(i,:);
end
end
A=C(:,1:na);
b=C(:,na+1);
%回带得方程的解
b1=b;
for i=na:-1:1
x(i)=b1(i)/A(i,i);
if i1
for j=i-1:-1:1
b1(j)=b1(j)-A(j,i)*x(i);
end
end
end
有x=Gauss（A，V）得方程组的解为：
X1= [2.5000 2.1512 2.5000 2.8488 2.5000 2.8488 3.3721 2.2093 2.5000 2.7907 2.5000 2.8488 3.3721 2.7907 3.3721]’
（3）、假设5-14,13-14,15_14,23-14的4条线路上的点阻有4%的增加，
则A矩阵变为：
此时
V= [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 5 0 5]’
由x=Gauss（A，V）
则此时方程组的解为：
X2= [2.5900 2.2481 2.6662 2.9458 2.5900 2.9458 3.4413 2.2924 2.5900 2.8738
2.5900 2.9458 3.4413 2.8738 3.4413]’
由X2-X1=[ 0.0900 0.0970 0.1662 0.0970 0.0900 0.0970 0.0693 0.0831 0.0900 0.0831 0.0900 0.0970 0.0693 0.0831 0.0693]’
由此可知稀疏矩阵的几个元素很小的变化，带来解的很大的变化。
实验2.2 时间序列模型的求解
实验目的 求解特殊结构矩阵的方程，要尽量利用矩阵的结构特征，使计算的效率和精度都能提高。
本实验室讨论时间序列模型的辨识问题，得到Toeplitz方程组
（1）、该矩阵是对称矩阵。同时其特征值都为正数，即其为正定矩阵
（2）、用计算机随机生成数 ，以及计算系数矩阵的条件数的程序如下：
clear;
%产生随机数的程序
a=rand(1,30);
for i=1:1:29
for j=1:1:(30-i)
if(a(j)a(j+1))
temp=a(j+1);
a(j+1)=a(j);
a(j)=temp;
end
end
end
%求系数矩阵条件数的程序
a1=cond(a);
由于是随机产生的，所以每次运行程序的结果都是不同的。
（3）、顺序Gauss消元法的求解程序如下：
clear all;
close all;
a=rand(1,30);
%对a进行从大到小排序
for i=1:1:29
for j=1:1:(30-i)
if(a(j)a(j+1))
temp=a(j+1);
a(j+1)=a(j);
a(j)=temp;
end
end
end
%构造矩阵B
for i=1:1:30
A(i,i)=1;
if i1
for k=1:1:i-1
A(i,k)=a(i-k)
end
end
for j=
文档评论（0）