2012_ch2 MATLAB.docVIP

一些简单的数学模型： 例一：研究传染病的经典模型 1.SI 模型 2，SIS模型 3.SIR模型 例二：求阿基米德（Archimedes）螺线在到间的曲线长度函数，并求出时的曲线长度。 据数学分析知，，为弧长元素，而曲线长度。 符号计算 符号对象和符号表达式 数学表达式和方程的基本组成是：数字、参数、变量；运算符号（加减乘除等）；数学函数（如三角函数、指数函数等）。 MATLAB作为面向对象的科学计算语言，也是依靠基本符号对象（包括数字、参数、变量）、运算符及一些预定义函数来构造和衍生符号表达式、符号方程的。 符号对象的创建和衍生 生成符号对象的基本规则 符号数字 完全准确 的数字是所谓的符号（类）数字。它的定义格式如下 〖说明〗 具体数字Num必须处于单引号内。 【例2.1-1】符号（类）数字与数值（类）数字之间的差异。本例演示：它们在创建方式、显示形式不同、类别判断、具体数值上的不同。 a=pi+sqrt(5) %创建数值类常数 sa=sym(pi+sqrt(5)) %创建符号类常数 sa1=sym(pi+sqrt(5)) Ca=class(a) %判断a的数据类别 Csa=class(sa) %判别sa的数据类别 vpa(sa-a) %在32位精度意义上计算两类数字之间的差 a = 5.3777 sa = pi+sqrt(5) sa1 = 6054707603575008*2^(-50) Ca = double Csa = sym ans = .138223758410852e-16 符号参数 表达式中的称谓参数，它们事先并不指定具体数值。它们是符号参数，是符号运算的基本组成元素。定义格式如下： 〖说明〗 参数名Para, Para1, Para2, ParaN分别代表符号的命名规则没有什么特殊。建议：符号参数名不用“x两侧的字母”开头。 Flag表示参数属性。它可具体取以下词条： positive 限定取正实数； real 限定为实数； unreal 不限定的复数。 符号变量 谁是待解的变量？ 可实现独立自变量的自动认定的指令如下 〖说明〗 EXPR可以是符号矩阵。确认对整个矩阵进行。 x是首选符号变量，其后的次序规则是：与x的ASCII码值之差的绝对值小的字母优先；差绝对值相同时，ASCII码值大的字母优先。例如，符号变量字母的优先次序为x, y, w, z, v。 【例2.1-2】用符号计算研究方程的解。本例演示：自动识别变量的能力；指定变量不同，解也不同；熟悉符号法解题基本指令；专用指令findsym。 （1）不指定变量情况 syms u v w z %定义符号参数和变量 Eq=u*z^2+v*z+w; %构成符号表达式 result_1=solve(Eq) %采用自动识别变量解方程u*z^2+v*z+w=0 result_1 = -u*z^2-v*z 问题：这是关于哪个变量的解？w，这可有以下指令运行后看出。 findsym(Eq,1) %自动识别出的符号变量 ans = w （2）把z指定为变量的情况 result_2=solve(Eq,z) %对于指定变量z解方程 result_2 = 1/2/u*(-v+(v^2-4*u*w)^(1/2)) 1/2/u*(-v-(v^2-4*u*w)^(1/2)) 【例2.1-3】对独立自由符号变量的自动辨认。本例演示：什么是独立自由符号变量；findsym的使用；符号表达式的生成方法。 （1）生成符号变量和表达式 syms a b x X Y %定义符号参数/变量 k=sym(3); %符号常数 z=sym(c*sqrt(delta)+y*sin(theta)); %直接定义符号表达式 EXPR=a*z*X+(b*x^2+k)*Y; %构成衍生符号表达式 （2）找出EXPR中的全部自由符号变量 findsym(EXPR) findsym(EXPR,1) ans = X, Y, a, b, c, delta, theta, x, y ans = x 问题：k和z为什么没列出？k 是常数；z是衍生符号对象，不独立。 （3）在EXPR中确定一个自由符号变量 findsym(EXPR,1) （4）在EXPR中确定2个和3个自由变量时的执行情况 findsym(EXPR,2),findsym(EXPR,3) 注意：第3个是theta，不是X! 后者A