原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
一、实验目的掌握 MATLAB 的数值运算及其运算中所用到的函数，掌握结构数组和细胞数组的操作。二、实验内容：(1) 多项式运算。(2) 多项式插值和拟合。(3) 数值微积分。(4) 结构数组和细胞数组。三、实验步骤：1. 多项式运算(1) 多项式表示。在MATLAB 中，多项式表示成向量的形式。如：在MATLAB 中表示为S=[ 1 3 -5 0 9] S=[1 3 -5 0 9]S =1 3 -5 0 9(2) 多项式的加减法相当于向量的加减法，但须注意阶次要相同。如不同，低阶的要补0。如多项式S1=[0 0 2 3 11 ]S2=[1 3 -5 4 7 ]S3=S1+S2 S1=[0 0 2 3 11]S1 = 0 0 2 3 1S2 = 1 3 -5 4 7 S3=S1+S2S3 = 1 3 -3 7 18(3) 多项式的乘、除法分别用函数conv 和deconv 实现S1=[ 2 3 11 ]S2=[1 3 -5 4 7 ]S3=conv(S1,S2)S4=deconv(S3,S1) S1=[2 3 11]S1 =2 3 11 S2=[1 3 -5 4 7]S2 =1 3 -5 4 7 S3=conv(S1,S2)S3 =2 9 10 26 -29 65 77 S4=deconv(S3,S1)S4 = 1 3 -5 4 7(4) 多项式求根用函数roots S1=[ 2 4 2 ] roots(S1)S1 = 2 4 2 roots(S1)ans = -1 -1(5) 多项式求值用函数polyval S1=[ 2 4 1 -3 ] S1 = 2 4 1 -3 polyval(S1,3) %计算x＝3 时多项式的值 x＝1:10 polyval(S1,3)ans = 90 x=1:10x = 1 2 3 4 5 6 7 8 9 10 y=polyval(S1,x) %计算x 向量对应的值得到y 向量 polyval(S1,x)ans = Columns 1 through 5 4 31 90 193 352 Columns 6 through 10 579 886 1285 1788 24072. 多项式插值和拟合有一组实验数据如附表1-1 所示。请分别用拟合(二阶至三阶)和插值(线性和三次样条)的方法来估测 X=9.5 时Y 的值。以下是实现一阶拟合的语句。x=1:10 x = 1 2 3 4 5 6 7 8 9 10y=[16 32 70 142 260 436 682 1010 1432 1960] y =Columns 1 through 5 16 32 70 142 260 Columns 6 through 10 436 682 1010 1432 1960p1=polyfit(x,y,1) %一阶拟合 pl = 204.8000 -522.4000y1=polyval(p1,9.5) %计算多项式p1 在x＝9.5 时的值 y1 = 1.3208e+0033. 数值微积分(1) 差分使用diff 函数实现。x=1:2:9 x = 1 3 5 7 9diff(x) ans = 2 2 2 2 (2) 可以用因变量和自变量差分的结果相除得到数值微分。x=linspace(0,2*pi,100);y=sin(x);plot(x,y)y1=diff(y)./diff(x);plot(x(1:en
文档评论（0）