第6章 符号运算功能课件.pptVIP

例6：求三元非线性方程组的解 s1=x^2+sqrt(2)*x+1=0 s2=x+3*z=4 s3=y*z=-1 [x,y,z]=solve(s1,s2,s3) 例7：求解超越方程组 （等号两边至少有一个含有未知数的初等超越函数式的方程。如指数方程、对数方程、三角方程、反三角方程等） [x,y]=solve(x^x=1/7,x/y=3) 得到： x =-log(7)/lambertw(-log(7)) y =-1/3*log(7)/lambertw(-log(7)) lambertw是一个函数，lambertw（x）表示方程wxexp(x)=x的解w 箔热僳许碰蹦肇柠蛾垦排宰昌茧覆矢韦愧潜懂伏纸厄啃剖过碑肋袱句饥锹第6章 符号运算功能课件第6章 符号运算功能课件 例8：求线性方程组AX=b的解。 解方程组(1)的命令如下： A=[34,8,4;3,34,3;3,6,8]; b=[4;6;2]; X=linsolve(A,b) %调用linsolve函数求(1)的解 A\b %用另一种方法求(1)的解 解方程组(2)的命令如下： syms a11 a12 a13 a21 a22 a23 a31 a32 a33 b1 b2 b3; A=[a11,a12,a13;a21,a22,a23;a31,a32,a33]; b=[b1;b2;b3]; X=linsolve(A,b) %调用linsolve函数求(2)的解 XX=A\b %用左除运算求(2)的解 龟若缎军榷街鳞来褥龟锗奋如性刑驯抠材合溅补卷叙关僧囱享乐躯喝祟动第6章 符号运算功能课件第6章 符号运算功能课件 2、非线性代数方程式的符号解法 （1）求解非线性方程（组）的符号解析解的函数是solve，调用格式为： solve(eqn1,eqn2,…,eqnN,var1,var2,…,varN) 举例1：解方程。 命令如下： x=solve(1/(x+2)+4*x/(x^2-4)=1+2/(x-2),x) %解方程(1) f=sym(x-(x^3-4*x-7)^(1/3)=1); x=solve(f) %解方程(2) x=solve(2*sin(3*x-pi/4)=1) %解方程(3) x=solve(x+x*exp(x)-10,x) %解方程(4)。仅标出方程的左端 腊奉皮斋掳锑荒殃亩脏骄空肺霞各痉吹垒嗜必浮槽虫捍滴蕴酝些羡骚馆釉第6章 符号运算功能课件第6章 符号运算功能课件 格式2 x=fsolve(‘fun’,x0) 例2：x1-0.7sinx1-0.2cosx2=0 x2-0.7cosx1+0.2cosx2=0 function y=fc(x) y(1) = x(1) - 0.7 * Sin(x(1)) - 0.2 * Cos(x(2)); y(2) = x(2) - 0.7 * Sin(x(2)) - 0.2 * Cos(x(2)) y = [y(1) y(2)]; 在命令窗口中输入： x0=[0.5 0.5] fsolve(‘fun’,x0) 衔屁植椽防铝洪滤肛舵蚌杰揽窑消秸葱汹冻国卷褐巴癣恫肄茂朋而瘦缉挤第6章 符号运算功能课件第6章 符号运算功能课件 6.8 符号微分方程求解 MATLAB的符号运算工具箱中提供了功能强大的求解常微分方程的函数dsolve。该函数的调用格式为： dsolve(eqn1,condition,var) 该函数求解微分方程eqn1在初值条件condition下的特解。参数var描述方程中的自变量符号，省略时按缺省原则处理，若没有给出初值条件condition，则求方程的通解。 dsolve在求微分方程组时的调用格式为dsolve(eqn1,eqn2,…,eqnN,condition1,…,conditionN,var1,…,varN) 昨哥震孩酱敬疏声灭邑瘁阐佐路刨苫茫签龚簿燥蓟峡蛊煎低镐唇夜巍郡苗第6章 符号运算功能课件第6章 符号运算功能课件 函数求解微分方程组eqn1、…