untitled3.m

利用矩阵分析的方法对线性方程组进行求解分析

%探索svd算法
w=10;
B=1;
row1=20;%采样区间采样点数
t1=linspace(-pi/10,pi/10,row1);
A1=matrixA(row1,5,t1);% 定义一个15*7大的A矩阵，用于解Ax=G
[U S V]=svd(A1)% 对A进行svd分解，得到U S V
%G1=matrixG(row1,t1)  %定义一个15*1大的G矩阵，
%S(1,1)=0;
%S(2,2)=0;
%S(3,3)=0;
%S(4,4)=0;
%S(5,5)=0;
%S(6,6)=0;
%S(7,7)=0;%matlab中有此函数（对奇异值进行取舍的函数），博库书城查《matlab数值计算范例教程》
%t=linspace(pi/(60*B),pi/2,15);
%t=linspace(-pi,pi,20);
%X=inv(V')*pinv(S)*inv(U)*G1 %svd中已去掉两个小的奇异值。分别为S(6,6)=0,S(7,7)=0
for i=1:5
    if S(i,i)>0
        S(i,i)=1/S(i,i);
    end
end
 S