untitled22.m

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

%用来画图，SVD值，与W,B都有关，并且与采样点，采样范围有关\
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=pinv(A1)*G1 %svd中已去掉两个小的奇异值。分别为S(6,6)=0,S(7,7)=0
row2=300;%比较范围所取点数
t2=linspace(-2*pi,2*pi,row2)
A2=matrixA(row2,5,t2)
p=A2*X
%subplot(1,2,1);
%t=linspace(-pi,pi,20);
plot(t2,p,'r')
hold on
%title('经过SVD后得到的图形')
t=linspace(-2*pi,2*pi,3000);
y=(2/pi)*((sin((B*t)/2)./t).^2).*cos(w*t)
plot(t,y)
%G2=matrixG(500,t2)
%plot(t2,G2)

%subplot(1,2,2);
%title('不经过SVD后的原来的图形')

%plot(t,p,'r',t,y)
xlabel('t');
ylabel('y值');
title('在(-pi/60,pi/60)之间采样60个点,在(-2*pi,2*pi)之间比较')
hold off

%legend('红：SVD值,绿：原方程值')