MATLAB科学计算7.ppt

（7）output —求解信息，可选 iterations: 迭代计算次数 algorithm: 所用算法 message: Exit message …… （8）jacobian—所得解处的Jacobian矩阵 [x,fval,exitflag,output]=fsolve(fun,x0,options) 例6　Kepler方程 [x,fval,exitflag,output]=fsolve(fun,x0) 已知 e=0.01，M=32*pi/180，求E。 其中 E-偏近点角，弧度 e-偏心率 M-平近点角，弧度 X(0)=(100，100)T要求精度满足　 例7 [x,fval,exitflag,output]=fsolve(fun,x0,options) 练习不同方式定义fun [x,fval,exitflag,output]=fsolve(fun,x0,options) 例8 　 X(0)=(1,1,1)T ，要求精度满足　 功能 求多项式方程(组)解析解。 对非多项式方程(组)，有解析解则给出全部解，否则调用数值求解器返回一个数值解。 可解线性系统 5、solve 符号运算 5、solve 解方程 f(x)=0 格式 syms x a b%定义所有符号变量包括自变量和参变量 eq=f(x)%待解方程f(x)=0左端表达式s=solve(eq,var) 优选格式 用户指定自变量 解向量,一个解时 1 维, n 个解时 n 维 eq ――待解方程 此时可以写为p(x)==q(x)， 如果q(x)是0， ==可以省略不写 可直接写在solve中 syms x a b s=solve(eq, var) 说明 f(x)=0 var缺省格式 s = solve(eq) 通过symvar函数确定自变量 排在字母表中前面的字母作为参变量，排在后面的字母表示自变量。 f(x)=0 解方程组 5、solve syms x1 x2 a b… eq1=f1 (…) ; eq2=f2 (…); … S=solve(eq1,eq2,..., var1, var2, … )%方程 之间逗号间隔 5、solve 格式 F(X)=0 varn――用户指定待求变量 S=solve(eq1,eq2,..., var1, var2, … ) eqn —— 左端表达式 S――方程组 解的结构数组 只有一组解时，S是个解向量 多个解时，S是个矩阵 输入　 S.var1　 得到自变量var1的解 S=solve(eq1,eq2,..., var1, var2, … ) eqn ――待解方程组 此时eqn还可以写为p(x)==q(x)， 如果q(x)是0， ==可以省略不写 也可直接写在solve中： syms x a b S=solve(eq1, eq2,…var1,var2…) 说明 其他输出格式 [x1,x2,…]=solve(eq1,eq2,..., var1, var2, … ) var1,var2,...输出的解的顺序 x1,x2,...对应于var1,var2,...的顺序，给输出的解命名 说明 还可添加特殊要求，对解进行选择 s = solve(eq,var, Name,Value) 例如, 将Name ,value设置为 ‘Real,true 说明 例8 　 S=solve(eq1,eq2,..., var1, var2, … ) 例9　求P=15 atm, T=460 K 时的Vm 已知 B= -0.265×10-3 m3·mol-1 C=0.03025×10-3 m6·mol-2 s=solve(eq, var) 作业 读书报告 解非线性方程（组）的功能函数 (已介绍的除外,如vpasolve)内容: 函数功能,原理, 使用方法,举例提交:纸版准备:ppt 1、roots 功能求多项式方程 f(x)=0 的全部根，包括实根和复根 f(x)=0 格式 roots(c) c：待解方程中f(x)的系数向量 降幂排列 缺少的幂次系数为零，要补齐 1、roots f(x)=0 例1 解方程 roots 例2 求A的特征值 poly roots 例3