标准实用文案大全线性代数应用实例求插值多项式右表给出函数()ft上4个点的值,试求三次插值多项式230123()ptaatatat,并求(1.5)f的近似值。解:令三次多项式函数230123()ptaatatat过表中已知的4点,可以得到四元线性方程组:627931842033210321032100aaaaaaaaaaaaa对于四元方程组,笔算就很费事了。应该用计算机求解了,键入:>>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=100030100-20010-200011得到01233,2,2,1aaaa,三次多项函数为23()322ptttt,故(1.5)f近似等于23(1.5)32(1.5)2(1.5)(1.5)1.125p。在一般情况下,当给出函数()ft在n+1个点(1,2,,1)itin上的值()ift时,就可以用n次多项式2012()nnptaatatat对()ft进行插值。在数字信号处理中的应用-----数字滤波器系统函数数字滤波器的网络结构图实际上也是一种信号流图。它的特点在于所有的相加节点都限定为双输入相加器;另外,数字滤波器器件有一个迟延一个节拍的运算,它也是一个线性算子,它的标注符号为z1。根据这样的结构图,也可以用类似于例7.4的方法,求它的输入输出之间的传递函数,在数字信号处理中称为系统函数。图1表示了某个数字滤波器的结构图,现在要求出它的系统函数,即输出y与输入u之比。先在它的三个中间节点上标注信号的名称x1,x2,x3,以便对每个节点列写方程。ti0123f(ti)30-16图1某数字滤波器结构图u2x1y1/41/4z1x3x2z13/8标准实用文案大全由于迟延算子z1不是数,要用符号代替,所以取qz1,按照图示情况,可以写出:1223312311844xqxuxqxuxx写成矩阵形式为112233002311008440100qxxxqxxxxux=Qx-Pu经过移项后,系统函数W可以写成:W=x/u=inv(I-Q)*P现在可以列写计算系统函数的MATLAB程序ea705,symsq%规定符号变量Q(1,2)q;Q(2,3)=3/8*q1/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/(83*q^22*q)2*q/(83*q^22*q)][2*(3*q2)/(83*q^22*q)2/(83*q^22*q)][16/(83*q^22*q)2*q/(83*q^22*q)]我们关心的是以yx3作为输出的系统函数,故再键入pretty(W(3))整理后得到1222116288(3)8321.541.54yqqzWuqqqqzz用线性代数方法的好处是适用于任何复杂系统,并能用计算机解决问题。信号与系统课程中的应用-----线性时不变系统的零输入响应描述n阶线性时不变(LTI)连续系统的微分方程为,dddddddddd111121ubtubtubyatyatyatyammmmnnnnnn≥m已知y及其各阶导数的初始值为y(0),y(1)(0),⋯,y(n-1)(0),求系统的零输入响应。解:当LTI系统的输入为零时,其零输入响应为微分方程的齐次解(即令微分方程等号右端为0),其形式为(设特征根均为单根)tpntptpnCCCtyeee)(2121其中p1,p2,⋯,pn是特征方程a1n+a2n-1+⋯+an+an+1=0的根,它们可用roots(a)语句求得。各系数C1,⋯,Cn由y及其各阶导数的初始值来确定。对此有C1+C2+⋯+Cn=y0y0=y(0)p1C1+p2C2+⋯+pnCn=Dy0(Dy0表示y的导数的初始值y(1)(0))⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯标准实用文案大全011212111DyCpCpCpnnnnnn写成矩阵形式为0100211121121DD111yyyCCCppppppnnnnnnn即V·C=Y0,其解为C=V\Y0式中112000[,,,];[,D,,D]nnCCCyyyTT0CY1121121111nnnnnppppppVV为范德蒙矩阵,在MATLAB的特殊矩阵库中有vander函数可直接生成。MATLAB程序ea703.ma=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=')t=0:dt:tf;y=zeros(1,length(t));fork=1:ny=y+c(k)*exp(p(k)*t);endplot(t,y),grid程序运行结果用这个通用程序来解一个三阶系统,运行此程序并输入a=[3,5,7,1];dt=0.2;tf=8;而Y0取[1,0,0];[0,1,0];[0,0,1]三种情况,用holdon语句使三次运行生成的图形画在一幅图上,得到图2。减肥配方的实现设三种食物每100克中蛋白质、碳水化合物和脂肪的含量如下表,表中还给出了80年代美国流行的剑桥大学医学院的简捷营养处方。现在的问题是:如果用这三种食物作为每天的主要...