《Experiment with MATLAB》读书笔记（四）

读书笔记（四）

这是第四部分数组与矩阵 将代码复制到m文件即可运行 函数部分需新建m文件保存 %% 向量与矩阵 x = [2; 4] % 向量 A = [4 -3; -2 1] % 矩阵 A*x A'*A % 转置 A*A' %% 随机矩阵 R = 2*rand(2,2)-1 %% 连线画图 X = [ -6 -6 -7 0 7 6 6 -3 -3 0 0 -7 2 1 8 1 2 -7 -7 -2 -2 -7 ]; X(:,end+1) = X(:,1); % 所有元素后移一位 plot(X(1,:),X(2,:),'.-','markersize',18,'linewidth',2) axis(10*[-1 1 -1 1]) axis square % 从头连线到尾

%% 旋转矩阵 theta = pi/6 % 方便的弧度制 G = [cos(theta) -sin(theta); sin(theta) cos(theta)] theta = 30 % 多一个d的角度 G = [cosd(theta) -sind(theta); sind(theta) cosd(theta)] subplot(1,2,1) dot2dot(G*X) subplot(1,2,2) dot2dot(G'*X) % 以原点为中心旋转

%% 更多表达形式 % 向量用方括号，空格隔开 v = [0 1/4 1/2 3/4 1] % 矩阵不同行之间用分号隔开，或者直接换行 A = [8 1 6; 3 5 7; 4 9 2] A = [8 1 6 3 5 7 4 9 2] %% 直接产生特殊矩阵 Z = zeros(3,4) % 三行四列都是0 E = ones(4,3) % 四行三列都是1 I = eye(4,4) % 单位矩阵，因为i已经被用又eye和I同音 M = magic(3) % 幻方矩阵（阶数） R = rand(2,4) % 两行四列随机矩阵 [K,J] = ndgrid(1:4) % 按行与列递增的矩阵 %% 冒号与分号 % 冒号规则产生向量 v = 0:0.25:1 n = 10 y = 1:n % 句尾加分号这行结果只计算不显示 n = 1000; y = 1:n; %% 矩阵乘法 KJ = K*J JK = J*K %% 数组乘法 % 数组运算是每个元素的运算 % 使用符号+ , - , .* , ./ , .\ 和 .^注意有点 K.*J v.^2 %% 转置 % 数组就单纯转置 % 复数矩阵变成共轭复数矩阵，再转置 % 加点就是直接转置 v = v' inner_prod = v'*v outer_prod = v*v' Z = [1 2; 3+4i 5]' Z = [1 2; 3+4i 5].'