一、matlab在高等数学中的应用(《数学建模算法与应用》P453)
1、求极限
syms x b=limit((sqrt(1+x^2)-1)/(1-cos(x)))
syms x a b=limit((1+a/x)^2,x,inf)
2、求导数
syms x dy=diff(log((x+2)/(1-x)),3) dy=simplify(dy) %对符号函数进行化简 pretty(dy) %分数线居中显示,这部一般不需要
差分:
a=[0,0.5,2,4]; da=diff(a)
3、求极值
syms x y=x^3+6*x^2+8*x-1; dy=diff(y); dy_zero=solve(dy) %表达式类型 dy_zero_num=double(dy_zero) %变成数值类型 ezplot(y)
syms x y=x^3+6*x^2+8*x-1; dy=diff(y); dy_zero=solve(dy) %表达式类型 dy_zero_num=double(dy_zero) %变成数值类型 ezplot(y) %符号曲线画图
4、求积分
不定积分:
syms x I=int(1/(1+sqrt(1-x^2))) pretty(I)
定积分 (符号解、数值解):
符号解:
syms x I=int(cos(x)*cos(2*x),-pi/2,pi/2)
数值解:
I=quadl(@(t)(t-3*t.^2+2*t.^3).^(-1/3),eps,0.5)
(积分上下限不能为inf,可以不变换而直接设一个很大的数
syms t %不加也可以 I=quadl(‘(t-3*t.^2+2*t.^3).^(-1/3)‘,eps,0.5) %后面两项分别为积分上下限,eps认为是0
二重数值积分:
I=dblquad(@(x,y)sqrt(1-x.^2-y.^2).*(x.^2+y.^2<=x),0,1,-0.5,0.5)
(matlab中的log就是ln
三重数值积分:
fun3=@(x,y,z)z.^2*log(x.^2+y.^2+z.^2+1)./(x.^2+y.^2+z.^2+1).*(z>=0&... z<=sqrt(1-x.^2-y.^2)); %...是续航符 I=triplequad(fun3,-1,1,-1,1,0,1)
注意把范围不等式点乘一下,最后的上下限只能为数字
5、级数求和
级数的上下限可以带inf
syms n f1=(2*n-1)/2^n; s1=symsum(f1,n,1,inf)
syms n f2=1/n^2 s2=symsum(f2,n,1,inf)