matlab程序设计(第二版)课后作业答案.docVIP

实验结果： 1、 利用函数eye，ones，zeros，rand，randn函数分别产生 首先编写脚本函数如下: a1=eye(3) a2=eye(15,8) b1=ones(3) b2=ones(15,8) c1=zeros(3) c2=zeros(15,8) d1=2*rand(3,3)-1 d2=2*randn(3,3)+1 r1=2*rand(15,8)-1 r2=2*randn(15,8)+1 2、使用函数diag和reshape函数： 编写脚本： a1=[8 -7 2]; a2=diag(a1); a3=fliplr(a2); a4=[5 3]; a5=diag(a4,-1); a6=fliplr(a5); a=a3+a6 a1=reshape(a,1,9) b1=[4 5 7]; b2=diag(b1); b3=fliplr(b2); b4=[2 0 8]; b5=diag(b4); b=b3+b5 b6=reshape(b,1,9) 3、直接在工作空间里输入： a=5-round(100*rand(50,2))/10 4、脚本文件： a=[-pi:2*pi/360:pi]; y1=sin(a) y2=cos(a) 5.、编写脚本： A=rand(10) b=A0.6 [i,j,v]=find(b) c=sum(v) 6、编写脚本： A=randn(10) b=-0.5A0.5 [i,j,v]=find(b); c=sum(v) 7、编写实现表中功能的语句： if and(a1,b=0.5) 语句1 elseif and(a1,b0.5) 语句2 elseif and(a=1,b=0.5) 语句3 elseif and(a=1,b0.5) 语句4 else disp(error) end 8、编写脚本如下： A=2*rand(3) b=A1 c=A.*b; [i,j,v]=find(c); v 9、编写脚本函数： A=101*rand(100,3) b=A100; c=A.*b; [i,j,d]=find(c); d 10、为节约空间，用5*5的矩阵来模拟100*100的矩阵：编写脚本函数如下： A=[0 2 0 3 0;0 0 0 0 0;0 8 0 5 0;0 4 0 5 0;0 4 0 2 0] b=any(A); [i,j,v]=find(1-b); A(:,j)=[]; B=A; c=any(B); [i,j,v]=find(1-c); B(:,j)=[]; A=B 实验结果： 1、输入脚本函数： x=0:pi/50:2*pi; k=[1 26 51 76 101]; x(k)=[]; figure(1) subplot(2,2,1) plot(x,sin(x)),grid on subplot(2,2,2) plot(x,cos(x)),grid on subplot(2,2,3) plot(x,tan(x)),grid on subplot(2,2,4) plot(x,cot(x)),grid on 2、分别利用title函数给图形加上标题，用text函数给坐标系中的建立文本对象，可将字符串放在图形的指定位置，用xlabel、ylabel函数添加x、y轴的标记，用legend函数给图形加上插图说明。写成脚本函数如下： x=0:0.1:2; y1=exp(x); figure(1); plot(x,y1); title(\itx); xlabel(x); ylabel(\ity=exp(x)); text(0.6,-0.28,exp(x)\rightarrow); legend(\ity=exp(x)); y2=log10(x); figure(2); plot(x,y2); title(\itx); xlabel(x); ylabel(\ity=log10(x)); text(0.6,-0.28,log10(x)\rightarrow); legend(\ity=log10(x)); 将脚本函数保存为xiti1.m，运行后产生图形如下： 3、编写脚本函数如下： x=0:pi/40:4*pi; figure(1); subplot(2,1,1); plot(x,sec(x)),grid on; title(正割函数图形); ylabel(y),xlabel(x); subplot(2,1,2); plot(x,csc(x)),grid on; title(余割函数图形); ylabel(y),xlabel(x); 4、编写脚本函数如下： x=(1:.1:10); y=exp(x+5)+x.^3; figure(1); semilogx(x,y,k-); title