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

《matlab程序设计(第二版)课后作业答案》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、Matlab程序设计实验报告一、实验结果：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)-1d2=2*randn(3,3)+1r1=2*rand(15,8)-1r2=2*randn(15,8)+12、使用函数diag和reshape函数：编写脚本：a1=[8-72];a2=diag(a1);a3=fliplr(a2);a4=[53];a5=diag(a4,-1);

2、a6=fliplr(a5);a=a3+a6a1=reshape(a,1,9)b1=[457];b2=diag(b1);b3=fliplr(b2);b4=[208];b5=diag(b4);b=b3+b5b6=reshape(b,1,9)3、直接在工作空间里输入：a=5-round(100*rand(50,2))/104、脚本文件：a=[-pi:2*pi/360:pi];y1=sin(a)y2=cos(a)5.、编写脚本：A=rand(10)b=A>0.6[i,j,v]=find(b)c=sum(v)6、编写脚本：A=randn(10)b=-0.5

3、ind(b);c=sum(v)7、编写实现表中功能的语句：ifand(a<1,b<=0.5)语句1elseifand(a<1,b>0.5)语句2elseifand(a>=1,b<=0.5)语句3elseifand(a>=1,b>0.5)语句4elsedisp('error')end8、编写脚本如下：A=2*rand(3)b=A>1c=A.*b;[i,j,v]=find(c);v9、编写脚本函数：A=101*rand(100,3)b=A>100;c=A.*b;[i,j,d]=find(c);d10、为节约空间，用5*5的矩阵来模拟100*100的矩阵：编写脚本函数如下：A=[0

4、2030;00000;08050;04050;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)=[];9Matlab程序设计实验报告figure(1)subplot(2,2,1)plot(x,sin(x)),gridonsubplot(2,2,2)plot(x,cos(x)),gridonsubplot(2,2,3)plot(x,tan(x)),gr

5、idonsubplot(2,2,4)plot(x,cot(x)),gridon2、分别利用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(

6、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)),gridon;title('正割函数图形');ylabel('y'),xlabel('x');su

7、bplot(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^)+x.^3');%此处加‘^’表示是写在指数部分的字符ylabel('y'),xlabel('x');5、绘制多峰函数[x,y]=meshgrid(-3:.125:3);z=peak