矫同浩上机报告.docVIP

（威海） MATLAB上机报告 学 号：110610207 姓 名：矫同浩 教 师：曲荣宁 日 期：2013年06月日 首先列举matlab中的常用符号功能 函数名 功能描述 函数名 功能描述 + 加 ... 续行标志 - 减 , 分行符(该行结果不显示) * 矩阵乘 ; 分行符(该行结果显示) .* 向量乘 % 注释标志 ^ 矩阵乘方 ‘ 转置 .^ 向量乘方 =赋值运算 \ 矩阵左除 == 关系运算之相等 / 矩阵右除 ~= 关系运算之不等 .\ 向量左除 关系运算之小于 ./ 向量右除 = 关系运算之小于等于 : 向量生成或子阵提取 关系运算之大于 () 下标运算或参数定义 = 关系运算之大于等于 [] 矩阵生成 逻辑运算之与 . 结构字段获取符 | 逻辑运算之或 . 点乘运算,常与其他运算符联合使用(如.\) ~ 逻辑运算之非 ”和“：”可以实现矩阵的变换。代码如下： A=[1 2 3; 4 5 6; 7 8 0]%建立矩阵 B1=A(1:2:end,:)%提取A矩阵全部奇数行、所有列 B2=A([3,2,1],[1,1,1])%提取A矩阵3,2,1行、反复三次由由首列构成子矩阵 B3=A(:,end:-1:1)%将A矩阵左右翻转，即最后一列排在最前前面 生成矩阵如下： 部分课后题 第三章 1.试求出如下极限。 ① ， ② 【求解】极限问题由下面的语句可以直接求出。 syms x; f=(3^x+9^x)^(1/x); limit(f,x,inf) ans = 9 syms x; f=(x+2)^(x+2)*(x+3)^(x+3)/(x+5)^(2*x+5); limit(f,x,inf) ans = exp(-5) 2.求导数并作图 syms x; f=sin(x)/x^2+4*x+3); f1=diff(f) f1 = cos(x)/(x^2 + 4*x + 3) - (sin(x)*(2*x + 4))/(x^2 + 4*x + 3)^2 x1=0:.01:5; y=subs(f,x,x1); y1=subs(f1,x,x1); plot(x1,y,x1,y1,:) 3.试求下面的双重极限①, ② 【求解】双重极限问题可以由下面语句直接求解。 syms x y fa=(x^2*y+x*y^3)/(x+y)^3; limit(limit(fa,x,-1),y,2) ans = -6 fb=x*y/(sqrt(x*y+1)-1); limit(limit(fb,x,0),y,0) ans = 2 fc=(1-cos(x^2+y^2))*exp(x^2+y^2)/(x^2+y^2); limit(limit(fc,x,0),y,0) ans = 0 第四章 写出多项式的伴随矩阵 P=[2 7 0 3 4 7 1]; A=compan(P) A = -3.5000 0 -1.5000 -2.0000 -3.5000 -0.5000 1.0000 0 0 0 0 0 0 1.0000 0 0 0 0 0 0 1.0000 0 0 0 0 0 0 1.0000 0 0 0 0 0 0 1.0000 0 2.试判定下面矩阵是否为正定矩阵，如是，则得出其 Cholesky分解矩阵。, B = 【求解】矩阵能进行 Cholesky分解，故为正定矩阵。 A=[9,2,1,2,2; 2,4,3,3,3; 1,3,7,3,4; 2,3,3,5,4; 2,3,4,4,5]; chol(A) ans = B矩阵进行 Cholesky分解失败，故为非正定矩阵。 B=[16,17,9,12,12; 17,12,12,2,18; 9,12,18,7,13; 12,2,7,18,12; 12,18,13,12,10]; chol(B) Error using == chol Matrix must be positive definite. 第五章 1. 已知下述各个Z变换表达式，试对它们分别进行Z反变换①， ②， ③，④ ⑤