matlab数值计算.ppt

1、MATLAB数值修正，MATLAB具有很好的数值修正能力，占据世界数值修正软件的主导地位，有数值运算的功能，矩阵运算多项式运算线性方程的数值统一修正线性插值法函数最佳差分方程的数值解，一、如何编写命令行基本操作沉积基质方法直接输入法的规则：沉积基质要素， 必须用逗号或空格分隔沉积基质元素的行和行之间必须用分号分隔。 符号角色，逗号和分号角色逗号和分号如果命令后面显示可用作指令之间分隔符号的分号，则结果不会显示在屏幕上。 注意：任何具有指定值的变量都保存在工作区中，无论是否显示在屏幕上，以后可以随时查看或调用。 变量名称尽量不要重复。 不重复的话会被复盖。 如果命令或矩阵太长，可以使用连续的行。

2、 id角色用于生成相等间隔的向量，而差动奥尔特的间隔为1。 选择矩阵以指定行、列和元素。 循环文，a=1，2，3； 四、五、六七、八、九或a=十二三； 四五六七八九、a=一、二、三； 四、五、六十七、八、九或a=十二三。 按四十五六七八九的顺序喀呖声。 对于相对较大的矩阵，可以使用ENTER代替分号，对于同一行的内容，可以使用继续符号(.)将一行的内容分成两行输入。 注意：此时换行字符为shift enter、矩阵的维数、a=1、2、3 4、5、6、7、8、9 size (a ) ans=3、矩阵的矩阵最大值、a=1、2、3； 可以将矩阵元素输入为4、5、6、7、8，9长度(a ) ans=3

3、，并且可以将矩阵元素输入为任意matlab方程式、实数、复数、复数输入为特殊函数I，并且可以将j输入为a=1 2 3。 4 5 6 x=2 pi/2。 sqrt(3) 3 5i、矩阵元素、x=2 pi/2； sqrt(3) 3 5i，x=2. 0001.57081.7321.00005.0000 I，MATLAB行政许可以矩阵为元素创建函数，其中C=a； 10 11 12，c=12346791012，MATLAB行政许可对矩阵的单个元素指派值，并且如果给定矩阵大于原始矩阵范围，则自动进行扩展，其中，未指派值的元素是0 .a (2，3 )=100，a=。 a=1234510079010，矩阵元素

4、可以由以下语句生成： S1:S2:S3、y=0:pi/4:pi、x=1:5、y=0. 78541.5770456789 a (4，5 )=10、a (133605 ) 使用a(:1:233605)ans=13046079、a(:)ans=12304060790、a (: ) ans=1470250300000010，2.MATLAB函数来表示矩阵rand随机矩阵eye单位矩阵zeros的所有元素为0、矩阵ones的所有元素为1的矩阵，b=eye(3) b=1 0 0 0 1 0 0 0 1、c=eye (size (b ) ) c=100000001d=ones (3) d=11111、ran

5、d 注意：因为matlab严格区分大小写，所以a和a是两个不同的变量。 matlab函数名称必须是小写字母。 在matlab工作区中长期存储有用数据的方法是生成mat数据文件。 save将工作区中的所有变量保存到matlab.mat文件中。 二、数据的保存和检索、差动奥尔特文件名、save data工作区内的所有变量保存到data.mat文件。 save data a b将工作区中的a变量和b变量保存到data.mat文件中。 下次运行matlab时，可以使用load命令调用生成的mat文件。 load load data load data a b mat文件是标准的二进制文件，也可以存储为

6、ASCII查询密码格式。 可以通过喀呖声恢复所有存储的变量。 矩阵加法、减法(，)运算规则：加法、减法的2个矩阵必须加上同一行和列的2个矩阵对应元素。 可以参加运算的两个矩阵之一是标量。 标量和矩阵的所有元素分别进行加减运算。三、矩阵运算，a=1 2 3； 4 5 6 7 8 9 b=a1b=234567910，a=1 2 3。 456789 b=a-1 b=012345678，2 .矩阵乘法()运算规则： a矩阵的列数必须等于b矩阵的行数标量，可以乘以任意矩阵。 a=1 2 3； 4 5 6 7 8 0 b=1； 2； 3； c=a*b c=14 32 23，d=-1。 0； 2； f=pi

