【大学MATLAB课件】第三讲 matlab的科学计算.ppt

第四章 MATLAB的科学计算 本章主要内容如下： 4.1 概述 4.2 MATLAB的多项式计算 4.3 MATLAB的符号变量 4.3 符号运算在高等数学中的应用 4.5 符号运算在工程数学中的应用 4.1.1 概述 MATLAB具有用法简单、灵活、结构性强、延展性好等优点，逐渐成为科技计算、视图交互系统和程序中的首选语言工具。 功能强大的数值运算功能 强大的符号运算功能 强大的图形处理能力 高级但简单的程序环境 丰富的工具箱与模块集 数学实验？？？ 4.1.1 概述 4.1.1 概述 4.1.1 概述 4.2 多项式计算 4.2 多项式计算 4.2 多项式计算 4.2 多项式计算 4.2 多项式计算 4.2 多项式计算 4.2 多项式计算 4.2.5 多项式相除—deconv a=[1 2 3]; c = [4.00 14.00 28.00 27.00 18.00]; d=deconv(c,a) d =4.00 5.00 6.00 4.2 多项式计算 4.2 多项式计算 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号微积分 符号代数方程求解 符号微分方程 符号积分变换 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; findsym(EXPR) ans = X, Y, a, b, c, delta, theta, x, y ◆在MATLAB早期的版本里,符号的四则运算有相应的函数 ◆在6.0以后的版本里，MATLAB符号运算里的基本运算符、关系运算符、三角函数、双曲函数及其反函数、指数对数函数和矩阵代数几乎与数值运算中的一样 ◆略有不同之处请查相关帮助 ◆代数方程的求解由函数solve实现： S=solve(‘exp1’,’exp2’,--,’expn’,’V1’,’V2’--’Vn’) ● ‘exp1’,’exp2’,--:表示方程的字符串或字符串 表达式 ● ‘expn’中不含等号，默认表示‘expn＝0’方程 ●’V1’,’V2’,--:表示求解的符号变量 ●S是一个构架数组，显示结果，用S.Vn引用 求下列非线性方程组关于y，z的解。 S=solve(u*y^2+v*z+w=0,y+z+w=0,y,z) disp(‘S.y’),disp(S.y),disp(‘S.z’),disp(S.z) S = y: [2x1 sym] z: [2x1 sym] S.y [ -1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))-w] [ -1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))-w] S.z [ 1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))] [ 1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))] Matlab 解初值问题 dsolve 求微分方程 dsolve 的使用 dsolve 举例 dsolve 举例 dsolve 举例 g=diff(f,v)：求符号表达式 f 关于 v 的导数 g=diff(f)：求符号表达式 f 关于默认变量的导数 g=diff(f,v,n)：求 f 关于 v 的 n 阶导数 diff syms x; f=sin(x)+3*x^2; g=diff(f,x) g=diff(f,x,2) g =cos(x)+6*x g =-sin(x)+6 4.3 高数应用 4.4.2 计算导数 int(f,v,a,b): 计算定积分 int(f,a,b): 计算关于默认变量的定积分 int(f,v): 计算不定积分 int(f): 计算关于默认变量的不定积分 例：计算 syms x; f=(x^2+1)/(x^2-2*x+2)^2; I=int(f,x); I=int(f,x,0,1) I = 1/4*(2*x-6)/(x^2-2*x+2)+3/2*atan(x) I = -1/4+3/8*pi 4.4 高数应用 4.4.3 计算积分 K=int(exp(-x^2),x,0,inf) K1=int(exp(-x^2),x,0,2);double(K1) 例：计算