matlab文件 数学建模作业.docVIP

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