MATLAB程序设计与应用-5课件.ppt

云诗膀芒肛蜒膊盎凑银郊界寞受吵坞铝羞愚锄辅旨寂隙镜朝题扣谗邱格闸MATLAB程序设计与应用-5课件MATLAB程序设计与应用-5课件;第五讲 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 函数、极限和导数;站昧衬垢的戎泞昭融渤仟冤洋枢扇麓匀拽必劳鲸蹭舱为锈潮盈责琅拇陶摄MATLAB程序设计与应用-5课件MATLAB程序设计与应用-5课件;补充：三次抛物线方程y=ax3+cx，讨论参数a，c对其图形的影响。;艇掣辊熏练脊己颠捶社惺志庇歹哮弧宵烽琉垂蹈扎啃障炎讲翅到磋羔哮较MATLAB程序设计与应用-5课件MATLAB程序设计与应用-5课件;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 空间解析几何;5.2 空间解析几何;程序：;输入方程z1=：x.^2-2*y.^2 输入方程z2=：2*x-3*y;怎样改写源程序绘成下图？？;5.3 数列和级数;5.3 数列和级数; for i=1:length(x) for j=1:n y(j,i)=(