线性的代数应用实例.docVIP

实用标准文案 精彩文档 线性代数应用实例 求插值多项式 右表给出函数上4个点的值，试求三次插值多项式，并求的近似值。 ti 0 1 2 3 f(ti) 3 0 -1 6 解：令三次多项式函数过表中已知的4点，可以得到四元线性方程组： 对于四元方程组，笔算就很费事了。应该用计算机求解了，键入： A=[1,0,0,0;1,1,1,1;1,2,4,8;1,3,9,27], b=[3;0;-1;6], s=rref([A,b]) 得到x = 1 0 0 0 3 0 1 0 0 -2 0 0 1 0 -2 0 0 0 1 1 得到，三次多项函数为，故近似等于。 在一般情况下，当给出函数在n+1个点上的值时，就可以用n次多项式对进行插值。 在数字信号处理中的应用----- 数字滤波器系统函数 u2x1y1/4? u 2 x1 y 1/4 ?1/4 z?1 x3 x2 z?1 3/8 图1 某数字滤波器结构图 图1表示了某个数字滤波器的结构图，现在要求出它的系统函数，即输出y与输入u之比。先在它的三个中间节点上标注信号的名称x1,x2,x3，以便对每个节点列写方程。由于迟延算子z ?1不是数，要用符号代替，所以取q? z ?1，按照图示情况，可以写出： 写成矩阵形式为 经过移项后，系统函数W可以写成： 现在可以列写计算系统函数的MATLAB程序ea705， syms q % 规定符号变量 Q(1,2)?q; Q(2,3)=3/8*q?1/4; Q(3,1)=1; % 给非零元素赋值 Q(3,3)=0; % 给右下角元素Q（3,3）赋值后，矩阵中未赋值元素都自动置零 P=[2;1/4;0] % 给P赋值 W=inv(eye(3)?Q)*P % 用信号流图求传递函数的公式 程序运行的结果为 W = [?16/(?8?3*q^2?2*q)?2*q/(?8?3*q^2?2*q) ] [ ?2*(3*q?2)/(?8?3*q^2?2*q)?2/(?8?3*q^2?2*q)] [?16/(?8?3*q^2?2*q)?2*q/(?8?3*q^2?2*q)] 我们关心的是以y?x3作为输出的系统函数，故再键入 pretty(W(3)) 整理后得到 用线性代数方法的好处是适用于任何复杂系统，并能用计算机解决问题。 信号与系统课程中的应用-----线性时不变系统的零输入响应 描述n阶线性时不变（LTI）连续系统的微分方程为 n≥m 已知y及其各阶导数的初始值为y(0)，y(1)(0)，…，y(n-1)(0)，求系统的零输入响应。 解：当LTI系统的输入为零时，其零输入响应为微分方程的齐次解（即令微分方程等号右端为0），其形式为（设特征根均为单根） 其中p1，p2，…，pn是特征方程a1?n+a2?n-1+…+ an?+ an+1 =0的根，它们可用roots(a)语句求得。各系数C1，…，Cn由y及其各阶导数的初始值来确定。对此有 C1+ C2+…+Cn = y0 y0 = y(0) p1C1+ p2C2+…+ pnCn=Dy0 (Dy0表示y的导数的初始值y ………………………………… 写成矩阵形式为 即 V·C = Y0 ， 其解为 C =V \ Y0 式中 V为范德蒙矩阵，在MATLAB的特殊矩阵库中有vander函数可直接生成。 MATLAB程序ea703.m a=input(输入分母系数向量a=[a1，a2，...]= ); n=length(a)-1; Y0=input(输入初始条件向量 Y0=[y0，Dy0，D2y0，...]= ); p=roots(a);V=rot90(vander(p));c= V\Y0; dt=input(dt=); tf=input(tf= ) 图2 三阶系统的零输入响应t=0:dt:tf; y=zeros(1，length(t)); 图2 三阶系统的零输入响应 for k=1:n y= y+c(k)*exp(p(k)*t);end plot(t，y)，grid ? 程序运行结果 用这个通用程序来解一个三阶系统，运行此程序并输入 a=[3，5，7，1]; dt=0.2; tf=8; 而Y0取 [1，0，0];[0，1，0];[0，0，1] 三种情况，用hold on语句使三次运行生成的图形画在一幅图上，得到图2。 减肥配方的实现 设三种食物每100克中蛋白质、碳水化合物和脂肪的含量如下