matlab程序题复习.doc

PAGE  PAGE 9 2008matlab程序题复习 1.编写m文件要求利用matlab数值绘图功能画出衰减振荡曲线及其它的包络线。的取值范围是，数值采样间隔点为pi/50，在图中用红色是线段绘出，用蓝色虚线绘出，且两函数作在同一张图上。 1.t=0:pi/50:4*pi; y0=exp(-t/3); y=exp(-t/3).*sin(3*t); figure; plot(t,y,-r) hold on; plot(t,y0,:b); plot(t,-y0,:b); hold off; title(‘y=exp(-t/3)sin(3t)’); xlabel(‘x’) ylabel(‘y’) 2.利用matlab数值绘图功能，画出所表示的三维曲面。的取值范围是，要求去除当x,y均等于0时，函数z的不连续点，数值坐标采样间隔为0.5. 2.clear;x=-8:0.5:8; y=x; ［X Y］=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps; Z=sin(R)./R; surf(X,Y,Z); (mesh(X,Y,Z)) colormap(cool) xlabel(x),ylabel(y),zlabel(z) Title(‘z=sin(sqrt(x^2+y^2)/sqrt(x^2+y^2)); 3.利用matlab数值绘图功能，绘制t在区间［-2pi 2pi］时函数y=sin(t)/t的图形，采样间隔为pi/10，要求利用逻辑运算求出x=0时函数近似极限，并修补图形缺口。 3.t=-2*pi:pi/10:2*pi; y=sin(t)./t; tt=t+(t==0)*eps; yy=sin(tt)./tt; subplot(1,2,1),plot(t,y),axis([-7,7,-0.5,1.2]), xlabel(t),ylabel(y)； subplot(1,2,2),plot(tt,yy),axis([-7,7,-0.5,1.2]) 4.编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。 4.function c=collatz(n) c=n; while n1 if rem(n,2)==0 n=n/2; else n=3*n+1; end c=[c n]; end 5.有一周期为4p的正弦波上叠加了方差为0.1的正态分布的随机噪声的信号，用循环结构编制一个三点线性滑动平均的程序。（提示：①用0.1*randn(1,n)产生方差为0.1的正态分布的随机噪声；②三点线性滑动平均就是依次取每三个相邻数的平均值作为新的数据，如x1(2)=(x(1)+x(2)+x(3))/3，x1(3)=( x(2)+x(3)+x(4))/3……）。 5. t=0:pi/50:4*pi; n=length(t); y=sin(t)+0.1*randn(1,n); ya(1)=y(1); for i=2:n-1 ya(i)=sum(y(i-1:i+1))/3; end ya(n)=y(n); plot(t,y,c,t,ya,r,linewidth,2) 6.有一组测量数据如下表所示，数据具有y=x2的变化趋势，用最小二乘法求解y。并要求图示之。 x11.522.533.544.55y-1.42.735.98.412.216.618.826.2 6. x=[1 1.5 2 2.5 3 3.5 4 4.5 5] y=[-1.4 2.7 3 5.9 8.4 12.2 16.6 18.8 26.2] e=[ones(size(x)) x.^2] c=e\y x1=[1:0.1:5]; y1=[ones(size(x1)),x1.^2]*c; plot(x,y,ro) hold on; plot(x1,y1,k); xlabel(x); ylabel(y); title(拟合曲线); legend; 7.，当x和y的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口中绘制出三维线图、网线图、表面图和带渲染效果的表面图。