412符号矩阵的输入 将数值矩阵A转换成符号矩阵,以便求出解 析解 →转换成符号矩阵的语法格式 B=sym(A) 些特殊矩阵,需要编写重载的函数,并置 于目录asym下 e compan, hankel, vander 2/20/2021星期六, 2008-9-6,13:08:30 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:08:30 4.1.2 符号矩阵的输入 将数值矩阵A转换成符号矩阵，以便求出解 析解 转换成符号矩阵的语法格式 一些特殊矩阵，需要编写重载的函数，并置 于目录@sym下 compan, hankel, vander
→构造符号伴随矩阵 function A=compan(c) if min(size(c))>1 error(,Input argument must be a vector. 2) end n =length(c) ifn<=1,A=[] elseif n==2, A=-C(2)/C(1) else, c=c(: A=sym(diag ones(1, n-2), -1)); A(1,:)=-c(2:n)./c(1); end 2/20/2021星期六, 2008-9-6,13:08:30 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:08:30 构造符号伴随矩阵
例46 给定多项式如下 P(入)=a1)9+a23+a3X7+…+a8A2+a9+a10 ↓用解析方法建立其伴随矩阵 MATLAB求解语句 E >>syms a1 a2 a3 a4 a5 a6 a7 a8 a9 a10 A=compan(La1, a2, a3, a4, a5 a6,a7,a8,a9,a10]) 2/20/2021星期六, 2008-9-6,13:08:30 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:08:30 例 4.6 给定多项式如下 用解析方法建立其伴随矩阵 MATLAB求解语句
构造符号 Hankel矩阵 function H=hankel(c,r) c=c(: ) nc=length(c); if nargin==1, r=zeros(size(c));end r=r(: nr=length(r) x=[c;r((2:nr)’)];cidx=(1:nc)); ridx=0:(nr-1) H1=cidx(:, ones (nr, 1))+ridx(ones(nc, 1),:) H=x(H1); 2/20/2021星期六, 2008-9-6,13:08:30 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:08:30 构造符号Hankel矩阵
构造符号 Vandermonde矩阵 function A=vander(v) n=length(v); V=v(:); A=sym (ones (n)) f Or1=n-1: A(:;j)=v.*A(∷,j+1); end 以上三个函数应当存在asm目录下 避免与现有的函数冲突 BACK) 2/20/2021星期六, 2008-9-6,13:08:30 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:08:30 构造符号Vandermonde 矩阵 以上三个函数应当存在@sym目录下 避免与现有的函数冲突