matlab在微积分中的应用

matlab在微积分中的应用

语法知识
8.1差分与近似微分
差分函数diff(X)用来计算向量x在相邻二元素间差分。diff(X)当x为向量时，则它返回一个x的差分向量，他比x少一个元素。diff(X)=[x(2)-x(1),x(3)-x(2),…,x(n)-x(n-1)]，每一列是原来列向量的差分。diff(x,n)返回x的n次差分。diff(x,n,dim)沿着dim所指定的维数进行计算。
8.2微分运算
diff(s,‘v’,n)s为符号数学表达式，v表达式中符号变量，n为正整数。
diff(s)对s表达式中缺省变量求微商，缺省变量由函数findsym确定。
diff(s,n)对s中的符号变量v求微商。
diff(s,‘v’,n)用eval将微分表达式转换成数值表达式。
8.3不定积分与定积分运算
R=int(s)s为数学符号表达式，返回s的不定积分。
R=int(s,v)v为积分的变量，
R=int(s,a,b)a，b为积分的上下限。
R=int(s,v,a,b)
8.4数值积分
1.梯形积分法： trapz(x,y)x为输入向量，y为被积函数的输出向量。
2.递推的辛普森积分法：q=quad(fun,a,b,tol,trace,p1,p2,…)fun为被积函数和表达式，a，b为积分的上限和下限，tol为积分允许误差。trace取1显示递推过程数据表，为0不显示。p1，p2用来对函数直接传递参数。
fun的三种表达式
（1）用字符串数字表达式。
（2）用内联函数inline。
（3）用函数句柄用m文件。
3.数值二重积分:q=dblquad(fun,xmin,xmax,ymin,ymax,tol)fun为被积函数，tol为允许误差。
4.三重数值积分：q=triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol)。
实战演练
8-1

X=1:6;
y1=diff(X)
y1 =
1 1 1 1 1

y2=diff(X,2)
y2 =
0 0 0 0

8-2

x=magic(4);
y1=diff(x)
y1 =
-11 9 7 -5
4 -4 -4 4
-5 7 9 -11

y2=diff(x,1,2)
y2 =
-14 1 10
6 -1 -2
-2 -1 6
10 1 -14

8-3

X=1:6;
Y=X.^2;
deltaY=diff(Y)./diff(X)%近似微商；
deltaY =
3 5 7 9 11

8-4

t=1:10;
s=[2 6 12 20 30 42 56 72 90 110];
v=diff(s)./diff(t)
v =
4 6 8 10 12 14 16 18 20

a=diff(v)./diff(t(1:9))
a =
2 2 2 2 2 2 2 2

8-5

syms x
y=xsin(x)+15;
diff(y)
ans =
sin(x) + xcos(x)

x=pi/4;
eval(ans)
ans =
1.2625

8-6

syms t
y=t^6;
diff(y,6)
ans =
720

8-7

syms x t
y=exp(-t)sin(xt);
yt=diff(y,‘t’)
yt =
x*exp(-t)cos(tx) - exp(-t)sin(tx)

yx=diff(y,‘x’)
yx =
t*exp(-t)cos(tx)

8-8

syms a t
x=acos(t)^3;
y=asin(t)^3;
x1=diff(x)
x1 =
-3acos(t)^2*sin(t)

y1=diff(y)
y1 =
3acos(t)*sin(t)^2

yx=x1/y1
yx =
-cos(t)/sin(t)

8-9

syms x y
z=x^2+2*x*y-y2-2x;
zx=diff(z,‘x’)
zx =
2x + 2*y - 2

zy=diff(z,‘y’)
zy =
2x - 2y

yx=-zx/zy
yx =
-(2x + 2y - 2)/(2x - 2y)

8-10

syms x y a
z=x^(2/3)+y(2/3)-a^(2/3);
zx=diff(z,‘x’)
zx =
2/(3*x^(1/3))

zy=diff(z,‘y’)
zy =
2/(3*y^(1/3))

yx=-zx/zy
yx =
-y^(1/3)/x(1/3)

8-11

syms x
y=sin(x)sin(2x)sin(3x);
y10=diff(y,10)
y10 =
15378432cos(2x)cos(3x)sin(x) + 14854144cos(2x)sin(3x)cos(x) + 14854656cos(3x)sin(2x)cos(x) - 15378944sin(2*x)sin(3x)*sin(x)

factor(262144)
ans =
1 至 10 列
2 2 2 2 2 2 2 2 2 2
11 至 18 列
2 2 2 2 2 2 2 2

8-12

syms x
y=x^2/(1-x);
y8=diff(y,8)
y8 =
(80640x)/(x - 1)^8 - 40320/(x - 1)^7 - (40320x^2)/(x - 1)^9

8-13

syms x t
y=exp(t)*(sin(t)+cos(t));
y1=diff(t)
y1 =
1

y1=diff(y)
y1 =
exp(t)(cos(t) + sin(t)) + exp(t)(cos(t) - sin(t))

t=log(x);
t1=diff(t)
t1 =
1/x

yx=y1t1
yx =
(exp(t)(cos(t) + sin(t)) + exp(t)*(cos(t) - sin(t)))/x

8-14

syms x
y=1/(1+x^2);
R=int(y)
R =
atan(x)

8-15

syms x
y=cos(x)cos(2x)cos(3x);
R=int(y)
R =
x/4 + sin(2x)/8 + sin(4x)/16 + sin(6*x)/24

8-16

syms a b x
y=exp(ax)sin(bx);
R=int(y)
R =
-(exp(ax)(bcos(bx) - asin(b*x)))/(a^2 + b^2)

8-17

syms x
y=(1-sin(x2))^(1/2);
R=int(y,0,pi/2)
R =
22^(1/2) - 2

eval®
ans =
0.8284

8-18

syms x y
V=int(int(10-x.^2/4-y.2/9,x,-4,4),y,-4,4)
V =
13952/27

eval(V)
ans =
516.7407

8-19

syms a b t
s=-absin(t)^2;
A=int(s,2pi,0)
A =
pia*b

8-20

y=[0 1 4 9 16 25 36];
z=trapz(y)
z =
73

cumtrapz(y)
ans =
1 至 4 列
0 0.5000 3.0000 9.5000
5 至 7 列
22.0000 42.5000 73.0000

8-21

x=0:pi/100:pi;
y1=sin(x).^2;
delta=pi/100;
z=(1/pideltatrapz(y1))^(1/2)
z =
0.7071

8-22

syms x y
v=int(int(6-x.^2/4-y.2/9,y,-4,4),x,-4,4)
v =
7040/27

eval(v)
ans =
260.7407

[x,y]=meshgrid(-4:0.2:4);
z=6-(x/2).^2-(y/3).2;
surfc(x,y,z)

学习心得
1.学1习了积分和微分