matlab正比符号,第3章-MATLAB符号计算.ppt

3.6.1 代数方程 当方程不存在解析解又无其他自由参数时，MATLAB可以用solve命令给出方程的数值解。 语法： solve('eq', 'v') %求方程关于指定变量的解 solve('eq1', 'eq2', 'v1', 'v2',…) %求方程组关于指定变量的解 【例3.22】 求方程ax2+bx+c=0和sinx=0的解。 >> f1=sym('a*x^2+b*x+c')%无等号 f1 = a*x^2+b*x+c >> solve(f1) %求方程的解x ans = -(b + (b^2 - 4*a*c)^(1/2))/(2*a) -(b - (b^2 - 4*a*c)^(1/2))/(2*a) >> f2=sym('sin(x)') f2 = sin(x) >> solve(f2,'x') ans = 0 【例3.23】 求三元非线性方程组 的解。 >> eq1=sym('x^2+2*x+1'); >> eq2=sym('x+3*z=4'); >> eq3=sym('y*z=?1'); >> [x,y,z]=solve(eq1,eq2,eq3) %解方程组并赋值给x、y、z x = ?1 y = ?3/5 z = 5/3 程序分析：输出结果为“结构对象”，如果最后一句为“S=solve(eq1,eq2,eq3) ”，则结果如下。 S = x: [1x1 sym] y: [1x1 sym] z: [1x1 sym] MATLAB提供了dsolve命令，可以用于对符号常微分方程进行求解。 语法： dsolve('eq', 'con', 'v') %求解微分方程 dsolve('eq1,eq2…', 'con1,con2…', 'v1,v2…') %求解微分方程组 (1)当y是因变量时，微分方程'eq'的表述规定为： y的一阶导数 或 表示为Dy； y的n阶导数 或 表示为Dny。 (2)微分初始条件'con'应写成'y(a)=b，Dy(c)=d'的格式；当初始条件少于微分方程数时，在所得解中将出现任意常数符C1、C2……，解中任意常数符的数目等于所缺少的初始条件数。 【例3.24】 求微分方程 ，y(1)=0，y(0)=0的解。 >> y=dsolve('x*D2y?3*Dy=x^2','x') %求微分方程的通解 y = C3*x^4 - x^3/3 + C2 >> y=dsolve('x*D2y?3*Dy=x^2','y(1)=0,y(0)=0','x') %求微分方程的特解 y = x^4/3 - x^3/3 【例3.25】 求微分方程组 的解。 >> [x,y]=dsolve('Dx=y,Dy=?x') x = C8*cos(t) + C7*sin(t) y = C7*cos(t) - C8*sin(t) 程序分析：默认的自由变量是t，C1、C2为任意常数，程序也可指定自由变量，结果相同： >> [x,y]=dsolve('Dx=y,Dy=?x','t') 3.7.1 符号函数的绘图命令 1．ezplot和 ezplot 3命令 ezplot命令是绘制符号表达式的自变量和对应各函数值的二维曲线，ezplot3命令用于绘制三维曲线。 语法： ezplot(F,[xmin,xmax],fig) %画符号表达式的图形 【例3.26】 绘制线性系统 的相平面图，如图3.4所示。 相平面图以x和导数 为坐标轴，在每个时刻t，状态x和导数 按时间先后顺序绘制系统状态运动轨迹。 >> syms x; >> x1=dsolve('D2x+0.25*x=0','x(0)=1','Dx(0)=0')%初始值为1 x1 = cos(t/2) >> dx1=diff(x1) dx1 = -sin(t/2)/2 >> ezplot(x1,dx1,[-100,100]);hold on%绘制相平面图 >> x2=dsolve('D2x+0.25*x=0','x(0)=2','Dx(0)=0');%初始值为2 >> dx2=diff(x2); >> ezplot(x2,dx2,[-100,100]); >> x3=dsolve('D2x+0.25*x=0','x(0)=3','Dx(0)=0');%初始值为3 >> dx3=diff(x3); >> ezplot(x3,dx3,[-100,100]); 语法： ezplot3(x,y,z,[tmin,tmax],'animate') %绘制三维曲线 【例3.27】 用ezplot