《数值计算》.ppt

该【《数值计算》 】是由【相惜】上传分享，文档一共【50】页，该文档可以免费在线阅读，需要了解更多关于【《数值计算》 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。第四章数值计算1精选课件数值计算内容:数值计算包括:多项式运算,线性方程组求解,矩阵特征值问题的解,卷积,数据分析,泛函指令的使用,:f(x)=a0xn+a1xn-1+a2xn-2+…+an-1x+an对于这种表示形式,很容易用一个行向量表示,即: T=[a0,a1,a2,…an-1,an]在MATLAB中,多项式正是用这样一个行向量表示的,,因此,-1直接输入向量构造f(x)=2x5+5x4+4x2+x+4T=[2,5,0,4,1,4];fx=poly2sym(T)函数poly2sym是符号工具箱中的函数,在用此种方式构造多项式时,必须把多项式各项的系数写完整,而不管此项的系数是否为0fx=2*x^5+5*x^4+4*x^2+x+43精选课件r=[1,2,3,4];T1=poly(r);y=poly2sym(T1)y_class=class(y)-2用多项式的根构造多项式,根为r=[1,2,3,4]T1=1-1035-5024y=x^4-10*x^3+35*x^2-50*x+24y_class=,导数运算,估值运算,>多项式的四那么运算多项式四那么运算主要是多项式的加,减,乘,,多项式的加减运算要求两个相加,减的多项式的大小必须相等,此时加,,---+/-乘/除---conv,deconvT=deconv(T1,T3)[T,r]=deconv(T1,T3)商多项式余式T3为分母5精选课件T1=[2,5,0,4,1,4]; T2=[0,0,5,1,3,2]; T3=[5,1,3,2];%除法运算中分母多项式第一个系数不能为0 T=T1+T2; %必须是同维的才能相加 T_add=poly2sym(T) T=T1-T2; T_sub=poly2sym(T) T=conv(T1,T2); %乘法不要求同维 T_mul=poly2sym(T) [A_coe,A_r]=deconv(T1,T3); T_coe=poly2sym(A_coe) T_rem=poly2sym(A_r)-3多项式的加减乘除运算f1(x)=2x5+5x4+4x2+x+4,f2(x)=5x3+x2+3x+-4多项式求值,求上式f1(x)在x==[2,5,0,4,1,4]; x=; y=polyval(T1,x)y=>多项式的导数运算----6求多项式f1(x)=2x5+5x4+4x2+x+4的导数T1=[2,5,0,4,1,4];h=polyder(T1);poly2sym(h)2>多项式求根----5求多项式f1(x)=2x5+5x4+4x2+x+4的根T1=[2,5,0,4,1,4];root=roots(T1);root =-+--+--*x^4+20*x^3+8*x+18精选课件4>拟合和插值---polyfit,-7对以下数据对(x0,y0)=0::1;y0=[-,,,,,,,,,,];p=polyfit(x0,y0,n);[p,s]=polyfit(x0,y0,n);x0,y0---给定数据对n----拟合出的多项式次数p---多项式向量s---偏差信息yi=interp1(x0,y0,xi,'cubic');xi,yi---得到的新的插值点cubic---插值方式9精选课件x0=0::1;y0=[-,,,,,,,,,,];n=3;%设定拟合次数为3[p,s]=polyfit(x0,y0,n);%得到拟合多项式向量和相关偏差信息xx=0::1;yy=polyval(p,xx);%按拟合曲线计算采样值n1=6;%设定拟合次数为6[p1,s1]=polyfit(x0,y0,n1);yy1=polyval(p1,xx);plot(xx,yy,'-b',xx,yy1,'-m',x0,y0,'.r','MarkerSize',20);y3=polyval(p,);y6=polyval(p1,);text(,y3-,'n=3');text(,y6+,'n=6');10精选课件