第九讲 符号运算(续).pptVIP

函数的傅立叶级数 MATLAB 5.x版中，尚未提供求函数傅立叶级数的内部函数。下面我们自己设计一个简化的求任意函数的傅立叶级数的函数文件。 function mfourier mfourier f,n syms x a b c; mfourier int f,-pi,pi /2; %计算a0 for i 1:n a i int f*cos i*x ,-pi,pi ; b i int f*sin i*x ,-pi,pi ; mfourier mfourier+a i *cos i*x +b i *sin i*x ; end Return 调用该函数时，需给出被展开的符号函数f和展开项数n，不可缺省。 例6.26在[-π，π]区间展开函数为傅立叶级数。 命令如下： x sym x ;a sym a ; f x; mfourier f,5 %求f x x的傅立叶级数的前5项 f abs x ; mfourier f,5 %求f x |x|的傅立叶级数的前5项 syms a; f cos a*x ; mfourier f,6 %求f x cos ax 的傅立叶级数的前6项 f sin a*x ; mfourier f,4 %求f x sin ax 的傅立叶级数的前4项 线性方程组求解 非线性方程的根 非线性方程的根 例： Matlab 符号方程求解器 Matlab 符号方程求解器 其它运算 本文观看结束！！！ 需要特别注意的是： ① solver 可以取以上五个函数之一，不同的函数代表不同的内部算法：ode23 运用组合的 2/3 阶龙格—库塔—费尔贝算法，ode45 运用组合的 4/5 阶龙格—库塔—费尔贝算法。通常使用函数 ode45； ② f 是由待解方程写成的m文件的文件名； ③ ts [t0, tf]，t0、tf为自变量的初值和终值； ④ x0为函数的初值； ⑤ options 用于设定误差限（可以缺省，缺省时设定为相对误差 10?3，绝对误差 10?6），程序为 options odeset ‘reltol’, rt, ‘abstol’, at 其中rt和at分别为设定的相对误差和绝对误差； ⑥ 在解 n 个未知函数的方程组时，x0、x 均为 n 维向量，m 文件中待解方程组应以 x 的分量形式写成； ⑦ 使用 Matlab 软件求数值解时，高阶微分方程必须等价地变换成一阶微分方程组。 例 8.5.4 求解下列微分方程 解：令 y1 x，y2 y1?，则微分方程变为一阶微分方程组： 1 建立 m 文件 vdp1000.m 如下： function dy vdp1000 t,y dy zeros 2,1 ; dy 1 y 2 ; dy 2 1000* 1-y 1 ^2 *y 2 -y 1 ; 2 取 t0 0，tf 3000，输入命令： [T,Y] ode15s vdp1000,[0 3000],[2 0] ; plot T,Y :,1 ,- 运行程序，得到如图的结果。 例 8.5.5 求解下列微分方程组 1 建立 m 文件 rigid.m 如下： function dy rigid t,y dy zeros 3,1 ; dy 1 y 2 *y 3 ; dy 2 -y 1 *y 3 ; dy 3 -0.51*y 1 *y 2 ; 2 取 t0 0，tf 12，输入命令： [T,Y] ode45 rigid,[0 12],[0 1 1] ; plot T,Y :,1 ,-,T,Y :,2 ,*,T,Y :,3 ,+ 运行程序，得到如图的结果。图中，y1 的图形为实线，y2 的图形为“*”线，y3 的图形为“+”线。 反函数 finverse f,v ：求 f 关于指定变量 v 的反函数 finverse f ：求 f 关于默认变量的反函数 syms x t; f x^2+2*t; g1 finverse f,x g2 finverse f,t 例：计算函数 的反函数 * * 在前面讨论过求和函数sum，sum处理的级数是有穷级数。对于无穷级数求和，sum是无能无力的。求无穷级数的和需要使用符合表达式求和函数symsum。 1. 级数的符号求和 级数符号求和函数symsum，调用格式为： symsum a,n,n0,nn 级 数 例1 求级数之和。 1 常数项级数 n sym n ; s1 symsum 1/n^2,n,1,inf %求s1 s1 pi^2/6 2 s2 symsum -1 ^ n+1 /n,1,inf %求s2。未指定求和变量，缺省为n s2 log 2 3 函数项级数 s3 syms