11、944374842393221638v3 =0.0000000000000002660111416629094726767991785515 n /5求符号矩阵的行列式值和逆 ， 所得结果应采用“子表达式置换”简洁化 。
目的l 理解subexpr指令 。
解答A=sym(a11 a12 a13;
a21 a22 a23;
a31 a32 a33)DA=det(A)IA=inv(A);
IAs,d=subexpr(IA,d) A = a11, a12, a13 a21, a22, a23 a31, a32, a33DA =a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12 。

12、*a23*a31 + a13*a21*a32 - a13*a22*a31IAs = d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22) -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21) d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)d =1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12 。

13、*a23*a31 + a13*a21*a32 - a13*a22*a31) n /6求的符号解 ， 并进而用该符号解求 ， 的准确值 。
目的l symsum, subs的应用 。
l 从实例中 ， 感受指令所给出的关于符号解的含义 。
解答syms x kf=x(k);
Z1=symsum(f,k,0,inf)Z1 =piecewise(1 0 此式总成立 ， 说明“无约束” 。
n 情况二：-(x-1)0此为“约束” ， 满足题意 。
n /8（1）通过符号计算求的导数 。
（2）然后根据此结果 ， 求和 。
目的l diff, limit指令的应用 。
l 如何理解运行结果 。
解答syms ty=abs(sin(t)d=diff(y) %求dy/ 。

14、dtd0_=limit(d,t,0,left) %求dy/dt|t=0-dpi_2=limit(d,t,pi/2) %求dy/dt|t=pi/2 y =abs(sin(t)d =sign(sin(t)*cos(t)d0_ =-1dpi_2 =0 n /9求出的具有64位有效数字的积分值 。
目的l 符号积分的解析解和符号数值解 。
l 符号计算和数值计算的相互校验 。
解答（1）符号积分syms x clearsyms xy=exp(-abs(x)*abs(sin(x)si=vpa(int(y,-10*pi,1.7*pi),64) %vpa 指定精确位数y =abs(sin(x)/exp(abs(x)si 。

15、 =1.087849499412904913166671875948174520895458535212845987519414166 （2）数值计算复验xx=-10*pi:pi/100:1.7*pi;
sn=trapz(exp(-abs(xx).*abs(sin(xx)*pi/100 %trapz梯形法求积分sn =1.0877 n /10计算二重积分 。
目的l 变上限二重积分的符号计算法 。
解答syms x yf=x2+y2;
r=int(int(f,y,1,x2),x,1,2) r =1006/105 n /11在区间 ， 画出曲线 ， 并计算 。
目的l 在符号计算中 ， 经常遇到计算结果是特殊经典函数的情况 。

16、 。
l 如何应用subs获得超过16位有效数字的符号数值结果 。
l 初步尝试ezplot指令的简便 。
解答（1）符号计算syms t x;
f=sin(t)/t;
y=int(f,t,0,x)% 将得到一个特殊经典函数y5=subs(y,x,sym(4.5)ezplot(y,0,2*pi) y =sinint(x)y5 =1.6541404143792439835039224868515（2）数值计算复验tt=0:0.001:4.5;
tt(1)=eps;
yn=trapz(sin(tt)./tt)*0.001 yn =1.6541 n /12在的限制下 ， 求的一般积分表达式 ， 并计算的32位有效数字表达 。
目的 。

17、l 一般符号解与高精度符号数值解 。
解答syms xsyms n positivef=sin(x)n;
yn=int(f,x,0,pi/2) y3s=vpa(subs(yn,n,sym(1/3)y3d=vpa(subs(yn,n,1/3) yn =beta(1/2, n/2 + 1/2)/2y3s =1.2935547796148952674767575125656y3d =1.2935547796148951782413405453553 n 13.有序列 ， （在此 ， ） ， 求这两个序列的卷积 。
目的l 符号离散卷积直接法和变换法 。
解答（1）直接法syms a b k nx=ak;
h=bk;
w=syms 。

18、um(subs(h,k,n)*subs(x,k,k-n),n,0,k)%据定义y1=simple(w)w =piecewise(a = b, bk + bk*k, a b, (a*ak - b*bk)/(a - b) （2）变换法（复验）syms zX=ztrans(ak,k,z);
H=ztrans(bk,k,z);
y2=iztrans(H*X,z,k)%通过Z变换及反变换 y2 =piecewise(b 0.5);
id=sub2ind(size(A),ri,cj);
ri=ri;
cj=cj;
disp( )disp(大于0.5的元素的全下标)disp(行号 ,int2str(ri)disp(列号 。

19、 ,int2str(cj)disp( )disp(大于0.5的元素的单下标)disp(id) A =0.8147 0.9134 0.2785 0.9649 0.95720.9058 0.6324 0.5469 0.1576 0.48540.1270 0.0975 0.9575 0.9706 0.8003大于0.5的元素的全下标行号 1 2 1 2 2 3 1 3 1 3列号 1 1 2 2 3 3 4 4 5 5大于0.5的元素的单下标1 2 4 5 8 9 10 12 13 15 n 7.采用默认全局随机流 ， 写出产生长度为1000的“等概率双位（即取-1 ， +1）取值的随机码”程序指令 ， 并给出 。

20、 -1码的数目 。
目的l 两种基本随机发生器的使用 。
l 关系运算产生逻辑数组可用于数组的元素的标识和寻访 。
l 逻辑数组的应用 。
l 如何判断两个整数数组是否相等 。
解答（1）运用均匀随机数解题法解法1rng default%为以下结果重现而设；产生默认随机流 。
详见第4.3.2节A=rand(1,1000);
a=2*(A0.5)-1;
Na=sum(a=-1) Na =512 （2）运用正态随机数解题法解法2 randn(state,123)B=randn(1,1000);
b=2*(B0)-1;
Nb=sum(b=-1) Nb =462 （3）直接发生法解法3c=randsrc(1,1000,-1,1);
。

21、Nc=sum(c=-1) Nc =482 n 2.已知矩阵 ， 运行指令B1=A.(0.5), B2=A(0.5), 可以观察到不同运算方法所得结果不同 。

来源：(未知)

【学习资料】网址：/a/2021/0121/0021123973.html

标题：MATLAB|MATLAB教程2012a习题解答1-7章完整版-张志涌-北航( 二 )