Matlab与线性代数

文章目录

  • 多项式求解
    • 1.2 多项式四则运算
    • 1.3 多项式的分解与合并
    • 行列式求解
  • 3、矩阵基本运算
    • ➢ 3.2 矩阵的取块和变换
    • ➢ 3.3 矩阵的基本运算
    • 4、求解线性方程组

多项式求解

➢ 1.1 多项式表达式与根
有关多项式函数表达式与根的Matlab命令:
poly2sym§ 返回由多项式系数转为多项式函数p(x)表达式
polyval(p,a) 返回多项式函数p(x)当x=a时的值
roots§ 返回多项式函数p(x)的所有复数根
poly® 返回由根组成的向量r创建的多项式函数p(x)的系

p=[1 0 -2 3]
px=poly2sym(p) 

px=x^3-2*x+3

1.2 多项式四则运算

有关多项式函数四则运算的Matlab命令:
conv(p1,p2) 返回多项式函数p1(x)和p2(x)的乘积
[q r]=deconv(p1,p2) 返回多项式函数p1(x)除了p2(x)的商式q(x)和余式r(x)

clear
p1=[1 0 -2 1];
p2=[-1 0 1 0 -2 3];
m= length(p1); n=length(p2);
t= max(m, n);
p1=[zeros(1,t-m),p1]
p2=[zeros(1,t-n),p2]
p=p1+p2

结果为
p= -1 0 2 0 -4 4
p(x)=-x5+2x3-4x+4

例1-3: 求多项式𝑝 𝑥 = −𝑥8 + 3𝑥6 − 𝑥5 − 4𝑥4 + 4𝑥3 + 4𝑥2 − 8𝑥 + 3分别被多项式p1 𝑥 = 𝑥3 − 2𝑥 + 1, 𝑝2 𝑥 = 𝑥3 − 2𝑥 + 3相clear

p=[-1 0 3 -1 -4 4 4 -8 3];
p1=[1 0 -2 1];
p2=[1 0 -2 3];
[q1,r1]=deconv(p,p1)
[q2,r2]=deconv(p,p2)

结果为
q1= -1 0 1 0 -2 3
r1= 0 0 0 0 0 0 0 0 0
q2= 1 0 1 2 -2 5
r2= 0 0 0 0 0 0 -6 8 -12除后的结果

1.3 多项式的分解与合并

有关多项式的分解与合并的Matlab命令:

syms x
collect(f) 对符号多项式f进行合并同类项
expand(f) 对符号多项式f进行展开
horner(f) 对符号多项式f进行嵌套分解
factor(f) 对符号多项式f进行因式分解

例1-4 合并同类项𝑓1 = 𝑥 − 1 𝑥 − 2 𝑥 − 3 , 𝑓2 = 1 + 𝑥 𝑡 + 𝑡𝑥。

syms x t
f1=(x-1)*(x-2)*(x-3);f2=(1+x)*t+t*x;
p1=collect(f1)
p2=collect(f2)

结果为
p1=-6+x ^ 3 - 6 * x ^ 2 +11 * x
p2=2 * t * x +t

行列式求解

求解行列式的Matlab命令:
det(A) 计算A的行列式值,A为数值或符号方阵

D=[2 -3 -1 2; 1 -5 3 -4; 0 2 1 -1; -5 1 3 -3];
det(D)

结果为
-75

3、矩阵基本运算

