MATLAB绘图及数值计算.docVIP

信管102班 蓝梦源 学号：1011200240 实验三 MATLAB绘图及数值计算 实验目的 1. 熟悉二维和三维绘图函数命令，了解图形修饰方法。 2. 掌握数值导数、数值积分、代数方程求解、常微分方程数值求解的方法。 实验内容 1. 在同一图形窗口绘制，的曲线，要求用不同的颜色和线型、数据点标记与 x轴范围为，y轴范围为； clc x=[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]内的图形。 clc x=-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. 绘制三维网格曲面：，，。 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）求其根。 clc p=[1 0 1 -2 0 1 0 0 1] d=polyder(p) y=polyval(p,1.32) r=roots(p) p = 1 0 1 -2 0 1 0 0 1 d = 8 0 6 -10 0 3 0 0 y = 9.7919 r = -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.7341i 求函数在[1，3]的1阶差分，并计算在的导数的近似值。 x=linspace(1:0.5:2); y=sqrt(x^3+2*x^2-x+12)+sqrt(x+5)+5*x+2 dy=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.01 ans = 55 7. 计算数值积分： 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.0982 8. 求解线性方程组： A=[2 1 -5 1;1 -3 0 -6;0 2 -1 2;1 4 -7 6] B=[8 9 -5 0] X=A/B A = 2 1 -5 1 1 -3 0 -6 0 2 -1 2 1 4 -7 6 B = 8 9 -5 0 X = 0.2941 -0.1118 0.1353 0.4647 求微分方程的数值解：，其中，。 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/7] function fy=vdp(x,t) fy=[t(2);exp^(2*t)-9*t(1)+10*t(2)]; y0=[33/7;6/7] [x,t]=ode45(@vdp,[0,40],y0] y=t(:,1); dy=t(:,2) ans = 0