7、*d f=-3.1416 0 6.2832，矩阵除法运算不在线性代数，有矩阵逆运算，matlab有两种矩阵除法运算：的左除法和右除法运算。 通常，x=AB是A*x=B的解，其中x=B/A是x*A=B的解，且通常两个这些个不相等a，p都是行列，ap没有意义。 a=1，2，3； 在4、5、6、7、8、9 a2 ans=30 36 42 66 81 96 102 126 150，如果在一个方阵中存在复特征值或负实际特征值，则非整数的平方是复矩阵。 a0.5ans=0. 44980.7623 I0. 55260.2068 I0. 6555-0.3487 I1. 0185 I1. 0842 I1. 25

8、15.0228 I1. 4844-0.0385 I1. 5873-0.5940 I1. 9500 340. 2717 I， 求inv矩阵逆det行列式的值的eig矩阵的特征值diag对角阵矩阵转置sqrt矩阵的开方、4 .矩阵的其他运算、5 .矩阵的特殊操作、矩阵的变量维a=1:12； b=形状(a，3，4 ) c=零(3，4 )； c(:)=a(: )矩阵的变换方向rot90:旋转90度。 fliplr:的上升flipud:提取下行翻转矩阵诊断3360提取主对角线tril:提取主下三角tril(A，k) triu:提取主上三角。 排队的扩张、关系运算，A=8 1 6 3 5 7 4 9 2

9、B=3 4 5。 6 9 7 1 8 2 A=B A=B，关系运算子对比较的两个序列a和b进行关系运算。 返还值与a和b的位数是相同的排列。 在a和b的对应位置进行了关系运算的结果为真时，结果数组的对应位置为1，否则为0 .数组运算指的是元素对元素的算术运算，与通常意义上的符号所表示的线性代数矩阵运算不同的数组加减(，-) a b a- b，5 .矩阵的数组运算，对应元素相加减() a=1 2 3； 四五六七八九b=二四六； 1 3 5 7 9 10 a.*b ans=2 8 18 4 15 30 49 72 90，a=1 2 3。 四五六七八九b=二四六； 13五七九十a * b ans=2

10、5 37 46 55 85 109 85 133 172，a./b=b.a.b=b./a./b=a，所有a的元素都由于b的对应元素而除去a.b。 b=4 5 6； c1=a.b； c2=b./AC1=4. 00002.0000 c2=4. 00002.50002.0000、a、b对应元素间的商.3 .对于排列幂(.) 元素的元素的应排列幂以求a的a的元素为底，b的对应位置的元素进行幂运算z=a.2 z=1. 004.009.00 z=a.BZ=1. 0032.00729.00，自动修正下：A=1 1 1。 1 2 3 1 3 6 B=2 x=A.2 y=2.A Z=X.B，结果： x=1119

11、36y=2224864z=111681811296，matlab f(x)=anxn an-1xn-1 a0可以用行向量p=an an-1 a1 a0来表示。 poly产生特征多项式系数向量特征多项式必须是n 1维的特征多项式的最初的要素是1、4，多项式运算，例如：a=1 2 3。 456780 p=多边形(a ) p=1. 00-6.00-72.00-27.00 p是多项式p(x)=x3-6x2-72x-27的矩阵描述方法，我们可以利用p 1，456780 p=多边形(a ) p=1. 00-6.00 -。3.conv、convs多项式乘法，例如：a(x)=x2 2x 3； b(x)=4x2 5x 6。 c=(x2 2x 3)(4x2 5x 6) a=1 2 3。 b=4 5 6； c=conv(a，b )=conv (123，456 )。 c=4. 0013.0028.0027.0018.00 p=多边形2天(c，x ) p=4x 413 x 328 x 227 x 18，4.deconv多项式除法，a=1 2 3。 c=4. 0013.0028.0027.0018.00 d=去除(c，a ) d=4. 005.006.00，5，5 .多项式导数，matlab提供了对polyder函