➢ 3.1 矩阵的生成
矩阵除了直接在命令窗口键入,也可以从外部数据文件(*.mat)和自编M文件
导入,还可以利用Matlab系统内部提供的一些函数生成,常见的函数如下:
zeros(m,n)函数 生成m行n列全部元素为0的矩阵
ones(m,n)函数 生成m行n列全部元素为1的矩阵
rand(m,n)函数 生成m行n列全部在0到1的均匀分布随机元素的矩阵
randn(m,n)函数 生成m行n列全部为标准正态分布随机元素的矩阵
magic(n)函数 生成n阶幻方方阵(魔阵),即每行每列和对角线上元素之和
相等
diag(M)函数 从矩阵M对角线元素生成一列数组或将一列数组M生成一个
对角阵
triu(M)函数 取矩阵M的对应元素生成上三角矩阵
tril(M)函数 取矩阵M的对应元素生成下三角矩阵
length(M)函数 返回向量M的长度
size(M)函数 返回矩阵M的行数和列数
eye(n)函数 生成n阶单位阵
hilb(n)函数 生成n阶希尔伯特病态矩阵
pascal(n)函数 生成n阶的帕斯卡矩阵

➢ 3.2 矩阵的取块和变换

在矩阵运算中,有时需要提取其中的一部分元素参与运算,比如提取某个元素,
某行,某列甚至某一个子阵,还有将原来的矩阵改变形式,常见的命令如下:
A(i ,:) 提取矩阵A的第i行 A(:, j) 提取矩阵A的第j列 A(:) 将矩阵A的各列元素依次排成一列向量
A(i;j) 将矩阵A(:)中的第i个到第j个的元素依次排成一行向量
A(i:j,:) 提取矩阵A的第i行到第j行的所有元素所成的矩阵
A(:,i:j) 提取矩阵A的第i列到第j列的所有元素所成的矩阵
A(i:j,k:l) 提取矩阵A的第i行到第j行以及第k列到第l列的所有元素所成的子阵
B=reshape(A,m,n) 将矩阵A的元素依次组成一个m行n列的矩阵B
B=rot90(A) 将矩阵A逆时针旋转90度得到新矩阵B
B=fliplr(A) 将矩阵A 中的元素左右对称得到新矩阵B
3.2 矩阵的取块和变换
B=flipud(A) 将矩阵A 中的元素上下对称得到新矩阵B
flipdim(A,1)=flipud(A) flipdim(A,2)=fliplr(A)
A(2,:)=[] 将矩阵A的第二行删除
A(A>5)=1 将矩阵中的所有大于6的元素赋值为

➢ 3.3 矩阵的基本运算

在矩阵运算中常见的命令如下:
𝐴 ± 𝐵 矩阵A加减矩阵B
A+k 矩阵A的所有元素加上数k
AB 矩阵A乘以矩阵B
K
A,A*k 矩阵A的所有元素乘以数k
A.*B 将矩阵A的各个元素对应和矩阵B中的元素相乘
A. /B 将矩阵A的各个元素对应和矩阵B中的元素相除
A/B 右除
A\B 左除
A.’ 得到矩阵A的转置
A’ 得到矩阵A共轭转置,在实数域内就是转置
inv(A)或A^(-1) 得到矩阵A的逆矩阵,A应该为可逆阵
A^k 矩阵A的k次幂
sqrtm(A)或A^(1/2) 矩阵A的开方,其平方为A
sqrt(A) 矩阵A的对应元素开方

4、求解线性方程组

rank(A) 得到矩阵A的秩
rref(A) 得到矩阵A的行最简形
null(A) 得到系数矩阵为A的齐次方程组基础解系
Null(A, ‘r’) 得到系数矩阵为A的齐次方程组有理数形式的基础解系
trace(A) 得到矩阵A的迹
poly(A) 得到矩阵A的特征多项式系数
[a,b]=eig(A) 得到矩阵A的特征列向量矩阵a和对应特征值组成的对角阵b
B=orth(A) 正交化空间,即矩阵B的列向量正交且生成的线性空间与矩阵A的列向量
生成的线性空间等价
r = roots(p)返回由p表示的多项式的根作为列向量。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318356.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【 CF1186D,E,F】Vus the Cossack and Numbers/Vus the Cossack and a Field/Vus the Cossack and a Graph

太ex了,哭了哭了orz 后面两道平均一道花了我一天啊! 文章目录D:Vus the Cossack and Numbers题意翻译题解代码实现E:Vus the Cossack and a Field题意翻译题解代码实现F:Vus the Cossack and a Graph题目暴力题解代码实现官方题解…

