matrixg.m

来自「利用矩阵分析的方法对线性方程组进行求解分析」· M 代码 · 共 14 行

function[G]=matrixG(m,t);   
%矩阵G，与W,B都有关，并且与才采点,采样范围有关

w=10;
B=1;%所给的仿真条件
%t=[pi/(60*B) pi/(50*B) pi/(45*B) pi/(40*B) pi/(30*B) pi/(28*B) pi/(25*B) pi/(20*B) pi/(15*B) pi/(13*B) pi/(10*B) pi/(8*B) pi/(5*B) pi/(4*B) pi/(3*B)];
%在pi/(60*B)和pi/(3*B)之间取15个值。
%t=linspace(pi/(60*B),pi/2,15);
%t=linspace(-pi/60,pi/60,20);
G=zeros(m,1);%先行分配一个m*1大的0矩阵，有利于提高运算速度。
for i=1:m,
    G(i,1)=(2/pi)*((sin((B*t(i))/2)/t(i)).^2).*cos(w*t(i));  %利用for循环来构造这样一个m*1的矩阵
end