![](https://img-blog.csdnimg.cn/20201014180756724.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
MATLAB
MATLAB 学习笔记
Encarta1993
Talk is cheap. Show me the code
展开
-
MATLAB学习笔记:总结
因为不太了解概率,这一课程学习便结束了总结这一段时间的学习,有很大的收获最后,感谢中国大学MOOC这一学习平台,北京理工大学 数学实验 课程 和 各位老师的认真讲授课程链接如下:https://www.icourse163.org/learn/BIT-1001756010?tid=1002116020原创 2018-01-27 19:50:23 · 381 阅读 · 0 评论 -
MATLAB学习笔记:区间估计
若已经知道了一组数据来自正太分布总体,但是不知道正态分布总体的函数。这时可以利用normfit()命令来完成对参数的点估计和区间估计。调用格式:[muhat,sigmahat,muci,sigmaci]=normfit(X,alpla)此命令以alpha为显著性水平,在数据X下,对参数进行估计(alpha缺省时设定为0.05).返回值中muhat是正态分布的均值的点估计值s原创 2018-01-27 16:23:01 · 9171 阅读 · 1 评论 -
MATLAB学习笔记:极大似然估计
极大似然估计的步骤:1、写出似然函数2、对似然函数取对数,并整理3、求导数4、解似然方程极大似然估计的Matlab命令mle调用格式:phat=mle(data)返回服从正态分布的数据参数的极大似然估计。phat=mle(data,'distribution','dist')返回服从dist确定的数据参数的极大似然估计。原创 2018-01-27 16:02:51 · 40553 阅读 · 3 评论 -
MATLAB学习笔记:矩估计
参数估计:1、点估计2、区间估计点估计:点估计问题就是要构造一个只依赖于样本的量,作为未知参数或未知参数的函数的估计值。点估计的主要方法有:矩估计法、极大似然估计法、最小二乘法、贝叶斯估计法等。准则特点:无偏性、有效性、一致性。矩估计:>> x=[1.2 3.5 4.2 0.8 1.4 3.1 4.8 0.9];>> u=mean(x);原创 2018-01-27 15:49:58 · 6092 阅读 · 0 评论 -
MATLAB学习笔记:常用统计量2
峰度:kurtosis()>> x=normrnd(0,1,1000,1);>> y=exprnd(5,1000,1);>> z=unifrnd(-1,1,1000,1);>> kurtosis(x)ans = 3.0214>> kurtosis(y)ans = 8.8943>> kurtosis(z)ans = 1.7777原创 2018-01-26 19:49:57 · 763 阅读 · 0 评论 -
MATLAB学习笔记:常用统计量1
常用统计量:平均数、中位数、方差、标准差、极差、偏度、峰度等。平均值与中位数:平均值 mean()中位数 median()>> a=[1 5 3 6 2];>> junzhi=mean(a)junzhi = 3.4000>> zhongweishu=median(a)zhongweishu = 3>> A=[7 4 3 7;1 4 3 7原创 2018-01-26 19:31:28 · 1786 阅读 · 0 评论 -
MATLAB学习笔记:随机模拟法计算数值积分
>> x=unifrnd(-2*pi,2*pi,10000,1);>> f=(1-x.^2).*sin(5*x);>> fmax=max(f)fmax = 34.711289618412117随机模拟法求定积分:数值积分:>> f=inline('exp(x.^2+1)','x');>> quad(f,0,1)ans原创 2018-01-25 19:40:32 · 7245 阅读 · 1 评论 -
MATLAB学习笔记:随机数
产生随机数:random('name',A1,A2,A3,m,k)name是相应分布的名称,比如poisson,normalA1,A2,A3是该分布中的参数m是随机数的行k是随机数的列>> x=random('poisson',3,1,10)x = 3 2 4 3 5 1 3原创 2018-01-24 18:54:10 · 424 阅读 · 0 评论 -
MATLAB学习笔记:随机变量的逆累计分布函数
正态分布:>> y=[0.1 0.3 0.5 0.9];>> norminv(y,0,2)ans = -2.5631 -1.0488 0 2.5631>> ezplot('norminv(y,0,2)')>> y=[0.1 0.2 0.3 0.9];>> binoinv(y,原创 2018-01-24 18:27:10 · 4998 阅读 · 1 评论 -
MATLAB学习笔记:连续型随机变量
均匀分布:x=0:0.1:10;f=unifpdf(x,3,7);F=unifcdf(x,3,7);plot(x,f,x,F)legend('概率密度函数','分布函数')指数分布:x=-1:0.1:30;f=exppdf(x,2);F=expcdf(x,2);plot(x,f,x,F原创 2018-01-24 18:13:20 · 4607 阅读 · 0 评论 -
MATLAB学习笔记:离散型随机变量
排列与组合function y=pailie(n,k)y=factorial(n)/factorial(n-k);function y=zuhe(n,k)y=pailie(n,k)/factorial(k);>> p=zuhe(5,2)*z原创 2018-01-23 20:18:51 · 5114 阅读 · 0 评论 -
MATLAB学习笔记:二次型
化二次型为标准型可以用schur和eig实现:[Q,D]=schur(A)[Q,D]=eig(A)其实,A为二次型的矩阵,D为A的特征值构成的对角矩阵,Q为正交矩阵。>> A=[-2 2 2 2;2 -6 0 0;2 0 -9 3;2 0 3 -9];>> [Q,D]=schur(A)Q = 0.0000 -0.5原创 2018-01-23 10:51:32 · 16204 阅读 · 1 评论 -
MATLAB学习笔记:方阵的相似对角化
>> A=[11 -6 4 -10 -4;-3 5 -2 4 1;-8 12 -3 12 4;1 6 -2 3 -1;8 -18 8 -14 -1];>> [V D]=eig(A)V = -0.3244 -0.4983 -0.7759 -0.2343 -0.5752 0.1622 0.1878 0.0887 -0.0186 0.3890原创 2018-01-21 19:07:01 · 21646 阅读 · 1 评论 -
MATLAB学习笔记:方阵的特征值与特征向量
>> A=[1 2;2 1];>> eigshow(A)MATLAB计算特征值和特征向量的命令:d=eig(A) 仅计算A的特征值(以向量方式d存放)[V,D]=eig(A) D是由A的特征值构成的对角矩阵,V是由特征向量按列排成的矩阵,AV=VD。>> A=[11 -6 4 -10 -4;-3 5原创 2018-01-21 15:55:31 · 8160 阅读 · 1 评论 -
MATLAB学习笔记:行列式及其应用
行列式的计算det(A)>> A=[7 -3 2 -3;11 -2 3 -4;-5 4 -6 3;3 -9 2 -5];>> D=det(A)D = 138.0000>> syms a b c d;>> A=[a 2*b a+2*b;a a+3*c a+3*c;a a+4*原创 2018-01-20 19:37:10 · 1701 阅读 · 0 评论 -
MATLAB学习笔记:线性方程组的几何意义与应用
线性方程组的几何意义:>> B1=[2 1 3;3 -1 2];>> B2=[1 2 -5;3 6 2];>> B3=[1 2 3;2 4 6];>> B4=[1 2 -5;3 6 2;1 -2 5];>> U1=rref(B1)U1 = 1 0 1 0 1原创 2018-01-20 12:54:38 · 1770 阅读 · 0 评论 -
MATLAB学习笔记:非齐次线性方程组的求解
>> A=[1 -1 1 -2;2 0 -1 4;3 2 1 0;-1 2 -1 2];>> b=[2 4 -1 -4]';>> r=rank(A)r = 4>> r1=rank([A,b])r1 = 4>> X=inv(A)*bX = 1.0000 -2.0000 0 0.5000方程组有唯一解。原创 2018-01-19 19:42:55 · 35597 阅读 · 2 评论 -
MATLAB学习笔记:齐次线性方程组的基础解系
求齐次线性方程组的基础解系的命令:B=null(A) B的列向量是AX=0的规范正交的基础解系B=null(A,'r') B的列向量是AX=0的有理数形式的基础解系>> A=[1 -1 -1 0 3;2 -2 -1 2 4;3 -3 -1 4 5;1 -1 1 4 -1];>> B1原创 2018-01-19 18:21:17 · 44479 阅读 · 1 评论 -
MATLAB学习笔记:向量空间基础
1、求向量组的秩和极大线性无关组>> a1=[1;2;2;3];>> a2=[1;4;-3;6];>> a3=[-2;-6;1;-9];>> a4=[1;4;-1;7];>> a5=[4;8;2;9];>> A=[a1,a2,a3,a4,a5];>> r=rank(A)r = 3向量组的秩为3,因为r>> [R,j]=rref(A)R = 1 0原创 2018-01-18 20:31:16 · 5714 阅读 · 0 评论 -
MATLAB学习笔记:矩阵基础
矩阵的秩和简化阶梯型:求矩阵A的秩: rank(A)求A的简化阶梯型: rref(A)[R,j]=rref(A) R是A的简化阶梯型,j是一个向量,表示主元所在的列的标号。解1:>> A=[1 1 -2 1 4;2 4 -6 4 8;2 -3 1 -1 2;3 6 -9 7 9];>> r=rank(A)r = 3>> a=rref(A)a = 1 0原创 2018-01-18 18:35:43 · 866 阅读 · 0 评论 -
MATLAB学习笔记:插值方法
分段线性插值:interp1(x0,y0,x)三次样条插值:interp1(x0,y0,x,'spline') 或者 spline(x0,y0,x)分段三次插值:interp1(x0,y0,x,'cubic')其中,x0和y0是已知的节点数组,x是待计算的插值点数组>> x0=-5:0.5:5;>> y0=1./(1+x0.^2);>> x=-5:0.1:5;>> y1=interp1(x原创 2018-01-17 20:45:28 · 1285 阅读 · 0 评论 -
MATLAB学习笔记:数列求和与级数
数列求和命令一:sum(x)例1:x为向量:>> a=[1 2 3]a = 1 2 3>> sum(a)ans = 6例2:x为矩阵:>> b=[1 2 3;4 5 6;7 8 9]b = 1 2 3 4 5 6 7原创 2018-01-16 20:26:38 · 30176 阅读 · 0 评论 -
MATLAB学习笔记:曲面积分
>> syms x y;>> s=sqrt(3)*int(int(2-y,y,0,1-x),x,0,1) s = (5*3^(1/2))/6>> syms y z;>> s=int(int(sqrt(1-y^2),y,0,1),z,0,3) s = (3*pi)/4原创 2018-01-13 19:59:13 · 5187 阅读 · 0 评论 -
MATLAB学习笔记:曲线积分
>> syms theta a k;>> z=a^2*k*sqrt(a^2+k^2)/2*int(theta*sin(2*theta),0,pi) z = -(pi*a^2*k*(a^2 + k^2)^(1/2))/4>> syms t a;>> x=a*(t-sin(t));>> y=a*(1-cos(t));>原创 2018-01-13 19:50:43 · 9539 阅读 · 0 评论 -
MATLAB学习笔记:重积分
>> syms x y;>> z=int(int(y*sqrt(1+x^2-y^2),y,x,1),x,-1,1) z = 1/2 >> >> syms x y;>> z=int(int(sin(y^2),x,0,y),y,0,1) z = sin(1/2)^2如果先对y求积分:>> syms x y;原创 2018-01-12 20:23:10 · 7115 阅读 · 2 评论 -
MATLAB学习笔记:数值积分
当(1)被积函数的原函数不能用初等函数表示。(2)被积函数难以用公式表示,而是用图形或表格给出的。就应该建立定积分的近似计算方法:数值积分方法。梯形法:z=trapz(x,y)>> x=0:0.5:1;>> y=exp(-x.^2);>> z=trapz(x,y)z = 0.7314>> x=0:0.05:1;原创 2018-01-11 20:37:41 · 68449 阅读 · 7 评论 -
MATLAB学习笔记:定积分与广义积分
计算定积分和广义积分的命令int调用格式:一元函数定积分:int(f(x),a,b)二元函数定积分:int(f(x,y),x,a,b)>> syms x;>> f=sin(x);>> s=int(f,0,pi)s = 2>> syms x;>> f=abs(x-1);>> s=int(f,0,2)s = 1原创 2018-01-10 20:04:11 · 8833 阅读 · 2 评论 -
MATLAB学习笔记:不定积分
不定积分的定义:若函数F(x)在区间 I 上满足F'(x)=f(x),则称F(x)是f(x)在 I 上的一个原函数。函数f(x)在区间 I 上的原函数全体,称为f(x)在区间 I上的不定积分。即:不定积分:int(f(x))二元不定积分:int(f(x,y),x)输入命令:>原创 2018-01-09 21:24:20 · 21140 阅读 · 3 评论 -
MATLAB学习笔记:求数值微分
数值微分是用离散方法近似计算函数的导数值或偏导数值。向前差商公式:向后差商公式:中心差商公式:二阶导数的中心差商公式:例1:利用数值微分求y=4*x^2+3*sin(x)在x=1处的近似导数。先利用导函数求出其一阶和二阶原创 2018-01-08 16:14:00 · 20373 阅读 · 1 评论 -
MATLAB学习笔记:常微分方程的数值解
常微分方程数值求解的命令:求常微分方程的数值解,MATLAB的命令格式为:[t,y]=solver('odefun',tspan,y0,options)其中solver选择ode45等函数名,odefun为根据待解方程或方程组编写的m文件名,tspan为自变量的区间[t0,tf],即准备在那个区间上求解,y0表示初始值,options用于设定误差限制。命令格式为:optio原创 2018-01-04 21:47:54 · 23463 阅读 · 0 评论 -
MATALB学习笔记:读取和写入Excel
读取Excel:xlsread('Excel路径+Excel名称')如果Excel文件和编写程序的文件在同一文件夹,不需要填写路径。例1:a=xlsread('工作簿1')例2:读取指定的sheet b=xlsread('工作簿1','sheet2')例3:读取指定的sheet中的指定的行或列c1=xlsread('工作簿1','sheet2','a1:a3原创 2018-01-03 20:18:37 · 1296 阅读 · 0 评论 -
MATLAB学习笔记:解常微分方程
求常微分方程(组)的通解:例1:dsolve('Dx+4*x=sin(2*t)')自变量默认为t。例2:dsolve('Dx=t')和syms x;dsolve('Dy=x',x)结果相同。例3:dsolve('D2x+2*Dx-3*x=exp(t)')例4:方程组syms t;[x,y]=dsolve('Dx=-3*y','Dy原创 2017-12-29 20:37:58 · 9408 阅读 · 1 评论 -
MATLAB学习笔记:导数的应用2
验证洛必达法则:例1: syms x;f=3^x-2^x;g=x;limit(f/g,x,0)syms x;f=3^x-2^x;g=x;limit(diff(f,x)/diff(g,x),x,0)其结果相同,都为:log(3) - log(2)验证泰勒展开式:例2:syms x f;f=exp(x)*sin(x);t2=taylor(f,x,原创 2017-12-28 19:38:40 · 427 阅读 · 0 评论 -
MATLAB学习笔记:导数的应用1
求方程(组)的根:例1:syms x;f=x^2+3*x+2;solve(f,x)syms x y b;f=x^2+6*y+2;g=x+y-b;[x,y]=solve(f,g,x,y)例2:syms x;f=5*sin(2*x)-exp(x);solve(f,x)当不能求得符号解时,返回数值近似解代替。求函数在一定范围内的零点:求函数在点x原创 2017-12-27 20:23:54 · 577 阅读 · 0 评论 -
MATLAB学习笔记:求偏导
求多元函数的偏导数计算多元函数关于某个变量的偏导数:例1:u=x^2+ln(y)+sqrt(z),求du/dx,du/dz。syms x y z;diff(x^2+log(y)+sqrt(z),x) syms x y z;diff(x^2+log(y)+sqrt(z),z)计算多元函数关于某个变量的n阶偏导数:例2:syms x y z;diff(x^2+原创 2017-12-26 19:32:52 · 39819 阅读 · 1 评论 -
MATLAB学习笔记:求导数
求一元函数的导数:例1:求y=3*sin(x)+4*x^2的导数syms x;diff(3*sin(x)+4*x^2)例2:求y=ln(x+cos(x))syms x;diff(log(x+cos(x)))求一组一元函数的导数:例3:求f(x)={sin(x)+1; cos(x)+x;原创 2017-12-25 20:16:49 · 29762 阅读 · 1 评论 -
MATLAB学习笔记:求极限
(1)函数趋于a点的极限: limit(f,x,a) 或 limit(f,a)(2)函数趋于a点的右极限: limit(f,x,a,'right')(3)函数趋于a点的左极限: limit(f,x,a,'left')(4)函数趋于无穷的极限: limit(f,x,inf)注:1.limit(f)=limit(f,x,0) 2.求数列极限可以视为求函数极限例1:函数f(原创 2017-12-20 20:52:46 · 11095 阅读 · 1 评论 -
MATLAB Study
开始学习了原创 2017-11-23 16:09:32 · 252 阅读 · 0 评论