n=4;%确定需要LU分解的矩阵维数

%A=zeros(n,n);

L=eye(n,n);P=eye(n,n);U=zeros(n,n);%初始化矩阵

tempU=zeros(1,n);tempP=zeros(1,n);%初始化中间变量矩阵

A=[1 2 -3 4;4 8 12 -8;2 3 2 1;-3 -1 1 -4];%需要LU分解矩阵赋值

for p=1:n %将A矩阵赋值给U

for q=1:n

U(p,q)=A(p,q);

end

end

jt=1;kt=0;

for i=1:n-1

jt=jt+1;

kt=kt+1;

ii=U(i,i);

if ii==0 %主元为零，进行行变换

for m=i:n

if U(m,i)~=0

tempU=U(i,:);

U(i,:)=U(m,:);

U(m,:)=tempU;

ii=U(i,i);

%%

tempP=P(i,:); %行变换结果存储在P中

P(i,:)=P(m,:);

P(m,:)=tempP;

break;

end

end

%disp(ii);

end

disp(ii);

for j=jt:n %%两重循环，完成高斯消元

perj=U(j,i)/ii;

L(j,i)=perj;

for k=kt:n

U(j,k)=U(j,k)-perj*U(i,k);

end

end

end

savefile='LUdapart';

save(savefile)

矩阵LU分解的MATLAB与C++实现

一:矩阵LU分解 矩阵的LU分解目的是将一个非奇异矩阵\(A\)分解成\(A=LU\)的形式,其中\(L\)是一个主对角线为\(1\)的下三角矩阵:\(U\)是一个上三角矩阵. 比如\(A= \beg ...

矩阵LU分解分块算法实现

本文主要描述实现LU分解算法过程中遇到的问题及解决方案,并给出了全部源代码. 1. 什么是LU分解? 矩阵的LU分解源于线性方程组的高斯消元过程.对于一个含有N个变量的N个线性方程组,总可以用高斯消去 ...

计算方法 -- 解线性方程组直接法(LU分解、列主元高斯消元、追赶法)

#include #include #include #include

MATLAB矩阵的LU分解及在解线性方程组中的应用

作者:凯鲁嘎吉 - 博客园http://www.cnblogs.com/kailugaji/ 三.实验程序 五.解答(按如下顺序提交电子版) 1.(程序) (1)LU分解源程序: function [ ...

matlab 求解线性方程组之LU分解

线性代数中的一个核心思想就是矩阵分解,既将一个复杂的矩阵分解为更简单的矩阵的乘积.常见的有如下分解: LU分解:A=LU,A是m×n矩阵,L是m×m下三角矩阵,U是m×n阶梯形矩阵 QR分解: 秩分解 ...

matlab实现高斯消去法、LU分解

朴素高斯消去法: function x = GauElim(n, A, b) if nargin < 2 for i = 1 : 1 : n for j = 1 : 1 : n A(i, j) ...

线性代数笔记10——矩阵的LU分解

在线性代数中, LU分解(LU Decomposition)是矩阵分解的一种,可以将一个矩阵分解为一个单位下三角矩阵和一个上三角矩阵的乘积(有时是它们和一个置换矩阵的乘积).LU分解主要应用在数值分析 ...

矩阵分解---QR正交分解，LU分解

相关概念: 正交矩阵:若一个方阵其行与列皆为正交的单位向量,则该矩阵为正交矩阵,且该矩阵的转置和其逆相等.两个向量正交的意思是两个向量的内积为 0 正定矩阵:如果对于所有的非零实系数向量x ,都有 x ...

第五节、矩阵分解之LU分解

一.A的LU分解:A=LU 我们之前探讨过矩阵消元,当时我们通过EA=U将A消元得到了U,这一节,我们从另一个角度分析A与U的关系 假设A是非奇异矩阵且消元过程中没有行交换,我们便可以将矩阵消元的EA ...

随机推荐

bzoj 4695&colon; 最假女选手

……一道丧病线段树膜板题…… 被常数卡的死去活来……QAQ 学到了些奇技淫巧:把取min标记 和 区间最小值 合并 可以快很多…… #include #de ...

远方的塔--Pylons

转自:https://en.wikipedia.org/wiki/Pylons_project#Pylons_Framework Pylons

NSBundle

属性: .使用类方法创建一个NSBundler对象 + (NSBundle *)mainBundle; eg:[NSBundle mailBundle]; .使用路径获取一个NSBundle 对象,这 ...

FIM相关报错汇总

1.错误1:FIM在修改MA名字的时候报错: The management agent cannot be deleted or renamed because the working directo ...

Greedy&colon;Packets&lpar;POJ 1017&rpar;

装箱问题1.0 题目大意:就是一个工厂制造的产品都是正方形的,有1*1,2*2,3*3,4*4,5*5,6*6,高度都是h,现在要包装这些物品,只能用6*6*h的包装去装,问你怎么装才能使箱子打到最小 ...

常用jQuery代码02

一.each函数拿到每个元素的宽度 setTimeout(function () { $(".sticker_list img").each(function () { var W ...

LevelDB的源码阅读(二) Open操作

在Linux上leveldb的安装和使用中我们写了一个测试代码,内容如下: #include "leveldb/db.h" #include #in ...

mybatis源码解析之Configuration加载(二)

概述 上一篇我们讲了configuation.xml中几个标签的解析,例如,,等,今天我们来介绍 ...

腾讯云cos封装

public class CosUtil { int _appId = xxxxx; string _secretId = "xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx ...

PID控制器(比例-积分-微分控制器)- II

Table of Contents Practical Process Control Proven Methods and Best Practices for Automatic PID Cont ...