用LU分解求解线性方程组数值 课程设计毕业设计（论文）.docVIP

目 录 实验一：用LU分解求解线性方程组 1 1.1实验目的与要求 1 1.2原理与程序 1 1.2.1算法原理 1 1.2.2程序 2 1.3实例 3 1.4算法总结 5 实验二：Lagrange插值 5 2.1实验目的与要求 5 2.2原理与程序 5 2.2.1算法原理 5 2.2.2程序 6 2.3实例 6 2.4算法总结 8 实验三：变步长的梯形积分公式 8 3.1实验目的与要求 8 3.2原理与程序 9 3.2.1算法原理 9 3.2.1程序 9 3.3实例 10 3.4算法总结 12 实验四：常微分方程数值解(Runge-Kutta方法) 12 4.1实验目的与要求 13 4.2原理与程序 13 4.2.1算法原理 13 4.2.2程序 13 4.3实例 14 4.4算法总结 15 实验五：数值方法实际应用 16 5.1实验目的与要求 16 5.2实例 16 5.2.1具体问题 16 5.2.2问题假设 16 5.2.3建立模型 17 5.2.4模型求解 17 5.2.5模型分析 21 5.3模型总结 22 参考文献 22 实验一：用LU分解求解线性方程组 1.1实验目的与要求 1.熟悉LU分解求解线性方程组的基本原理 2.了解LU分解求解线性方程组的计算流程 3.能编程实现LU分解并求解线性方程组 由《数值分析与实验学习指导》第31页定理2.1有：只有矩阵A的各阶顺序主子行列式都不等于零时才能用LU分解。 LU分解思想如下： 所以： a11= u11,a12=u12···, a1n= u1n ; u11=a11,u12=a12, ···,u1n=a1n a21 = l21u11, ··· , an1 = ln1u11 ; l21=a21/u11,··· , ln1=an1/u11 a22=l21u12+ u22, ··· , a2n =l21u1n+ u2n ; u22=a22 - l21u12,··· , u2n=a2n- l21u1n a32=l31u12+ l32u22, ··· , an2=ln1u12+ ln2u22 ; l32=(a32- l31u12)/u22, ··· , ln2=(an2- ln1u12)/u22 对A的元素aij ,当 j≥k 和 i≥k时 矩阵L和矩阵U的元素计算公式 function [L,U,x]=lux(A,b) [n,n]=size(A); %LU 分解法解线性方程组（列主元LU分解） for k=1:n-1 for i=k+1:n m=A(i,k)/A(k,k); for j=k+1:n A(i,j)=A(i,j)-m*A(k,j); end A(i,k)=m; end end A %输出紧凑型LU矩阵 %解下三角矩阵 Ly=b for i=2:n s=0; for j=1:i-1 s=s+A(i,j)*b(j); end b(i)=b(i)-s; end b %输出中间向量y %解上三角方程组 Ux=y b(n)=b(n)/A(n,n); for i=n-1:-1:1 s=0; for j=i+1:n s=s+A(i,j)*b(j); end b(i)=(b(i)-s)/A(i,i); end b例1：用LU分解法求解性线方程组： A=[2,3,4;3,5,2;4,3,30]; b=[6;5;32]; lux(A,b) 回车得到： A = 2.0000 3.0000 4.0000 1.5000 0.5000 -4.0000 2.0000 -6.0000 -2.0000 b = 6 -4 -4 b = -13 8 2 故原方程组的解为：x=[-13,8,2]T 解：在matlab中新建m文件，输入如下代码，并保存为LU1.m function [L,U,x]=lU(A,b) A=[10,-7,0;-3,2,6;5,-1,5]; %系数矩阵 b=[7;4;6];% [n,r]=size(A); for k=1:n-1 %LU 分