实验八MATLAB在数学中的应用.docVIP

实验八 MATLAB在数学中的应用 （范数特征值，多项式，矩阵统计及函数分析） 1、实验目的 1）掌握矩阵的分析方法 2）掌握多项式的运算规则 3）掌握数据的分析与统计 4）掌握函数分析与数值积分 注意：每个题目都需要创建独立的m文件，每个完成后命名为“e学号_题目编号”，最后打包成“学号_实验八”,发送到邮箱“ HYPERLINK mailto:2085565215@ 2085565215 HYPERLINK mailto:2085565215@ @”。邮件标题“班级+姓名+实验八”，附件的形式上传“学号_实验八zip”。 2、实验内容 2.1 矩阵分析 2.1.1 范数 a=[1, 2, 3; 4, 3, 1; 5, 3, 8] n=norm(a) % n = 11.2269 n=norm(a,2) % n = 11.2269 b=[1 2 4 3 5] % 1 2 4 3 5 m=norm(b,2) % m = 7.4162 2.1.2 矩阵求逆inv(A)，行列式值：det(A) a=[1,2,3;4,3,1;5,3,8] inv(a) det(a) a*inv(a) 2.1.3 线性方程组的解 例1. m=n det(A)≠0 A = [6 3 4;-2 5 7;8 -4 -3] B = [3 -4 -7] X = A\B det(A) 例2. mn, 超定方程组 A(4,:) = [1 5 -7] B(4,:) = [9] X = A\B 例3. mn, 解不唯一 A=[2,4,2,1;-1,2,0,2;3,5,2,1] B=[1;4;6] X = A\B 2.1.4 矩阵分解-lu分解 A=[2,4,1; -1,2,2;3,5,1] [L,U,P]=lu(A) [L,U]=lu(A) 试比较结果的不同及原因 2.1.5 矩阵特征值 A=[0,-6,-1;6,2,-16;-5,20,-10] lanmuda=eig(A) [V,D]=eig(A) 2.1.6 矩阵的幂次运算 A=[1 1 1;1 2 3;2 3 5] X=A^3 Y=A^(-3) Z=A^(2/3) 2^A exp(A) 2.2 多项式运算 2.2.1 多项式表示及其四则运算 a=[1 2 3 4] b=[2 3] b=[0 0 b] c=a+b a=[1 2]; b=[1 3]; c=conv(a,b) [q,r]=deconv(c,a) 2.2.2 多项式求导函数polyder 例1 a=[1,4,3], b=[1,3] c=polyder(a) d=polyder(b) 例2 a=[1,1] b=[1,2] c=polyder(a,b) conv(a,b) 2.2.3 多项式求根 p =[ 1 -12 0 25 116] r=roots(p) 由根求多项式：poly r = [-1,-2]’ poly(r) 2.2.4 多项式求值 例1：利用m文件运行下列程序 % 求H(jw)并画图 figure(1) a1=[2,4,6,8];b1=[3,6,9]; w=linspace(0,10); A=polyval(a1,j*w); B=polyval(b1,j*w); subplot(2,1,1), plot(w,abs(B./A)) subplot(2,1,2) plot(w,angle(B./A)) 例2：利用m文件运行下列程序 figure(2) w1=logspace(-1,1); H=polyval(b1,j*w1)./polyval(a1,j*w1); subplot(2,1,1) semilogx(w1,abs(H)) subplot(2,1,2) semilogx(w1,angle(H)) 例3：求的时域响应并画图。 num=[3,0.5,4]; den=[1,5,4,7]; [r,p,k]=residue(num,den) clf t=[0:0.2:10]; y=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(3)*t); plot(t,y) 2.2.5 多项式拟合 线性拟和 x=[0:10] y=exp(-0.1*x)+0.5 p=polyfit(x,y,1) z=p(1)*x+p(2); plot(x,y,x,z) 二次拟和 x=[0:10] y=exp(-0.1*x)+0.5 p=polyfit(x,y,3) z1=