原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
1. 表1 用三次样条方法插值计算0-90 度内整数点的sin 值和0-75 度内整数点的
正切值，然后用5 次多项式拟合方法计算相同的函数。
a（度） 0 15 30 45 60 75 90
Sin（a） 0 0.2588 0.5000 0.7071 0.8660 0.9659 1.0000
tan（a） 0 0.2679 0.5774 1.0000 1.7320 3.732
解：分别对应的程序如下：
Sin x函数和tan x函数x1= pi*(0:90)/180;
y1 = sin(x1);
xx1 = pi*(0:.25:90)/180;
x2 = pi*(0:75)/180;
y2 = tan(x2);
xx2 = pi*(0:.25:75)/180;
yy2 = spline(x2,y2,xx2);
yy1 = spline(x1,y1,xx1);
plot(x1,y1,mo,xx1,yy1,x2,y2,ro,xx2,yy2)；
Sin x与figure
x1=pi*(0:15:90)/180;
y1=[0,0.2588,0.5,0.7071,0.866,0.9659,1.0];
x2=pi*(0:15:75)/180;
y2=[0,0.2679,0.5774,1,1.732,3.732];
xx1=pi*(1:90)/180;
pp1=polyfit(x1,y1,5);
yy1=polyval(pp1,xx1);
xx2=pi*(1:75)/180;
pp2=polyfit(x2,y2,5);
yy2=polyval(pp2,xx2);
plot(x1,y1,-mo,xx1,yy1,x2,y2,-ro,xx2,yy2);
legend(描点显示,五次拟合)
2. 采用最近点法、线性法和3 次样条法插值计算1-100 整数间平方根
n 1 4 9 16 25 36 49 64 81 100
Sqtr（n）1 2 3 4 5 6 7 8 9 10
解：程序如下：
x=[1,4,9,16,25,36,49,64,81,100];y=x.^0.5;
xx=0:100;
subplot(1,3,1);
y1=interp1(x,y,xx,nearest);
plot(x,y,-*,xx,y2,xr);
title(最近点法)
subplot(1,3,2);
y2=interp1(x,y,xx);
plot(x,y,-*,xx,yy,rx);
title(线性法);
subplot(1,3,3);
y3=interp1(x,y,xx,spline);
plot(x,y,-*,xx,y3,xr);
title(3次样条法);
仿真的结果：
3. 已知p(x)=2x^4-3x^3+5x+13，求p（x）的全部根，由方程p（x）=0 的根构造
一个多项式f（x），并和p（x）比较。
解：
程序如下：
P=[2,-3,0,5,13];
K=roots(P);
P=poly(K);
F=poly2str(P,x);
4. 有两个多项式p(x)=2x^4-3x^3+5x+13;q(x)=x^2+5x+8, 求p(x),p(x)q(x) 和
p(x)/q(x)的导数。
解：对应的程序如下：
P=[2,-3,0,5,13];
Q=[1,5,8];
format rat
K1=polyder(P);
K2=polyder(P,Q);
[p,q]=polyder(P,Q);
p=poly2str(p,x);
q=poly2str(q,x);
D1=poly2str(K1,x);
D2=poly2str(K2,x);
D3=poly2str([p,q],x)
运行的结果如下：
D1 =
8 x^3 - 9 x^2 + 5
D2 =
12 x^5 + 35 x^4 + 4 x^3 - 57 x^2 + 76 x + 105
D3 =
x^82 + x^81 + x^80 + 4 x^79 + x^78 + x x^77 + ^ x^76 + 5 x^75 + x^74 + + x^73 + x^72 + 2 x^71 + 7 x^70 + x^69 + x x^68 + ^ x^67 + 4 x^66 + x^65 + + x^64 + x^63 + 3 x^62 + 4 x^61 + x^60 + x x^59 + ^ x^58 + 3 x^57 + x^56 + - x^55 + x^54+
文档评论（0）