数值分析LU分解 要求matlab编程 求L 和U detA

网上很多高手编的程序，随便找找就是。

找了一个，你看看。

1.lu分解程序

function

[l,u,detA]

=

mylu(

a

)

%

该程式主是对系数矩阵a进行LU分解,具体算法可参见数值分析.

n=size(a,2)

u=zeros(size(a))

l=eye(size(a))

u(1,:)=a(1,:)

l(2:end,1)=a(2:end,1)/a(1,1)

for

r=2:n

for

j=r:n

u(r,j)=a(r,j)-l(r,1:r-1)*u(1:r-1,j)

end

for

i=r+1:n

l(i,r)=(a(i,r)-l(i,1:r-1)*u(1:r-1,r))/u(r,r)

end

end

detA=det(a)

2.调用

>>

A=magic(4)

>>

[L,U,detA]=mylu(A)

L

=

1.0000

0

0

0

0.3125

1.0000

0

0

0.5625

0.5663

1.0000

0

0.2500

1.3012

-3.0000

1.0000

U

=

16.0000

2.0000

3.0000

13.0000

0

10.3750

9.0625

3.9375

0

0

-0.8193

2.4578

0

0

0

0

detA

=

0

>>

在不同的城市代码里有不同的解释如：

1I是代表圆弧的直径 是X方向 平面为 G17

2 J 也是是代表圆弧的直径 是Y方向 平面为 G18

3 R 也是是代表圆弧的直径 是Z方向 平面为 G19 在钻孔中R代表的是离工件的距离

4 K 在钻孔中代表的是次数

5 P在G04P1000里代表的是暂停时间为10秒 M98P1 代表的是呼叫O1号程序

注：以上是发那科系统

---