matlab练习题及答案.docVIP

一、填空题(本题满分30分，每空3分) 1.设有程序 A=[2,0,1,0;1,0,1,5;0, 6 ,2, 3];B=[1,2,-1;1,3,1];C=A(1:3,[1,2]); D=size(C)- length(B(2,:));F= C.*B’ 将下列命令的运行结果填在横线上 D = ; F= . 2. A=[2,3,2,1;1,2,2,1; 1,2,3,-2]; A1=sum(AA(2,1));A(2,:)=[],A2=A.^2-4 A1= ; A2= ; 3.P=[1,2,5,4;2,2,1,3];max(P,[ ],1)= ;mean(P,2)= . 二、(本题满分10分)试编写计算程序画出曲线与曲面图形 %1 %2 t=0:0.1: 2*pi; [x,y]=meshgrid(-2:0.1:2); x= sin(t).^3; z=(x.^2+1).* exp(-x.^2- y.^2+)+x.^2.*y; y= cos (t).^3; mesh(x,y,z) plot(x,y) 三、（本题满分12分）编写程序解决下列问题： （1）试产生100(4阶矩阵使其元素在0—100之间的随机整数,可将此矩阵理解为100名同学的四门课程(课程1, 课程2, 课程3, 课程4)的成绩. （2）计算出每门课程和每位同学的平均成绩; （3）统计出总的优秀人次(90分及以上为优秀); （4）统计出四门课程成绩至少有两门课程不及格的人数. 解：(1) a=fix((100-0+1)*rand(100,4)); (2) mean(a) mean(a,2) (3)sum(sum(a=90,2)==4) (4)sum(sum(a60,2)=2) 四、（本题满10分） (1). function y=fun1(x,n) (2). x=0:0.1:2*pi; y=0; y1= fun1 (x,20); for k=1:n y2= fun1 (x,50); y=y+ ((-1)^(k+1)*sin(k*x))/k; subplot(2,1,1),plot(x,y1) end subplot(2,1,2),plot(x,y2) y 五、(本题满分10分) 试利用微分方程的数值解法写出下列微分方程的求解步骤和MATLAB程序. .(t的变化区间为[0,10]). 第一步：先将高阶微分方程转化为一阶微分方程组： 选择状态变量，则原方程化为 第二步：建立函数文件 function dxdt = fun2(t,x) dxdt = [-3*t*x(1)+ x(2)+t; x(1)]; 第三步：求解微分方程，命令如下： [t,x]=ode45(@fun2,[0,10],[0;2]) 5.求级数之和 六、参考程序： solve(exp(x)=x^2+2) 2．syms x s A=[x*sin(x),log(x);exp(-x),1/(s+x)] diff(A) int(A) A*A 3．syms x limit((cos(x)-cos(x)^(1/2))/x/(exp(sin(x))-1)) 4．dsolve(D4x+x=2*exp(t),x(0)=1,Dx(0)=1,D2x(0)=1,D3x(0)=1) 5.syms n; s=symsum(1/(n^2+1),1,inf) 解一：function y=fc(x) y(1)=4*x(1).^2+x(2).^2+2*x(1).*x(2)-x(2)-2 y(2)=2*x(1).^2+x(2).^2+3*x(1).*x(2)-3 y=[y(1) y(2)]; x0=[1 1]; fsolve(@fc,x0) 解二: syms x1 x2 [x1 x2]=solve(‘4*x1^2+ x2^2+2* x1* x2- x2-2=0’,’ 2*x1^2+ x2^2+3* x1* x2-3=0’,’ x1’,’ x2’) 八．已知，