matlab概率统计ppt课件
10MATLAB符号运算与概率统计.ppt
26.03.2019
MATLAB符号运算
dsolve('D2y = -a^2*y', 'y(0) = 1', 'Dy(pi/a) = 0') ans= cos(a*t) dsolve('Dx = y', 'Dy = -x') x= cos(t)*C1+sin(t)*C2 y= -sin(t)*C1+cos(t)*C2
MATLAB符号运算与概率统计
26.03.2019
1
MATLAB符号运算
1.设置符号变量(sym, syms) S = sym(A):将已知量声明成符号变量。 x = sym(‘x’):声明x为符号型变量。 x = sym('x','real') x = sym('x','unreal') S = sym(A,flag) flag 可以是 ‘r’, ‘d’, ‘e’, 或 'f'.
12
26.03.2019
MATLAB符号运算
A = solve('a*u^2 + v^2', 'u - v = 1', 'a^2 - 5*a + 6') A= a: [1x4 sym] u: [1x4 sym] v: [1x4 sym] 这里A.a = [ 2, 2, 3, 3]
26.03.2019
8
MATLAB符号运算
v = [(1 + a/x)^x, exp(-x)]; limit(v,x,inf,'left') => [exp(a), 0] 2.4.级数和 Symsum(s,t,a,b) s中t从a到b求和。
概率统计与MATLAB精品PPT课件
功能:产生M lambda)
功能:计算分布密度p(x)在x的值
21.10.2020
x0 x0
7
§1 随机变量及其分布
均匀分布X~U(a,b) 命令1:Fx=unifcdf(x, a,b) 功能:计算累积概率Fx=P{X≤x}=F(x) 命令2:x=unifinv(p, a,b) 功能:计算随机量x,使得p=P{X≤x} 命令3:X=unifrnd(a,b,M,N) 功能:产生M*N维随机数矩阵X 命令4:Px=unifpdf(x, a,b) 功能:计算分布密度p(x)在x的值 补充:rand()---(0,1)均匀分布随机数
21.10.2020
12
§1 随机变量及其分布
例1某人向空中抛硬币100次,落下为正面的概率 为0.5。这100次中正面向上的次数记为X: (1)试计算x=45的概率和x≤45的概率; (2)绘制分布函数图象和分布列图象。
程序:》clear;
px=binopdf(45,100,0.5) % 计算x=45的概率
命令2:x=hygeinv(p,M, N,K)
功能:在已知参数M、N 、 K和p的情况下计算随 机量x,使得p=P{0≤次品数X≤x}
命令3:X=hygernd(M,N,K,m,n)
功能:在已知参数M,N ,K的情况下产生m*n维符合
超几何分布的随机数矩阵X
21.10.2020
2
§1 随机变量及其分布
21.10.2020
6
§1 随机变量及其分布
指数分布X~exp(λ)
1ex
P{Xx}
0
命令1:Fx=expcdf(x, lambda)
功能:计算累积概率Fx=P{X≤x}=F(x)
10MATLAB符号运算与概率统计.ppt
26.03.2019
8
MATLAB符号运算
v = [(1 + a/x)^x, exp(-x)]; limit(v,x,inf,'left') => [exp(a), 0] 2.4.级数和 Symsum(s,t,a,b) s中t从a到b求和。
26.03.2019
9
MATLAB符号运算
int([exp(t),exp(alpha*t)]) ans= [exp(t), 1/alpha*exp(alpha*t)] 2.3.极限 limit(F,x,a) limit(F,a)(x->a时F的极限) limit(F) (x->0时F的极限) limit(F,x,a,'right') limit(F,x,a,'left')
26.03.2019
4
MATLAB符号运算
2.2积分 R = int(S) R = int(S,v) R = int(S,a,b) R = int(S,v,a,b) 例如: int(-2*x/(1+x^2)^2) ans= 1/(1+x^2)
26.03.2019 5
26.03.2019
11
MATLAB符号运算
如: solve('a*x^2 + b*x + c') ans= [ 1/2/a*(-b+(b^2-4*a*c)^(1/2)), 1/2/a*(-b-(b^2-4*a*c)^(1/2))] solve('a*x^2 + b*x + c','b') ans= -(a*x^2+c)/x solve('x + y = 1','x - 11*y = 5') y = -1/3, x = 4/3
概率统计的数值实验MATLAB在概率统计教学中的应用-PPT精选
P A iA jA k P A k P A iA jA k n n 1 1 n 2 1 i j k n
PA1A2A3Ann1!
于是
n
P Ai i1
n
1, P Ai A j
1i jn
n 2 15
k 1
模拟Galton钉板试验的步骤: (1) 确定钉子的位置:将钉子的横、纵坐标存储在两个矩阵X和
Y中。 (2) 在Galton钉板试验中,小球每碰到钉子下落时都具有两种
4/5 0.5134
4/5 0.5086
4/5 0.5093
4/5 0.5093
π的近似值 3.1116 3.1165 3.1460 3.1418 3.1418
试验次数n
5千
1万
10万 100万 1000万
针长l/平行间 距d
相交频率
17/20 0.5432
17/20 0.5452
17/20 0.5420
概率
生日各不相同的概率 至少两人生日相同的概率 1
0.8
0.6
0.4
0.2
0 0 10 20 30 40 50 60 70 80 90 100 人数
• p1(30)=0.7063, p1(60)= 0.9941
分析:在30名学生中至少两人生日相同的概率为70.63%。 下面进行计算机仿真。
随机产生30个正整数,代表一个班30名学生的生日,然后观
解 记事件 A i 为第i个人拿到自已枪,事件 A i 为第i个人 没拿到自己枪,易知:
PAi
1 n
P Ai
n1 n
i1,2, ,n
又记 p 0 为没有一个人拿到自己枪的概率。
《概率统计》学习课件
134 第4章概率统计本章介绍MATLAB在概率统计中的若干命令和使用格式,这些命令存放于MatlabR12\Toolbox\Stats中。
4.1 随机数的产生4.1.1 二项分布的随机数据的产生命令参数为N,P的二项随机数据函数binornd格式R = binornd(N,P) %N、P为二项分布的两个参数,返回服从参数为N、P的二项分布的随机数,N、P大小相同。
R = binornd(N,P,m) %m指定随机数的个数,与R同维数。
R = binornd(N,P,m,n) %m,n分别表示R的行数和列数例4-1>> R=binornd(10,0.5)R =3>> R=binornd(10,0.5,1,6)R =8 1 3 7 6 4>> R=binornd(10,0.5,[1,10])R =6 8 4 67 5 3 5 6 2>> R=binornd(10,0.5,[2,3])R =7 5 86 5 6>>n = 10:10:60;>>r1 = binornd(n,1./n)r1 =2 1 0 1 1 2>>r2 = binornd(n,1./n,[1 6])r2 =0 1 2 1 3 14.1.2 正态分布的随机数据的产生命令参数为μ、σ的正态分布的随机数据函数normrnd格式R = normrnd(MU,SIGMA) %返回均值为MU,标准差为SIGMA的正态分布的随机数据,R可以是向量或矩阵。
R = normrnd(MU,SIGMA,m) %m指定随机数的个数,与R同维数。
R = normrnd(MU,SIGMA,m,n) %m,n分别表示R的行数和列数例4-2>>n1 = normrnd(1:6,1./(1:6))n1 =2.1650 2.31343.02504.0879 4.8607 6.2827>>n2 = normrnd(0,1,[1 5])n2 =0.0591 1.7971 0.2641 0.8717 -1.4462>>n3 = normrnd([1 2 3;4 5 6],0.1,2,3) %mu为均值矩阵n3 =0.9299 1.9361 2.96404.12465.0577 5.9864>> R=normrnd(10,0.5,[2,3]) %mu为10,sigma为0.5的2行3列个正态随机数R =9.7837 10.0627 9.42689.1672 10.1438 10.59554.1.3 常见分布的随机数产生常见分布的随机数的使用格式与上面相同表4-1 随机数产生函数表函数名调用形式注释Unifrnd unifrnd ( A,B,m,n) [A,B]上均匀分布(连续) 随机数Unidrnd unidrnd(N,m,n) 均匀分布(离散)随机数Exprnd exprnd(Lambda,m,n) 参数为Lambda的指数分布随机数Normrnd normrnd(MU,SIGMA,m,n) 参数为MU,SIGMA的正态分布随机数chi2rnd chi2rnd(N,m,n) 自由度为N的卡方分布随机数Trnd trnd(N,m,n) 自由度为N的t分布随机数Frnd frnd(N1, N2,m,n) 第一自由度为N1,第二自由度为N2的F分布随机数gamrnd gamrnd(A, B,m,n) 参数为A, B的γ分布随机数betarnd betarnd(A, B,m,n) 参数为A, B的β分布随机数lognrnd lognrnd(MU, SIGMA,m,n) 参数为MU, SIGMA的对数正态分布随机数nbinrnd nbinrnd(R, P,m,n) 参数为R,P的负二项式分布随机数ncfrnd ncfrnd(N1, N2, delta,m,n) 参数为N1,N2,delta的非中心F分布随机数nctrnd nctrnd(N, delta,m,n) 参数为N,delta的非中心t分布随机数ncx2rnd ncx2rnd(N, delta,m,n) 参数为N,delta的非中心卡方分布随机数raylrnd raylrnd(B,m,n) 参数为B的瑞利分布随机数weibrnd weibrnd(A, B,m,n) 参数为A, B的韦伯分布随机数binornd binornd(N,P,m,n) 参数为N, p的二项分布随机数geornd geornd(P,m,n) 参数为p的几何分布随机数hygernd hygernd(M,K,N,m,n) 参数为M,K,N的超几何分布随机数Poissrnd poissrnd(Lambda,m,n) 参数为Lambda的泊松分布随机数4.1.4 通用函数求各分布的随机数据命令求指定分布的随机数135函数random格式y = random('name',A1,A2,A3,m,n) %name的取值见表4-2;A1,A2,A3为分布的参数;m,n指定随机数的行和列例4-3 产生12(3行4列)个均值为2,标准差为0.3的正态分布随机数>> y=random('norm',2,0.3,3,4)y =2.3567 2.0524 1.8235 2.03421.9887 1.94402.6550 2.32002.0982 2.2177 1.9591 2.01784.2 随机变量的概率密度计算4.2.1 通用函数计算概率密度函数值命令通用函数计算概率密度函数值函数pdf格式Y=pdf(name,K,A)Y=pdf(name,K,A,B)Y=pdf(name,K,A,B,C)说明返回在X=K处、参数为A、B、C的概率密度值,对于不同的分布,参数个数是不同;name为分布函数名,其取值如表4-2。
matlab概率统计第2讲ppt课件
A(i:k , j:l) 返回由二维矩阵A中的第i到k行行向量和
第j到l列列向量组成的子阵。
17
e. 特殊矩阵
命令 hilb(n) toeplitz(k,r) pascal(n) rosser gallery(n) wilkinson(n) magic(n)
……
运行结果 生成n×n的希尔波特矩阵 生成非对称的托普利兹矩阵 pascal矩阵(帕斯卡矩阵) rosser矩阵 数字分析中有名的n×n试验矩阵 返回wilkinson特征值测试矩阵 魔方矩阵
阵,B为n×p矩阵,则C=A*B为m×p矩阵。
39
(3) 矩阵除法 在MATLAB中,有两种矩阵除法运算:\和/,分别表
示左除和右除.如果A矩阵是非奇异方阵,则A\B和B/A 运算可以实现.A\B等效于A的逆左乘B矩阵,也就是 inv(A)*B,而B/A等效于A矩阵的逆右乘B矩阵,也就是 B*inv(A).
1 4
7
%A的2,3行 %A的第一列
33
键入: A=[1 2 3;4 5 6;7 8 9]; C=A(1:2, [1 3]) %A的第1,2行, %A的第1,3列
输出:C= 13 46
还有A(1:2:3, 3:-1:1)
34
5) 矩阵的拼接
将几个矩阵接在一起称为拼接,左右拼 接行数要相同,上下拼接列数要相同。
15
c. 三角矩阵
命令
triu(A) triu(A,k) tril(A) tril(A,k)
运行结果
生成一个和A维数相同的上三角矩阵。该矩 阵主对角线及以上元素取自A中相应元素。 其余元素为0。
生成一个和A维数相同的上三角矩阵。该矩 阵第k条对角线及以上元素取自A中相应元 素。其余元素为0。
matlab课件--第5讲 概率统计实验
2. exprnd函数
例7: 产生4行5列的指数分布的随机数. 程序如下: y=exprnd(3,4,5) %参数=3
Matlab 软件实习
三、随机变量与概率分布密度
1. 几个常用的离散型分布密度函数(…pdf )
(1)均匀分布 P(X=xn)=1/n 密度函数调用格式:y=unidpdf(X,N) 例8:求X取值为1,2,3,4,5,6,7,8时服从均匀分布的概率值. 程序如下: X=1:8,N=8; Y=unidpdf(X,N)
% 参数SIGMA为正数
Matlab 软件实习
四、随机变量与概率分布函数
累积分布函数(…cdf )—在工具箱中分布函数亦称累积分布函 数,即表示事件的概率P{Xx}。
累积分布函数表 分布类型名称 函数名称 函数调用格式
离散均匀分布
二项分布 泊松分布 几何分布
unidcdf
binocdf poisscdf geocdf
程序如下: Y=[1500 2000 2500 3000]; P=[0.0952 0.0861 0.0779 0.7408]; EX=Y*P’
Matlab 软件实习 (2) 连续型 EX=int(x*f(x),-inf,inf)
例2:
1 , f (v ) a 0,
0va 其它
EV=int(v*1/a,0,a)
DV=int(v^2*1/a,0,a)-EV^2
Matlab 软件实习
3. 常见分布的期望与方差函数
分布类型名称 离散均匀分布 二项分布 几何分布 超几何分布 函数名称 unidstat binostat geostat hygestat 函数调用格式 [E,D]=unidstat(N) [E,D]=binostat(N,P) [E,D]=geostat(P) [E,D]=hygestat(M,K,N)
基于MATLAB的概率统计数值实验ppt课件
快捷的学习可借助MATLAB的系统帮助,通过指令doc 获得具体函数的详细信息,语法是 doc <函数名>
5/60
2. 二项分布实验
已知Y~b(20, 0.3)求Y分布率的值,并划出图形
在Matlab中输入以下命令:
binopdf(10,20,0.2) x=0:1:20; y=binopdf(x,20,0.2) plot(x,y,’r.’)
例9 某种重大疾病的医疗险种,每份每年需交保险费100元,若在 这一年中,投保人得了这种疾病,则每份可以得到索赔额10000元, 假设该地区这种疾病的患病率为0.0002,现该险种共有10000份保 单,问: (1)保险公司亏本的概率是多少? (2)保险公司获利不少于80万元的概率是多少?
15/60
(1) 每小时恰有4次呼叫的概率
(2) 一小时内呼叫不超过5次的概率 (3) 画出分布律图像
(1)
( 2)
P ( X 4)
4
4!
5
e
3k 3 P ( X 5) P ( X k ) e k 0 k 0 k!
5
34 3 e 4!
在Matlab中输入以下命令: (1)p1= poisspdf(4,3) (2)p2= poisscdf(5,3) (3)x=0:1:20;y=poisspdf(x,3);plot(x,y)
(2) σ=0.5, μ=1,2,3,4
(1)命令: x=-6:0.1:6; y1=normpdf(x,3,0.5); y2=normpdf(x,3,0.7); y3=normpdf(x,3,1); y4=normpdf(x,3,1.5); y5=normpdf(x,3,2); plot(x,y1,'.',x,y2,'+',x,y3,'*',x,y4,'d',x,y5)
matlab概率统计ppt课件
例7: 参数方程绘圆心在原点,半径为2的圆.
程序如下: t=linspace(-2,2,30)*pi; x=2*cos(t);y=2*sin(t); plot(x,y,’r’); text(-0.25,0,’x^2+y^2=4’); % 在图中标记曲线方程
2. 一元符号表达式函数绘图法
① 首先定义x是符号变量,再定义y是x的符号表达式函数; ② 用绘图命令ezplot或fplot绘图.
程序如下: t= -pi:pi/200:pi; comet(t,tan(sin(t))-sin(tan(t)))
(2) 三维动态轨线图
程序如下: t= 0:0.05:100; x=t;y=sin(t);z=sin(2*t); comet3(x,y,z)
4. 统计图形绘制
(1) 条形图
bar(x,y) %竖直条形图,其中x是横坐标向量,y是向量或矩阵 barh(x,y) %水平条形图 bar3(x,y) %三维条形图 bar3h(x,y) %三维水平条形图
例6 :在同一坐标内,分别用不同线型和颜色绘制曲线 y1=0.2e-0.5xcos(4x) 和y2=2e-0.5xcos(x),标记两曲线 交叉点.
程序如下:
x=linspace(0,2*pi,1000);
y1=0.2*exp(-0.5*x).*cos(4*pi*x);
y2=2*exp(-0.5*x).*cos(pi*x);
② ezplot(f,[a,b])表示在给定区间a<x<b上绘制函数f(x)的图形.
③ ezplot(f(x,y))表示在默认区域-2*pi<x<2*pi, -2*pi<y<2*pi上 绘制隐函数f(x,y)=0的函数图形.
用matlab计算各种概率分布(ppt)
x=0.01:0.1:8.01; y=fpdf(x,4,10); plot(x,y)
抽样分布: t 分布
设随机变量 X ~ N (0,1), Y ~ χ 2 ( n) ,且 X 与 Y 相 互独立,则称随机变量
wblplot(x)
统计绘图函数,进行 Weibull 分布检验。
Matlab相关命令介绍
其它函数
cdf 系列函数:累积分布函数 inv 系列函数:逆累积分布函数 rnd 系列函数:随机数发生函数 stat 系列函数:均值与方差函数
例: p=normcdf(-2:2,0,1)
离散分布: Poisson 分布
泊松分布也属于离散分布,是1837年由发个数 学家 Poisson 首次提出,其概率分布列为:
k! 记做:X ~ P ( λ )
) P( X = k =
λk
e −λ
k (=
0, 1, 2, , λ > 0 )
泊松分布是一种常用的离散分布,它与单位时间(或单 位面积、单位产品等)上的计数过程相联系。如:单位时 间内,电话总机接到用户呼唤次数;1 平方米内,玻璃上的 气泡数等。
指数分布举例
例: λ=4 时的指数分布密度函数图
x=0:0.1:30; y=exppdf(x,4); plot(x,y)
离散分布:几何分布
几何分布是一种常见的离散分布
在贝努里实验中,每次试验成功的概率为 p,设试验进行 到第 ξ 次才出现成功,则 ξ 的分布满足:
P (ξ = k = ) pq k −1
X T= Y /n
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
⑥ ezplot(x,y,[a,b])表示在给定区间a<t<b上绘制由参数方程x=x(t), y=y(t)确定的一元函数图形.
例9: 画出参数方程 x=sin3tcost , y=sin3tsint 在[0,]上 的图形.
例2: 绘制曲线. 程序如下: t=0:0.1:2*pi; x=t.*sin(3*t); y=t.*sin(t).*sin(t); plot(x,y)
(3) plot(x1,y1,x2,y2,…,xn,yn)
用这种形式可以在同一个窗口绘制多条曲线.
plot(x1,y1,‘选项1’,x2,y2,‘选项2’,…,xn,yn,‘选项n’)
(4) 具有两个纵坐标标度的图形 在MATLAB中,如果需要绘制出具有不同纵坐标
标度的两个图形,可以使用plotyy绘图函数.调用格式 为:
plotyy(x1,y1,x2,y2) 其中x1,y1对应一条曲线,x2,y2对应另一条曲线.横坐 标的标度相同,纵坐标有两个,左纵坐标用于x1,y1数 据对,右纵坐标用于x2,y2数据对.
例5: 用不同标度在同一坐标内绘制曲线 y1=0.2e-0.5xcos(4x) 和 y2=2e-0.5xcos(x).
程序如下: x=0:pi/100:2*pi; y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=2*exp(-0.5*x).*cos(pi*x); plotyy(x,y1,x,y2);
例14: 画三维曲线
x t
y
sin
t
.
z sin 2 t
程序如下: t=0:0.05:100; x=t; y=sin(t); z=sin(2t); plot3(x,y,z)
2. 空间曲面绘图法
(1) 一般曲面的绘制
首先将自变量(x,y)由向量给出,然后由[x,y]=meshgrid(x,y)构 成网格矩阵.
颜色 黄色 紫色 青色 红色 绿色 蓝色 白色 黑色
线型符号 : -. --
线型 实线 点线 点划线 虚线
数据点形式
.
o x + * ∨ ^
标记符号 点
圆圈 x标记 加号 星号 向下的三角形 向上的三角形
数据点形式 < > s d h p
标记符号 小于号
大于号 正方形
菱形 六角星 五角星
平面曲线的图形绘制
程序如下: syms t ezplot(sin(3*t)*cos(t),sin(3*t)*sin(t),[0,pi])
3. 一元函数极坐标绘图法
polar(t,r),其中t和r分别为角度向量和幅值向量, 且向量的长度相同.
例10: 绘图程序如下.
t=0:0.1:4*pi; r=(cos(t/6)+0.5) polar(t,r)
多元函数绘图
1. 空间曲线绘图法
自变量t为数组,x、y、z是以t为参数的函数数组: plot3(x,y,z)
x cos( t )
例13: 画空间螺旋线
Hale Waihona Puke ysin(t)
.
z t
程序如下: t=0:pi/50:10*pi; x=cos(t); y=sin(t); z=t; plot3(x,y,z) xlabel(‘x’),ylabel(‘y’)
4. 复函数绘图法
plot(z),其中z是复数,绘图是以z的实部为横坐标, 虚部为纵坐标绘制的关系曲线.
例11: 绘制 z=sin(t)+cos(t)i 的关系图形.
程序如下: t=0:0.1:2*pi; x=sin(t); y=cos(t); z=x+y*i; plot(z)
5. 分块绘图
分块绘图函数: subplot(m,n,p) 其作用为将图形窗口分割为m行n列的子窗口,然后选定 第p号子窗口为当前窗口.
程序如下: syms x y=f(x) ezplot(y)
或 fplot(‘y’,[a,b],’r*’)
例8: 画曲线 y e x .
程序如下:
syms x y=exp(x) ezplot(y)
ezplot的几种格式:
① ezplot(f)表示在默认区间-2*pi<x<2*pi上绘制函数f(x)的图形.
例6 :在同一坐标内,分别用不同线型和颜色绘制曲线 y1=0.2e-0.5xcos(4x) 和y2=2e-0.5xcos(x),标记两曲线 交叉点.
程序如下:
x=linspace(0,2*pi,1000);
y1=0.2*exp(-0.5*x).*cos(4*pi*x);
y2=2*exp(-0.5*x).*cos(pi*x);
1. 离散数据绘图法
① 首先定义自变量X 的取值向量; ② 再定义函数Y 的取值向量; ③ 用plot(x,y)命令给出平面曲线图.
在绘图参数中可以给出绘制图形的线型和颜色的参数.
基本的绘图命令
plot(x, y, ‘color-linestyle-marker’)
转折点
连线颜色、线型、数据点图标
程序如下: t= -pi:pi/200:pi; comet(t,tan(sin(t))-sin(tan(t)))
(2) 三维动态轨线图
程序如下: t= 0:0.05:100; x=t;y=sin(t);z=sin(2*t); comet3(x,y,z)
4. 统计图形绘制
(1) 条形图
bar(x,y) %竖直条形图,其中x是横坐标向量,y是向量或矩阵 barh(x,y) %水平条形图 bar3(x,y) %三维条形图 bar3h(x,y) %三维水平条形图
② 当输入参数有矩阵形式时,配对的x,y按对应列元 素为横、纵坐标分别绘制曲线,曲线条数等于矩阵的列 数.
例3: 绘制曲线. 程序如下: x1=0:0.1:2*pi; x2=1:0.1:3*pi; plot(x1,sin(x1),x2,cos(x2));
例4: 分析下列程序绘制的曲线. x1=linspace(0,2*pi,100); x2=linspace(0,3*pi,100); x3=linspace(0,4*pi,100); y1=sin(x1); y2=1+sin(x2); y3=2+sin(x3); x=[x1;x2;x3]'; y=[y1;y2;y3]'; plot(x,y,x1,y1-1)
y
Sine and COsine Curves 1
0.8 ? sinx
0.6
0.4
0.2
0
cosx ?
-0.2
-0.4
-0.6
-0.8
-1
0
1
2
3
4
5
6
7
x
① 当输入参数都为向量时,x1和y1,x2和y2,…,xn和 yn分别组成一组向量对,每一组向量对的长度可以不同. 每一向量对可以绘制出一条曲线,这样可以在同一坐标 内绘制出多条曲线.
例7: 参数方程绘圆心在原点,半径为2的圆.
程序如下: t=linspace(-2,2,30)*pi; x=2*cos(t);y=2*sin(t); plot(x,y,’r’); text(-0.25,0,’x^2+y^2=4’); % 在图中标记曲线方程
2. 一元符号表达式函数绘图法
① 首先定义x是符号变量,再定义y是x的符号表达式函数; ② 用绘图命令ezplot或fplot绘图.
(2) 球面的绘制
调用格式: [x,y,z]=sphere(r,n); mesh(x,y,z)
例17: 绘一球面.
程序如下: [x,y,z]=sphere(30); surf(x,y,z) axis(‘equal’) %调整纵横坐标比为1:1
3.动态图形
(1) 二维动态轨线图 调用格式: comet(x,y,p) 说明: 平面曲线y=y(x),其中p为尾长参数,缺省值为0.1. 例18: 画二维动态图.
(2) 柱面的绘制
柱面采用“母线”旋转生成.“母线”用向量r定义,旋转 圆周上的分格线条数用n定义,默认值n=20.
调用格式: [x,y,z]=cylinder(r,n); mesh(x,y,z)
例16: 旋转图.
程序如下: t=0:pi/2:3*pi; r=abs(5*t).*sin(t); [x,y,z]=cylinder(r,30); mesh(x,y,z) colormap([1 0 0]) %装入色图矩阵
k=find(abs(y1-y2)<1e-2); %查找y1与y2相等点 (近似相等)的下标
x1=x(k);
%取y1与y2相等点的x坐标
y3=0.2*exp(-0.5*x1).*cos(4*pi*x1); %求y1与y2值 相等点的y坐标
plot(x,y1,x,y2,'k:',x1,y3,'bp');
给出曲面函数矩阵:z=f(x,y)
有如下的绘制曲面函数: plot3(x,y,z) %线条图 mesh(x,y,z) %网格图 meshc(x,y,z) %具有基本等高线的网格图 meshz(x,y,z) %具有基准平面的网格图 surf(x,y,z) %表面图 surf(x,y,z) %表面图 surfc(x,y,z) %具有基本等高线的表面图 surfl(x,y,z) %具有关照效果的表面图 surface(x,y,z) %具有表面图到xoy面的投影图
例19: 画竖直条形图.
程序如下: x=1:12; y=[2 3.5 5 7 6 5 7.5 8 4.3 3 2.1 1.2]; bar(x,y)
(2) 直方图