高等数值分析次实验作业副本.docVIP

高等数值分析第一次实验 （一） 矩阵生成 1.生成对角阵 （1）第一种方法产生等差特征值，最大值为10，最小值为1，均匀分布； （2）第二种方法产生随机正特征值，最大值为100，最小值为1，随机均匀分布； （3）第三中方法产生正特征值，最大为100，最小为1，较多的特征值靠近100； （4）第四种方法产生正特征值，最大为100，最小为1，较多的特征值靠近1； （5）第五种方法产生特征值，最大为10，最小为-10，不包括[-1,1]中的值，随机分布； %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %按不同方法产生5类不同的对角矩阵 D=cell(5,1); d=cell(5,1); %%第一种方法产生随机正特征值，在区间[1,10]中选取等差特征值 d1=linspace(10,1,100); d{1}=d1; D{1}=diag(d1); %%第二种方法产生随机正特征值，条件数较大：区间[1,100]中随机选取特征值 d2=sort(rand(100,1),descend); d2=1+(d2-min(d2))/(max(d2)-min(d2))*99; d{2}=d2; D{2}=diag(d2); %%第三种方法产生随机正特征值分布，较多的特征值靠近100； d3=sort(normrnd(100,10,200,1),descend); d3=d3(101:200); d3=1+(d3-min(d3))/(max(d3)-min(d3))*99; d{3}=d3; D{3}=diag(d3); %%第四种方法产生随机数，较多的特征值靠近1； d4=sort(normrnd(100,10,200,1),descend); d4=d4(1:100); d4=1+(d4-min(d4))/(max(d4)-min(d4))*99; d{4}=d4; D{4}=diag(d4); %%第五方法产生随机特征值，生成非正定阵，但A非奇异随机数的分布为[-10，-1],[1,10]各50个数的均匀分布 d51=sort(rand(50,1),descend); d51=-10+(d51-min(d51))/(max(d51)-min(d51))*9; d52=sort(rand(50,1),descend); d52=1+(d52-min(d52))/(max(d52)-min(d52))*9; d5=[d52;d51]; d{5}=d5; D{5}=diag(d5); 2. 生成正交阵 正交的生成为由生成的矩阵的分解得到，首先，在矩阵的选择方面，使用随机分布和典型的矩阵来产生，[0,1]中均匀分布的矩阵，魔方分布矩阵， 分布矩阵，泊松分布矩阵，Pascal矩阵，希尔伯特矩阵，F分布矩阵，卡方分布矩阵，hankel矩阵，使用 分解可以到的由这些典型矩阵分解得到的正交阵 ，其生成程序如下所示： %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%产生正交阵， M=cell(10,1); Q=cell(10,1); R=cell(10,1); %产生矩阵M M{1}=rand(100); M{2}=magic(100); M{3}=hilb(100); M{4}=pascal(100); M{5}=poissrnd(50,100,100); M{6}=normrnd(50,10,100,100); M{7}=frnd(10,10,100,100); M{8}=trnd(10,100,100); M{9}=chi2rnd(5,100,100); M{10}=hankel(1:100); for n=1:10; [Q{n},R{n}]=qr(M{n}); end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 3. 生成对称矩阵 因为 是对角阵，所以 是对称矩阵，并且的特征值即为 的对角元，所以对角阵可以按式来生成，生成程序如下： %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% A=cell(5,10); for i=1:5 for j=1:10 A{i,j}=Q{j}*D{i}*(Q{j}); End； end ； %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 通过以上三步可以生成矩阵 ，