matlab2013实验1参考答案.docVIP

MATLAB实验一 MATLAB数值计算 试验报告说明： 1 做试验前请先预习，并独立完成试验和试验报告。 2 报告解答方式：将MATLAB执行命令和最后运行结果从命令窗口拷贝到每题的题目下面，请将报告解答部分的底纹设置为灰色，以便于批阅。 3 在页眉上写清报告名称，学生姓名，学号，专业以及班级。 3 报告以Word文档书写。 文档命名方式: 学号+姓名+_（下划线）+2010试验几.doc 如：110400220张三_试验1.doc 4 试验报告doc文档以附件形式发送到maya_email@。凡文档命名不符合规范，或者发送方式不正确，不予登记。 一 目的和要求 1 熟练掌握MATLAB变量的使用 2 熟练掌握矩阵的创建 3 熟练掌握MATLAB 的矩阵和数组的运算 4 掌握稀疏矩阵的运算 5 熟悉MATLAB多项式的运用 6 使用元胞数组和结构数组 二 试验内容 1 创建矩阵（必做） 1.1使用直接输入，from:step:to，linspace，logspace等方式创建矩阵。 1.2 输入矩阵 1.2-1）分别使用全下标和单下标达方式取出元素“6” a=[1 2 3 4;2 4 6 8;3 6 9 12] a(2,3) a(8) 1.2-2）分别用不同的方式从矩阵a中取出子矩阵 a([1,3],[1 3 4]) %双下标 a([1 7 10;3 9 12]) %单下标 a(logical([1 0 1]),logical([1 0 1 1])) %利用逻辑向量 1.2-3）分别写出以下的执行结果： a ( :, 3 ) a (:,end) a (:, 2 : 3) a ([1 3],: ) a (:) a (2 : 3) a([2 3 4]) 1.2-4）找出矩阵a中所有大于3并且小于9的元素和下标： ind=find((a3)(a9)) %返回下标 a(ind) %返回元素 2 矩阵的运算（必做） 2.1使用format命令，将 pi分别用short和long格式显示其值。 format short pi ans = 3.1416 format long pi ans = 3.14159265358979 2.2-1) 设, 。利用MATLAB分别计算下面的表达式：: 1) 2) x=3; v=4; log(x^2+v^2) ans = 3.21887582486820 sqrt(x-3)/(x-2*v)^2 %此处，注意 / 和 \ 的使用区别，有同学混淆了 ans = 0 2.2-2) 定义矩阵a=[2.3 ,5.4, 9]，b=[6 ,3.14,2] 1) 计算a的正弦值 a=[2.3,5.4,9],b=[6,3.14,2] sin(a) 2）对a的每个元素加2 a+2 3）a的每个元素求平方 a.^2 3）a和b的对应元素相乘 a.*b 4）计算a和b的转置的乘积 a*b’ 2.3 计算函数的值，其中t 的范围从0到pi，步长0.1pi；Z为0.707。为的部分，计算的值。 t=0:0.1*pi:pi; z=0.707; f=10*exp(-2*t).*sin(4*t)/sqrt(1-z^2) %或f=10/sqrt(1-z^2)*exp(-2*t).*sin(4*t) l=(f=0); % 关系运算，返回l 为逻辑向量 f1=f(l) % 利用逻辑向量l 取出f中满足要求的元素 %也可以用此指令代替： % li=find(f0) %找出大于0元素的下标 % f1=f(li)?; 2.4 设A=[1 2 1; 2 4 2; 2 1 1]; B=[3 1 2; 1 2 2; 3 1 4]; C=A+i*B. 利用MATLAB观察下列表达式的结果：1) C’ 和 C.’ 2) A’ 和 A.’（注意比较’ 和.’的区别） 求解过程略。注意区别：’ 是共轭转置，.’是非共轭转置。 2.5利用MATLAB求解线性方程组。（思考：可以用几种方法实现） a=[1 3 1;2 1 3; 1 1 0] b=[15 20 4] %注意符号 ，也可以写成b=[15;20;4]， 注意符号 ；漏写则不是列向量 a\b %此处注意 / 和\ 的区别。 %解法二： a=[1 3 1;2 1 3; 1 1 0] b=[15 20 4] inv(a)*b %此处，注意顺序，有同学写b*inv(a) 思考：错误在哪里？ %解法三 a=[1 3 1;2 1 3; 1 1 0] b=[15 20 4] a^-1*b % 此题考用数值方法求解方程组。采用符号方法，解