matlab文件数学建模作业.docVIP

下载本文档

4
0
约1.63万字
约 10页
2017-02-24 发布于浙江
举报
版权申诉

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文档。上传文档

1.1 代码：syms x;m=433;limit((m*x-sin(m*x))/(x^3),x,0) 结果：ans6 代码:syms x;m=433;limit((m*x-sin(m*x))/(x^3),x,Inf) 结果：ans=0 1.2 代码： syms x;m=433; y=exp(x)*cos(m*x/1000); diff(y,x,2) 结果：(812511*cos((433*x)/1000)*exp(x))/1000000 - (433*sin((433*x)/1000)*exp(x))/500 1.3 代码： dblquad(exp(x.*x+y.*y),0,1,0,1) 结果： ans = 2.1394 1.4 代码： syms x;m=433; int(x^4/(m^2+4*x*x),x) 结果： atan((2*x)/433))/32 - (187489*x)/16 + x^3/12 1.5 代码： syms x;m=433; diff(exp(x)*cos(m*x),10) 结果： - 231619285105554448813874688*cos(433*x)*exp(x) - 5350116883610368850333216*sin(433*x)*exp(x) 1.6 代码： syms x;m=433; simple(taylor(sqrt(m/1000.0+x),x,0,Order,5)) 结果： 4330^(1/2)*(- (390625000*x^4)/35152125121 + (625000*x^3)- (1250*x^2)/187489 + (5*x)/433 + 1/100) 1.7 代码： n=2;A=[1,1]; while n20;A(n+1)=A(n)+A(n-1);n=n+1;end; fprintf(A=(); for i=1:size(A,2);fprintf(%u,A(i));if(i~=size(A,2)) fprintf(,);end;end; fprintf()); 结果： A=(1,1,2,3,5,8,13,21,34,55,89,144,233,377,610,987,1597,2584,4181,6765) 1.8 代码： m=433;A=[-2,1,1;0,2,0;-4,1,m/1000]; inv(A) %逆矩阵 [P,D]=eig(A) %求P,D和特征值特征向量 det(A) %求行列式 A^6 %求A6 结果：逆矩阵为： P： 0.1382 0.0905 -0.3191 0.2720 - 0.3550i 0.2720 + 0.3550i 0.2425 0 0.5000 0 0 0 0.9701 1.2763 -0.3191 -0.6382 0.8944 0.8944 0 D： -0.7835 + 1.5875i 0 0 0 -0.7835 - 1.5875i 0 0 0 2.0000 行列式：6.2680 A6： 37.4463 6.6384 -7.3857 0 64.0000 0 29.5429 -7.3857 19.4769 特征向量：P的列向量特征值：与P的列向量同列的D的对角线的值为对应的特征值 1.9 -------------------------------------------- f.m -------------------------------------------- function y=f(x) if(x=0 x=0.5);y=2*x; elseif(x0.5 x=1);y=2*(1-x); else;y=NaN;end; -------------------------------------------- fplot(@f,[0,1]); 1.1