MATLAB-上机实验练习题.docVIP

下载本文档

66
0
约5.54千字
约 15页
2017-09-21 发布于江西
举报
版权申诉

MATLAB-上机实验练习题.doc

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.分别使用poly和polyfit命令，将下列多项式表示为幂级数形式：解：P1： 5 -15 -65 135 180；代码与计算结果： r=[-1,-3,4,3]; p1=poly(r) p1 = 1 -3 -13 27 36 p1=5*poly(r) p1 = 5 -15 -65 135 180 x=[0,1,3,5,7.9]; y=polyval(p1,x) y = 1.0e+003 * 0.1800 0.2400 0 0.4800 9.2693 a=polyfit(x,y,length(x)-1) a = 5.0000 -15.0000 -65.0000 135.0000 180.0000 （2）P2： 4 20 -28 -116 120 0 代码与计算结果： r=[0,1,2,-3,-5]; p2=4*poly(r); x=[0,2,4,6,8,10]; y=polyval(p1,x) y = 0 0 6048 47520 192192 561600 a=polyfit(x,y,length(x)-1) a = 4.0000 20.0000 -28.0000 -116.0000 120.0000 0.0000 2.分别求出经过下述数据点集多项式的幂级数形式： (-1,1), (1,4) (-2,2), (0,-1), (2,1) (-1,-1), (0,0), (1,2), (2,5) 先用解方程的方法求解，再用polyfit命令验证结果. 答案：（1）： x=[-1,1]; y=[1,4]; a=polyfit(x,y,length(x)-1) a = 1.5000 2.5000 （2）： x=[-2,0,2]; y=[2,-1,1]; a=polyfit(x,y,length(x)-1) a = 0.6250 -0.2500 -1.0000 （3）： x=[-1,0,1,2]; y=[-1,0,2,5]; a=polyfit(x,y,length(x)-1) a = 0.0000 0.5000 1.5000 0.0000 3.分别在指定定义域内绘制下列函数图形: 将三个函数曲线话在同一幅图中,图例指出每条曲线代表的函数. 解： x1=0:pi/20:4*pi; y1=sin(x1)./(1+cos(x1)+0/ (0+eps)); Warning: Divide by zero. plot(x1,y1,b); hold on x2=0:0.001:4; y2=1./(1+(x2-2).^2); plot(x2,y2,r) x3=0:0.001:10; y3=exp(-x3).*x3.^2; plot(x3,y3,g) grid on 蓝色代表（1）、红色代表（2）、绿色代表（3）。（1）： x=0:pi/20:4*pi; y=sin(x)./(1+cos(x)+0/ (0+eps)); Warning: Divide by zero. plot(x,y,b) grid on （2）： x=0:0.01:4; y=1./(1+(x-2).^2); plot(x,y,r) grid on 图：（3）： x=0:0.001:10; y=exp(-x).*x.^2; plot(x,y,g) grid on 4.用LU分解求解下列方程：答案： A=[2,-1,0;-1,2,-1;0,-1,2]; [L,A]=lu(A); B=[1;2;3]; x=A\(L\B) x = 2.5000 4.0000 3.5000 5.设，（1）用基于点的二次拉格朗日多项式，求和的近似值.（2）用基于点，，和的三次拉格朗日多项式，求和的近似值. 解：拉格朗日函数定义： function [C,L]=lagran(x,y) w=length(x); n=w-1; L=zeros(w,w); for k=1:n+1 V=1; for j=1:n+1 if k~=j V=conv(V,poly(x(j)))/(x(k)-x(j)); end end L(k,:)=V; end C=y*L （1）： clear,clf x=[1,2,2.5]; y=x+2./x; [C,L]=lagran(x,y); C = 0.4000 -1.2000 3.