matlab数学实验复习题有答案Word文件下载.docx

1、axis equal纵、横坐标轴采用等长刻度text（1，2，y=sin（x）在x=1,y=2处加上字符串y=sin（x）；hold on把新的plot产生的图形画在原来的图形上。title（y=sin（x）在图形正上方加上字符串y=sin（x）16、Matlab中自定义函数M文件的第一行必须以function开头；17、二种数值积分的库函数名为：quad;quadl18、unifrnd（1，2，3，4）的功能是：随机生成3行4列均匀分布，每个元素服从（1，2）的矩阵19、binornd（20，0.3，3，4）的功能是随机生成3行4列服从（20，0.3）的二项分布的矩阵20、eig（A）的功能

2、是矩阵A的特征值21、设x是一向量，则hist（x）的功能是作出将X十等分的直方图22、interp1（1,2,3,3,4,5,2.5）Ans=4.523、建立一阶微分方程组的函数M文件。（做不出来）二、写出运行结果：1、eye（3，4）=2、size（1,2,3）=1；33、设b=round（unifrnd（-5，5，1，4），则=3 5 2 -5x,m=min（b）；x=-5；m=4,x,n=sort（b）-5 2 3 54 3 1 2mean（b）=1.25，median（b）=2.5，range（b）=104、向量b如上题，则any（b）,all（b2）,all（b5 6;7 87 8

3、;5 6=6、若，则7、diag（diag（B）=8、4:-2:1.*-1,6=-4 129、acos（0.5）,atan（1）ans=.0471*0.78539816339744810、norm（1,2,3）Ans=3.74165738677394111、length（1，3，-1）=312、x=0:0.4:2;plot（x,2*x,k*）13、zeros（3,1）;14、ones（3）=，vander（2,3,5）=16、floor（1：0.3：3）=1 1 1 1 2 2 218、subplot（2,2,1）; fplot（sin,0,2*pi）;subplot（2,2,2）;plot（

4、1,2,-1）;x=linspace（0,6*pi）;subplot（2,2,3）;plot3（cos（x）,sin（x）,x）;subplot（2,2,4）;polar（x,5*sin（4*x/3）;19、t=linespace（0，2，11）0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.020、a,b=binostat（15,0.2）a=3 b=2.4y1=binopdf（5,10,0.7）=0.1029,y2=binocdf（5,10,0.7）=0.150321、log10（1,10,100）=0 1 222、p=1;for k=2:3:9 p=p*

5、k;end;p p=8023、s=0;9 s=s+k;s s=1524、Ans=3.864725、a1=norminv（0.6,3,4）a1=4.013426、unifinv（0.4,1,5）,unifpdf（0.4,1,5）,unifpdf（2,1,5）Ans=2.6 0 0.2527、A=0 1-1;2 1 0;1-1 1;0 1 -11 -1 1A（1,3,:）A（3,1,:）=1 -1 1 0 1 -1A（2,:）=2 1 0-2*A（1,:）= 0 -2 228、quad（sin（x）,0,pi/2）=1.000029、trapz（3,4,6,1,2,3）=6.500030、 int

6、（x-sin（x）,0,1）Ans=cos（1） - 1/231、round（3:5）,ceil（3:5）;floor（3:5）3 3 4 4 5 53 3 3 4 4 5limit（1+1/（3*x）x,inf）=1diff（sin（3*x）+x3,2）=6*x-9*sin（3*x） taylor（exp（3*x）,5,1）:命令输入： y=taylor（exp（3*x）,x,1,Order,5） exp（3） + 3*exp（3）*（x - 1） + （9*exp（3）*（x - 1）2）/2 + （9*exp（3）*（x - 1）3）/2 + （27*exp（3）*（x - 1）4）/8

7、a1=mod（15,4）,b1=rem（15,4）=3,3a2=mod（-15,-4）,b2=rem（-15,-4）=-3,-3a3=mod（15,-4）,b3=rem（15,-4）=-1,-3a4=mod（-15,4）,b4=rem（-15,4）=1,-334、x=binornd（20,0.4,2,4）8 7 10 810 7 9 12sign（x）,1 1 1 1y=-poissrnd（8,2,4）-16 -10 8 -7-7 -8 -6 -9sign（y）-1 -1 -1 -135、a1,b1=binostat（20,0.4） a1=8 b1=4.8a2,b2=poisstat（8）an

8、s=8,8a3,b3=chi2stat（15）ans=15 3036、运行M文件：chi2fign=5;a=0.9;xa=chi2inv（a,n）;0.1:15;y=chi2pdf（x,n）;plot（x,y,b）;hold on;xf=0:xa;yf=chi2pdf（xf,n）;fill（xf,xa,yf,0,gtext（xa*1.01,0.005,num2str（xa）;text（2.5,0.05,alpha=0.9,fontsize,20）;text（9,0.09,Xchi2（4）,16）;37、t=linspace（0,2*pi）;polar（t,3*t,g*）38、quadl（exp

9、（2*x）.*log（3*x）,1,3）ans = 398.635239、x0=0:2*pi/6:2*pi;y0=sin（x0）.*cos（x0）;x=linspace（0,2*pi,100）;y=sin（x）.*cos（x）;y1=spline（x0,y0,x）;x;y;y1k,x,y1,b-注：此处省略100组数据40、A=round（unifrnd（0,100,3,3）;L,U=lu（A）L = 0.9897 0.4699 1.0000 0.1649 1.0000 0 1.0000 0 0U = 97.0000 80.0000 92.0000 0 35.8041 26.8247 0 0

10、-89.656841、a=sparse（1 3 3,2 3 5,1 2 3,4,5）;s=full（a）s = 0 1 0 0 0 0 0 0 0 0 0 0 2 0 3三、编程1、分别用矩形公式、梯形公式、辛普森公式、Gauss-Lobatto公式及随机模拟方法计算数值积分，并与符号运算计算的结果进行比较。format long0.01:pi/2;y=exp（3*x）.*sin（2*x）;s1=sum（y）*0.01;s2=trapz（x,y）;s3=quad（exp（3*x）.*sin（2*x）,0,pi/2）;s4=quadl（n=10000;x=unifrnd（0,pi/2,1,n）;

11、y=unifrnd（0,exp（5.5）,1,n）;k=0;for i=1:n if y（i）=exp（3*x（i）.*sin（2*x（i） k=k+1; endends5=k/n*pi/2*exp（5.5）;syms xs=int（exp（3*x）.*sin（2*x）,0,pi/2）;s6=double（s）;s1,s2,s3,s4,s5,s6 输出结果： Columns 1 through 3 17.278609048277868 17.277724710546092 17.279658142557587 Columns 4 through 6 17.279658229217087 17.219381240184841 17.2796582292086502、用雅可比迭代求解线性方程组，其中随机取。要求使用函数型M文件，并有对其迭代格式的收敛性进行判断的功能。雅可比迭代M文件；function x,m=yakebi（A,b,x0,to