方程求解和微积分.ppt

第4讲 方程求解和微积分 1.求解单个代数方程 MATLAB具有求解符号表达式的工具,如果表 达式不是一个方程式(不含等号)，则在求解之前函数 solve将表达式置成等于0。 一、方程求解 格式：solve( ‘f(x) ’ )，缺省变量为“x” 即：解方程 f(x)=0 例：1.解方程 2.解方程 3.解方程 solve(x^2+2*x+1) solve(x^2+3*x+2) solve(x^2+2*x+5) 例：4.解方程 eval(字符串) 求出字符串的值 方法一：solve(x^3+2*x^2-2*x+1) 方法二：p=[1 2 -2 1]; roots(p) 带有等号的符号方程也可以求解。 例： f=solve( cos(x)=sin(x) ) 注意：在求解周期函数方程时，有无穷多的解。 在这种情况下，solve对解的搜索范围限制在 接近于零的有限范围，并返回非唯一的解的子集。 如果不能求得符号解，就计算可变精度解。 x=solve( exp(x)=tan(x) ) x= 1.306326940423079 2.代数方程组求解 格式：[x,y, …]=solve(‘eq1’,‘eq2’,…,‘eqn’) 解方程组 未知变量 方程组 “=0”可省略 例：解方程组 [x,y]=solve(‘2*x+y-3’,’3*x-y+2’) 练习：解方程组 [x,y,z]=solve(‘3*x+2*y-z-10’,’-x+3*y+2*z-5’,’x-y-z+1’) 二、极限 Syms x 定义自变量为x limit(F,a) 当x→a时,计算符号表达式F=F(x)的极限值; limit(F) 当x→0时,计算符号表达式F=F(x)的极限值; limit(F,a,‘left’) 计算符号函数F的左极限x→a-; limit(F,a,‘ right ’) 计算符号函数F的右极限x→a+ 。 函数 limit limit((sqrt(1+x)-sqrt(1-x))/sin(3*x)) syms x; limit((x^2+5)/(x-3),2) limit((1/(1-x)-3/(1-x^3)),1) 例： 练习： limit(((2*x+1)/(2*x-1))^(3*x),inf) 三、导数 函数 diff 可用来求一元函数导数与二元函数的偏导数 syms x diff(f(x)) 对表达式f(x)中指定符号变量x的1阶导数； diff(f(x),n) 对表达式f(x)中指定的符号变量x的n阶导数。 diff(f(x,y),x) 对表达式f(x,y)中指定的符号变量x的1阶偏导数。 diff(f(x,y),x,n) 对表达式f(x,y)中指定的符号变量x的n阶偏导数。 格式： 例：求 的一阶、二阶导数 diff(x^3+4*x^2+1) diff(x^3+4*x^2+1,2) 例：求 对x的一阶偏导 diff(y^2*sin(x^2),x) 练习： 四、积分 int(f(x)) 对符号表达式f(x)中指定的符号变量x计算 不定积分，后面没有带任意常数C； int(f(x),a,b) 对表达式f(x)中指定的符号变量x计算从 a到b的定积分。 定义变量x syms x 例：求 int(tan(x)) 例：求 int(x^3) 例：求 int(x*log(1+x),0,1) log24 log2(4) lg100 log10(100) ln3 log(3) 对数的输入法则：