原创力文档- 1、本文档共21页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
实验一 牛顿下山法
一、 实验目的：
1、 掌握牛顿下山法求解方程根的推导原理。
2、 理解牛顿下山法的具体算法与相应程序的编写。
二、 实验内容：
采用牛顿下山法求方程2x3-5x-17=0在2附近的一个根。
三、 实验实现：
1、 算法：
下山因子从开始，逐次将减半进行试算，直到能使下降条件成立为止。再将得到的循环求得方程根近似值。
2、 程序代码如下：
function [p,k]=NewtonDownHill(f,df,p0)
N=2000;Tol=10^(-5);e=10^(-8);
for k=1:N
lamda=1;
p1=p0-lamda*f(p0)/df(p0);
while (abs(f(p1)) = abs(f(p0)) lamdae)
lamda=lamda/2;
p1=p0-lamda*f(p0)/df(p0);
end
if abs(p1-p0)Tol break
end
p0=p1;
end
ans=p1
3、 运行结果：
四、 实验体会：
牛顿下山法可以较快求的方程结果，对于该题，只需要5步。运用计算机的数值迭代法可以很快求得满足精度要求的结果。
实验二 矩阵的列主元三角分解
一、 实验目的：
学会矩阵的三角分解，并且能够用MATLAB编写相关程序，实现矩阵的三角分解，解方程组。
二、 实验内容：
用列主元消去法求解方程组（实现PA=LU) 要求输出：
（1）计算解X;
（2）L,U;
（3）正整型数组IP(i),(i=1,···,n) (记录主行信息)。
三、 实验实现：
1、算法：
列主元三角分解和普通Dooliitle分解不同，第k步分解时为了避免用绝对值很小的数作除数，
引进量
若，则将矩阵的第t行与第k行元素互换，再进行正常的Doolittle分解。
2、程序代码如下：
clear all;clc;
A=[1 1 1 1 1 1 1;
2 1 1 1 1 1 1;
3 2 1 1 1 1 1;
4 3 2 1 1 1 1;
5 4 3 2 1 1 1;
6 5 4 3 2 1 1;
7 6 5 4 3 2 1];
b=[7;8;10;13;17;22;28];
n=length(A);
IP=eye(n);
U=zeros(n);
L=eye(n);
[m,p]=max(A(:,1));
C1=b(1);b(1)=b(p);b(p)=C1;
C2(1:n)=IP(1,:);IP(1,:)=IP(p,:);IP(p,:)=C2(1:n);
C2(1:n)=A(1,:);A(1,:)=A(p,:);A(p,:)=C2(1:n);
U(1,:)=A(1,:);
for k=2:n
L(k:n,k-1)=A(k:n,k-1)/U(k-1,k-1);
j=1;
for i=k:n
S(j)=A(i,k)-L(i,1:k-1)*U(1:k-1,k);
j=j+1;
end
[m,p]=max(abs(S(:))) ;
C2(1:n)=A(k,1:n);A(k,1:n)=A(k+p-1,1:n);A(k+p-1,1:7)=C2(1:n);
C2(1:n)=IP(k,1:n);IP(k,1:n)=IP(k+p-1,1:n);IP(k+p-1,1:n)=C2(1:n);
U(k,k:n)=A(k,k:n)-L(k,1:k-1)*U(1:k-1,k:n);
if kn
A(k+1:n,k)=A(k+1:n,k)-L(k+1:n,1:k-1)*U(1:k-1,k);
end
end
disp(L=);disp(L);
disp(U=);disp(U);
disp(IP=);disp(IP);
y=L\b;
x=U\y;
disp(y=);disp(y);
disp(x=);disp(x);
4、 运行结果：
四、 实验体会：
列主元三角分解原理很简单，用MATLAB实现起来却不容易。对MATLAB还需要深入学习。
实验三 SOR迭代法
一、 实验内容：
采用SOR方法求解方程
采用0向量为初始向量，迭代以结束，但迭代次数不操作1000次，松弛因子使用 其中，给出使用的各松弛因子及对应迭代次数。
二、 实验实现：
1、程序代码如下：
A=[ 10 -1 0;
-1 10 -2;
0 -2 10];
b=[9;7;6];
x0=[0;0;0];e=10^(-6);
文档评论（0）