数学软件报告.pptVIP

数学建模常用数学软件： 1、Matlab 2、Lindo/Lingo 例9：已知检测参数f随时间t的采样结果，用数值插值法计算t=2，7，12，17，22，17，32，37，42，47，52，57时f的值。 命令如下： T=0:5:65; X=2:5:57; F=[3.2015,2.2560,879.5,1835.9,2968.8,4136.2,5237.9,6152.7,... 6725.3,6848.3,6403.5,6824.7,7328.5,7857.6]; F1=interp1(T,F,X) %用线性方法插值 F1=interp1(T,F,X,nearest) %用最近方法插值 F1=interp1(T,F,X,spline) %用三次样条方法插值 F1=interp1(T,F,X,cubic) %用三次多项式方法插值 所谓符号计算是指在运算时,无须事先对变量赋值,而将所得到结果以标准的符号形式来表示。 MathWorks公司以Maple的内核作为符号计算引擎（Engine），依赖Maple已有的函数库，开发了实现符号计算的两个工具箱：基本符号工具箱和扩展符号工具箱。 常见符号计算：解方程、微积分、简化方程表达式等。 * * 六、Matlab符号计算功能 * * 求极限 函数limit用于求符号函数f的极限。系统可以根据用户要求，计算变量从不同方向趋近于指定值的极限值。该函数的格式及功能： 1、微积分 * * limit(f,x,a)：求符号函数f（x）的极限值。即计算当变量x趋近于常数a时，f（x）函数的极限值。 limit(f,a)：求符号函数f（x）的极限值。由于没有指定符号函数f（x）的自变量，则使用该格式时，符号函数f（x）的变量为函数findsym(f)确定的默认自变量，既变量x趋近于a。 limit(f)：求符号函数f（x）的极限值。符号函数f（x）的变量为函数findsym(f)确定的默认变量；没有指定变量的目标值时，系统默认变量趋近于0，即a=0的情况。 limit(f,x,a,right)：求符号函数f的极限值。right表示变量x从右边趋近于a。 limit(f,x,a,left)：求符号函数f的极限值。left表示变量x从左边趋近于a。 * * 【例1】求极限 syms x; %定义符号变量 f=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/sin(x)^3; %确定符号表达式 w=limit(f) %求函数的极限 运行结果： w = -1/2 * * 【例2】求导数： x = sym(x); %定义符号变量 diff(sin(x^2)) %求导运算 运行结果： ans = 2*cos(x^2)*x * * 【例3】求下述积分。 求积分： syms x int(1/(1+x^2)) ans = atan(x) * * 【例4】求级数的和:键入：1/12+1/22+1/32+1/42+ …… syms k symsum(1/k^2,1,Inf) %k值为1到无穷大 ans = 1/6*pi^2 其结果为：1/12+1/22+1/32+1/42+ ……=π2/6 * * 【例5】将表达式(x^9-1)分解为多个因式。 syms x factor(x^9-1) ans = (x-1)*(x^2+x+1)*(x^6+x^3+1) 2、简化方程表达式 * * 代数方程的求解 求解代数方程组 微分方程的求解 微分方程组的求解 复合方程 反方程 3、解方程 * * 【例6】解代数方程：a*x2-b*x-6=0 syms a b x solve(a*x^2-b*x-6) ans = [ 1/2/a*(b+(b^2+24*a)^(1/2))] [ 1/2/a*(b-(b^2+24*a)^(1/2))] 即该方程有两个根: x1=1/2/a*(b+(b^2+24*a)^(1/2))； x2=1/2/a*(b-(b^2+24*a)^(1/2)) * * 4、解方程 【例7】作符号计算： a,b,x,y均为符号运算量。在符号运算前，应先将a,b,x,y定义为符号运算量 * * a=sym(a); %定义‘a’为符号运算量，输出变量名为a b=sym(b); x=sym(x); y=sym(y); [x,y]=solve(a*x-b*y-1,a*x+b*y-5,x,y) %以a,b为符号常数，x,y为