matla 数值计算 (1.docVIP

第三章课程复习 要求在闭区间上产生具有10个等距采样点的一维数组。试用两种不同的指令实现。 linspace(0,2*pi,10) 0:2*pi/9:2*pi 在使用123作为rand随机数发生器的初始化状态的情况下，写出产生长度为1000的“等概率双位（即取-1，+1）取值的随机码”程序指令，并给出 -1码的数目。 〖答案〗 A=randsrc(1000,1,[-1,1],123); B=sum(A==-1) 已知矩阵，运行指令B1=A.^(0.5), B2=A^(0.5), 可以观察到不同运算方法所得结果不同。（1）请分别写出根据B1, B2恢复原矩阵A的程序。（2）用指令检验所得的两个恢复矩阵是否相等。 A=[1 2;3 4] B1=A.^(0.5), B2=A^(0.5) B1.^2 B2^2 先运行clear,format long,rand(state,1),A=rand(3,3)，然后根据A写出两个矩阵：一个对角阵B，其相应元素由A的对角元素构成；另一个矩阵C，其对角元素全为0，而其余元素与对应的A阵元素相同。 clear,format long,rand(state,1),A=rand(3,3) B=diag(diag(A)) C=A-B 先运行指令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）请写出绘制相应的“无裂缝”图形的全部指令。 x=-3*pi:pi/15:3*pi; y=x; [X,Y]=meshgrid(x,y); warning off; Z=sin(X).*sin(Y)./X./Y; A=sum(sum(isnan(Z))) surf(X,Y,Z); shading interp xx=x+(x==0)*eps; yy=y+(y==0)*eps; [X,Y]=meshgrid(xx,yy); Z=sin(X).*sin(Y)./X./Y; surf(X,Y,Z); shading interp; 〖答案〗 下面有一段程序，企图用来解决如下计算任务：有矩阵，当依次取10, 9, 8, 7, 6, 5, 4, 3, 2, 1时，计算矩阵“各列元素的和”，并把此求和结果存放为矩阵Sa的第k行。例如时，A阵为，此时它各列元素 的和是一个行数组，并把它保存为Sa的第3行。问题：该段程序的计算结果对吗？假如计算结果不正确，请指出错误发生的根源，并改正之。 for k=10:-1:1 A=reshape(1:10*k,k,10); Sa(k,:)=sum(A); end Sa 数值计算 与符号计算相比，数值计算在科研和工程中的应用更为广泛。 和十年、二十年前相比较，在当今计算机软硬件的支持下，人们所能拥有的计算能力已经得到了巨大的提升。。 数值微积分 近似数值极限及导数 MATLAB数值计算中，没有专门的求极限和导数的指令。 原因：数值精度有限，不能表示无穷小量，不能准确描述一个数的邻域。 【例4.1-1】设，，试用机器零阈值eps替代理论0计算极限，。 x=eps; L1=(1-cos(2*x))/(x*sin(x)), L2=sin(x)/x, syms t f1=(1-cos(2*t))/(t*sin(t)); f2=sin(t)/t; Ls1=limit(f1,t,0) Ls2=limit(f2,t,0) clear k=logspace(0,14,15); x=k*eps; f1=(1-cos(2.*x))./x./sin(x)/2; f2=sin(x)./x; format short e disp([blanks(7),x,blanks(11),f1,blanks(11),f2]) disp([x,f1,f2]) clf n1=sum(f110^-10); subplot(2,1,1) semilogx(x(1:n1),f1(1:n1),ro,LineWidth,3) hold on semilogx(x(n1+1:end),f1(n1+1:end),b.-,LineWidth,3) hold off text(10^-14,0.2,错误区段),text(10^-6,0.8,正确区段) axis([10^-20,1,0,1.1]) grid on