matlab实验四【DOC精选】.doc

实验四　曲线绘图 【实验目的】 1．了解曲线的几种表示方法。 2．学习掌握MATLAB软件有关的命令。 【实验内容】 绘制下列三种曲线： 以直角坐标方程表示的正、余弦线。 t=0:0.1:10; y1=sin(t);y2=cos(t);plot(t,y1,r,t,y2,b--); x=[1.7*pi;1.6*pi]; y=[-0.3;0.8]; s=[sin(t);cos(t)]; text(x,y,s); title(正弦和余弦曲线); legend(正弦,余弦) xlabel(时间t),ylabel(正弦、余弦) grid axis square 以参数方程表示的平面曲线（单位圆）。 t=0:0.1:2*pi; x=cos(t); y=sin(t); plot(x,y) title(x=cos(t),y=sin(t)); axis square 以参数方程表示的空间曲线。 t=0:0.1:20; x=exp(-0.2*t).*cos(pi/2*t); y=pi/2*exp(-0.2*t).*sin(t); z=t; plot3(x,y,z); xlabel(x轴) ylabel(y轴) zlabel(z轴) title( x=exp(-0.2*t)*cos(pi/2*t),y=pi/2*exp(-0.2*t)*sin(t),z=t, t∈[0：20]) 以极坐标方程表示的心脏线。 theta=0:0.1:2*pi; polar(theta,1-cos(theta)) title(‘r=a(1+cosθ)，θ∈[0,2π]’) 或者 clear; t=0:2*pi/40:2*pi; r=1+cos(t); x=r.*cos(t); y=r.*sin(t); plot(x,y) axis square 【实验准备】 1．平面、空间曲线的表示形式 对于平面曲线，常见的有三种表示形式，即以直角坐标方程，以参数方程，和以极坐标表示等三种形式。而对于空间曲线，常见的是用参数方程表示。 2．曲线绘图的MATLAB命令 MATLAB中主要用plot,fplot,plot3三种命令绘制不同的曲线。 plot(x,y) 作出以数据(x(i),y(i))为节点的折线图，其中x,y为同维数的向量。 plot(x1,y1,x2,y2,…)　作出多组数据折线图 fplot(‘fun’,[a,b])　作出函数fun在区间[a,b]上的函数图。 plot3(x,y,z)　空间曲线图，其中x,y,z为同维数的向量。 可以用help plot, help fplot, help plot3查阅有关这些命令的详细信息 【实验方法与步骤】 练习1 作出函数的图形，并观测它们的周期性。先作函数在上的图形，用MATLAB作图的程序代码为： x=linspace(-4*pi,4*pi,300); %产生300维向量x y=sin(x); plot(x,y) %二维图形绘图命令 结果如图4.1，上述语句中%后面如“%产生300维向量x”是说明性语句，无需键入。 图4.1 的图形 此图也可用fplot命令，相应的MATLAB程序代码为： clear; close; %clear清理内存；close关闭已有窗口。 fplot(sin(x),[-4*pi,4*pi]) 结果如图4.2. 图1.2 的图形 如果在同一坐标系下作出两条曲线和在上的图形，相应的MATLAB程序代码为： x=-2*pi:2*pi/30:2*pi; %产生向量x y1=sin(x); y2=cos(x); plot(x,y1,x,y2,’:’) %’:’表示绘出的图形是点线 结果如图4.３其中实线是的图形，点线是的图形。 图4.３ 的图形 练习２ 设，要求以0.01秒为间隔，求出y的151个点，绘出y及其导数的图形。相应的MATLAB程序代码为： dt=0.01; %设定时间间隔 t=0:0.01:1.5; %设定自变量t w=4*sqrt(3); %固定频率 y=sqrt(3)/2*exp(-4*t).*sin(w*t+pi/3); %注意用数组运算式 subplot(2,1,1),plot(t,y),grid %绘制曲线并加上坐标网格 title(绘图示例),xlabel(时间t),ylabel(y(t)) %加标注 %求导数并绘制导数曲线，注意数组求导后长度减少1 for i=1:length(t)-1 t1(i)=t(i); end Dy=diff(y)/dt; subplot(2,1,2), plot(t1,Dy), grid, ylabel(Dy(t)) %加标注 曲