绝对值或复数模 abs
开平方 sqt
四舍五入到最接近的整数 round
除后余数 rem
以e为底的指数 exp
自然对数 log
以10为底的对数 log10
正弦 sin
余弦 cos
正切 tan
余切 cot
反正弦 asin
反正切 atan
反余切 acot
双曲正弦 sinh
双曲余弦 cosh
π pi
描点画图 plot(x,y,’r’)
绿 b
红 r
定义变量 syms sym
清除变量 clear(all)
显示变量 whos
因式分解 factor
多项式展开 expand
将式子化简 simplify
求方程的根 solve
求极限 limit(f(x),x,a,’right’)
无穷大 inf
求导 diff(f,x,n)
带入求值 subs(y,x1,x2)
泰勒展开 taylor(y,10,1)
求函数局部极小值 fminbnd(f,x1,x2)
求不定积分 int(f)
求定积分 int(f,x1,x2)
定积分的数值解法 quadl,quad(f,x1,x2)
这里的f定义上与数组类似。
给出固定位数的数值 vap(s,8)
绘制三维曲线 plot3(x,y,z)
绘制网格图 mesh(x,y,z)
绘制表格图 surf(x,y,z)
构造x-y平面上自变量的格点矩阵 [x,y]=meshgrid(x,y)
绘制出瀑布图形 waterfall(x,y,z)
绘制三维等高线图形 contour3(x,y,z)
多元函数求极限 limit(limit(f,x,x0),y,y0)
多元函数的偏导数 diff(diff(f,x,2),y,3)
矩形区域二重积分的数值解 dblquad(f,x0,x1,y0,y1)
离散函数定义 (与数组运算类似,需加点)
inline(‘f’,’x’,’y’)
一般区域二重积分的数值解
quad2dggen(f,x1,x2,y1,y2)
quad2dggen(f,y1,y2,x1,x2,eps)
三重积分 triplequad
第一类线积分
作变化ds=sqrt(diff(x,t)^2+ diff(y,t)^2+ diff(z,t)^2)dt
变为一重积分(这里t可以为x,y,z或其它)
第二类线积分
第一类面积分
第二类面积分
数组求和 sum
符号序列求和 symsum(f,v,a,b)
正无穷 inf
求微分方程的解析解 dsolve(‘方程1’,’方程2’…’方程n’,’初始条件’,’自变量’)
例:输入命令dsolve(’Du=1+u^2’,’t’)
高阶微分 Dn(n代表阶数)
[x,y,z]=dsolve('Dx=2*x-3*y+3*z','Dy=4*x-5*y+3*z','Dz=4*x-4*y+2*z', 't');
x=simple(x) % 将x化简
y=simple(y)
z=simple(z)
用Matlab求常微分方程的数值解
例:求解常微分方程,,
x=0:0.01:0.5;
fun=inline('-2*y+2*x*x+2*x');
[x,y]=ode15s(fun,[0:0.01:0.5],1)
隐函数的导数 maple('implicitdiff(f(u,x,y,z,…,)=0,u,x)')
多项式求根 roots(p)
也可以 f = @(x)x.^3-2*x-5;
z = fzero(f,2)
e exp(1)
erf(x)表示函数:erf(x)=(2/根号下派)*(exp(-z方)对z积分,积分下限是0,上限是x),误差函数从形式上很像正态分布的分布函数Φ(x),是对一个形如正态分布的概率密度函数做变上限积分的结果
conj 共轭
一次多项式拟合: p1 = polyfit(x,y,1)
三次多项式拟合: p3 = polyfit(x,y,3)
带入求多项式的值 y=polyval(p,x)
这里的p为多项式的系数向量,x为给定的x向量
一维数组的创建 方法1 a=1:2:10
方法2 a=[1,2,3,4,5]
方法3 x=linspace(x1, x2, n)
方法3 x=linspace(x1, x2, n)
这里的x1为第一个元素,x2为第二个元素
n为元素的个数
ones 生成全1矩阵 例ones(3) ones(3,4)
zeros 生成全0矩阵
大小 size
Length 列数
例3-9 计算
解:先定义函数,文件名:f.m
function y=f(x)
y=1/sqrt(2*pi)*exp(-x.^2/2);
保存后,在命令窗口键入
format long
>> [y,n]=quad('f',-3,3)
则显示结果为:
y =
0.99729991863154matlab等高线填充颜
n =
57 (表示被积函数f的调用次数)
例3-10 计算
其中D为直线y = x和抛物线y = x2所围部分。
解:由数学方法可得:
(1)画出积分区域示意图
>>syms x y
>> f=(2-x-y)/2;
>> y1=x;
>> y2=x^2;
>> ezplot(y1);hold on
>> ezplot(y2);
>> axis([0,2,0,2])
(2)确定积分限
a=fzero('x-x^2',0)
a =
0
>> b=fzero('x-x^2',1)
b =
1
这是第2次积分的上下限。
(3)积分运算
f_dy=int(f,y,x^2,x) %先对y积分
f_dy =
x-5/4*x^2-1/2*x*(x-x^2)+1/4*x^4
>> I=int(f_dy,a,b) %再对x积分
I =
11/120
例3-11 试求下面的二次积分
解:先定义函数,函数名为:my2dfun.m
function z=my2dfun(x,y)
z=exp(-x.^2/2).*sin(x.^2+y);
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论