计算机编程与数字信号处理实习.doc

计算机编程与数字信号处理实习 能够利用Matlab熟悉地画图，内容包括：X、Y坐标轴上的label，每幅图上的title，绘画多条曲线时的legend，对图形进行适当的标注等。 （1）在一副图上画出多幅小图； clc;clear;close; x=linspace(0,5,200); y1=sin(3*x); y2=sin(x.^3); y3=exp(x-2) y4=sin(2*x)+cos(3*x); for k=1:4; if k==1 subplot(2,2,k) plot(x,y1) hold on; title(x与y1) elseif k==2 subplot(2,2,k) plot(x,y2) hold on; title(x与y2) elseif k==3 subplot(2,2,k) plot(x,y3) hold on; title(x与y3) elseif k==4 subplot(2,2,k) plot(x,y4) hold on; title(x与y4) end end print -djpeg -r0 xinhao1.jpeg （2）画出一组二维图形；（以分段函数为例） clc;clear;close; x=-3:0.01:3; y1=zeros(size(x)); y2=zeros(size(x)); y3=zeros(size(x)); N=length(x); for k=1:N if x(k)-1x(k)=-3; y1(k)=(-x(k).^2-4*x(k)-3)/2; elseif x(k)=-1x(k)1 ; y2(k)=-x(k).^2+1; else x(k)=3x(k)=1 ; y3(k)=(-x(k).^2+4*x(k)-3)/2; end end y=y1+y2+y3; plot(x,y) xlabel(x);ylabel(y); print -djpeg -r0 xinhao2.jpeg （3）画出一组三维图形； clear x=0:0.2:3*pi; y=0:0.25:5*pi; [xx,yy]=meshgrid(x,y); z1=sin(xx).*sin(yy); x=-3:0.25:3; y=x; [xx,yy]=meshgrid(x,y); z2=xx-0.5*xx.^3+0.2*yy.^2+1; x=-8:0.5:8; y=x; [xx,yy]=meshgrid(x,y); r=sqrt(xx.^2+yy.^2)+eps; z3=sin(r)./r; subplot(2,2,1),mesh(z1); title(z1) subplot(2,2,2),mesh(z2); title(z2) subplot(2,2,3),waterfall(z2); title(z2) subplot(2,2,4),mesh(z3); title(z3) print -djpeg -r0 xinhao3.jpeg 说明：第一个图和第二个图中对x1赋值方式不同，结果相同。可以看出，无论是第一段中对X1的定义，还是第二段中对X1的定义。只要X1的维度为10，就符合surf的输入条件。且第一段中X1的范围就是常数9，X2中同样不变，所以画出的图像当然是一样的。 （4）画出复数的实部与虚部。 clc;clear;close; x=20; y=rand(x,1)+rand(x,1)*i; m=real(y); n=imag(y); subplot(2,2,1) plot(y,.);title(复数); hold on;grid on; subplot(2,2,2) plot(m,n,.);title(复数（实部为横，虚部为纵)); xlabel(实部);ylabel(虚部); hold on;grid on; subplot(2,2,3) plot(m,.);title(复数实部); hold on;grid on; subplot(2,2,4) plot(n,.);title(复数虚部); hold on;grid on; print -djpeg -r0 xinhao4.jpeg （5）完成对一个源程序进行详细注释。（仍以分段函数为例） clc;clear x=-3:0.01:3; %从-3到+3以0.01为步长给x赋值 y1=zeros(size(x)); %初始化y1的值 y2=zeros(size(x)); %初始化y1的值 y3=zeros(size(