MATLAB程序设计与应用期末课程作业.doc

1、MATLAB程序设计与应用 期末课程作业 学院：信息科学与技术学院专业和班级：09计算机班姓名：学号：老师：我做的是MATLAB 7.0上的题。由于你说用matlab notebook做的那节课我正好请假了，当我知道是已经做完了！ 第一张 MATLAB语言入门1：简述MATLAB语言的特点。 答：MATLAB语言是一种高级编程语言，它的主要特点是语言简单易学，代码简短高效，计算功能强大，绘图非常方便，扩充功能强大，帮助功能完整。3：MATLAB 7.0的帮助文件怎么使用？ 答：（1），使用函数在线帮助，即使用help （2），是要全部帮助，在MATLAB 7.0的自述文件中的全部帮助，包括MA

2、TLAB 7.0语言介绍，函数说明，工具箱说明等。 第二章 MATLAB 7.0 程序设计基础1：简述MATLAB中矩阵的建立和引用的方法. 答：矩阵建立的方法有，逐个输入法，冒号法，特殊方法，即调用函数linspace和logscape等。引用同样是用冒号的方法。3:建立一个 A(1,1)=中国，A(1,2)=上海，A(2,1)=北京，A(2,2)=1 2 3;的细胞数组。答： A(1,1)=中国;A(1,2)=上海;A(2,1)=北京; A(2,2)=1 2 3; AA = 中国 上海 北京 1 2 35：用MATLAB语言的inline语言实现函数，这里u,v是自变量。答： f=inli

3、ne(sin(u)+cos(v),u,v)f = Inline function: f(u,v) = sin(u)+cos(v) 第三章 矩阵线性代数算法实现1：试在MATLAB中生成20 30内均匀分布的4介随即矩阵。 m=20; n=30; x=m+(n-m)*rand(4); x=m+(n-m)*rand(4)x = 24.2176 26.5574 26.7874 26.5548 29.1574 20.3571 27.5774 21.7119 27.9221 28.4913 27.4313 27.0605 29.5949 29.3399 23.9223 20.31833：已知A=1 2;

4、3 4,B=2 3;4 5,试用MATLAB分别计算A和B的乘积，点乘，点积.答： A=1 2;3 4; B=2 3;4 5; X=A*BX = 10 13 22 29 C=A.*BC = 2 612 205：已知A=1 2;3 4,试用MATLAB分别计算A的秩，转置。答： A=1 2;3 4; rank(A)ans = 2 Aans = 1 3 2 47：求解下列方程组；答案如下：A=1 -2 3 -4;-2 3 -4 5;-3 4 5 -6;4 -5 -6 7;b=1 2 3 4;C=A b; n=rref(C)n = 1 0 0 0 12 0 1 0 0 10 0 0 1 0 -29

5、0 0 0 1 -24 第五章 符号计算1：求在（0,2）上的极值。 clearf=inline(x.3-2*x-5); x1=fminbnd(f,0,2)x1 =0.81653：求答案如下：q1=quad(sqrt(4*cos(2*t).2+sin(t).2)+1,0,3*pi)q1 = 23.56895：求答案如下：q=dblquad(y*sin(x)+x*cos(y),3*pi,4*pi,pi,2*pi)q = -29.60887：计算积分，答：代码和答案如下 syms x y z; r1=int(x/(1+z.2),z)r1 = x*atan(z) r2=int(x*log(1+x),

6、0.1)? Error using = mapcheckstat at 29Error: Second argument must be either of form x or x=a.b intError in = map at 15 checkstat(result,status,nargout);Error in = sym.int at 40 r = reshape(map(f(:),int,x),size(f); r2=int(x*log(1+x),0,1)r2 =1/4r3=int(int(x.2+y.2+1,y,x,x+1),x,0,1)r3 =5/2第七章 图形处理1：试用MA

7、TLAB中在同一坐标系中画出，在 x=0:2*pi;y1=cos(x);y2=cos(2*x);y3=cos(3*x);y4=cos(4*x);subplot(2,2,1);plot(x,y1) subplot(2,2,2);plot(x,y2) subplot(2,2,3);plot(x,y3) subplot(2,2,4);plot(x,y4) grid on? Error using = grid at 104Unknown command option. grid on3：在同一只脚坐标系中画出，代码和图形如下； clearx=0:2*pi;y=sin(tan(x);plot(x,y,

8、b-)hold ony2=sin(5*x);plot(x,y2,r-) grid on legend(曲线一,曲线二,);? legend(曲线一,曲线二,); |Error: Unbalanced or unexpected parenthesis or bracket. legend(曲线一,曲线二,1); title(函数曲线绘制,显示比例) xabel(横轴)? Undefined function or method xabel for input arguments of type char. xlabel(横轴) ylabel(纵轴) 5：试在，上的网格表面图像代码如下 x=-5

9、:0.1:5; y=-4:0.1:4; x,y=meshgrid(x,y); z=0.5*x.3+y.3; mesh(x,y,z) colormap(hsv) brighten(0.60)画出的图像算法及图像如下；x=-5:0.01:5;y=-5:0.01:5;x y=meshgrid(x,y);z=0.5*x3+y2; subplot(2,2,1);mesh(x,y,z) subplot(2,2,2);surf(x,y,z); subplot(2,2,3);meshc(x,y,z); subplot(2,2,4);surfc(x,y,z); x=-5:0.01:5;y=-5:0.01:5;x

10、 y=meshgrid(x,y);z=+(x2+y2)*cot(x);Warning: Divide by zero. In cot at 10 plot3(x,y,z,b); 画出函数的图像及写出它的程序 f=x2+2.5siny+z2*x2*z2+sin(tan(tan（x); x=-5:0.01:5;y=-5:0.01:5; z=-5:0.01:5; x y=meshgrid(x,y); plot3(x,y,z,b); 画出的图形代码和图形如下 x=-2*pi:0.001:12*pi;h=figureh = 3 y=2*cos(x)+sin(cot(x); plot(y,xg); 7：动画制作：制作正弦曲线在上从起点延伸的图像 clear s=0.2; x1=0; nframes=50; for k=1:nframes x1=x1+s; x=0:0.001:x1; y=sin(x); plot(x,y); axis(0 2*pi -1 1) grid off M(k)=getframe;endmovie(M,5)