数学实验考试题精选.docVIP

数学实验考试 1、求微分方程在初始条件y（0）=0，下的特解，并画出解函数的图形。 程序： dydx?=?@(x,y)?[y(2);?sin(2*x)?-?y(1)?]; xspan?=?[0,?10];y0?=?[0;?1];[x,y]?=?ode45(dydx,?xspan,?y0);%?Plot?of?the?solutionplot(x,y(:,1))xlabel(x)ylabel(solution?y 运行结果： 2、求代数方程组的解。 程序： syms x y; eq1=2x+3y; eq2=4x^2+3y-1; [x,y]=solve(eq1,eq2) x=eval(x) y=eval(y) 运行结果： x = 1/4-1/4*5^(1/2) 1/4*5^(1/2)+1/4 y = -1/6+1/6*5^(1/2) -1/6-1/6*5^(1/2) x = -0.3090 0.8090 y = 0.2060 -0.5393 3、.用ode23(或ode45)求解处置问题 ，y(0)=0，的数值解。 分析：令，，则，。 编写函数文件fun.m，作为ode45输入函数。 程序： 编写M文件.m function?dy=cdq(x,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=-2*y(2)-y(1)+cos(x); 编写M文件cdq1.m function?dY=cdq1(x,Y) dY=[0?1;-1?-2]*Y+[0;1]*cos(x); 命令窗口 ?subplot(1,2,1),ode23(@?cdq,[0,2*pi],[0,3/2]),grid ?subplot(1,2,2),ode45(@?cdq,[0,2*pi],[0,3/2]),grid ?[x,Y]=ode23(@?cdq1,[0,2*pi],[0,3/2]),grid 求微分方程的数值解，并画出图形。 程序： 运行结果： 4、求平面3x+5y+6z=0和平面x-3y-6z-1=0的交线与球面的交点。 分析：要求两平面的交线与球面的交点，将它们联立成方程组求解即可。 编写函数文件fun2.m，建立方程左边的向量值函数。 程序： r=3; the=0:pi/100:pi; phi=0:pi/100:2*pi; [the,phi]=meshgrid(the,phi); z1=r.*sin(phi); x1=r.*cos(phi).*cos(the); y1=r.*cos(phi).*sin(the); mesh(x1,y1,z1); axis equal xlabel(x);ylabel(y);zlabel(z); view(123,40); set(gcf,units,normalized,position,[0,0,1,1]); hold on; syms x y z fx1=3*x+5*y+6*z; fx2=x-3*y-6*z-1; fx11=subs(fx1,x,-1.5); fx22=subs(fx2,x,-1.5); [y11,z22]=solve(fx11,fx22) fx111=subs(fx1,x,2); fx222=subs(fx2,x,2); [y111,z222]=solve(fx111,fx222) plot3([-1.5,2],[eval(y11),eval(y111)],[eval(z22),eval(z222)],r,linewidth,3); 运行结果： y11 = 7/2 z22 = -13/6 y111 = -7/2 z222 = 23/12 求解方程组，得到交点。 程序： 运行结果： 5、已知数据，见表3-11-1. 表3-11-1 X 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 y 0.3 0.5 1 1.4 1.6 1.9 0.6 0.4 0.8 1.5 2 画出用线性、三次样条和三次多项式插值所得[0，1]区间内的曲线图，并求当=0.25,0.35,0.45时的的值。 画曲线图。 程序： 代码如下：x=0:0.1:1; y=[0.3 0.5 1 1.4 1.6 1.9 0.6 0.4 0.8 1.5 2]; ` subplot(2,2,1) ` plot(x,y,‘b+’) ` title(‘散点图’) 线性插值： 三次样条插值：