第3章 线性代数应用.pptVIP

线性微分方程的解已知传递函数，将其分解为部分分式Num=[1,2];Den=[143];[r,p,k]=residue(Num,Den)r=0.50000.5000p=-3-1k=[]第29页，共47页，星期日，2025年，2月5日求解线性常微分方程在输入u（t）为单位脉冲和单位阶跃信号时的解析解Num1=[30.54];Den1=[1547];[r1,p1,k1]=residue(Num1,Den1)t=0:0.2:10;yi=r1(1)*exp(p1(1)*t)+r1(2)*exp(p1(2)*t)…+r1(3)*exp(p1(3)*t);subplot(1,2,1)plot(t,yi)第30页，共47页，星期日，2025年，2月5日Den2=[15470];[r,p,k]=residue(Num1,Den2)ys=r(1)*exp2(p(1)*t)+r(2)*exp(p(2)*t)…+r(3)*exp(p(3)*t)+r(4)*exp(p(4)*t);subplot(1,2,2)plot(t,ys)r1=3.2288-0.1144+0.0730i-0.1144-0.0730ip1=-4.4548-0.2726+1.2235i-0.2726-1.2235ik1=[]r=-0.72480.0767+0.0764i0.0767-0.0764i0.5714p=-4.4548-0.2726+1.2235i-0.2726-1.2235i0k=[]第31页，共47页，星期日，2025年，2月5日第32页，共47页，星期日，2025年，2月5日多项式拟合polyfit(x,y,n)x,y是已知的N个数据点坐标向量，n是拟合多项式的次数用于在最小均方误差准则下对数据进行n阶多项式曲线拟合。其中x，y是已知的N个数据点坐标向量，其长度均为N，现欲用一个n阶多项式来表示x，y之间的函数关系，待定的是该n阶多项式的系数。多项式阶数n越大，拟合程度越好，设p为用于拟合的多项式系数个数，p=n+1。第33页，共47页，星期日，2025年，2月5日第1页，共47页，星期日，2025年，2月5日数据分析和傅立叶变换基本数据分析场论数据分析随机数据分析相关分析和傅立叶分析3.1数据分析(datafun)第2页，共47页，星期日，2025年，2月5日基本数据分析函数Max(data)Min(data)mean(data)sum(data)第3页，共47页，星期日，2025年，2月5日场论的数据分析函数Gradient()二维场和三维场的近似梯度参照helpgradient。Cross()向量的矢量积Dot()向量的数量积第4页，共47页，星期日，2025年，2月5日随机过程的数据分析函数rand(m,n)均匀分布随机矩阵参照rand.m文件，熟悉rand的各种用法。randn(m,n)正态分布随机矩阵第5页，共47页，星期日，2025年，2月5日相关分析和傅立叶分析函数corrcoef两个同长信号的相关系数con协方差矩阵conv向量的卷积，参照conv.m。fft傅立叶变换，参照fft.m。ifft傅立叶反变换sound声音输出函数第6页，共47页，星期日，2025年，2月5日3.2矩阵分析和线性代数(matfun)线性方程组的系数矩阵矩阵的分解矩阵的特征值特殊矩阵库第7页，共47页，星期日，2025年，2月5日1.逆、秩与行列式计算求逆：inv(A)求行列式：det(A)求秩：rank（A）要求矩阵必须为方阵A\B=inv(A)*BA/B=A*inv(B)第8页，共47页，星期日，2025年，2月5日【例3.1】a=[123;456;235];b=inv(a)%求方阵a的逆c=det(a)%求方阵a的行列式d=rank(a)%求a中线性无关的行数第9页，共47页，星期日，2025年，2月5日b=-2.33330.33331.00002.66670.3333-2.0000-0.6667