文档介绍:学生实验报告
实验课程名称
开课实验室
学 院_年级专业班
学生姓名 学 号
开课时间_至—学年—学期
总成绩
教师签名
精品文档
课程
名称
数值分析
实验项目
名 称
Gauss数
[m, n]=size(A);
% 先检查系数正确性
if m~=n
error( '矩阵A的行数和列数必须相同');
return ;
end
if m~=size(b)
error( ' b的大小必须和A的行数或A的列数相同');
return ;
end
%再检查方程是否存在唯一解
if ran k(A)~=ra nk([A,b])
error( ' A矩阵的秩和增广矩阵的秩不相同 ,方程不存在唯一解’);
return ;
end
c=n+1;
A(:,c)=b; %(增广)
for k=1: n-1
精品文档
[r,m]=max(abs(A(k: n,k))); %选主元
m=m+k-1; %修正操作行的值
if (A(m,k)〜=0)
if (m~=k)
A([k m],:)=A([m k],:); %换行
end
A(k+1: n, k:c)=A(k+1: n, k:c)-(A(k+1: n,k)/ A(k,k))*A(k, k:c); %消去
end
end
x=zeros(le ngth(b),1); %回代求解
x( n)=A( n,c)/A( n,n);
for k=n-1:-1:1
x(k)=(A(k,c)-A(k,k+1: n)*x(k+1: n))/A(k,k);
end
y=x;
format short ; %设置为默认格式显示 ,显示5位
建立MATLAB界面
利用MATLAB的GUI建立如下界面求解线性方程组:
详见程序
五、 计算实例、数据、结果、分析
下面我们对以上的结果进行测试,求解:
精品文档
2 10 0
3 x1
10
3 4 12
13 x2
5
1 2 3
4 X3
2
4 14 9
13 x4
7
输入数据后点击1斥E 和
GL求罢方程
,得到如下结果:
Gauss直接法求解结果
1
2
3
4
T-
Gauss列主元法求解结果
1
2
3
•L求毎忑忌
斗
Copyright: yfg Tel:**********
更改以上数据进行测试,求解如下1
方程组:
4
3
2
1
X1
1
3
4
3
2
X2
1
2
3
4
3
X3
1
1
2
3
4
X4
1
得到如下结果:
求解线性方程组Ax=b
精品文档
精品文档
请按m atlab中输入矩阵方式输入矩阵A
[4321;3432;2343;1 2 34]
精品文档
精品文档
请按matlab中输入矩阵方弍输入b
[11-1 -1]
Gauss直接法求輕结果
2r56e-0l7
1
-1
-1 WSe-016
创
2775fie-(M7
Gauss列主元法求解结果
1
-1
-13323e-016
六、 实验中遇到的问题及解决办法
在本实验中,遇到的问题主要有两个:
精品文档
(1) 如何将上述的Gauss消元法的算法在 MATLAB中实现
针对此问题我借鉴了网上以及 课本上的算法的 MATLAB实现的程序;
(2) 如何将建立界面使得可以随意输入想要求解的相关矩阵后就可以直接求解
针对此问题,我通过网上的一些关于 MATLAB的GUI设计的相关资料,总结经验完
成了此项任务。
七、 实验结论
通过以上的测试,我们发现以上算法和程序能够求出线性方程组的比较精确解。
八、 参考文献
杨大地,:科学出版社
百度文库,百度知道
教师签名
年 月曰
课程
名称
数值分析
实验项目
名 称
插值方法
实验项目类型
验证
演示
综合
设计
其他
指导
教师
何光辉
成 绩
实验目的:
(1) 学会拉格朗日插值、牛顿插值等基本方法
(2) 设计出相应的算法,编制相应的函数子程序
(3) 会用这些函数解决实际问题
实验内容
(1) 设计拉格朗日插值算法,编制并调试相应的函数子程序