资源描述：
《MATLAB绘图及数值计算》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、信管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*s
2、in(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.绘制三维网格曲面：，，。clc[x,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）求其
3、根。clcp=[101-201001]d=polyder(p)y=polyval(p,1.32)r=roots(p)p=101-201001d=806-100300y=9.7919r=-0.5159+1.4187i-0.5159-1.4187i0.9776+0.3710i0.9776-0.3710i-0.7267+0.3617i-0.7267-0.3617i0.2650+0.7341i0.2650-0.7341i6.求函数在[1，3]的1阶差分，并计算在的导数的近似值。x=linspace(1:0.5:2);y=sqrt(x^3+2*x^2-x+12)+sq
4、rt(x+5)+5*x+2dy=diff(y)x=[2,2.1]y[=[16+sqrt7,sqrt(2.1^3+2*2.1^2-2.1+12)+sqrt(7.1)+12.5]dy=diff(y)d=dy/0.01ans=557.计算数值积分：f=inline('sin(x).^5.*sin(5*x)')I1=quad(f,0,pi)f=Inlinefunction:f(x)=sin(x).^5.*sin(5*x)I1=0.09828.求解线性方程组：A=[21-51;1-30-6;02-12;14-76]B=[89-50]X=A/BA=21-511-30-6
5、02-1214-76B=89-50X=0.2941-0.11180.13530.46479.求微分方程的数值解：，其中，。y=dsolve('d2y-10*y+9*y=exp^(2*x)','y(0)=33/7,dy(0)=6/7','x')functiondy=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/7]functionfy=vdp(x,t)fy=[t(2);exp^(2*t)-9*t(1)+10*t
6、(2)];y0=[33/7;6/7][x,t]=ode45(@vdp,[0,40],y0]y=t(:,1);dy=t(:,2)ans=0
