第二章符号计算(新)讲述.ppt

2.7 代数状态方程求符号传递函数： 4)、编写M代码 新建一个M文件 代码： syms G1 G2 G3 G4 H1 H2 H3; A=[0, 0, 0, 0, 0, 0, -G1;G2, 0, 0, 0, 0, -G2, 0; 0, G3, 0, 0, G3, 0, 0; 0, 0, G4, 0, 0, 0, 0; 0, 0, 0, H2, 0, 0, 0; 0, 0, 0, H1, 0, 0, 0; 0, 0, 0, H3, 0, 0, 0]; b=[ G1; 0; 0; 0; 0; 0; 0]; c=[ 0, 0, 0, 1, 0, 0, 0]; Y2Ua=c*((eye(size(A))-A)\b); %利用“左除”取代“求逆”，计算传递函数 disp([blanks(5),传递函数 Y2Ua 为]) pretty(Y2Ua) 传递函数 Y2Ua 为 G1 G2 G3 G4 ------------------------------------------------------------- G2 G3 G4 H1 - G3 G4 H2 + G1 G2 G3 G4 H3 + 1 1 2 3、4 5 6 7 2.8 符号计算结果的可视化： 2.8.1 直接可视化符号表达式 指令名 含义 可执行示例 ezcontour 画等位线 见表2.8-1 P96 ezcontourf 画填色等位线 ezmesh 画网格线 ezmeshc 画带等位线的网格线 ezplot 画二维曲线 ezplot3 画三维曲线 ezpolar 画极坐标曲线 ezsurf 画曲面图 ezsurfc 画带等位线的曲面图 2.8 符号计算结果的可视化： 单独立变量符号函数的可视化 ezplot(Fx, [xmin, xmax, ymin, ymax]) %绘制y=f(x)的二维图形 ezplot(Fxy, [xmin, xmax, ymin, ymax]) %绘制f(x,y)=0的二维图形 ezplot(xt, yt, [tmin, tmax]) %绘制x=x(t)，y=y(t)的二维图形 ezplot3(xt, yt, zt, [tmin, tmax]) %绘制[x(t), y(t), z(t)]的三维图形 例2.8-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) 1 2 3 4 5 6 7 8 9 10 2.8 符号计算结果的可视化： 双独立变量符号函数的可视化 ezsurf(Fxyz, dom_f) %绘制F(x,y,z)=0的曲面图形， dom_f为指定的矩形区域 ezsurf(x, y, z, dom_f) %绘制[x(t), y(t), z(t)]的曲面图形 例2.8-2 使用球坐标参量画部分球壳。 clf x=‘cos(s)*cos(t)’; %注意：x是char；也可以将s、t定义为符号变量 y=cos(s)*sin(t); z=sin(s); ezsurf(x,y,z,[0,pi/2,0,3*pi/2]) %注意：ezsurf( )参数中函数部分支持char类型 view(17,40) shading interp colormap(spring) light(position,[0,0,-10],style,local) light(position,[-1,-0.5,2],style,local) material([0.5,0.5,0.5,10,0.3]) 2.8 符号计算结果的可视化： P101-104（自学） 注意以下语句的使用： 1)、figure(1)、figure(2)、……