实验一Matlab的矩阵生成和运算

1、matlab课内实验报告姓名：袁清晨学号：1009300131学院：理 学 院班级：数学101课程名称：matlab语言程序设计实验题目：matlab的矩阵生成和运算指导教师：王冠 2012年03月09日一、实验目的1. 熟悉matlab的工作环境；2. 掌握matlab中创建矩阵和向量的方法；3. 掌握matlab中的逻辑运算和关系运算；4. 掌握matlab中的矩阵的代数运算。二、实验内容和结果1.计算；保存；清除工作空间中的变量；将保存的读取到工作空间；清空命令窗口。>> format short>> s=(2*sin(pi/5+2)-0.4)/(cos(6)+4

2、)+4s = 4.11732.将format设为rat;输入h=hilb(4);查找函数hilb的帮助文件>> format rat>> h=hilb(4)h = 1 1/2 1/3 1/4 1/2 1/3 1/4 1/5 1/3 1/4 1/5 1/6 1/4 1/5 1/6 1/7 >> help hilb hilb hilbert matrix. hilb(n) is the n by n matrix with elements 1/(i+j-1), which is a famous example of a badly conditioned m

3、atrix. see invhilb for the exact inverse. hilb(n,classname) produces a matrix of class classname. classname must be either 'single' or 'double' (the default). this is also a good example of efficient matlab programming style where conventional for or do loops are replaced by vectoriz

4、ed statements. this approach is faster, but uses more storage. see also invhilb. reference page in help browser doc hilb3.直接输入生成矩阵；取矩阵的第1、3行元素生成矩阵；将矩阵按列合成一个列向量。>> a=8 7 9;3 5 6; 9 2 4a = 8 7 9 3 5 6 9 2 4>> b=a(1 3,:)b = 8 7 9 9 2 4>> c=b(:)c = 8 9 7 2 9 44.以为首项，为步长，为末项生成向量；计算向量的正弦

5、值，赋值到变量；以向量为对角线生成对角矩阵。>> x=-pi:0.2*pi:pix = columns 1 through 8 -3.1416 -2.5133 -1.8850 -1.2566 -0.6283 0 0.6283 1.2566 columns 9 through 11 1.8850 2.5133 3.1416>> y=sin(x)y = columns 1 through 8 -0.0000 -0.5878 -0.9511 -0.9511 -0.5878 0 0.5878 0.9511 columns 9 through 11 0.9511 0.5878 0

6、.0000>> a=diag(y)a = columns 1 through 8 -0.0000 0 0 0 0 0 0 0 0 -0.5878 0 0 0 0 0 0 0 0 -0.9511 0 0 0 0 0 0 0 0 -0.9511 0 0 0 0 0 0 0 0 -0.5878 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.5878 0 0 0 0 0 0 0 0 0.9511 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 columns 9 through 11 0 0 0 0 0 0 0 0

7、0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.9511 0 0 0 0.5878 0 0 0 0.00005.生成一个5阶希尔伯特矩阵；取矩阵的第1、3、5行，第2到4列元素生成矩阵；左右翻转矩阵得矩阵；对矩阵逆时针旋转90度得矩阵；取矩阵的主对角线向量；以向量为第-3对角线生成对角矩阵（矩阵得阶数是多少？）。>> h=hilb(5)h = 1.0000 0.5000 0.3333 0.2500 0.2000 0.5000 0.3333 0.2500 0.2000 0.1667 0.3333 0.2500 0.2000 0.1667 0.1429 0.2500

8、 0.2000 0.1667 0.1429 0.1250 0.2000 0.1667 0.1429 0.1250 0.1111>> m=h(1 3 5,2 3 4)m = 0.5000 0.3333 0.2500 0.2500 0.2000 0.1667 0.1667 0.1429 0.1250>> f=fliplr? f=fliplr |error: incomplete or misformed expression or statement.>> f=fliplr(m)f = 0.2500 0.3333 0.5000 0.1667 0.2000 0.2

9、500 0.1250 0.1429 0.1667>> r=rot90(f)r = 0.5000 0.2500 0.1667 0.3333 0.2000 0.1429 0.2500 0.1667 0.1250>> v=diag(r)v = 0.5000 0.2000 0.1250>> s=diag(v,-3)s = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.5000 0 0 0 0 0 0 0.2000 0 0 0 0 0 0 0.1250 0 0 0>> size(s)ans = 6 66.生成随机逻辑矩阵和；对

10、矩阵和作与、或、非以及异或运算，并将结果分别赋给四个变量名；选择其中两个测试逻辑函数any和all。>> a=round(rand(6)a = 1 0 1 0 0 0 0 0 1 1 0 1 1 1 0 0 0 0 0 0 0 0 1 1 1 1 1 1 0 0 1 1 1 0 0 0>> b=round(rand(6)b = 1 1 0 0 0 0 1 0 0 1 1 0 0 1 0 0 1 0 1 1 1 1 1 1 1 1 0 1 1 1 0 0 1 1 1 0>> c=a&bc = 1 0 0 0 0 0 0 0 0 1 0 0 0 1 0

11、0 0 0 0 0 0 0 1 1 1 1 0 1 0 0 0 0 1 0 0 0>> d=a|bd = 1 1 1 0 0 0 1 0 1 1 1 1 1 1 0 0 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0>> e=ae = 0 1 0 1 1 1 1 1 0 0 1 0 0 0 1 1 1 1 1 1 1 1 0 0 0 0 0 0 1 1 0 0 0 1 1 1>> f=bf = 0 0 1 1 1 1 0 1 1 0 0 1 1 0 1 1 0 1 0 0 0 0 0 0 0 0 1 0 0 0 1 1 0 0 0 1>> d=any(d)d = 1 1 1 1 1 1>> c=all(c)c = 0 0 0 0 0 0>>7.生成三阶魔方矩阵和三阶pascal矩阵；计算和的加减乘除运算以及按位乘除运算。>> a=magic(3)a = 8 1 6 3 5 7 4 9 2>> b=pascal(3)b = 1 1 1 1 2 3 1 3