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

精品文档 一、实验结果： 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-72]; a2=diag(a1); a3=fliplr(a2); a4=[53]; a5=diag(a4,-1); a6=fliplr(a5); a=a3+a6 a1=reshape(a,1,9) b1=[457]; b2=diag(b1); b3=fliplr(b2); b4=[208]; 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、编写实现表中功能的语句： ifand(a1,b=0.5) 语句1 elseifand(a1,b0.5) 语句2 elseifand(a=1,b=0.5) 语句3 elseifand(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=[02030;00000;08050;0 4050;04020] 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=[1265176101]; x(k)=[]; 精品文档 精品文档 figure(1) subplot(2,2,1) plot(x,sin(x)),gridon subplot(2,2,2) plot(x,cos(x)),gridon subplot(2,2,3) plot(x,tan(x)),gridon subplot(2,2,4) plot(x,cot(x)),gridon 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)\rightarro w); 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)\rightar row); legend( \ity=log10(x) ); 将脚本函数保存为xiti1.m，运行后产生图形如下： 3、编写脚本函数如下： x=0:pi/40:4*pi; figure(1); subplot(2,1,1); plot(x,sec(x)),gridon; title(正割函数图形); ylabel(y),xlabel(x); subplot(2,1,2);  plot(x,csc(x)),gridon; title(余割函数图形); ylabel(y),xlabel(x); 4、编写脚本函数如下： x=(1:.1:10); y=exp(x+5)+x.^3; figure(1); semilogx(x,y,k-); title(y=e^(^x^+^5^