（1）请分别写出根据B1, B2恢复原矩阵A的程序 。
（2）用指令检验所得的两个恢复矩阵是否相等 。
目的l 数组运算和矩阵运算的不同 。
l 如何判断两个双精度数组是否相等 。
l norm指令的应用 。
解答A=1,2;
3,4;
B1=A.0.5B2=A0.5A1=B1.*B1;
A2=B2*B2;
norm(A1-A2,fro)% 求误差矩阵的F-范数 ， 当接近eps量级时 ， 就认为实际相等B1 =1.0000 1.41421.7321 2.0000B2 =0.553 。

22、7 + 0.4644i 0.8070 - 0.2124i1.2104 - 0.3186i 1.7641 + 0.1458ians =8.4961e-016 n 4.在时间区间 0,10中 ， 绘制曲线 。
要求分别采取“标量循环运算法”和“数组运算法”编写两段程序绘图 。
目的l 加强理解数组运算的机理和应用 。
l 初步使用subplot, plot, xlabel, ylabel等指令绘图 。
解答%标量循环运算法t=linspace(0,10,200);
N=length(t);
y1=zeros(size(t);
for k=1:Ny1(k)=1-exp(-0.5*t(k)*cos(2*t(k);
endsub 。

23、plot(1,2,1),plot(t,y1),xlabel(t),ylabel(y1),grid on%数组运算法y2=1-exp(-0.5*t).*cos(2*t);
subplot(1,2,2),plot(t,y2),xlabel(t),ylabel(y2),grid on n 8.先运行clear,format long,rng(default),A=rand(3,3) ， 然后根据A写出两个矩阵：一个对角阵B ， 其相应元素由A的对角元素构成；另一个矩阵C ， 其对角元素全为0 ， 而其余元素与对应的A阵元素相同 。
目的l 常用指令diag的使用场合 。
解答clear,format longrng(defa 。

24、ult)A=rand(3,3)B=diag(diag(A)C=A-B A =0.814723686393179 0.913375856139019 0.2784982188670480.905791937075619 0.632359246225410 0.5468815192049840.126986816293506 0.097540404999410 0.957506835434298B =0.814723686393179 0 00 0.632359246225410 00 0 0.957506835434298C =0 0.913375856139019 0.278498218867 。

25、0480.905791937075619 0 0.5468815192049840.126986816293506 0.097540404999410 0 n 9.先运行指令x=-3*pi:pi/15:3*pi;
y=x;
X,Y=meshgrid(x,y);
warning off;
Z=sin(X).*sin(Y)./X./Y;
产生矩阵Z 。
（1）请问矩阵Z中有多少个“非数”数据？（2）用指令surf(X,Y,Z);
shading interp观察所绘的图形 。
（3）请写出绘制相应的“无裂缝”图形的全部指令 。
目的l 初步感受三维曲面的绘制方法 。
l 非数NaN的产生 ， 非数的检测 ， 和对图形的影 。

26、响 。
l sum的应用 。
l eps如何克服“被零除”的尴尬 。
解答x=-3*pi:pi/15:3*pi;
y=x;
X,Y=meshgrid(x,y);
warning offZ=sin(X).*sin(Y)./X./Y;
NumOfNaN=sum(sum(isnan(Z)%计算“非数”数目subplot(1,2,1),surf(X,Y,Z),shading interp,title(有缝图)%产生无缝图XX=X+(X=0)*eps;
YY=Y+(Y=0)*eps;
ZZ=sin(XX).*sin(YY)./XX./YY;
subplot(1,2,2),surf(XX,YY,ZZ),shading interp 。

27、,title(无缝图) NumOfNaN =181n 10.下面有一段程序 ， 企图用来解决如下计算任务：有矩阵 ， 当依次取10, 9, 8, 7, 6, 5, 4, 3, 2, 1时 ， 计算矩阵“各列元素的和” ， 并把此求和结果存放为矩阵Sa的第k行 。
例如时 ， A阵为 ， 此时它各列元素 的和是一个行数组 ， 并把它保存为Sa的第3行 。
问题：该段程序的计算结果对吗？假如计算结果不正确 ， 请指出错误发生的根源 ， 并改正之 。
目的l 正确理解sum的工作机理 。
l reshape的应用 。
解答（1）企图用以下程序完成题目要求 。
for k=10:-1:1A=reshape(1:10*k,k,10);
Sa(k,:)=sum(A) 。

28、;
endSa Sa =55 55 55 55 55 55 55 55 55 553 7 11 15 19 23 27 31 35 396 15 24 33 42 51 60 69 78 8710 26 42 58 74 90 106 122 138 15415 40 65 90 115 140 165 190 215 24021 57 93 129 165 201 237 273 309 34528 77 126 175 224 273 322 371 420 46936 100 164 228 292 356 420 484 548 61245 126 207 288 369 450 531 。

29、 612 693 77455 155 255 355 455 555 655 755 855 955 （2）正确性分析除k=1外 ， 计算所得Sa所有行的结果都正确 。
但k=1时 ， Sa的第1行应该与相同 。
上述程序的错误是对sum理解不正确 。
sum对二维数组 ， 求和按列施行；而对一维数组 ， 不管行数组或列数组 ， 总是求那数组所有元素的和 。
正确的程序应该写成for k=10:-1:1A=reshape(1:10*k,k,10);
Sa(k,:)=sum(A);
if k=1Sa(k,:)=A;
endendSa Sa =1 2 3 4 5 6 7 8 9 103 7 11 15 19 23 27 31 35 396 。

30、 15 24 33 42 51 60 69 78 8710 26 42 58 74 90 106 122 138 15415 40 65 90 115 140 165 190 215 24021 57 93 129 165 201 237 273 309 34528 77 126 175 224 273 322 371 420 46936 100 164 228 292 356 420 484 548 61245 126 207 288 369 450 531 612 693 77455 155 255 355 455 555 655 755 855 955 第4章 数值运算习题 4 及解答n 。

31、 1.根据题给的模拟实际测量数据的一组和 试用数值差分diff或数值梯度gradient指令计算 ， 然后把和曲线绘制在同一张图上 ， 观察数值求导的后果 。

来源：(未知)

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

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