- 1、本文档共70页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
[研究生入学考试]矩阵理论线性代数与计算1
线性代数与相关计算的Matlab应用
行列式(determinant)的定义与理论计算
行列式相应具有若干性质
行列式的理论计算
利用行列式的性质将行列式化成上(下)三角形行列式
行列式在解析几何中的两个基本应用
求平面三角形的面积
,其中为三角形三个顶点坐标,且按逆时针排列。
求空间四面体的体积
,其中为四面体四个顶点坐标,且按右手系法则排列。
行列式的Matlab计算
在Matlab命令窗口输入矩阵A,然后用det(A)
例如:
a=[2 3 4;6 7 8;9 10 11];
d=det(a)
再如
d=det(magic(5))
Matlab中如何计算行列式呢?
可在命令窗口键入doc det 。
解线性方程组
,
其矩阵形式为,
其中
2.1)按照克莱姆法则,可知当系数行列式时,方程组有解且是唯一的。解的形式
其中克莱姆法则的一个应用为
问题:过平面上的个不同观测点的次数不超过的多项
式存在且唯一。
假设过这些观测点的多项式为,于是有
由于
(范得蒙行列式)
所以得到问题的结论。但是如何来获得这个多项式呢?用Matlab做实验。
x=linspace(0,2*pi,5); %均匀生成[0,2*pi]区间上100个点
y=cos(x);
p2=polyfit(x,y,4); %用4阶多项式拟合cos(x)
xx=linspace(0,2*pi,100);
y2=polyval(p2,xx); %求4阶拟合多项式在生成点x的值
yy=cos(xx);
plot(x,y,go,xx,y2,b-,xx,yy,r-)
title(用4阶多项式拟合cos(x))
2.2)用求逆的办法,求解线性方程组(理论求法)
(行初等变换),然后有,
但是对于高阶的线性方程组,用理论求法无法实现,为此需要借助计算机来实现。
2.3)用 Matlab求解线性方程组以及方阵求逆
输入A,b,x=A\b即可,
例1
A=magic(5),b=[1,2,3,4,5];x=A\b
Ainv=inv(A),如果det(A)=0,可以求逆吗?例如A=magic(4),a=inv(A)会得到什么结果呢?
例2
求的解。
a=[1 1 0;-3 0 3;1 -3 4];b=[1 4 2];
x=a\b
此例也可以用符号运算求解:
syms x1 x2 x3; %定义符号变量
[x1,x2,x3]=solve(x1+x2=1,-3*x1+3*x3=4,x1-3*x2+4*x3=2)%solve函数还可以用于非线性解。
vap([x1,x2,x3],5) % 将解用5位有效数字表示
或
a=sym([1 1 0;-3 0 3;1,-3 4]) %定义符号变量
b=sym([1; 4 ;2]);% 定义符号变量
x=linsolve(a,b);x1=a\b; %线性方程组求解
vpa([x,x1],5)
向量、矩阵数据的生成,矩阵的代数运算以及与方阵的相关数值计算
3.1)向量数据生成
a)手工输入,例a=[1 3 7 8 10 -1];
b)命令输入,例x=linspace(-2,3,10);,x=-2:0.01:3;,x=logspace(-2,3,20);,
[x,y]=meshgrid(-1:0.01:,1,-2:0.01:2);
此种数据生成在作图时非常有用,例如:
x=linspace(-pi,pi,100);
y=sin(x);
plot(x,y,r-)
title(sin曲线)
另外,
[x,y]=meshgrid(-1:0.1:1);
z=x.^2+y.^2;
figure(1)
mesh(x,y,z);
figure(2)
surf(x,y,z)
c) 特殊矩阵生成
单位矩阵eye(n),零矩阵zeros(m,n),全一矩阵ones(m,n),魔方矩阵magic(n),pascal矩阵pascal(n),随机矩阵rand(m,n)和randn(m,n)。
d) 子矩阵的提取与整合
Matlab实现,例给定
a=magic(5);b=pascal(5);
al=a(2,:);%提取a的第二行;
ac=a(:,4);%提取a的第四列;
alc=a(2:3,1:4); %提取a的第二行到第三行以及第一列到第四列的子矩阵。
ab=[a,b];ac=[a;b];%组装两个矩阵,注意两者的区别。
ab=horzcat(a,b);ac=vertcat(a,b)%和上面两个运算相同。
A = magic(3);
B = [-5 -6 -9; -4 -4 -2];
C = eye(2
您可能关注的文档
- [研究生入学考试]02 第二章 GIS设计思想与方法.ppt
- [研究生入学考试]02-2光谱法仪器与光学仪器.ppt
- [畜牧兽医]第四章 蛋白质与氨基酸营养.ppt
- [畜牧兽医]鸡鸭猪畜牧养殖草种知识.doc
- [畜牧兽医]鸡新城疫.ppt
- [研究生入学考试]02-数据通信的基础知识.ppt
- [研究生入学考试]03—12英语考研真题完美打印.doc
- [研究生入学考试]02第一章古代建筑发展概况.ppt
- [研究生入学考试]03-07东财投资经济学考研真题.pdf
- [研究生入学考试]05第九章 第5节 隐函数的求导公式.ppt
- 《GB/T 32879-2025电动汽车更换用电池箱连接器》.pdf
- 中国国家标准 GB/T 21649.2-2025粒度分析 图像分析法 第2部分: 动态图像分析法.pdf
- 中国国家标准 GB/T 20899.9-2025金矿石化学分析方法 第9部分:碳量的测定.pdf
- 《GB/T 20899.9-2025金矿石化学分析方法 第9部分:碳量的测定》.pdf
- GB/T 20899.9-2025金矿石化学分析方法 第9部分:碳量的测定.pdf
- 《GB/T 33820-2025金属材料 延性试验 多孔状和蜂窝状金属高速压缩试验方法》.pdf
- GB/T 33820-2025金属材料 延性试验 多孔状和蜂窝状金属高速压缩试验方法.pdf
- 中国国家标准 GB/T 33820-2025金属材料 延性试验 多孔状和蜂窝状金属高速压缩试验方法.pdf
- GB/T 45910-2025信息技术 生物特征识别模板保护方案的性能测试.pdf
- 《GB/T 45910-2025信息技术 生物特征识别模板保护方案的性能测试》.pdf
文档评论(0)