MATLAB程序设计与应用_5.ppt

Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;第五讲 MATLAB数值计算 符号计算 授课教师：田 鹏;数值计算+符号计算;复习 分子物理学绘图;程序： 1）子程序（mksw.m）： function f=mksw(T,mu,v) R=8.31; %气体常数 k=1.381*10^(-23); %玻尔茨曼常数 NA=6.022*10^23; %阿伏伽德罗数 m=mu/NA; %分子质量 f=4*pi*(m/(2*pi*k*T)).^(3/2) ... .*exp(-m*v.^2./(2*k*T)).*V.^2; %速度分布率;2）主程序：;结果：;复习．我国古代数学家张丘在“算经”里提出一个世界数学史上有名的百鸡问题：鸡翁一，值钱五，鸡母一，值钱三，幼鸡三，值钱一，百钱买百鸡，问各几何？ ;方法1-1：;方法1-1：运用符号函数 syms+solve;方法1-2：;方法1-3：;方法2-1：;一.单变量函数的计算和绘图 例5.1：已知 要求以0.01秒为间隔，绘出y及其导数的曲线. 分析：用diff(y,n)求Dy ， 每求导一次，y的维数减一。 Dy=diff(y)结果为Dy=y1-y2, 故 y’=Dy/Dx= diff(y)/Dx;b=0.1;t=0:b:1.5;w= 4*sqrt(3); y=5*sqrt(3)*exp(-2*t).*sin(w*t+pi/3); plot(t,y); title(‘单变量绘图’); xlabel(x);ylabel(y(t)); grid on;hold on; Dy=diff(y)/b; plot(t(1:length(t)-1),Dy,’*’);plot(t(2:length(t)),Dy,’p’) y1=-10*sqrt(3)*exp(-2*t)… .*sin(w*t+pi/3)+60*exp(-2*t)… .*cos(w*t+pi/3); hold on ; plot(t,y1,‘r) legend(y,Dy1, Dy2, y1) ;Dy和y1不重合呢？;例5.2.求点u=(1,2,3)到平面 3x-2y+z=4的距离;u=input(点的坐标u[]=); v=input(直线方程系数v[]=); r=abs(u*v(1:3)+v(4))/sqrt(sum(v(1:3).^2)) ;5.1 函数、极限和导数;5.1 函数、极限和导数;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;补充：三次抛物线方程y=ax3+cx，讨论参数a，c对其图形的影响。;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;5.1 函数、极限和导数;5.1 函数、极限和导数;5.2 空间解析几何;程序——求切点 x=sym(x); y=x^3+3*x-2; f=diff(y); x1=solve(f-100); y1=x1.^3+3*x1-2; c=y1-100*x1;hold on; plot(eval(x1),eval(y1),*); x=-100:0.1:100; y1=100*x+ eval(c(1)); y2= 100*x+ eval(c(2)); plot(x,y1,x,y2,r) y3=x.^3+3*x-2; y4=100*x-1; plot(x,y3,x,y4,r) axis([-10,10,-500,500]);5.2 空间解析几何;5.2 空间解析几何;5.2 空间解析几何;5.2 空间解析几何;5.2 空间解