简单物理实验模拟.ppt

第一章 简单物理实验模拟 矩阵操作 连接: A=[1 2]; B=[3 4]; C=[A;B] C = 1 2 3 4 D=[A B] D = 1 2 3 4 转置: A=[1 2;3 4] A = 1 2 3 4 A ans = 1 3 2 4 求逆：inv（） 求逆：inv（） inv(A) ans = -2.0000 1.0000 1.5000 -0.5000 翻转 左右翻转fliplr(): A=[1 2 3;4 5 6] A = 1 2 3 4 5 6 fliplr(A) ans = 3 2 1 6 5 4 上下翻转flipud(): flipud(A) ans = 4 5 6 1 2 3 旋转90度rot90(): B = rot90(A,k) 逆时针旋转k*90度,k默认为1 rot90(A) ans = 3 6 2 5 1 4 rot90(A,2) ans = 6 5 4 3 2 1 魔方矩阵magic(): magic(3) ans = 8 1 6 3 5 7 4 9 2 抽取 对角元素抽取diag(X，k):抽取矩阵X的第k条对角线的元素向量。 A = 1 2 3 4 5 6 7 8 9 x=diag(A) x = 1 5 9 生成对角矩阵diag(v,k):使得向量v为所得矩阵的第k条对角元素。 v=[1 2 3]; X=diag(v) X = 1 0 0 0 2 0 0 0 3 X=diag(v,1) X = 0 1 0 0 0 0 2 0 0 0 0 3 0 0 0 0 X=diag(v,-1) X = 0 0 0 0 1 0 0 0 0 2 0 0 0 0 3 0 上三角矩阵和下三角矩阵的抽取 tril(X,k):提取矩阵X的第k条对角线下面的部分（包括第k条对角线） triu(X,k):提取矩阵X的第k条对角线上面的部分（包括第k条对角线） v = 1 2 3 4 5 6 7 8 9 a1=tril(v,-1) a1 = 0 0 0 4 0 0 7 8 0 a1=tril(v,1) a1 = 1 2 0 4 5 6 7 8 9 a1=triu(v,1) a1 = 0 2 3 0 0 6 0 0 0 a1=triu(v,-1) a1 = 1 2 3 4 5 6 0 8 9 a1=triu(v) a1 = 1 2 3 0 5 6 0 0 9 函数定义 生成*.m文件 函数包括:定义行, H1行,帮助文本,函数体,注释. 定义行: function [y1,y2,…ym]=ff(x1,x2,…,xn) 帮助文本: % 注释，介于函数声明名行和第一行命令之间的若干注释， H1: % 注释的第一行 函数体: (可包含子函数) function y = f(x) z=g(l); function z=g(l)