IdentityServer4与ocelot实现认证与客户端统一入口

关于IdentityServer4与ocelot博客园里已经有很多介绍我这里就不再重复了。ocelot与IdentityServer4组合认证博客园里也有很多,但大多使用ocelot内置的认证,而且大多都是用来认证API的,查找了很多资料也没看到如何认证oidc,所以这里…

CF1556D-Take a Guess【交互】

正题 题目链接:https://codeforces.com/contest/1556/problem/D 题目大意 现在有nnn个你不知道的数字,你有两种询问操作 询问两个下标的数字的andandand询问两个下标的数字的ororor 要求在2n2n2n次操作以内求出第kkk小的数字 1≤n≤104,0≤ai≤1091\leq n\leq 1…

YBTOJ:彩球抽取(期望)

文章目录题目描述解析代码题目描述 解析 首先,可以使用dp解决本题 设fi,j,k:操作i轮之后编号j的小球有k个的概率 转移和统计答案就都不难了 但是还有一个问题 不难发现这个题循环下去是可以无穷无尽的 所以限定一个i的上界(如500000&#xf…

魔改森林

题意: 曾经有一道叫做迷雾森林的题目,然而牛牛认为地图中的障碍太多,实在是太难了,所以删去了很多点,出了这道题。 牛牛给出了一个n行m列的网格图 初始牛牛处在最左下角的格点上(n1,1),终点在右上角的格点…

基于IdentityServer4 实现.NET Core的认证授权

IdentityServer4是什么?IdentityServer4是基于ASP.NET Core实现的认证和授权框架,是对OpenID Connect和OAuth 2.0协议的实现。OpenID Connect 和 OAuth2.0是什么OpenID Connect:OpenID Connect由OpenID基金会于2014年发布的一个开放标准, 是建立在OAuth …

常见存储、查找算法

存储 散列存储&#xff1a;即哈希的存储方式。 顺序存储&#xff1a;数组的存储方式 链式存储&#xff1a;链式前向星、vector<> 压缩存储 索引存储 查找 常见查找算法 顺序查找 一个一个往下找&#xff0c;复杂度 \(O(\dfrac{n1}{2})\) 。 适合顺序存储&#xff0c…

CF1556E-Equilibrium【栈,树状数组】

正题 题目连接:https://codeforces.com/contest/1556/problem/E 题目大意 两个长度为nnn的序列a,ba,ba,b&#xff0c;qqq次询问一个区间[l,r][l,r][l,r]。 在这个区间中你每次可以选择一个长度为偶数的下标递增的序列&#xff0c;让奇数位置的aaa加一&#xff0c;偶数位置的…

[COCI 2018#5]Parametriziran

这道题呢&#xff01; 算了&#xff0c;不要让这玩意儿活着祸害众生吧&#xff01;让我们来拯救苍生于苦海之中&#xff01;&#xff01; 骚话连篇ing 题目 由小写英文字母和问号组成的字符串成为参数化单词&#xff08;例如&#xff1a;??cd,bcd,??&#xff09;。如果两…

P2324 [SCOI2005]骑士精神(迭代加深搜索,dfs)

传送门 文章目录解析解析 很显然&#xff0c;让马走的话状态记录和转移都会比较复杂 所以转化成让空位跳会更好做一点 但这不是重点 初看本题&#xff0c;其实第一感觉是bfs 但是状态数理论上最差可以达到815&#xff0c;&#xff08;当然基本不可能跑满&#xff09;&#xff…

AcWing 1402. 星空之夜 1月28

AcWing 1402. 星空之夜 1月28 题意&#xff1a; 一个星群是指一组非空的在水平&#xff0c;垂直或对角线方向相邻的星星的集合。 一个星群不能是一个更大星群的一部分。 星群可能是相似的。 如果两个星群的形状、包含星星的数目相同&#xff0c;那么无论它们的朝向如何&#…

RMQ、ST表

ST表 \(\text{ST}\) 表是用于解决可重复贡献问题的数据结构。 可重复贡献问题&#xff1a;区间按位和、区间按位或、区间 \(\gcd\) 、区间最大、区间最小等满足结合律且可重复统计的问题。 模板预处理&#xff1a;(以区间最大值为例) void pre_work() {for(int i2;i<n;i) lg…

NumSharp v0.6 科学计算库发布,新增 LAPACK 的线性库支持

NumSharp&#xff08;Numerical .NET&#xff09;可以说是C&#xff03;中的科学计算库。 它是用C&#xff03;编写的&#xff0c;符合.netstandard 2.0库标准。 它的目标是让.NET开发人员使用NumPy的语法编写机器学习代码&#xff0c;从而最大限度地借鉴现有大量在python代码的…

[COCI] Zamjena

连这种模拟题都能。。。orz ex&#xff0c;太恶心了&#xff01; 驰骋坑底这么久了&#xff0c;我明白了 开始吧&#xff01;我发誓&#xff0c;这个超级兵&#xff0c;我就算用小书包平A都要A了它 题目 Vlatko喜欢使用整数数组&#xff0c;他在一张纸上写下了两个数组&…

P7470-[NOI Online 2021 提高组]岛屿探险【Trie,CDQ分治】

正题 题目链接:https://www.luogu.com.cn/problem/P7470 题目大意 给出nnn个二元组(a,b)(a,b)(a,b)。 qqq次询问给出(l,r,c,d)(l,r,c,d)(l,r,c,d)表示询问[l,r][l,r][l,r]中有多少二元组满足cxora≤min(b,d)c\ xor\ a\leq min(b,d)c xor a≤min(b,d)。 1≤n,q≤1051\leq n,q\…

P2601 [ZJOI2009]对称的正方形(二维哈希)(二分)

洛谷传送门 文章目录题目描述解析代码题目描述 解析 做三个hash 分一下正方形边长的奇偶性 然后枚举中心点&#xff0c;二分边长即可 有点类似模拟赛那道红十字的题 我一开始觉得分奇偶好麻烦啊 为什么不直接枚举左上方的点二分呢&#xff1f;awa 很遗憾的是… 那样答案就没有…

初赛—错题集

计算机基础知识 LAN&#xff1a;局域网&#xff0c;WAN&#xff1a;广域网&#xff0c;MAN&#xff1a;城域网 汇编语言是(依赖于具体计算机)的低级程序设计语言 计算机操作的最小时间单位是(时钟周期)。 注意所需空间需要 \(\div 8\) &#xff01;&#xff01;&#xff01;…

.NET Core 和 DevOps

关键要点无论你目前使用什么样的技术栈&#xff0c;DevOps 都是值得一试的。闭源、专有软件和构建过程与 DevOps 实践不兼容。.NET Core 是开源的&#xff0c;是基于 DevOps 构思和构建的。.NET Core CLI 和 Roslyn API 让整个交付流程变得更加开放&#xff0c;且具有更强的适应…

2021 NOI游记

文章目录前言day1总结T1 轻重边(0/50)考后得分T2 路径交点(0/20)考后得分T3 庆典(28/44)赛后得分day2总结T1 量子通信(20/12)赛后得分T2 密码箱(0/35)赛后得分T3 机器人游戏(12/12)赛后得分总结前言 62pts 太惨了 &#xff08;只是参加了同步赛而已&#xff09; 感觉心态受到了…

【每日一题】1月29日题目 和与或

题意&#xff1a; 给你一个数组R&#xff0c;包含N个元素&#xff0c;求有多少满足条件的序列A使得 0 ≤ A[i] ≤ R [ i ] A[0]A[1]…A[N−1] A[0] | ]A[1]… | A [ N − 1 ] 输出答案对1e99取模 题解&#xff1a; 参考博客 数位dp问题 如果和等于或的话&#xff0c;说明两…