西工大matlab第五章精要.ppt

第五章 MATLAB图形处理 第五章 MATLAB图形处理 六、符号表达式绘图 例题：使用命令fplot作出函数y=sin（x）的图形，自变量区间分别为：x在[0,2*pi]之间；x在[pi,3*pi]；x在[-pi,pi],y在[-pi,pi]；x在[0 2*pi]（此时要求作图误差小于0.001）。 第五章 MATLAB图形处理 六、符号表达式绘图 clf subplot(2,2,1),fplot(sin(x),[0 2*pi]); subplot(2,2,2),fplot(sin(x),[pi 3*pi]); subplot(2,2,3),fplot(sin(x),pi*[-1 1 -1 1]); subplot(2,2,4),fplot(sin(x),[0 2*pi],1e-3); 第五章 MATLAB图形处理 六、符号表达式绘图 2 ezplot（fun，lims，fig） fun: 绘制函数fun的图形； lims：作图区间； fig：图形窗口句柄。 重要提示：ezplot可作多元函数的图形。默认作图区间为[-2*pi 2*pi]. 第五章 MATLAB图形处理 六、符号表达式绘图 例题：绘制函数u2-v2-1=0在u区间[-3，2]，v区间[-2,3]上的图形。 clf ezplot(u^2-v^2-1,[-3 2 -2 3]);? ezplot(u^2-v^2-1,[-2 3 -3 2]);? 第五章 MATLAB图形处理 六、符号表达式绘图 例：符号表达式绘图程序如下： subplot(2,2,1); ezplot(x^2+y^2-9);axis equal subplot(2,2,2); ezplot(x^3+y^3-5*x*y+1/5) subplot(2,2,3); ezplot(cos(tan(pi*x)),[ 0,1]) subplot(2,2,4); ezplot(8*cos(t),4*sqrt(2)*sin(t),[0,2*pi]) 对MATLAB初学者来说，易混淆plot和ezplot的不同适用场合： 前者只适于“离散数值数据”表达的曲线绘制； 后者则适于“函数解析式”表达的曲线绘制。 练习1：绘制 和 它的积分 在 间的图形 。 解： syms t tao y=2/3*exp(-t/2)*cos(sqrt(3)/2*t) s=subs(int(y,t,0,tao),tao,t) subplot(2,1,1) ezplot(y,[0,4*pi]),ylim([-0.2,0.7]) grid on subplot(2,1,2) ezplot(s,[0,4*pi]) grid on title(s = \int y(t)dt) 练习2：图示微分方程 的通解和奇解的关系。 （1） clear all y=dsolve((Dy)^2-x*Dy+y=0,x) y = x^2/4 C3*x - C3^2 （2） clf hold on hy1=ezplot(y(1),[-6,6,-4,8],1); set(hy1,Color,r,LineWidth,5) for k=-2:0.5:2 y2=subs(y(2),C3,k); ezplot(y2,[-6,6,-4,8],1) end hold off box on legend(奇解,通解,Location,Best) ylabel(y) title([\fontsize{14}微分方程, (y )^2 – xy + y = 0 ,的解]) 练习3：求解两点边值问题： 并画图表示。 y=dsolve(x*D2y-3*Dy=x^2,y(1)=0,y(5)=0,x) ezplot(y,[-1,6]) hold on plot([1,5],[0,0],.r,MarkerSize,20) text(1,1,y(1)=0) text(4,1,y(5)=0) title([x*D2y - 3*Dy = x^2,, y(1)=0,y(5)=0]) hold off 第五章 MATLAB图形处理 补充： 画一个circle。 t = linspace(0,2*pi); x = cos(t);y = sin(t); plot(x,y) ？发现什么问题？尝试以下命令观察变化。 axis equal axis square