MATLAB基础和应用教程(人民邮电出版社-蔡旭辉)第二章c.ppt

第2章 矩阵、数组和符号运算;第2章 矩阵、数组和符号运算;第2章 矩阵、数组和符号运算;第2章 矩阵、数组和符号运算;b. 常微分方程的符号解 函数 dsolve 用来求常微分方程的符号解。 在符号方程中，用符号表达式中包含的字母“ D”来代替微分运算，符号 D2、D3、…DN 分别对应于第二、第三、…第 N 阶导数。因变量是位于 D后面的变量，缺省的自变量为 t。 dsolve 的调用格式为： S=dsolve(’eq1’,’eq2’, ... ,cond1,cond2,...,v ) 输入参数包含三部分内容：微分方程、初始条件和指定独立变量。 初始条件或边界条件写成‘y(a)=b’或‘Dy(a)=b’等。a、b 可以是变量使用符以外的其他字符。 当初始条件少于微分方程数时，在所得解中将出现任意常数符 C1，C2，….，解中任意常数符的数目等于所缺少的初始条件数。;例1： 求微分方程dy/dx=ay的通解和当y(0)=b时的特解。 dsolve(‘Dy=a*y’) dsolve(Dy=a*y,x) dsolve(‘Dy=a*y’,’y(0)=b’,’x’) ans = C1*exp(a*t) 默认自变量为t ans = C1*exp(a*x) ans = b*exp(a*x) 例2： 求微分方程 =-a2y当y(0)=1及 时的特解。 dsolve(D2y=-a^2*y,y(0)=1,Dy(pi/a)=0,x) ans = cos(a*x);求 ， ， f(0)=1, g(0)=2 的解。 S = dsolve(’Df = f + g’,’Dg = -f + g’,’f(0) = 1’,’g(0) = 2’) S = f: [1x1 sym] g: [1x1 sym] S.f ans = exp(t)*(2*sin(t)+ cos(t)) S.g ans = -exp(t)*(-2*cos(t)+sin(t)) 函数 dsolve 命令求解微分方程时，如果得不到其解，则给出警告信息。 ;6、符号函数的二维图 *ezplot(f) 绘制表达式f（x）的二维图形，轴坐标的近似范围为[-2π, 2π]. *ezplot(f,[xmin,xmax]) 绘制表达式f（x）的二维图形，轴坐标的近似范围为[xmin, xmax]. *ezpolar(f) 在极坐标下绘制函数表达式f(x)的二维图形。 syms x f=(x^2)^(cos(x)^2); ezplot(f);例1：绘制函数表达式x2-y4的二维图形 syms x y ezplot(x^2-y^4) 例2：绘制误差函数f(x)= 的二维图形. syms x ezplot(erf(x)) %erf(x)：误差函数 grid;例3：在极坐标下绘制函数表达式1+cos(t)的二维图形。 syms t figure (1) ezplot(1+cos(t)) figure ( 2) ezpolar(1+cos(t));7.符号函数三维绘图函数(Three dimensional plotting function) *ezplot3(x,y,z) 绘制由表达式x=x(t), y=y(t) 和z=z(t)定义的三维曲线，自变量t的变化范围为 [0, 2π]。 *ezplot3(x,y,z，[tmin, tmax]) 绘制由表达式x=x(t), y=y(t) 和z=z(t)定义的三维曲线，自变量t的变化范围为 [tmin, tmax]。 *ezplot3(…,’animate’) 如果在函数中增加参数’animate’，则绘制三维动态轨迹图。;例；根据表达式x=sin(t), y=cos(t) 和 z=t, 绘制三维曲线. syms t; ezplot3(sin(t),cos(t),t, [0,6*pi]) 绘动态轨迹图 syms t; ezplot3(sin(t),cos(t),t, [0,6*pi],animate);8.符号函数等高线绘图函数 *ezcontour(f) 绘制由表达式f(x,y)定义的等高线，自变量 x和 y 的变化范围为[-2π, 2π]。 *ezcontour(f,domain) 绘制由表达式f(x,y)定义的等高线，自变量 x和 y 的变化范围由domain确定，domai