第5章MATLAB符号运算_PPT教案.pptVIP

MATLAB 7.X程序设计 */44 5.6 方程求解 1.代数方程 格式： solve(f,v) 解方程f，自变量为v solve(f1,…,fn) 解由f1,f2…,fn构成的 方程组的解 MATLAB 7.X程序设计 */44 f=sym(a*x*x+b*x+c=0) solve(f) ans = 1/2/a*(-b+(b^2-4*c*a)^(1/2)) 1/2/a*(-b-(b^2-4*c*a)^(1/2)) solve(f,a) ans =-(b*x+c)/x^2 例 代数方程 ax2+bx+c=0 MATLAB 7.X程序设计 */44 例 求方程组 的解 f1=sym(x+y+z=10); f2=sym(x+y-z=5); f3=sym(x+5*y-z=5); [x,y,z]=solve(f1,f2,f3) x=15/2 y=0 z=5/2 x+y+z=10 x+y-z =5 x+5y-z=5 MATLAB 7.X程序设计 */44 2.常微分方程 格式： dsolve(‘f1,f2,…’,’cond1,cond2,…’,’v’) f1,…fn 为常微分方程式 D表示求微分，后面数字表示几重微分 Dy y’ D2y y’’ cond1…condn 为初始条件 v为自变量 ，默认时对t求导 MATLAB 7.X程序设计 */44 例 常微分方程 dsolve(‘Dy=7’) %求y’=7的通解 ans =7*t+C1 dsolve(‘Dy=x^2’) %自变量为t ans =x^2*t+C1 dsolve(Dy=x^2,x) %自变量为x ans =1/3*x^3+C1 dsolve(D2y=1+Dy,y(0)=1,Dy(0)=0) ans =exp(t)-t MATLAB 7.X程序设计 */44 求解两点边值问题： y=dsolve(x*D2y-3*Dy=x^2,y(1)=0,y(5)=0,x) MATLAB 7.X程序设计 */44 5.7 符号表达式绘图 1.fplot函数 格式： fplot(fun,lims) fun是要绘制的函数（可以是定义函数的M文 件，也可是以X为变量的可计算字符串） lims=[XMIN XMAX YMIN YMAX] 限定了x,y轴上的绘图空间 MATLAB 7.X程序设计 */44 例 subplot(2,2,1) fplot(‘sin’,[0 10]) %sin是Matlab的内部函数 subplot(2,2,2) fplot(‘sin(1/x)’,[1 10]) %x为自变量 subplot(2,2,3) fplot(‘[sin(x),cos(x)]’,[0 2*pi]) %两条曲线 subplot(2,2,4) fplot(‘[sin(x),cos(x)]’,[0 2*pi 0 1]) %改变绘图空间 MATLAB 7.X程序设计 */44 MATLAB 7.X程序设计 */44 2.ezplot函数 格式： ezplot(f) (1) f 为包含一个变量x的字符串或数学函数f=f(x)，则默认x轴的绘图区域为 (2) f 为f=f(x,y)的包含x,y之间关系的隐函数,则默认绘图区域为 MATLAB 7.X程序设计 */44 格式： ezplot(f,[min,max]) 绘制f(x)=0在minxmax或f(x,y)=0在 minxmax和minymax之间的图形 ezplot(f,[xmin,xmax,ymin,ymax]) 绘制f(x,y)=0在 xminxxmax和yminyymax之间的图形 MATLAB 7.X程序设计 */44 例18 subplot(2,2,1) ezplot(‘cos(x)’) %x范围[-2*pi,2*pi] subplot(2,2,2) ezplot(cos(x), [0, pi]) %x范围[-pi,pi] subplot(2,2,3) ezplot(‘x^2+y^2 - 1’) %x,y范围[-2*pi,2*pi]