(1)**：syms x;m=433;limit((m*x-sin(m*x))/xa3),^0) 结果：ans=81182737/6
2 )**：syms x;m=433;limit((m*x-sin(m*x))/xa3)/x/inf) 结果：ans=0
**：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
**：dblquad(,exp(x.*x+y.*y)'/0/l/0/l)
结果：ans = 2.1394
**：syms x;m=433;
i nt(xa4/(ma 2+4*x*x),x)
结東：81182737*atan((2*x)/433))/32 - 187489*x)/16 + xa3/12
**：syms x;m=433;
diff(exp(x)*cos(m*x)/10)
结果：c231619285105554448813874688*cos(433*x)*exp(x) -5350116883610368850333216*sin(433*x)*exp(x)
**：syms x;m=433;
******(taylor(sqrt(m/1000.0+x)/x/0/,order,/5))
结果：4330a(1/2)*(390625000*xa4)/35152125121 + 625000*xa3)/81182737 (1250*xa2)/187489 + 5*x)/433 + 1/100)
**：n二 2;a二g];
while n<20;a(n+1)=a(n)+a(n-l);n=n+l；end;
fprintf(,a=(,
for i = l:size(a/2);fprintf( %u,/a(i));if(i~=size(a/2)) fprintf(7);end;end;
fprintfc)1);
结果：zv=(1j,2,3,5,8,13,21,34,55,89,144,233377,610,987,1597,2584,4181,6765)
function y二f(x)
if(x>=0 &&x<=0.5);y=2*x; elseif(x>0.5 &&x< =l);y=2*(l-x); else;y=nan;e nd;
**：t=0:0.01*pi:6*pi;
plot3(cos(t)zsinm/jhold on;grid on;
plot3(2*cos(t)/2*sin(t)/t/k,);
**：m=433;a=[4r2/2;-3a5;l/5*m/3]； b=[l/3a-2/0/-3;2rl/l];
det(a)/2*a-b/a*b/a.*b/a*ba(-l)/ inv(a)*b/aa2/a,结果：
x=-3:0.1:3;m=433;a = m/600;
f=@(x,a,b) l/sqrt(2*pi)/a*exp(-(x-b).a2./(2*a*a));plot(x，f(xao;hold on;
plot&,f&a-i));hold on;
plot(x,f(xai));
x 二二 433;
f=@(x/a/b) l/sqrt(2*pi)/a*exp(-(x-b).a2./(2*a*a));plot(x,f(x,l,0));hold on; plot(x,f&20;hold on; plot(x,f(xao));hold on;
p i ot(x,f (x, m/100r 0));
m=433;ezmesh('433*x.*x+ya4,,30); coiorm 叩([oqlj);
绘图。m=433;x=-2:0.01:2;
y= 5-m/200*x-0.1;
plot(x,y);grid on;
求根。vpa(solve(y),3)
得根为：1.22,-0.0462,-1.2
求单调区间：
vpa(solve(diff(y))z3)
得：0.811,-0.811
根据图象：单调增区间为(8 , 0.811)4 0.811, +8 )单调减区间为：(-0.811a811)
推理：当x<-0.811是"的导数大于0 ,故递增；-0.
8110.811时，y的导数大于0 ,故递增。x=-0.
811和0.811为函数的极大点和极小点，极小值小于0 ,极大值大于0 , x趋近于正无穷时，y趋近于正无穷，x趋近于负无穷时，y趋近于负无穷，根据零点定理，只有所求的3个解。
体会：利用计算机可以方便的分析函数。
画图：m二433;f二@(x) exp(x)-3*m*x./k2;
x=-10：0.1：14;plot(xzf(x));grid on;
求解：fsolve(f,[-5,20])得到解为：
**：m=433;f=@(x) (x+m/x)/2;
x=3;for k=l:12;x=f(x);fprintf('%8.8g ,,x);end;fprintf(',);结果：
所以｛x』收敛，收敛于20.808652
1)求不动点。
m二433;syms x;vpa(solve(x*x+m*xx+l),3)
得不动点为-0.0231和432
函数的收敛点。
function y=fun(x)
m=433;f=@(t) (t-l)./t+m);lastx=x;
for k=l：100;lastx=x;x=f(x);if(x==lastx);y=vpa(xz3);return;end;end; y=nan;
fun(0)=-0.00231 fun(10)=-0.00231
fun(-10) =0.00231 fun(-1000) =0.00231 fun(-432) =nan fun(-433) =nan
2)求不动点。
m二433;syms x;vpa(solve(x*x+m*x-x
得不动点为268.0和700
函数的收敛点。
function y=fun(x)
m=433;f=@(t) (t+m*m)./t+m);lastx=x;
for k= 1： 100;lastx=x;x=f(x);if(abs(x==lastx)<0.000001);y=vpa(x/3)；return;end;end;
y=nan;
fun(0)=268.0 fun(268.0)=268.0 fun(1000)=268.0
fun(-1000)=268.0 fun(-700)=268.0 fun(-699.8853169977664947509765625) =700
打印对应初值的x的迭代序列。
function fun(x)
f=@(m) (m>二0 & m<=0.5)*2*m+(m>0.5 & m<=l)*2*(l-m));for k=l：
100;x=f(x);fprintf('%7.4g'zx);
讦(mod(k/10)==0);f printf ('n' )end;end;
fun(pi/5):
绘制蛛网图：
fun function fun（xoq） f=@（x） a*x*（l-x）;
x=y=[j；
x ⑴=xo；x（2）=xo；y（i）=o；y（2）=f（xo）；
for k二 1:50
x(k*2+l)=y(k*2);x(k*2+2)=x(k*2+l);
y(k*2--l)=y(k*2)；y(k*2+2)=f(y(k*2));
endplot(x//)hold on;grid on;
数学建模matlab作业
数统应数 20121323003 王妍。数统应数 20121323022 胡可旺。1 生成5阶矩阵，使其元素满足均值为1，方差为4的正态分布 y 1 sqrt 4 randn 5 结果 2，生成一个20行5列矩阵a，其元素分别以概率 2取自 这三个整数，然后将 删除a的全零行 之后的结果赋给矩阵b ...
数学建模培训作业 MATLAB编程部分
数学建模培训作业。matlab编程部分 1 请使用switch语句将百分制的学生成绩转换为五级制的成绩输出。2 猜数游戏 首先由计算机随机产生一个 1,100 之间的一个整数，然后由用户猜测所产生的这个数。根据用户猜测的情况给出不同的提示，如果猜测的数大于产生的数，则显示 high 小于则显示 lo...
数学建模Matlab练习
matlab练习四。1.绘制两个一元函数的曲线图 例如y sin x y cos x 加上标注说明，包括坐标 抬头 曲线说明等。x 0 pi 30 2 pi y sin x plot x,y,ob markersize 10 hold on y1 cos x plot x,y1,markersize...