码迷,mamicode.com
首页 > 其他好文 > 详细

Matlab与微积分计算

时间:2015-06-30 21:46:44      阅读:170      评论:0      收藏:0      [点我收藏+]

标签:

一.极限问题的解析解

  1.1 单变量函数的极限

  格式1: L= limit( fun, x, x0)
  格式2: L= limit( fun, x, x0, ‘left’ 或 ‘right’)

>> syms x a b;
>> f=x*(1+a/x)^x*sin(b/x);
>> L=limit(f,x,inf) 
L = 
b*exp(a)

  可以结合函数图像来观察。

  1.2 多变量

  格式: L1=limit(limit(f,x,x0),y,y0)
  或 L1=limit(limit(f,y,y0), x,x0)

  如果x0 或y0不是确定的值,而是另一个变量的函数,如x->g(y),则上述的极限求取顺序不能交换。

相当于换元法,那肯定是先还原。

>> syms x y a; 
>> f=exp(-1/(y^2+x^2)) … *sin(x)^2/x^2*(1+1/y^2)^(x+a^2*y^2);
>> L=limit(limit(f,x,1/sqrt(y)),y,inf)
L =
exp(a^2)

 

二.函数导数的解析解  

  2.1 导数

  格式: y=diff(fun,x) %求导数
  y= diff(fun,x,n) %求n阶导数

  例:

>> syms x; f=sin(x)/(x^2+4*x+3);
>> f1=diff(f); pretty(f1)

pretty函数可以使多项式看起来更清晰。

cos(x) sin(x) (2 x + 4)
--------------- - -------------------
2 2 2
x + 4 x + 3 (x + 4 x + 3)

否则是 cos(x)/(x^2 + 4*x + 3) - (sin(x)*(2*x + 4))/(x^2 + 4*x + 3)^2

原函数及一阶导数图:
>> x1=0:.01:5; 
>> y=subs(f, x, x1);
>> y1=subs(f1, x, x1);
>> plot(x1,y,x1,y1,‘:’)

更高阶导数:
>> tic, diff(f,x,100); toc
elapsed_time =
4.6860

  2.2 多元函数偏导  

  和求多远极限的方法类似。
  格式: f=diff(diff(f,x,m),y,n)
  或 f=diff(diff(f,y,n),x,m)

  例:求其偏导数并用图表示。

>> syms x y; z=(x^2-2*x)*exp(-x^2-y^2-x*y);
>> zx=simple(diff(z,x))
zx =
-exp(-x^2-y^2-x*y)*(-2*x+2+2*x^3+x^2*y-4*x^2-2*x*y)

>> zy=diff(z,y)
zy =
(x^2-2*x)*(-2*y-x)*exp(-x^2-y^2-x*y)
直接绘制三维曲面
>> [x,y]=meshgrid(-3:.2:3,-2:.2:2);
>> z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y);
>> surf(x,y,z), axis([-3 3 -2 2 -0.7 1.5]) %直接surf(x,y,z)就行

axis([-3 3 -2 2 -0.7 1.5])产生三维坐标系。

  matlab中三维网格图(mesh)和三维曲面图(surf)区别?
  感觉没啥区别,颜色不同。
  http://www.cnblogs.com/hxsyl/archive/2012/10/10/2718380.html#3218180

三.积分问题解析解  

  格式: F=int(fun,x)
  定积分 格式: I=int(f,x,a,b)
  无穷积分 格式: I=int(f,x,a,inf)

  多重积分:顺序的改变使化简结果不同于原函数,但其误差为0,表明二者实际完全一致。这是由于积分顺序不同,得不出实际的最简形式。
用simple(f-f1),看结果是不是0

四.级数展开

  taylor(f,x,k,a)
  在x = a出展开的前k项。没a的话默认在0出。

例:对y=sinx进行Taylor幂级数展开,并观察不同阶次的近似效果。

>> x0=-2*pi:0.01:2*pi; y0=sin(x0); syms x; y=sin(x);
>> plot(x0,y0,‘r-.‘), axis([-2*pi,2*pi,-1.5,1.5]); hold on
>> for n=[8:2:16] p=taylor(y,x,n), y1=subs(p,x,x0); line(x0,y1) end % p=taylor(y,n) (2009b

 

  line([起点横坐标,终点横坐标],[起点纵坐标,终点纵坐标]),

  例line([1,2],[3,4])将画出(1,3)到(2,4)的一条直线,而不是(1,2)到(3,4)。

  当line(a,b)中,a b是相同大小的矩阵时,将会在对应的每一列做一条直线。

  当line(a,b,c)时,相应地会在三维图中画一条线。

  a b c 均为2XN矩阵。
  注意:实际上就是连接了几个点,说成折现比较好。

  其他的用到再说吧。

五.级数求和

 技术分享

技术分享

技术分享

 

>> x = [1 2 3]

x =

1 2 3

>> y = [4 5 6 ]

y =

4 5 6

>> x*y
错误使用 * 
内部矩阵维度必须一致。

>> x.*y

ans =

4 10 18

  别被.*理解错了,我理解的是对于两个行向量a和b,a*b会是向量乘法,就是的到了一个数,a.*b也会是向量乘法,其实这是因为对matlab不熟悉所致,想想对min/max函数的误解。其实*就是矩阵乘法,提示维度不同,.*就是对应位置相乘,且并没有自动连加操作。

技术分享

>> syms m,n
未定义函数或变量 ‘n‘。

>> syms m n
>> f = symsum(1/m,m,1,n)-log(n)

f =

eulergamma + psi(n + 1) - log(n)

>> limit(f,n,1,inf)
错误使用 mupadmex
Unexpected second parameter to mupadmex.

出错 sym/limit (line 50)
rSym = mupadmex(‘symobj::map‘, args{1}.s, ‘symobj::limit‘, args{2}.s, args{3}.s, dir);

>> limit(f,n,inf)

ans =

eulergamma

六.数值微分

  d(x+1)-d(x)/Δx或者d(x)-d(x-1)/Δx或者d(x+1)-d(x-1)/2*Δx。

Matlab与微积分计算

标签:

原文地址:http://www.cnblogs.com/hxsyl/p/4611009.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!