MATLAB绘图及数值计算.doc

信管102班 蓝梦源 学号：1011200240实验三 MATLAB绘图及数值计算实验目的1. 熟悉二维和三维绘图函数命令，了解图形修饰方法。2. 掌握数值导数、数值积分、代数方程求解、常微分方程数值求解的方法。实验内容1. 在同一图形窗口绘制，的曲线，要求用不同的颜色和线型、数据点标记与x轴范围为，y轴范围为；clcx=0:pi/100:2*pi;y=-2,2;y1=sin(x);y2=cos(x);plot(x,y1,k:,x,y2,r-)2. 在同一图形窗口绘制函数和在区间x=-20，20内的图形。clcx=-20:0.1:20;y1=200*sin(x)./x;y2=x.2;plot(x,y1,k:,x,y2,b-)3. 绘制三维曲线：，。t=0:pi/1000:2*pi;x=sqrt(t).*sin(2*t);y=sqrt(t).*cos(2*t);z=0.5*t;plot3(x,y,z)4. 绘制三维网格曲面：，。clcx,y=gmeshgrid(x,y);x=-10:0.1:10;y=-8:0.1:8;z=sin(sqrt(x.2+y.2)./sqrt(x.2+y.2);mesh(x,y,z)5. 对多项式进行下面计算：（1）求导数；（2）求在的值；（3）求其根。clcp=1 0 1 -2 0 1 0 0 1d=polyder(p)y=polyval(p,1.32)r=roots(p)p = 1 0 1 -2 0 1 0 0 1d = 8 0 6 -10 0 3 0 0y = 9.7919r = -0.5159 + 1.4187i -0.5159 - 1.4187i 0.9776 + 0.3710i 0.9776 - 0.3710i -0.7267 + 0.3617i -0.7267 - 0.3617i 0.2650 + 0.7341i 0.2650 - 0.7341i6. 求函数在1，3的1阶差分，并计算在的导数的近似值。 x=linspace(1:0.5:2);y=sqrt(x3+2*x2-x+12)+sqrt(x+5)+5*x+2dy=diff(y)x=2,2.1y=16+sqrt7,sqrt(2.13+2*2.12-2.1+12)+sqrt(7.1)+12.5dy=diff(y)d=dy/0.01ans = 557. 计算数值积分：f=inline(sin(x).5.*sin(5*x)I1=quad(f,0,pi)f = Inline function: f(x) = sin(x).5.*sin(5*x)I1 = 0.09828. 求解线性方程组：A=2 1 -5 1;1 -3 0 -6;0 2 -1 2;1 4 -7 6B=8 9 -5 0X=A/BA = 2 1 -5 1 1 -3 0 -6 0 2 -1 2 1 4 -7 6B = 8 9 -5 0X = 0.2941 -0.1118 0.1353 0.46479. 求微分方程的数值解：，其中，。y=dsolve(d2y-10*y+9*y=exp(2*x),y(0)=33/7,dy(0)=6/7,x) function dy=vdp(x,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=exp(2*x)-9*y(1)+10*y(2); x,y=ode15s(vdp,0,100,33/7,6/7function fy=vdp(x,t) fy=t(2);exp(2*t)-